本發(fā)明涉及圖像處理技術(shù)領(lǐng)域，具體涉及基于多尺度特征的紅外噪聲抑制和細(xì)節(jié)增強(qiáng)。

背景技術(shù)：

紅外圖像通過探測目標(biāo)物體表面和背景之間的熱輻射差而成像，其成像不受外界環(huán)境光明亮的影響，具備較強(qiáng)的檢測和識別偽裝能力；另外紅外熱像儀無需向外輻射能量，只需被動接受外界輻射，不容易被敵方探測器所偵查和監(jiān)視，具備較強(qiáng)的抗干擾和偽裝能力。因而紅外熱像儀被廣泛應(yīng)用在場所監(jiān)控^[1]、醫(yī)學(xué)成像^[2]和軍事遙感成像^[3]等方面。但同時受到應(yīng)用環(huán)境和采集設(shè)備的影響，紅外圖像往往存在：(1)灰度分布集中，目標(biāo)物體和背景環(huán)境對比度較低；(2)所探測目標(biāo)的輪廓和形狀不是很清晰，且物體的細(xì)節(jié)模糊；(3)同時伴隨著環(huán)境光線變換和熱像儀內(nèi)部電子元器件的運(yùn)行，時常會出現(xiàn)大量的散粒噪聲，使得場景中細(xì)節(jié)不易觀察，圖像整體視覺效果差^[4]。因此，研究一種既可有效抑制噪聲，又能很好保留細(xì)節(jié)和邊緣等重要信息的處理方法是一項非常有意義的工作，為后續(xù)的紅外目標(biāo)探測、視頻跟蹤和目標(biāo)識別等應(yīng)用提供重要的技術(shù)支撐。

為有效增強(qiáng)紅外圖像細(xì)節(jié)和提升紅外圖像對比度，目前最簡單有效的方法是全局直方圖均衡(he)?；叶戎狈綀D能直觀反映圖像灰度分布情況，he方法是從統(tǒng)計學(xué)的觀點(diǎn)為提高圖像對比度，改變圖像中灰度概率分布，使紅外圖像灰度層次分明^[5]。

給定一幅數(shù)字圖像，其離散的表達(dá)形式為:

式中：灰度級范圍為[0,l-1]，rk是第k級灰度值，nk是圖像中灰度為rk的像素個數(shù)，p(rk)為第k級灰度級上的概率密度。

數(shù)字圖像的直方圖均衡變換函數(shù)可表示為：

處理后，圖像像素灰度的動態(tài)范圍增加，圖像的對比度和視覺效果得到提高。

相對全局直方圖均衡，后續(xù)還有基于對比度限制自適應(yīng)局部直方圖均衡(clahe)^[6]。其中心思想不同于全局直方圖均衡，是對局部塊進(jìn)行直方圖均衡使變換后的直方圖分布匹配事先設(shè)計的直方圖分布，可以在一定程度上降低圖像平坦區(qū)域噪聲的放大。

全局直方圖均衡和基于對比度限制的自適應(yīng)局部直方圖均衡增強(qiáng)紅外圖像的效果如圖1所示。圖1(a)原始紅外圖像,(b)全局直方圖均衡結(jié)果,(c)為基于對比度限制的自適應(yīng)局部直方圖結(jié)果。

這兩種紅外圖像增強(qiáng)算法處理高信噪比或具有少量背景雜波紅外圖像都具有不錯的視覺效果，全局直方圖均衡能大大提升圖像的整體對比度，基于對比度限制的自適應(yīng)局部直方圖均衡能很好地凸顯紅外圖像局部細(xì)節(jié)。

然而，處理無背景噪聲紅外圖像時，全局直方圖均衡產(chǎn)生的過暗或過亮的區(qū)域會淹沒了目標(biāo)局部細(xì)節(jié)；對比度限制的自適應(yīng)直方圖容易在處理結(jié)果中產(chǎn)生一些令人意料之外的細(xì)節(jié)，例如圖1(c)中的草地區(qū)域產(chǎn)生類似于雜草的細(xì)節(jié)。同時，更重要的是，當(dāng)上述兩種增強(qiáng)算法直接處理背景噪聲較大的紅外圖像時，會顯著提高噪聲水平淹沒圖像的有用信息。如圖2所示，原始紅外圖像的背景噪聲被放大，處理結(jié)果的細(xì)節(jié)和結(jié)構(gòu)信息被破壞。圖2(a)有噪紅外圖像,(b)全局直方圖均衡放大噪聲,(c)對比度限制的自適應(yīng)局部直方圖放大噪聲。

