楊緒兵，覃欣怡，張福全

(南京林業(yè)大學(xué) 信息科學(xué)技術(shù)學(xué)院， 南京 210037)

基于樣條的林火圖像多閾值分割算法

楊緒兵*，覃欣怡，張福全

(南京林業(yè)大學(xué) 信息科學(xué)技術(shù)學(xué)院， 南京 210037)

基于光滑樣條原理，提出一種自適應(yīng)的多閾值分割算法HistSplineReg，即采用光滑樣條回歸圖像一維直方圖，再從回歸函數(shù)尋找極值從而實現(xiàn)圖像的多閾值自動分割。較之現(xiàn)有的閾值分割方法，HistSplineReg具有以下優(yōu)勢：1)設(shè)計方法符合人類直覺；2)基于光滑樣條設(shè)計算法，有理論依據(jù)；3)閾值個數(shù)及閾值位置可自動選擇；4)回歸函數(shù)可分析求解，計算規(guī)模主要集中在矩陣的Cholesky分解，矩陣大小由圖像像素水平級決定，而不是圖像尺寸; 5)只有一個待定參數(shù)，該參數(shù)用于平衡回歸經(jīng)驗誤差和回歸函數(shù)的光滑性。對林火識別問題，實驗提供一個經(jīng)驗參數(shù)供參考。最后，在紅綠藍(lán)顏色(RGB)模式的林火數(shù)字圖像上進(jìn)行實驗，從灰度圖像、多種顏色通道、各通道分割結(jié)果合成的彩色圖像等方面進(jìn)行驗證，與同樣采樣回歸思想的支持向量回歸(SVR)及多項式回歸(PolyFit)相比，HistSplineReg方法直觀分割效果更好，且三種方法都反映出紅色通道信息對林火圖像分割效果的影響更為顯著。

圖像分割；光滑樣條函數(shù)；林火識別問題；閾值

0 引言

森林火災(zāi)造成的損失巨大，為有效保護(hù)人們的生命財產(chǎn)安全和保護(hù)森林資源，世界各國不惜花費巨資，投入幾十億、上百億美元，用于研制森林防火預(yù)警系統(tǒng)。近年來已有很多火災(zāi)監(jiān)測系統(tǒng)，按研究歷史上可分為兩個階段：一是以硬件設(shè)備為主導(dǎo)的，如基于傳感器(如感溫、感煙報警器等)檢測煙或火情的預(yù)警方法[1-2]，此類方法主要弊病在于：火源與傳感器部署地之間總會存在距離，加之野外環(huán)境，極易造成誤報、漏報、延時等問題[3]。二是以軟件技術(shù)為主導(dǎo)，特別是計算機(jī)視覺技術(shù)的興起，為防火系統(tǒng)注入了新的生機(jī)[4]，基于數(shù)字圖像或視頻的可視化系統(tǒng)逐漸占據(jù)主導(dǎo)地位，該系統(tǒng)主要利用圖像的顏色、紋理、幾何特征、閃爍(flickering)、目標(biāo)運動軌跡等[5]?？v觀林火識別研究歷史，單純就靜態(tài)圖像而言(對應(yīng)稱視頻為動態(tài)圖像)，據(jù)作者所知，目前所能查到的文獻(xiàn)中，火焰識別的研究工作仍主要關(guān)注圖像的顏色信息[6]。一個可能的原因是圖像中目標(biāo)“火”的特殊性，如火焰的形狀、顏色、火勢走向等都是不確定的，使得已有的、性能卓越的目標(biāo)檢測方法失效，或至少無法直接應(yīng)用于林火識別任務(wù)。但對林火圖像，也存在如下先驗知識：如起火初期溫度低，顏色主要表現(xiàn)為紅色或黃色，隨著溫度升高，顏色逐漸變?yōu)榈S或白色；白天的火焰飽和度(saturation)高于晚上等。因此，基于顏色的火焰識別較之其他方法更有優(yōu)勢。常用的顏色空間有RGB(Red, Green, Blue)、HSI(Hue Saturation Intensity)[7]、YCbCr[8]等。

