資源描述：
《二維最大熵閾值分割法.doc》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫。
1、二維最大熵閾值分割算法[引用]杜峰，施文康，鄧勇等：《一種快速紅外圖像分割方法》1.二維最大熵閾值分割熵是平均信息量的表征。二維最大熵法是基于圖像二維直方圖。圖像二維直方圖定義如下：其中表示圖像大小，表示圖像灰度值為i，鄰域灰度平均值為j的像素個(gè)數(shù)。通常二維直方圖的平面示意圖可以用下圖1表示：圖1二維直方圖平面示意圖灰階L均值L0ts1234其中區(qū)域1和2表示背景和目標(biāo)像素，區(qū)域3和4通常表示邊界和噪聲信息。閾值向量(t,s)，t表示灰度值，s表示像素鄰域均值（通常是8鄰域）。對于L個(gè)灰度級(jí)的圖像，設(shè)在閾值(t,s
2、)定義區(qū)域1和2的概率P1，P2：，定義二維離散熵H的一般表示：對各區(qū)域概率進(jìn)行歸一化處理可得區(qū)域1的二維熵：同理區(qū)域2的二維熵：其中，H1，H2為：，那么整個(gè)圖像中目標(biāo)和背景熵之和的函數(shù)根據(jù)最大熵原則，存在最佳的閾值向量滿足條件：圖2顯示了一幅圖像的二維直方圖說明了背景和目標(biāo)的主要分布情況，其中圖2(b)橫坐標(biāo)表示鄰域的均值，縱坐標(biāo)表示灰度值分布：圖2目標(biāo)與背景的二維直方圖分布情況(a)原始紅外圖(b)二維直方圖的平面分布(c)二維直方圖的空間分布2.微粒群尋優(yōu)算法（PSO）PSO最早由Kenredy和Eberh
3、art于1995年提出。PSO把優(yōu)化問題的潛在解都當(dāng)做解空間的粒子，所有粒子都有一個(gè)適應(yīng)值（適應(yīng)值由被優(yōu)化函數(shù)決定），每個(gè)粒子還有一個(gè)速度決定它們飛翔的方向和距離。然后粒子們就追隨當(dāng)前的最優(yōu)粒子在解空間中搜索，初始化為一群隨機(jī)粒子（隨機(jī)解）然后通過迭代找到最優(yōu)解。最后在每一次迭代中粒子通過跟蹤兩個(gè)極值來更新自己，第一個(gè)就是粒子本身所找到的最優(yōu)解，稱為個(gè)體極值；另一個(gè)極值是整個(gè)種群目前找到的最優(yōu)解，稱為全局極值。本文的目標(biāo)是要找到滿足最大熵原則的最優(yōu)解，下面以圖文方式解釋PSO算法步驟原理：均值灰度值0圖3隨機(jī)初始化
4、的粒子群位置第1步：第1次迭代→如圖3在解空間有效范圍內(nèi)選定m個(gè)隨機(jī)解（即粒子）并初始化，如：、……其中為最優(yōu)解位置向量。①計(jì)算每個(gè)粒子的適應(yīng)度，即目標(biāo)函數(shù)的熵，。②計(jì)算當(dāng)前粒子群的全局最優(yōu)解(熵)及其對應(yīng)位置：③計(jì)算n次迭代后每個(gè)粒子自身找到的最優(yōu)解(熵)及其位置：；其中n表示迭代次數(shù)，Die表示最大迭代次數(shù)。首次迭代(n=1)時(shí)單個(gè)粒子最優(yōu)值即為其初始化時(shí)的隨機(jī)值。圖4粒子群的運(yùn)動(dòng)速度和更新位置均值灰度值0第2步：第n次迭代→如圖4更新粒子速度向量和位置，粒子運(yùn)動(dòng)服從如下方程：其中、為隨機(jī)數(shù)，服從(0,1)之間
5、的均分布，、為學(xué)習(xí)因子，通常，是慣性系數(shù)。表示第個(gè)粒子的位置向量(即)，表示第個(gè)粒子的運(yùn)動(dòng)速度，表示第個(gè)粒子自身的最優(yōu)位置。表示整個(gè)粒子群全局最優(yōu)位置。3.實(shí)驗(yàn)結(jié)果如圖5顯示了二維最大熵閾值分割的結(jié)果。(a)原始紅外圖圖5二維最大熵閾值分割結(jié)果(b)閾值分割后的二值化圖其中圖2(b)也就是圖5(a)對應(yīng)的二維直方圖分布。如何在圖2(b)找到最優(yōu)的閾值向量使得目標(biāo)圖像熵最大，一個(gè)最直接的方法就是窮盡搜索法。窮盡搜索法無目的性而且計(jì)算量大，需要進(jìn)行256×256次計(jì)算。本文采用PSO算法搜索最佳閾值，在實(shí)驗(yàn)中，令粒子群
6、為15個(gè)，迭代次數(shù)30，c1=c2=2，w=0.35。圖6顯示了粒子群在每次迭代中達(dá)到的局部最優(yōu)熵。完成整個(gè)迭代尋優(yōu)過程粒子群找到的全局最優(yōu)閾值向量為(105，103)全局最優(yōu)熵。從圖6可以看出：第14代的粒子群局部最優(yōu)熵就達(dá)到了5.1484，說明了迭代到第14代就至少有一個(gè)粒子尋找到了全局最優(yōu)位置。從第16代到30代之間，粒子群局部最優(yōu)熵一直保持5.1484，說明此時(shí)粒子群中至少且總有一個(gè)粒子到達(dá)了全局最優(yōu)位置。因此整個(gè)迭代過程中，尋找到全局最優(yōu)位置PSO的計(jì)算量為16×15次。圖6粒子群在各次迭代中的局部最優(yōu)熵
7、圖7首次迭代時(shí)初始化的隨機(jī)粒子分布圖7顯示了在首次迭代時(shí)初始化的15個(gè)隨機(jī)粒子位置分布圖，其中橫坐標(biāo)表示均值(s)，縱坐標(biāo)表示灰度值(t)。圖8顯示了在第14次迭代后粒子群的位置分布以及各個(gè)粒子的位置坐標(biāo)。從圖8可以看出第6個(gè)粒子首次尋找到全局最優(yōu)位置(105，103)。st10710211199104103105100110103105103106102106102110100111991109810899102111107102104102圖8第14次迭代粒子群的位置分布及其位置坐標(biāo)值圖9顯示了第30次迭代后粒
8、子群的分布情況。從圖中可以看出此時(shí)大部分粒子都收斂于全局最優(yōu)位置。st105103105103105103105103105102105103105103105103105103105103105103105103105103105103104103圖9第30次迭代粒子群的位置分布及其位置坐標(biāo)值總之，PSO算法中的粒子群從初始隨機(jī)位置經(jīng)過各次迭代過程遵照粒子