可知，對于存在一定噪聲的紅外圖像，對噪聲進(jìn)行抑制是十分有必要的。為有效濾除紅外圖像的背景噪聲，wavelet和curvelet將有噪圖像分解成多尺度具有時頻特性的成分，對每個尺度設(shè)定軟閾值降低噪聲成分達(dá)到抑制噪聲的目的^[7,8]。nlm認(rèn)為在自然圖像中往往存在一些空間位置不同但彼此相似的圖像片，由當(dāng)前像素點(diǎn)與其鄰域內(nèi)其它像素點(diǎn)為中心的圖像片之間的相似性來決定加權(quán)系數(shù)，能夠較好地消除像素的噪聲^[9]。bm3d把具有相似結(jié)構(gòu)的二維圖像片組合形成三維數(shù)組，然后這些三維數(shù)組進(jìn)行協(xié)同濾波，并把結(jié)果返回到圖像中，從而得到去噪后的圖像，方法不僅有較高的信噪比，而且視覺效果也很好^[10]。rudin提出噪聲圖像的總變分要明顯大于無噪聲圖像的總變分的先驗，提出了tv去噪模型。tv模型是一種基于最優(yōu)化理論有約束的擴(kuò)散濾波方法，具有全局最優(yōu)解和較好保持圖像結(jié)構(gòu)信息的優(yōu)點(diǎn)^[11]。

對于結(jié)構(gòu)信息不豐富的紅外圖像，只考慮單像素點(diǎn)的wavelet和curvelet，很容易將結(jié)構(gòu)和細(xì)節(jié)信息丟失?；趫D像結(jié)構(gòu)信息相似性的nlm和bm3d濾波方法忽略了圖像的不同性質(zhì)區(qū)域信息具有不同強(qiáng)度的相關(guān)性：在圖像平坦區(qū)域，像素間的相關(guān)性較強(qiáng)，而在細(xì)節(jié)區(qū)域，相關(guān)性則較弱。而tv但在噪聲環(huán)境下梯度值大小不足以作為區(qū)分噪聲和區(qū)域邊界的特征，去噪后的圖像存在平坦區(qū)域不夠光滑仍有明顯的噪聲和塊效應(yīng)?，F(xiàn)有的技術(shù)二對紅外圖像的背景噪聲抑制效果如圖3所示。圖3(a)有噪紅外圖像,(b)wt抑制噪聲結(jié)果,(c)bm3d抑制噪聲結(jié)果,(d)tv抑制噪聲結(jié)果。

從上述分析可知，解決邊緣過增強(qiáng)和抑制噪聲時的邊緣模糊問題，都需要對圖像的結(jié)構(gòu)信息正確地刻畫：一是在無噪聲干擾下能精細(xì)表征不同意義的圖像區(qū)域，二是在噪聲干擾下仍能準(zhǔn)確刻畫邊緣的位置和強(qiáng)度。

參考文獻(xiàn)：

[1]a.rossi,n.acito,m.diani,g.corsini,dynamicrangereductionandcontrastadjustmentofinfraredimagesinsurveillancescenarios,opt.eng.52(10)(2013).102002–102002.

[2]k.liang,y.ma,y.xie,etal.,anewadaptivecontrastenhancementalgorithmforinfraredimagesbasedondoubleplateaushistogramequalization,infraredphys.technol.55(4)(2012)309–315.

[3]x.bai,f.zhou,analysisofnewtop-hattransformationandtheapplicationforinfrareddimsmalltargetdetection,patternrecogn.43(2010)2145–2156.

[4]q.chen,l.f.bai,b.m.zhang,histogramdoubleequalizationininfraredimage,j.infraredmillimeterwaves22(6)(2003)428–430.

[5]k.s.sim,c.p.tso,y.y.ta,recursivesub-imagehistogramequalizationappliedtograyscaleimages,patterrecognit.lett.28(10)(2007)1209–1221.

[6]u.qidwai,infraredimageenhancementusingh∞boundsforsurveillanceapplications,ieeetrans.imageprocess.17(2008)1274–1282.