圖像分割的目的是將圖像的劃分成若干個區(qū)域，區(qū)域之間彼此不相交。該工作屬圖像預(yù)處理步驟，但分割質(zhì)量會直接影響后續(xù)的識別效果。目前，常用的圖像分割方法有：閾值分割、邊緣檢測、區(qū)域生長和基于其他理論的分割方法(如形態(tài)學(xué)方法)等[9]。本文僅關(guān)注圖像閾值分割方法。閾值分割方法很多，如基于聚類的k均值、FCM(Fuzzy C-means)、層次聚類[10-11]等方法，僅利用灰度直方圖求解閾值的方法就有Otsu、最大熵法[12]等。已有的研究成果告知：基于Otsu的分割方法對于目標(biāo)和背景所占面積接近效果較好，否則易失效，同時也存在分割不夠充分問題。最大熵法通過最大化目標(biāo)熵和背景熵之和來決定閾值，圖像的細(xì)節(jié)越豐富則分割效果越好，而且計算亦較為耗時[13]。這兩種方法的原型是基于單個閾值提出的，單閾值方法(亦稱雙層分割(bilevel segmentation))雖易于求解，但難以適用復(fù)雜場景下圖像分割。目前也有若干閾值的改良版本，限于篇幅在此不細(xì)述。但對多閾值分割問題，一個關(guān)鍵問題就是如何確定閾值個數(shù)，一個行之有效地辦法就是通過曲線擬合直方圖，再從擬合曲線獲得閾值個數(shù)。如文獻(xiàn)[14]采用支持向量回歸(Support Vector Regression, SVR)擬合直方圖，并通過邊界支持向量(Boundary Support Vector, BSV)來確定閾值，該方法雖然回避了求解極值，閾值選擇空間限制在邊界支持向量集中，但卻面臨著需要核及核參數(shù)的選擇問題；文獻(xiàn)[15-16]采用高斯過濾器、混合高斯函數(shù)回歸直方圖，但也只是把閾值選擇問題轉(zhuǎn)化為多個高斯函數(shù)的混合問題，而且同樣面臨多個參數(shù)的選擇問題；還有采用高階多項式實現(xiàn)擬線性回歸[17]，存在擬合經(jīng)驗誤差大、多項式階數(shù)選擇問題等，無法滿足自動化的要求。

考慮林火圖像顏色信息的重要性，如前文所述，本文采用光滑樣條進(jìn)行模型設(shè)計。本質(zhì)上說，光滑樣條采用的局部回歸的方法，每個子區(qū)間內(nèi)部通過多項式擬合，子區(qū)間端點處的光滑性由其左右高階導(dǎo)數(shù)相等來保證，有望解決回歸函數(shù)的經(jīng)驗誤差大、回歸曲線振蕩等問題。本文提出一種基于光滑樣條的回歸方法——HistSplineReg擬合直方圖，有以下幾點優(yōu)勢：1)方法簡單，易于理解，符合直方圖中尋找波谷(valley)的直覺；2)基于光滑樣條設(shè)計回歸算法，回歸函數(shù)的存在性、唯一性等有理論保證；3)閾值個數(shù)及閾值位置由回歸函數(shù)的極值確定，可實現(xiàn)自動化選擇；4)無需迭代，回歸函數(shù)可分析求解，且計算規(guī)模主要集中在矩陣Cholesky分解，矩陣大小由圖像的大小無關(guān)，僅取決于像素水平級；5)只有一個待定參數(shù)，該參數(shù)用于折中回歸誤差和函數(shù)的光滑性，而且本文的實驗部分對林火識別問題給出一個經(jīng)驗參數(shù)供讀者參考。

1 背景知識

本章將簡要介紹一下光滑樣條產(chǎn)生機(jī)理和直方圖的統(tǒng)計特性。

1.1 光滑樣條

設(shè)觀測集為{(xi,yi)}(i=1,2,…)，按式(1)所示的優(yōu)化目標(biāo)求解回歸函數(shù)g(·)。若g為光滑樣條函數(shù)[18]，需滿足以下優(yōu)化目標(biāo)：

(1)

其中：wi為權(quán)參數(shù)，g″為回歸函數(shù)g的二階導(dǎo)數(shù)。光滑參數(shù)η∈[0,1]，當(dāng)η→1時，模型關(guān)注于第一項的經(jīng)驗誤差部分，使得回歸曲線穿過盡可能多的觀測點；當(dāng)η→0時，模型關(guān)注第二項，顯然g趨近于直線。為照顧下文即將描述的一維直方圖擬合問題，本節(jié)簡單介紹三次樣條(Cubic Spline)，所闡述的問題均在一維空間上。

