丁帥帥，于 洋，馮迎賓

(沈陽理工大學(xué)自動化與電氣工程學(xué)院，沈陽 110159)

紅外相機通過捕獲電磁波進(jìn)行成像。 一般情況下，人體目標(biāo)的溫度會高于周圍環(huán)境，在紅外圖像上的顯示效果為人體目標(biāo)比周圍環(huán)境亮。 基于這種亮度差，提出一種基于箱形圖的自適應(yīng)非線性變換算法，可實現(xiàn)根據(jù)紅外圖像特點自動生成參數(shù)，并在保證圖像平滑過渡的前提下，增強紅外圖像[1]。

紅外圖像增強算法有很多，常用的有直方圖均衡化[2-3]、自適應(yīng)直方圖均衡化[4]、直方圖規(guī)定化[5]、對比度受限的自適應(yīng)直方圖均衡化[6]、灰度變換[7]以及空域濾波[8]等方法。 陳明亮等[9]針對灰度圖像分段線性變換[10]轉(zhuǎn)折點選取的問題，提出了在直方圖基礎(chǔ)上使用OTSU 法[11-12]和格拉布斯法[13]得到分段線性變換函數(shù)雙轉(zhuǎn)折點的方案，該方案可以自動選擇轉(zhuǎn)折點，確定變換曲線，但該方案僅適用于雙峰圖。 郭曉川[14]提出一種基于小波變換的紅外圖像增強算法，孫曉斐等[15]提出了基于特征融合的紅外圖像增強算法，陳龍等[16]提出一種改進(jìn)自適應(yīng)分段線性變換算法，通過尋找灰度集中區(qū)域最大灰度級，比較灰度直方圖簇的個數(shù)，在候選集中自動選取適合于當(dāng)前原始圖像的參數(shù)取值，但這種方法在灰度分布不均勻時會造成算法計算失誤。 郭鈺璐[17]提出一種融合邊緣信息的對比度增強算法，該方法將分層處理和增強濾波相結(jié)合實現(xiàn)紅外圖像的對比度增強。

本文在分段線性變換方法的基礎(chǔ)上，采用箱形圖思想，通過排序統(tǒng)計的方法，實現(xiàn)自動確定灰度集中分布區(qū)間的功能，一定程度上解決因圖片灰度分布不均造成的灰度分布區(qū)間選取錯誤問題，增強算法穩(wěn)定性，再結(jié)合引入分段點的三次樣條插值的非線性灰度變換，達(dá)到自適應(yīng)增強圖像的效果。

1 自適應(yīng)分段非線性變換算法

1.1 算法原理

本文采用一種基于箱形圖原理的自適應(yīng)分段非線性變換算法，該算法可以根據(jù)輸入紅外圖像的灰度特點，自動計算出紅外圖像集中分布區(qū)間，實現(xiàn)動態(tài)調(diào)整分段點，再經(jīng)過一種非線性灰度變換，達(dá)到圖像增強效果。 本文算法整體流程圖如圖1所示。

圖1 本文算法流程圖

使用紅外攝像頭采集8 位、256 灰度級的圖像。 首先采用箱形圖原理，將輸入的紅外圖像扁平化為一維數(shù)組形式，然后按照從小到大的順序排序分位，找出圖像灰度集中的灰度區(qū)間，獲得集中分段點，從灰度區(qū)間右側(cè)到灰度255 之間另取一點作為調(diào)整分段點，用以改變非線性曲線形狀，再加上灰度級0 和灰度級255 共5 個分段點，將這些分段點進(jìn)行三次樣條插值，通過改變調(diào)整分段點，可以生成不同的變換曲線。 對圖像像素進(jìn)行灰度變換，遍歷整個圖像得到的即是處理好的圖像。

1.2 求取自適應(yīng)分段點