[7]k.dabov,a.foi,v.katkovnik,k.egiazarian,imagedenoisingbysparse3dtransform-domaincollaborativefiltering,trans.imageprocess.16(8)(2007)2080–2095.

[8]l.rudin,s.osher,e.fatemi,nonlineartotalvariationbasednoiseremovalalgorithms,physicad,60(1992)259-268.

[9]l.tang,l.liu,j.h.su,modelingandsimulationresearchofinfraredimagenoise,infraredtechnol.36(2014)542-548.

[10]j.wu,f.liu,l.c.jiao,x.d.wang,compressivesensingsarimagereconstructionbasedonbayesianframeworkandevolutionarycomputation,trans.imageprocess.20(7)(2011):1904–1911.

[11]a.m.reza,realizationofthecontrastlimitedadaptivehistogramequalization(clahe)forreal-timeimageenhancement,j.vlsisignalprocess.syst.signal,imagevideotechnol.38(2004)35–44.

[12]x.h.liu,m.tanaka,m.okutomi,noiselevelestimationusingweaktexturedpatchesofasinglenoisyimage,ieeeint.conf.imageprocess.(2012)665-668.

技術(shù)實現(xiàn)要素：

針對現(xiàn)有技術(shù)的不足，本發(fā)明旨在提供一種基于多尺度特征的紅外噪聲抑制和細(xì)節(jié)增強(qiáng)方法，通過分析紅外圖像背景噪聲特性，提出了一種新的在噪聲干擾下仍能準(zhǔn)確刻畫圖像結(jié)構(gòu)信息的先驗?zāi)Ｐ?，?jīng)過連續(xù)多尺度提取優(yōu)化的先驗?zāi)Ｐ鸵肓素惾~斯和最大后驗估計準(zhǔn)則迭代平滑噪聲。同時為增強(qiáng)邊緣，利用在不同尺度下平滑提取的邊緣信息有針對性地增強(qiáng)圖像?？梢詮V泛地應(yīng)用于安防監(jiān)控、道路交通管理和紅外精準(zhǔn)制導(dǎo)。

為了實現(xiàn)上述目的，本發(fā)明采用如下技術(shù)方案：

一種基于多尺度特征的紅外噪聲抑制和細(xì)節(jié)增強(qiáng)方法，包括如下步驟：

s1搭建刻畫紅外圖像特征的先驗?zāi)Ｐ停?/p>

采用貝葉斯準(zhǔn)則作為抑制噪聲的模型，對于一幅清晰紅外圖像u受到噪聲n污染，退化圖像v表示為:

v＝u+n；

n是零均值的高斯白噪聲，對于退化圖像v和原始圖像u，貝葉斯準(zhǔn)則表述為：

p(u|v)∝p(v|u)p(u)；

其中，p(v|u)為似然密度函數(shù)，p(u)為原始圖像的先驗密度函數(shù)，p(u|v)稱為后驗密度函數(shù)；根據(jù)map估計準(zhǔn)則，即在給定觀測數(shù)據(jù)v的條件下，使后驗密度函數(shù)p(u|v)最大的估計量為：

對取負(fù)對數(shù)轉(zhuǎn)化為：

退化圖像v和原始圖像u之間的誤差滿足高斯白噪聲分布，因此似然密度函數(shù)p(v|u)表示為:

對于退化圖像v，像素點(diǎn)i,j∈p方塊內(nèi)的相似性取決于像素點(diǎn)之間的亮度和梯度差異，將亮度值差異的歐氏距離和梯度值差異的絕對距離之和的衰減函數(shù)定義為兩點(diǎn)之間的相似性r(i,j)：

其中，h和ρ均為大于0的常數(shù)，分別表示亮度值差異和梯度值差異的衰減系數(shù)；相似性0＜r(i,j)≤1，值越大，表示兩像素點(diǎn)間的相似程度越大；

在以i為中心像素點(diǎn)的鄰域內(nèi)，根據(jù)其他像素點(diǎn)與中心像素點(diǎn)之間相似性估計中心像素在鄰域p內(nèi)的隸屬性，并以此刻畫紅外圖像中的形狀先驗，定義中心像素點(diǎn)i在鄰域內(nèi)的隸屬度w(i)為：