在有界閉區(qū)間[a,b]上，存在該區(qū)間上的一個劃分，即存在實數(shù)序列t1,t2,…,tn，使得a

1)g為區(qū)間(a,t1)，(t1,t2)，…，(tn,b)上的三次多項式(此處的開區(qū)間也可用閉區(qū)間)；

2)g及其一階二階導(dǎo)數(shù)g′、g″在每個ti處連續(xù)，i=1,2,…,n。

若在條件1)、2)，再加上條件區(qū)間端點的二階、三階導(dǎo)數(shù)為0，即:

g″(a)=g″(b)=g?(a)=g?(b)=0

(2)

此時稱g為自然三次樣條函數(shù)(Natural Cubic Spline, NCS)。由式(2)知g在小區(qū)間[a,t1]、[tn,b]上為線性函數(shù)。

本質(zhì)上說，樣條回歸采用局部擬合思想，即先在每個小區(qū)間內(nèi)部作光滑曲線回歸，區(qū)間端點處加入如連續(xù)、可微等條件限制，使回歸曲線整體上滿足連續(xù)、光滑的要求。

1.2 直方圖的統(tǒng)計特性

接下來將要考慮采用光滑樣條回歸直方圖問題。

2 本文方法

對于直方圖而言，設(shè)待回歸的樣本集為{(i,p(i)}(i=0,1,…,L-1)。它具有兩個顯著特點：1)樣本數(shù)固定(由像素水平級決定)，與圖像尺寸無關(guān)；2)在[0,L-1]區(qū)間上劃分小區(qū)間盡可能不破壞直方圖的統(tǒng)計特性，可按水平級進(jìn)行劃分，小區(qū)間的長度可取固定值h，對直方圖，不妨取h=1。如此不僅降低了計算量，還兼顧了直方圖的實際意義。

2.1 直方圖光滑樣條回歸函數(shù)性質(zhì)

先給出光滑樣條在直方圖回歸問題中的幾個結(jié)論，以定理形式描述[19]。

記回歸函數(shù)為g(·)，灰度級i處的估計值g(i)=gi，g″(i)=γi。記回歸參數(shù)向量g=(g1,g2,…,gL-2)T，γ=(γ1,γ2,…,γL-2)T，由式(2)知，γ0=γL-1=0。在小區(qū)間[i,i+1]上，由前假設(shè)知:

(3)

在直方圖回歸問題中，主要結(jié)論如下。

引理1 在小區(qū)間[i,i+1]上,自然三次樣條函數(shù)g(·)滿足式(3)，則有

g″(t)=(t-i)γi+1+(i+1-t)γi;t∈[i,i+1]

(4)

成立。

證明 因三次樣條的二階導(dǎo)函數(shù)為線性函數(shù)，在[i,i+1]上，如圖1所示，(t,g″(t))為線性函數(shù)任一點，由式(3)及三點共線關(guān)系可得：

(5)

對式(4)化簡整理可立得結(jié)論。

圖1 自然三次樣條函數(shù)二階導(dǎo)函數(shù)構(gòu)造示意圖Fig. 1 Illustration for 2-order deviation function of natural cubic spline

引理2 在小區(qū)間[i,i+1](i=1,2,…,L-3)的自然三次樣條回歸函數(shù)，具有如下表達(dá)式：

t)[(1+t-i)γi+1+(2+i-t)γi];t∈[i,i+1]

(6)

對式(4)進(jìn)行關(guān)于變量t兩次積分，積分過程中產(chǎn)生兩個任意常數(shù)，將式(3)的條件代入即可消去。整理可得式(6)結(jié)果，在i=1,2,…,L-3時均成立，共產(chǎn)生L-3個未知數(shù)γi。證明略。

以上過程不包括兩個端點小區(qū)間[0, 1]和[L-2,L-1]。區(qū)間[0, 1]和[L-2,L-1]上，按三次樣條函數(shù)定義，端點處的二階導(dǎo)數(shù)為0，同樣是按線性函數(shù)處理。可得：

g(t)=

(7)

至此，樣條函數(shù)g(·)在區(qū)間[0,L-1]的表達(dá)式已全部給出，接下來就是如何求解參數(shù)向量g和γ了。

由式(6)、(7)，再加三次樣條的條件2)的一階、二階導(dǎo)數(shù)連續(xù)，可得如下定理：

定理1 設(shè)向量g、γ可確定光滑樣條g(·)當(dāng)且僅當(dāng)QTg=Rγ時成立，此時式(1)的第二項滿足下式：