能否準(zhǔn)確獲取灰度集中區(qū)域分段點關(guān)系到算法處理效果的好壞，過窄或者過寬的范圍都會影響最終效果。 通常情況下，紅外圖像大部分像素的灰度值集中在某一范圍內(nèi)，導(dǎo)致紅外圖像整體偏暗，因此找到該范圍區(qū)間成為解決圖像整體對比不明顯的關(guān)鍵所在。 一些自適應(yīng)求取灰度集中分布方法以灰度直方圖中灰度個數(shù)最大的灰度值作為灰度集中分布的中心，向兩邊搜索。 該方法在灰度圖中出現(xiàn)局部峰值(例如極白點處和極黑點處)時，可能會出現(xiàn)局部窄范圍峰值超過真正灰度分布區(qū)域峰值的情況，導(dǎo)致程序判斷錯誤。 本文獲取灰度集中分布區(qū)間的方法基于箱形圖，通過對圖像灰度值排序分位來判斷灰度分布情況，避免了因出現(xiàn)局部峰值導(dǎo)致判斷錯誤的情況，增強了算法穩(wěn)定性。 具體步驟如下。

(1)將紅外圖像數(shù)組扁平化為一維數(shù)組，按照從小到大的順序排列，記總長度為n；

(2)取排序后數(shù)組的第25%的元素值，記為Q1，則

(3)取排序后數(shù)組的第50%的元素值，即數(shù)組的中位數(shù)，記為Q2，則

(4)取排序后數(shù)組的第75%的元素值，記為Q3，則

(5)取Q1 與Q3 的差值為IQR，即

(6)記分段區(qū)間上限為MAX，則

(7)記分段區(qū)間下限為MIN，則

通過以上步驟，當(dāng)輸入紅外圖像時，能夠自動獲取灰度集中區(qū)域的分段點，輸入不同紅外圖像時，會針對不同的圖像得到不同的分段點，達(dá)到自適應(yīng)計算集中分段點的效果。

由于紅外設(shè)備參數(shù)的不同以及環(huán)境的復(fù)雜多樣，為使改進(jìn)算法有更好的適應(yīng)性，針對公式(5)和公式(6)，引入調(diào)整參數(shù)θ，取值范圍設(shè)為( -0.5，0.5)，則

θ值默認(rèn)為0，在圖像處理過程中用以對選定的范圍大小進(jìn)行補償，針對不同硬件設(shè)備屬性以及要達(dá)到的目標(biāo)效果，可以通過調(diào)整θ值實現(xiàn)。

1.3 非線性變換

傳統(tǒng)分段線性變換方式的優(yōu)點是算法原理簡單，比較容易實現(xiàn)，線性變換函數(shù)示意圖如圖2a 所示。 由于線性變換算法簡單，在復(fù)雜的圖像處理中，僅僅采用線性變換會導(dǎo)致圖像過渡非常生硬，尤其在灰度變化較大和集中分布區(qū)間端點附近的區(qū)域，會導(dǎo)致圖像失真。 變換結(jié)果對比如圖3所示，圖3a 經(jīng)線性變換獲得結(jié)果如圖3b 所示。 由圖3b 可見，圖像感興趣區(qū)域亮度得到了提升，背景區(qū)域灰度值得到了抑制，但是圖像邊緣過渡較為生硬，甚至出現(xiàn)了不連貫現(xiàn)象。 針對傳統(tǒng)線性變換的缺點，本文基于三次樣條插值方法，通過引入調(diào)整點參數(shù)，生成符合變換要求的三次樣條函數(shù)[18-20]進(jìn)行灰度變換，示意圖如圖2b 所示，變換結(jié)果見圖3c。

圖2 灰度變換函數(shù)示意圖

圖3 變換結(jié)果對比

灰度變換函數(shù)橫坐標(biāo)為輸入圖像灰度，縱坐標(biāo)為輸出圖像灰度，在未經(jīng)算法處理時，輸出圖像灰度等于輸入圖像灰度。 設(shè)灰度圖像灰度集中分布區(qū)間為(a，b)，則在灰度變換函數(shù)上集中分段點可設(shè)為(a，a)和(b，b)兩點，變換函數(shù)首尾兩端點為(0，0)和(255，255)，如果直接對這些分段點進(jìn)行三次樣條插值，則生成的變換曲線為一條斜率為1 的直線，圖像變換前后無變化，如圖4所示。

圖4 未引入調(diào)整分段點的灰度變換函數(shù)

為解決三次樣條插值后得出的變換函數(shù)無效果的問題，引入一調(diào)整分段點(c，d)，實現(xiàn)灰度變換函數(shù)形狀的改變，如圖5所示。