z(i)是規(guī)范化函數(shù)，σ是高斯模板的標(biāo)準(zhǔn)差，中心像素點(diǎn)的隸屬度0＜w(i)≤1，w(v)＝{w(i)丨i∈v}為圖像v的隸屬度運(yùn)算，構(gòu)成了一幅與原圖像尺寸一致的隸屬度圖；

定義基于隸屬度的約束控制函數(shù)如下：

c(w)＝c+(1-c)·sin(π·w)；

其中，c為約束控制函數(shù)的調(diào)整因子，其值介于0到1之間；w為由像素點(diǎn)隸屬度{w(i)}構(gòu)成對應(yīng)的隸屬度圖；

將原始圖像的先驗密度函數(shù)約束為:

后驗密度函數(shù)p(u|v)最大的估計量轉(zhuǎn)化為:

β為正則化系數(shù)；根據(jù)原始圖像的噪聲水平，噪聲水平越高正則化系數(shù)取值越??；

s2基于多尺度特征的紅外噪聲抑制：采用分裂bregman方法求解同時利用多尺度特征不斷更新迭代貝葉斯過程，獲取最優(yōu)估計；

s3基于多尺度特征的紅外細(xì)節(jié)增強(qiáng)：

將隸屬度圖和多尺度提取引入如下式所示的非銳化掩模算法，自適應(yīng)地控制增益系數(shù)：

其中，uen是最終輸出的增強(qiáng)圖像，un是步驟s2中獲得的最優(yōu)估計，un-j和un-j+1分別表示步驟s2中第n-j次和n-j+1次迭代的得到的估計，g(wn-j+1)是關(guān)于un-j+1的隸屬度圖wn-j+1的增益函數(shù)；g(w)其分段表示見下：

分段點(diǎn)

需要說明的是，步驟s1中，鄰域p尺寸為7。

需要說明的是，步驟s2中，利用分裂bregman方法來求解的方法為：

引入輔助量g，其中d是梯度算子；轉(zhuǎn)換為有關(guān)u和g的優(yōu)化約束模型:

添加關(guān)于u和g的二次懲罰項，得到與s.t.g＝c(w)d·u對應(yīng)的無約束優(yōu)化模型:

先固定gn利用關(guān)于u求導(dǎo)得到un+1，再固定un+1利用收縮技術(shù)得到gn+1，gn表示第n次迭代得到的g的更新值，un+1為第n+1次迭代得到的u的更新值；更新過程如下式所示：

m＝c(wn+1)d·un+1

其中，wn＝w(un)為n次迭代后圖像的隸屬度圖；迭代終止后得到的最終結(jié)果un即為求解結(jié)果。

進(jìn)一步需要說明的是，步驟s2中，需要利用多尺度特征不斷更新迭代過程中的隸屬度圖，具體方法為：

用最大值運(yùn)算更新顯著結(jié)構(gòu)特征，如下式所示：

wn＝max(w1,w2,...,wn-1)；

用前n-1次迭代的隸屬度圖更新第n次迭代過程所需的隸屬度圖。

進(jìn)一步需要說明的是，當(dāng)滿足下式時，停止迭代過程，其中為n次迭代中的隸屬度圖的均值：

本發(fā)明的有益效果在于：本發(fā)明針對所處理的紅外圖像存在背景噪聲大、對比度低以及細(xì)節(jié)模糊的問題，為達(dá)到噪聲抑制和細(xì)節(jié)增強(qiáng)的目的，采用貝葉斯最佳估計模型，提出了適用刻畫信噪比低的紅外圖像特征的先驗?zāi)Ｐ?，并?jīng)過多尺度結(jié)構(gòu)特征進(jìn)一步的優(yōu)化，最后基于多尺度結(jié)構(gòu)特征分別進(jìn)行噪聲抑制和細(xì)節(jié)增強(qiáng)。本發(fā)明在處理低信噪比的紅外圖像時能消除背景噪聲、保持原始圖像自然性且銳化目標(biāo)的邊緣輪廓；在處理高信噪比的圖像時，可以自適應(yīng)地增強(qiáng)目標(biāo)的細(xì)節(jié)和紋理同時能避免強(qiáng)邊緣的過增強(qiáng)。

附圖說明