其中:Q、R為帶狀常數(shù)矩陣，

且R對稱、正定(嚴(yán)格對角占優(yōu)矩陣)，K=QR-1QT。

因

γTQTg=γTRγ

證畢。

令模型式(1)的權(quán)值為1，有:

(8)

在式(8)目標(biāo)下考慮光滑樣條問題，懲罰因子α=1/η-1(>0)。問題求解用定理2描述。

定理2 式(8)的解存在且唯一，解為:

g=(I+αK)-1Y

(9)

其中Y=[y1,y2,…，yL-2]T，g、K定義同前，I是單位矩陣。

證明 用矩陣形式重寫目標(biāo)函數(shù)，并記為S(g)，由定理1知:

S(g)=(Y-g)T(Y-g)+αgTKg=

gT(I+αK)g-2YTg+YTY

(10)

對直方圖而言,有yi=p(i-1)(i=1,2,…,L-2)成立。關(guān)于g求導(dǎo)，得:

(I+αK)g=Y

(11)

矩陣I+αK是方陣且正定，展開式(11)，可得g=(I+αK)-1Y。

證畢。

以上結(jié)論均在h=1的條件下完成。為兼顧直方圖的幾何意義，建議h的取值能夠正好等分區(qū)間[0,L-1]。

至此，直方圖樣條回歸函數(shù)求解完畢，以算法形式總結(jié)如下。

2.2 直方圖樣條回歸函數(shù)求解算法HistSplineReg

2.1節(jié)介紹了光滑樣條回歸一維直方圖過程，并在自然三次樣條函數(shù)基礎(chǔ)上，提出了可用于回歸直方圖的光滑樣條函數(shù)求解方法HistSplineReg。為回避矩陣求逆，實際應(yīng)用中可采用Cholesky分解來進(jìn)一步減小計算量[20]，簡述如下。

將定理1中的QTg=Rγ代入式(11)并整理，有：

g=Y-αKg=Y-αQR-1QTg=Y-αQγ

用QT左乘g=Y-αQγ，可得:

QTY=(R+αQTQ)γ

(12)

顯然矩陣R+αQTQ對稱且嚴(yán)格正定，可對其進(jìn)行Cholesky分解，有R+αQTQ=LDLT，其中D是嚴(yán)格正定對角陣，L是下三角陣。式(12)可變形為：

QTY=(LDLT)γ

(13)

以上求解過程歸納為HistSplineReg算法描述如下。

算法1 HistSplineReg算法。

輸入 樣本集{(i,p(i))}(i=0,1,…,L-1)；常量矩陣Q,R(定理1)，正則化參數(shù)α或η。

輸出 回歸估計值g；

步驟1 計算Y，yi=p(i-1),i=1,…,L-2；

步驟2 計算QTY；

步驟3 計算R+αQTQ的Cholesky分解因子L和D；

步驟4 求解方程QTY=(LDLT)γ，獲得γ；

步驟5 計算g=Y-αQγ。

如需繪制回歸曲線，可將HistSplineReg計算得到的g和γ，代入式(6)～(7)中，即可獲得區(qū)間[0,L-1]上回歸函數(shù)g(t)。根據(jù)直方圖回歸算法，就可以實現(xiàn)多閾值的圖像分割。

2.3 基于HistSplineReg算法的圖像多閾值分割算法

因HistSplineReg算法返回的只是像素水平級的估計值，而不是回歸曲線，如果根據(jù)式(6)～(7)來求解極值，會導(dǎo)致分割速度慢。本節(jié)中求解閾值可通過估算gi的一階導(dǎo)數(shù)，即若(gi-gi-1)/h<0且(gi+1-gi)/h>0，則該水平級i為閾值，如算法2。

算法2 基于HistSplineReg的多閾值分割算法。

輸入 待分割圖像Im;

輸出 分割后的圖像F。

步驟1 獲得圖像直方圖，并以{(i,p(i))}(i=0,1,…,L-1)形式存儲；

步驟2 執(zhí)行HistSplineReg算法，獲得回歸估計向量g；

步驟3 按水平級順序檢查g中分量gi的一階導(dǎo)數(shù)，若導(dǎo)數(shù)符號從負(fù)變正，則該水平級為對應(yīng)回歸曲線的極小值點，其個數(shù)即為閾值個數(shù)；