圖5 引入調(diào)整分段點的灰度變換函數(shù)

點(c，d)位置的選取直接影響到曲線形狀，因此必須找到合適的位置引入調(diào)整分段點。 圖5中(a，b)區(qū)間為灰度圖集中分布區(qū)域，也是紅外圖像中背景色集中部分，因此可通過將(a，b)區(qū)間曲線下拉來達(dá)到抑制背景效果；(b，255)區(qū)間是感興趣區(qū)域分布點，可通過將(b，255)區(qū)間曲線上拉來突出目標(biāo)；(0，a)區(qū)間是暗色分布區(qū)，不做單獨處理。 因此，滿足以上要求的方法為在(a，b)區(qū)間做下拉點和(b，255)區(qū)間做上拉點。 由于(a，b)區(qū)間位于灰度級中間，與圖像暗色分布區(qū)相鄰，直接在該區(qū)間引入調(diào)整點下拉處理，對暗色區(qū)域影響較大，會導(dǎo)致圖像暗色過度拉亮，不符合增強效果，因此只能在(b，255)區(qū)間引入調(diào)整分段點(c，d)。

通過多次實驗，確定最佳c、d取值為

求出調(diào)整分段點(c，d)后，對(0，0)、(a，a)、(b，b)、(c，d)、(255，255)五點進(jìn)行三次樣條插值，即可獲得適合該圖片的非線性變換曲線，最后進(jìn)行灰度變換，獲得處理后結(jié)果。 對于不同的灰度圖像，算法自動計算出不同的(a，b)區(qū)間，同時根據(jù)不同的(a，b)區(qū)間引入不同的調(diào)整分段點，生成不同的變換曲線，實現(xiàn)各階段自適應(yīng)調(diào)整。

2 實驗結(jié)果

本文的實驗結(jié)果包括兩個方面，一是自適應(yīng)分段點的精確度，二是采用三次樣條曲線實現(xiàn)圖像灰度變換的效果。 實驗采用平臺為Inteli7 11800H、NVIDIAGeForce 3050Ti 的筆記本， 在Miniconda 環(huán)境下使用Python 語言編程實現(xiàn)。

2.1 求取自適應(yīng)分段點

由于本文研究背景為車載行人檢測，所以實驗選取的素材為紅外攝像頭在室外道路上拍攝的行人圖像。 圖6為四組行人紅外圖像，圖7為圖6中四組圖像的灰度直方圖。

圖6 行人紅外圖像

直方圖中峰形區(qū)域是圖像中灰度級聚集的地方，該區(qū)域所對應(yīng)的灰度范圍也是本文所求的灰度集中區(qū)域的分段點取值范圍。 從圖7直方圖中可直接得出圖6中各個紅外圖像的峰值范圍，如表1所示。

表1 直方圖峰值范圍

本文算法所得到的灰度集中區(qū)域分段區(qū)間如表2所示。

表2 自適應(yīng)算法分段區(qū)間

由表1和表2對比可知，本文所采用的基于箱形圖的自適應(yīng)分段非線性變換算法得出的結(jié)果與實際結(jié)果相差不大，且實際的峰值范圍并不是確定范圍，有一定的可擴(kuò)展性，因此采用該算法求取圖像的灰度集中區(qū)域分段區(qū)間是可行的。

2.2 非線性變換

在獲取圖像灰度集中分布區(qū)域后，即可引入調(diào)整分段點，將包含調(diào)整分段點在內(nèi)的五個點進(jìn)行三次樣條插值，生成對應(yīng)的變換曲線。 變換的曲線在各分段點過渡比較光滑，變換效果如圖8所示。

圖8 灰度變換前后對比

圖8a、8b、8c、8d 各組左邊是原圖，右邊是變換后圖像，通過對變換前后的圖像對比，可看出變換后圖像目標(biāo)主體更加突出，環(huán)境背景得到抑制。圖9a、9b、9c、9d 分別是圖8a、8b、8c、8d 采用的變換曲線，頂端為直線是因為灰度級最大為255，將圖像亮的地方灰度值拉高可以更加突出目標(biāo)。 從圖9中可以看出，對于不同的紅外圖像，算法生成的自適應(yīng)集中分段點不同，引入的調(diào)整分段點也不同，在目標(biāo)分段區(qū)間調(diào)整點上拉后，背景區(qū)間自動下拉，實現(xiàn)了通過一個調(diào)整點控制多個分段的效果，而且過渡區(qū)域比較平滑，為傳統(tǒng)分段變換算法所不具有的功能。