圖1為全局直方圖均衡和基于對比度限制的自適應(yīng)局部直方圖均衡增強(qiáng)紅外圖像的效果示意圖，其中圖1(a)原始紅外圖像,(b)全局直方圖均衡結(jié)果,(c)為基于對比度限制的自適應(yīng)局部直方圖結(jié)果；

圖2為全局直方圖均衡和基于對比度限制的自適應(yīng)局部直方圖均衡直接處理背景噪聲較大的紅外圖像的效果示意圖，其中圖2(a)有噪紅外圖像,(b)全局直方圖均衡放大噪聲,(c)對比度限制的自適應(yīng)局部直方圖放大噪聲；

圖3為基于圖像結(jié)構(gòu)信息相似性的nlm和bm3d濾波方法對紅外圖像的背景噪聲抑制效果示意圖，其中圖3(a)有噪紅外圖像,(b)wt抑制噪聲結(jié)果,(c)bm3d抑制噪聲結(jié)果,(d)tv抑制噪聲結(jié)果。

圖4為本發(fā)明的實施流程示意圖；

圖5為約束控制函數(shù)和隸屬度曲線圖；

圖6為清晰圖像及有噪圖像不同先驗擴(kuò)散圖，其中，圖6(a)是一幅清晰紅外圖像，圖6(b)是圖6(a)加入標(biāo)準(zhǔn)方差為25的零均值高斯噪聲的紅外圖像，圖6(c)為圖6(b)在gibbs先驗下的約束圖，圖6(d)為圖6(b)在本發(fā)明提出紅外圖像先驗下的約束圖。

圖7為迭代次數(shù)和多尺度特征圖均值關(guān)系示意圖；

圖8為自適應(yīng)增益函數(shù)和隸屬度曲線。

具體實施方式

以下將結(jié)合附圖對本發(fā)明作進(jìn)一步的描述，需要說明的是，本實施例以本技術(shù)方案為前提，給出了詳細(xì)的實施方式和具體的操作過程，但本發(fā)明的保護(hù)范圍并不限于本實施例。

本實施例以具備一定背景噪聲的紅外圖像作為實例，整個噪聲抑制和邊緣增強(qiáng)過程的技術(shù)方案包括：搭建適用刻畫紅外圖像特征的先驗?zāi)Ｐ?、基于多尺度特征的紅外噪聲抑制和基于多尺度特征的紅外細(xì)節(jié)增強(qiáng)等步驟。

如圖4所示，步驟一搭建適用刻畫紅外圖像特征的先驗?zāi)Ｐ停簻?zhǔn)確刻畫紅外圖像特征且易求解的先驗成為抑制紅外背景噪聲且能保護(hù)邊緣的關(guān)鍵。

本發(fā)明抑制噪聲所用的模型是貝葉斯準(zhǔn)則，該模型被廣泛的應(yīng)用于圖像處理計算機(jī)領(lǐng)域。一幅清晰紅外圖像u受到噪聲n污染：

v＝u+n(1)

根據(jù)上述背景噪聲的分析，n是零均值的高斯白噪聲。對于退化圖像v和原始圖像u，貝葉斯準(zhǔn)則表述為：

p(u|v)∝p(v|u)p(u)(2)

其中，p(v|u)為似然密度函數(shù)，p(u)為原始圖像的先驗密度函數(shù)，p(u|v)稱為后驗密度函數(shù)；根據(jù)map估計準(zhǔn)則，即在給定觀測數(shù)據(jù)v的條件下，使后驗密度函數(shù)p(u|v)最大的估計量為：

對式(3)取負(fù)對數(shù)可轉(zhuǎn)化為：

建立觀測模型求似然密度函數(shù)和建立紅外圖像先驗?zāi)Ｐ凸烙嬒闰灻芏群瘮?shù)成為預(yù)測最優(yōu)的估計量的關(guān)鍵。退化圖像和原始圖像之間的誤差滿足高斯白噪聲分布，因此似然密度函數(shù)p(v|u)可表示為:

這樣，原始圖像u的先驗?zāi)Ｐ蚿(u)的估計直接影響最終復(fù)原效果的好壞。