步驟4 根據(jù)獲得的多個閾值進(jìn)行圖像分割。

3 實驗驗證

本章中，將通過實驗來驗證上述結(jié)論。實驗將分為兩個部分：第一部分是直方圖擬合實驗；第二部分林火圖像上的分割實驗。由于該類問題缺乏標(biāo)準(zhǔn)數(shù)據(jù)庫，本文所用圖像采集自互聯(lián)網(wǎng)。實驗對比對象采用SVR[14]和多項式回歸方法(PolyFit)[21]。SVR選擇高斯核k(a,b)=exp((‖a-b‖2)/σ2)，正則化參數(shù)C和核參數(shù)σ選擇范圍為{10i|i=-4,-3,…,4}和{5×10i|i=-4,-3,…,4}；HistSplineReg只有一個正則化參數(shù)η，參數(shù)選擇范圍為(0,1)；多項式回歸需先指定多項式的最高階數(shù)。以下實驗環(huán)境是Window 7 SP1 64位操作系統(tǒng)，Inter Core i7-3632QM CPU, 2.2 GHz, 8 GB內(nèi)存, Matlab R2015b。SVR使用軟件包LIBSVM[22](序列最小化(Sequential Minimal Optimization, SMO)算法)。多項式回歸采用Matlab的Curve Fitting工具PolyFit函數(shù)。

3.1 直方圖擬合實驗

為方便可視化，避免直方圖p(i)(∈[0,1])與水平級i(i∈[0,255],i取整數(shù))之間的數(shù)值差距過大造成計算誤差，影響可視化效果，本節(jié)實驗中先進(jìn)行歸一化處理。如圖2所示。

圖2 直方圖及其擬合效果圖Fig. 2 Image histogram and its curve fitting results

圖2(b)中展示的SVR不敏感因子ε=0.005，高斯核參σ=50，平衡因子C=100。訓(xùn)練所得支持向量為167個，其中邊界支持向量BSV是142個。光滑樣條參數(shù)η=0.999 8。多項式回歸最高階為18。SVR、HistSplineReg和PolyFit的執(zhí)行時間分別為2.96 s，0.75 s和0.18 s(CPU時間，其中參數(shù)選擇時間未計算在內(nèi))。圖2可見，HistSplineReg方法的回歸效果更好。

3.2 林火圖像分割實驗

如何評價圖像閾值分割算法性能，目前國際上也一直沒有統(tǒng)一標(biāo)準(zhǔn)，現(xiàn)有文獻(xiàn)多是采用主觀視覺效果和算法的執(zhí)行效率作評估[12]。本文實驗對象是RGB林火圖像，為探究顏色信息對圖像分割效果的影響，分別從灰度圖像、RGB三通道亮度圖像及各通道分割完畢再次合成的彩色圖像來展示實驗效果。

圖4～5分別是SVR和多項式回歸的圖像分割效果表，這兩種方法對灰度圖像效果極差，圖像中除顯示幾個亮點外，看不到任何分割效果；SVR在綠色、藍(lán)色通道上分割區(qū)域較為清楚，但紅色通道效果不佳，因而造成了合成分割圖像中出現(xiàn)了詭異的黃紅色火焰，這是因為SVR在調(diào)參過程中，因參數(shù)空間的步長較大，未能選到更合適的參數(shù)，還有諸如核函數(shù)的選擇問題等。多項式擬合中，尤其是在直方圖尾部容易產(chǎn)生振蕩，階數(shù)越高，振蕩現(xiàn)象越嚴(yán)重。兩者對于林火的灰度圖像，可能因為在火焰的照明下，亮度值不同于其他類型圖像，如果閾值選擇不當(dāng)，如所選閾值比較接近，分割區(qū)域之間像素值差別不大，就會造成視覺感官上的全黑假象。

圖3 林火圖像HistSplineReg分割效果圖Fig. 3 Illustration for forest fire image segmentation results by HistSplineReg

圖4 林火圖像SVR分割效果圖Fig. 4 Illustration for forest fire image segmentation results by SVR

圖5 林火圖像多項式分割效果圖Fig. 5 Illustration for forest fire image segmentation results by Polynomial fitting

4 結(jié)語