圖9 灰度變換曲線

2.3 調(diào)整分段點對變換結(jié)果的影響

本文采用五點三次樣條插值方法生成灰度變換曲線，五點分別是灰度變換函數(shù)的首尾兩點(0，0)和(255，255)，通過箱形圖算法求取的自適應(yīng)分段點(a，a)和(b，b)以及引入的調(diào)整點(c，d)。針對圖8b 左側(cè)原圖，改變調(diào)整點的位置，其處理效果如圖10 所示。

圖10 調(diào)整分段點對變換結(jié)果的影響

由圖10 可以看出，不同的分段點選取會使得變換后的結(jié)果差別顯著：圖10a增強過度，背景被完全抑制，近似二值化處理，丟失信息過多；圖10b 目標(biāo)融入背景，無法突出目標(biāo)；圖10c 是根據(jù)本文算法生成調(diào)整點處理后的結(jié)果，在適當(dāng)保留背景信息的同時有效地突出目標(biāo)，取得了更加理想的效果。

2.4 算法對比

為了比較不同變換算法處理后的效果，本文使用同一張圖像，分別采用直方圖均衡化、自適應(yīng)直方圖均衡化、直方圖規(guī)定化、對比度受限的自適應(yīng)直方圖均衡化和本文算法，將5 種處理效果進(jìn)行對比，結(jié)果如圖11 所示。

直方圖均衡化用來增加全局對比度；自適應(yīng)直方圖均衡化用來改進(jìn)圖像的局部對比度；直方圖規(guī)定化產(chǎn)生處理后有特殊直方圖的圖像；對比度受限的自適應(yīng)直方圖均衡化能限制自適應(yīng)均衡化中噪聲的問題。 從圖11 中可以看出，圖11b、11c 和11e 對圖像背景抑制效果不好，圖11d 雖然背景抑制效果較好，但也增強了部分干擾信息。 綜合來看，在此環(huán)境中本文算法效果較佳。

圖11 算法對比

傳統(tǒng)的圖像質(zhì)量評價方法，例如峰值、信噪比等，僅僅是對像素點之間的純數(shù)學(xué)統(tǒng)計，沒有考慮到人眼視覺系統(tǒng)感知特性。 因此，本文采用更能貼近人主觀感受的結(jié)構(gòu)相似性圖像評價指標(biāo)SSIM 對圖像質(zhì)量進(jìn)行客觀評價。 針對圖11 中不同算法的評價結(jié)果如表3所示。

表3 SSIM 評價結(jié)果

SSIM 評價指標(biāo)結(jié)果范圍為0 至1，指標(biāo)結(jié)果越接近1，說明算法處理效果越好。 由表3可知，在五種圖像處理算法中，本文所采用的算法處理效果最好。

3 結(jié)論

針對紅外圖像成像效果不佳、目標(biāo)主體突出不明顯的問題，本文在研究分段線性變換算法的基礎(chǔ)上，提出一種基于箱形圖的自適應(yīng)分段非線性變換算法。 該算法能夠根據(jù)紅外圖像的不同，自動計算出不同圖像的灰度集中分布區(qū)間，無需人為輸入?yún)?shù)。 算法采用基于箱形圖的統(tǒng)計學(xué)思想，通過排序分位，避免了因出現(xiàn)局部少量灰度峰值使算法失效的問題，算法穩(wěn)定性更強。 在灰度非線性變換方面，引入調(diào)整分段點，生成不同分段變換曲線，從而保證了曲線過渡平滑，使經(jīng)過變換后的圖像過渡更自然。 通過實驗證明，本文提出的算法能夠針對不同圖像生成合適的自適應(yīng)參數(shù)，解決傳統(tǒng)算法無法自適應(yīng)的問題。 在變換結(jié)果方面，能夠提高圖像前景目標(biāo)的亮度，同時能抑制背景色，圖像過渡自然。 通過與其他幾種增強算法的對比，突顯了本文算法的優(yōu)勢，證明了本文算法取得的處理效果更好。