被噪聲污染的像素點(diǎn)在鄰域內(nèi)比較孤立，邊緣處的像素點(diǎn)與某一側(cè)鄰域相似，而平坦區(qū)的像素點(diǎn)則和鄰域內(nèi)其他像素點(diǎn)都具有較強(qiáng)相似性。對于噪聲污染的紅外圖像v＝{v(i)|i∈i}，像素點(diǎn)i,j∈p方塊內(nèi)的相似性取決于像素點(diǎn)之間的亮度和梯度差異，將亮度值差異的歐氏距離和梯度值差異的絕對距離之和的衰減函數(shù)定義為兩點(diǎn)之間的相似性r(i,j)：

式中h和ρ均為大于0的常數(shù)，定義分別為亮度值差異和梯度值差異的衰減系數(shù)。相似性0＜r(i,j)≤1，值越大，表示兩像素點(diǎn)間的相似程度越大。

在以i為中心像素點(diǎn)的鄰域內(nèi)，根據(jù)其他像素點(diǎn)與中心像素點(diǎn)之間相似性估計中心像素點(diǎn)在鄰域p內(nèi)的隸屬性，并以此刻畫紅外圖像中的形狀先驗。定義像素點(diǎn)i在領(lǐng)域內(nèi)的隸屬度w(i)為：

z(i)是規(guī)范化函數(shù)，σ是高斯模板的標(biāo)準(zhǔn)差，中心像素點(diǎn)的隸屬度0＜w(i)≤1，w(v)＝{w(i)丨i∈v}為圖像v的隸屬度運(yùn)算，構(gòu)成了一幅與原圖像尺寸一致的隸屬度圖，在本實施例中，取鄰域p尺寸為7。

為消除噪聲，對應(yīng)較小的隸屬度w(i)的像素點(diǎn)應(yīng)該具有較小的約束量，可以盡量平滑圖像的噪聲；在邊緣處，盡量迅速加大隸屬度適中的像素點(diǎn)的約束，保證邊緣在平滑過程中得以保持；在平坦區(qū)域內(nèi)，為了消除周圍噪聲像素點(diǎn)的殘余，應(yīng)該迅速減小對較大隸屬度w(i)像素點(diǎn)的約束。式(9)定義了基于隸屬度的約束控制函數(shù)，圖5描述了兩者之間關(guān)系。

c(w)＝c+(1-c)·sin(π·w)(9)

將gibbs先驗和擴(kuò)散控制函數(shù)結(jié)合的紅外圖像先驗約束為:

圖6(a)是一幅清晰紅外圖像，圖6(b)是圖6(a)加入標(biāo)準(zhǔn)方差為25的零均值高斯噪聲的紅外圖像，圖6(c)為圖6(b)在gibbs先驗下的約束圖，圖6(d)為圖6(b)在本發(fā)明提出紅外圖像先驗下的約束圖。與gibbs先驗相比，在本發(fā)明提出的先驗下的約束圖受噪聲影響較小，且基本能指示出正確的邊緣位置和強(qiáng)度。

基于map的貝葉斯估計轉(zhuǎn)化為:

β為正則化系數(shù)。根據(jù)原始圖像的噪聲水平，噪聲水平越高正則化系數(shù)取值越小^[12]。

步驟二基于多尺度特征的紅外噪聲抑制：采用分裂bregman方法求解模型(11)，同時利用多尺度特征不斷更新迭代貝葉斯過程，使得獲取最優(yōu)的濾波結(jié)果。

為利用分裂bregman方法來求解模型(11)，引入輔助量g，其中d是梯度算子。式(11)可轉(zhuǎn)換為有關(guān)g和u的優(yōu)化約束模型:

添加關(guān)于g和u的二次懲罰項，得到與式(14)對應(yīng)的無約束優(yōu)化模型:

利用收縮(shrinkage)技術(shù)，參數(shù)μ按2^k指數(shù)倍增加控制更新過程，得到關(guān)于g和u的更新公式，如式(14)和(15)所示：

m＝c(wn+1)d·un+1

其中wn＝w(un)，為n次迭代后圖像的隸屬度圖。從本質(zhì)上講，迭代的圖像也是一個鄰域濾波的過程，當(dāng)c(wn+1)＝1，轉(zhuǎn)化為gibbs的先驗，很明顯本實施例的約束先驗更加考慮了圖像原始結(jié)構(gòu)信息。