本文提出一種基于光滑樣條的直方圖回歸方法HistSplineReg，并在此基礎(chǔ)上提出一種多閾值圖像分割方法，在林火圖像上進(jìn)行了實驗。分析和實驗表明，HistSplineReg具有可分析求解回歸函數(shù)、擬合速度快、閾值個數(shù)及閾值位置選擇準(zhǔn)確等特點。實驗過程中發(fā)現(xiàn)林火類圖像的紅色信息更值得關(guān)注，但其他顏色通道也能提供更多的信息參考，此外，在參數(shù)選擇方面，HistSplineReg在實驗中建議更多的關(guān)注點應(yīng)該放在經(jīng)驗誤差項，而不是光滑性。

林火識別問題的研究目前仍處于起步階段，絕大多數(shù)的研究工作仍集中在顏色信息部分。加之國際上目前還沒有標(biāo)準(zhǔn)的林火類圖像標(biāo)準(zhǔn)庫，算法性能測試只能在研究者自己擁有的少數(shù)圖像上進(jìn)行實驗，性能比較缺乏公平性。據(jù)人類對火的認(rèn)知，除顏色以外，還會有溫度、灼熱感、甚至燃燒的聲音等，是一種綜合的體驗，而絕不只是依賴顏色來判斷。其他研究領(lǐng)域中有關(guān)目標(biāo)識別方法，能否對林火識別問題提供借鑒，還有圖像中的“煙”如何識別等問題，均是下一步的工作。

References)

[1] GOLLNER M J. Detection and suppression of fires: a cornerstone of fire protection engineering[J]. Fire Technology, 2016, 52(5): 1193-1196.

[2] KO B C, CHEONG K H, NAM J Y. Fire detection based on vision sensor and support vector machines[J]. Fire Safety Journal, 2009, 44(3): 322-329.

[3] SHEN P, ZHANG L, SONG J, et.al. DM642-based fire detection in video sequences using statistical color model[J]. International Journal of Digital Content Technology and Its Applications, 2013(7): 669-678.

[4] FLEURY A D T, TRIGO F C, MARTINS F P R. A new approach based on computer vision and non-linear Kalman filtering to monitor the nebulization quality of oil flames[J]. Expert Systems with Applications, 2013, 40(12): 4760-4769.

[5] YU C, MEI Z, ZHANG X. A real-time video fire flame and smoke detection algorithm[J]. Procedia Engineering, 2013, 62: 891-898.

[6] SHAO J, WANG G, GUO W, et.al. An image-based fire detection method using color analysis[C]// Proceedings of the 2012 International Conference on computer science and information processing. Piscataway, NJ: IEEE, 2012: 1008-1011.

[7] HORNG W, PENG J, CHEN C. A new image-based real time flame detection method using color analysis[C]// Proceedings of the 2005 IEEE International Conference on Networking, Sensing and Control. Piscataway, NJ: IEEE, 2005: 100-105.

[8] CELIK T, DEMIREL H. Fire detection in video sequences using a generic color model[J]. Fire Safety Journal, 2009, 44(2): 147-158.

[9] 蘇毅. 基于圖像處理的火焰識別與特征參數(shù)提取算法研究[D]. 淮南:安徽理工大學(xué)，2013. (SU Y. Research on the flame recognition and feature parameter extraction based on image processing[D]. Huainan: Anhui University of Science and Technology, 2013.)

[10] ZHAO Y, TANG G, XU M. Hierarchical detection of wildfire flame video from pixel level to semantic level[J]. Expert Systems with Applications, 2015, 42(8): 4097-4104.

[11] ARIFIN A Z, ASANO A. Image segmentation by histogram thresholding using hierarchical cluster analysis[J]. Pattern Recognition Letters, 2006, 27(13): 1515-1521.

[12] 吳一全，孟天亮，吳詩婳. 圖像閾值分割方法研究進(jìn)展20年(1994— 2014)[J]. 數(shù)據(jù)采集與處理，2015，30(1)：1-23.(WU Y Q, MENG T L, WU S H. Research progress of image thresholding methods in recent 20 years (1994— 2014)[J]. Journal of Data Acquisition and Processing, 2015, 30(1): 1-23.)

[13] ZHANG C, XIE Y, LIU D, et.al. Fast threshold image segmentation based on 2D fuzzy Fisher and random local optimized QPSO[J]. IEEE Transactions on Image Processing, 2017, 26(3): 1355-1362.

[14] CHEN S, WANG M. Seeking mutli-thresholds directly from support vectors for image segmentation[J]. Neurocomputing, 2005, 67(1): 335-344.