多尺度特征表達(dá)能夠在不同尺度上描述圖像局部結(jié)構(gòu)特征，多尺度可以獲取感興趣特征的最佳尺度，所以，使用多尺度平滑提取特征優(yōu)化約束先驗。改變貝葉斯平滑次數(shù)將會產(chǎn)生不同尺度的圖像，可以提取在不同尺度上圖像的顯著結(jié)構(gòu)特征，這些多尺度特征對特征的分析和應(yīng)用具有很大意義。定義原始圖像經(jīng)過n次貝葉斯平滑過程為:

[unwn]＝bays(un-1,wn-1)(16)

[u1w1]＝bays(u,w)(17)

用wn-1表示多尺度的隸屬度圖，由于重要的局部特征在每個尺度下都是最明顯的。所以，用最大值運(yùn)算更新顯著結(jié)構(gòu)特征，如式(18)所示：

wn＝max(w1,w2,...,wn-1)(18)

用前n-1次迭代的隸屬度圖更新第n次迭代過程所需的隸屬度圖。

由于正則化系數(shù)β是粗略估計，為了最佳保留邊緣和平滑噪聲，仍然需要精準(zhǔn)確定貝葉斯平滑過程。根據(jù)從圖7可以看出，圖6(b)隨著迭代次數(shù)的增加，越來越大的特征圖的均值給出了正確的變化趨勢。由于噪聲減少，而邊緣和紋理得到保持，平滑結(jié)果的結(jié)構(gòu)特征值會整體增加。采用相鄰兩幅平滑結(jié)果的隸屬度圖像來判斷圖像的平滑過程是否應(yīng)該停止。當(dāng)滿足式(21)時，停止平滑過程，其中為n次平滑結(jié)果特征圖的均值。

步驟三基于多尺度特征的紅外細(xì)節(jié)增強(qiáng)：根據(jù)不受噪聲干擾的多尺度特征圖能表征紅外圖像的結(jié)構(gòu)信息，本發(fā)明基于多尺度特征圖改進(jìn)傳統(tǒng)的非銳化掩膜增強(qiáng)方法來增強(qiáng)紅外圖像的邊緣。

um算法思想是先提取圖像的高頻分量加以放大，再疊加到原圖上，實現(xiàn)保持低頻信息基本不變，同時增強(qiáng)原圖像的邊緣和細(xì)節(jié)等高頻信息。根據(jù)前面分析可知，多尺度特征圖不受噪聲影響地指示邊緣和細(xì)節(jié)等形狀的位置和強(qiáng)度，而多尺度的結(jié)構(gòu)特征能提取更加詳細(xì)的結(jié)構(gòu)信息。所以，為有針對性地增強(qiáng)邊緣同時防止噪聲的放大，本發(fā)明將隸屬度圖和多尺度提取引入如式(22)所示的非銳化掩模算法，自適應(yīng)地控制增益系數(shù)：

其中，uen是最終輸出的增強(qiáng)圖像，un是貝葉斯模型的最優(yōu)估計，g(wn+j-1)是關(guān)于un+j-1的隸屬度圖wn+j-1的增益函數(shù)。自適應(yīng)的增益函數(shù)實現(xiàn)在不同區(qū)域進(jìn)行不同強(qiáng)度的增強(qiáng)。g(w)其分段表示見下：

分段點(diǎn)如圖8所示，以圖6(b)為例描述了結(jié)構(gòu)特征和自適應(yīng)增益函數(shù)的關(guān)系曲線。在其結(jié)構(gòu)特征值較小的區(qū)域，對應(yīng)圖像強(qiáng)邊緣，為了避免出現(xiàn)過沖現(xiàn)象，則進(jìn)行相對較弱的增強(qiáng)；隨著結(jié)構(gòu)特征值的增大，增強(qiáng)程度也應(yīng)緩慢提升；在結(jié)構(gòu)特征值較大的區(qū)域，對應(yīng)圖像平坦區(qū)域，由于人眼對噪聲較敏感，只做較小的增強(qiáng)，增益函數(shù)隨隸屬度迅速減小。

對于本領(lǐng)域的技術(shù)人員來說，可以根據(jù)以上的技術(shù)方案和構(gòu)思，作出各種相應(yīng)的改變和變形，而所有的這些改變和變形都應(yīng)該包括在本發(fā)明權(quán)利要求的保護(hù)范圍之內(nèi)。