[15] CHANG J, FAN K, CHANG Y. Multi-modal gray-level histogram modeling and decomposition[J]. Image and Vision Computing, 2002, 20(3): 203-216.

[16] BO C, LIU Z, WANG J. A new model of image segmentation with multi-threshold[J]. Journal of Chemical and Pharmaceutical Research, 2014, 6(7): 836-853.

[17] DIANAT R, KASAEI S. On automatic threshold selection in regression method for change detection in remote sensing images[EB/OL].[2016- 11- 20]. http://ce.sharif.edu/～ipl/Papers/Ist_08_RD.pdf.

[18] GREEN P, SILVERMAN B. Nonparametric Regression and Generalized Linear Models[M]. London: Chapman & Hall, 1994: 85-92.

[19] 韓圣同. 三次光滑樣本的一些應(yīng)用[D]. 上海: 華東師范大學(xué)，2009.(HAN S T. Some applications of cubic smoothing spline[D]. Shanghai: East China Normal University，2009.)

[20] REINSCH C. Smoothing by spline functions[J]. Numerische Machematik, 1967, 10(3):177-183.

[21] YUAN J, WANG D, LI R. Image segmentation using local spectral histograms and linear regression[J]. Pattern Recognition Letters, 2012, 33(5): 615-622.

[22] CHANG C C, LIN C J. LIBSVM: a library for support vector machines[J]. ACM Transactions on Intelligent Systems and Technology, 2011,2(3):1-27.

This work is partially supported by the National Natural Science Foundation of China (61472186，50375057), the Natural Science Foundation of Jiangsu Province (BK20161527).

YANGXubing, born in 1973, Ph.D., associate professor. His research interests include pattern recognition, neural computing.

TANXinyi, born in 1992, M.S. candidate. Her reseach interests include image processing, wireless sensor networks.

ZHANGFuquan, born in 1977, Ph.D., associate professor. His research interests include the Internet of things in forestry.

Forestfireimagesegmentationalgorithmwithadaptivethresholdbasedonsmoothsplinefunction

YANG Xubing*, TAN Xinyi, ZHANG Fuquan

(CollegeofInformationScienceandTechnology,NanjingForestryUniversity,NanjingJiangsu210037,China)

Based on smooth spline principle, a self-adaptive multi-threshold segmentation algorithm HistSplineReg (Spline Regression for Histogram) was proposed. HistSplineReg is a two-step method. Firstly, a smoothing spline function was regressed to fit the one-dimensional image histogram, and then the extreme value was found by the regression function to achieve multi-threshold automatic segmentation of the image. Compared to the existing multi-threshold methods, the advantages of HistSplineReg lie in 5 aspects: 1) it is quite consistent with the human intuition; 2) it is constructed on the smoothing spline, which is a solid mathematic basis; 3) both the number and the size of multiple thresholds can be automatically determined; 4) HistSplineReg can be analytically solved, and its computing burden is mainly concentrated on Cholesky decomposition of the matrix, while the size of matrix depends on the pixel level of the image, rather than the scale of the image; 5) it has only one trade-off parameter to balance the empirical error and regressor’s smoothness. Furthermore, for the forest fire recognition task, an experimental reference value was provided. Finally, experiments were conducted on some digital forest fire images in the RGB (Red, Green, Blue) mode. The experimental results show that the histSplineReg method is more effective than Support Vector Regression (SVR) and Polynomial Fitting (PolyFit), which is based on the grayscale image, the color channel, the color image synthesized by each channel segmentation. And the three methods all reflect the red channel information is most significant to the forest fire image segmentation effect.

image segmentation; smoothing spline function; forest fire recognition; threshold

2017- 05- 16;

2017- 07- 18。

國家自然科學(xué)基金資助項目(61472186，50375057);江蘇省自然科學(xué)基金資助項目(BK20161527)。

楊緒兵(1973—)，男，安徽六安人，副教授，博士，CCF會員，主要研究方向：模式識別、神經(jīng)計算； 覃欣怡(1992—)，女，山西交城人，碩士研究生，主要研究方向：圖像處理、無線傳感網(wǎng)絡(luò)； 張福全(1977—)，男，甘肅玉門人，副教授，博士，主要研究方向：林業(yè)物聯(lián)網(wǎng)。

1001- 9081(2017)11- 3157- 05

10.11772/j.issn.1001- 9081.2017.11.3157

(*通信作者電子郵箱xbyang@njfu.edu.cn)

TP391.41

A