成貴學(xué)，鄭曉楠

(上海電力大學(xué)計(jì)算機(jī)科學(xué)與技術(shù)學(xué)院，上海 200090)

1 引言

無人機(jī)在電力巡檢中的應(yīng)用日益廣泛，高質(zhì)量的數(shù)字圖像在電力設(shè)備故障診斷、電力系統(tǒng)監(jiān)控等領(lǐng)域發(fā)揮著關(guān)鍵作用。針對復(fù)雜拍攝條件如陰天、大霧、大風(fēng)等天氣，無人機(jī)拍攝的電力設(shè)備圖像存在細(xì)節(jié)模糊、亮度不均等問題，需對拍攝的彩色圖像進(jìn)行增強(qiáng)處理，提高圖像對比度，去除無人機(jī)抖動造成的噪聲，突出圖像中的有用信息。

常用的圖像增強(qiáng)方法有空域法，例如直方圖均衡化法、拉普拉斯算子法、中值濾波法等；還有頻域法，例如小波變換、傅里葉變換等。上述方法簡單易行，但存在圖像細(xì)節(jié)損失、邊緣模糊，參數(shù)無法自適應(yīng)選擇等問題。而脈沖耦合神經(jīng)網(wǎng)絡(luò)(Pulse coupled neural networks， PCNN)因其簡單、處理速度快、性能穩(wěn)定等優(yōu)勢被廣泛應(yīng)用于圖像增強(qiáng)。文獻(xiàn)[6]提出基于PCNN內(nèi)部活動項(xiàng)的彩色圖像增強(qiáng)；文獻(xiàn)[7]將改進(jìn)的PCNN 用于線陣CCD立靶測量系統(tǒng)的線陣圖像增強(qiáng)中；文獻(xiàn)[8]提出非下采樣Contourlet變換結(jié)合PCNN的圖像增強(qiáng)。但傳統(tǒng)的PCNN需要根據(jù)經(jīng)驗(yàn)選擇最佳參數(shù)組合，直接影響了圖像增強(qiáng)效果。因此，一些學(xué)者將參數(shù)選擇問題視為目標(biāo)函數(shù)優(yōu)化問題，利用啟發(fā)式智能優(yōu)化算法提高參數(shù)優(yōu)化效率。文獻(xiàn)[9]提出GA優(yōu)化PCNN模型的圖像增強(qiáng)方法。但由于GA算法存在易陷入局部最優(yōu)解，收斂精度低等缺陷，導(dǎo)致圖像增強(qiáng)效果不穩(wěn)定。

樽海鞘群算法(Salp Swarm Aalgorith， SSA)的尋優(yōu)能力主要是個體間的相互作用和影響，個體本身沒有變異機(jī)制，易早熟收斂。本文提出一種改進(jìn)樽海鞘群優(yōu)化算法(Improved Salp Swarm Aalgorith， ISSA)并將其應(yīng)用于PCNN圖像增強(qiáng)領(lǐng)域。解決樽海鞘群算法收斂精度差且易陷入局部最優(yōu)解的缺陷，在保證算法運(yùn)行效率的同時提高圖像增強(qiáng)的效果，對算法進(jìn)行了三部分改進(jìn)：在領(lǐng)導(dǎo)者位置更新中引入高斯變異，提高種群多樣性；在追隨者位置更新中引入自適應(yīng)慣性權(quán)重，平衡算法局部開發(fā)和全局搜索；設(shè)定早熟判別機(jī)制來判斷算法是否早熟收斂，并通過柯西變異跳出局部最優(yōu)解。最后，選取5個標(biāo)準(zhǔn)函數(shù)，對比ISSA、PSO、DE、BA、SSA算法，并在電力巡檢圖像上驗(yàn)證算法的有效性。

2 圖像增強(qiáng)模型

將電力巡檢彩色圖像從RGB顏色空間轉(zhuǎn)換到HSI顏色空間，得到圖像亮度分量I。對圖像亮度分量I進(jìn)行非下采樣剪切波變換(Non-subsampled shearlet transform， NSST)分解，得到低頻子帶圖像和高頻子帶圖像。

高頻子帶圖像包含原圖像的邊緣、輪廓和噪聲等信息。采取自適應(yīng)閾值對多尺度多方向的高頻子帶系數(shù)去噪，去噪過程如下

(1)

(2)

m

為

C

，，的均值。

低頻子帶圖像包含了原圖像的主體信息，利用PCNN模型增強(qiáng)低頻子帶圖像，PCNN神經(jīng)元與圖像像素點(diǎn)一一對應(yīng)，由接收器，調(diào)制器和脈沖發(fā)生器組成。

簡化的PCNN數(shù)學(xué)模型為

F

(

n

)=

S

(3)

L

(

n

)=

V

∑

W

Y

(4)

U

(

n

)=

F

(1+

βL

(

n

))

(5)

θ

(

n

)=

e

-

θ

(

n

-1)+

V

Y

(

n

-1)

(6)

(7)

其中，

S

為第(

i

，

j

)個神經(jīng)元的刺激信號作用于反饋輸入

F

(

n

)；

L

(

n

)為第(

i

，

j

)個神經(jīng)元的鏈接輸入；

Y

(

n

)為輸出值。

V

為通道幅值；

W

為

L

(

n

)的加權(quán)系數(shù)；

β

為鏈接系數(shù)；

V

、

α

分別為

θ

(

n

)的固有電勢和時間衰減常數(shù)。圖1為PCNN神經(jīng)元數(shù)學(xué)模型：

圖1 PCNN神經(jīng)元的簡化模型

利用PCNN增強(qiáng)低頻子帶圖像的過程中需確定參數(shù)

α

、

β

、

V

的取值，且參數(shù)取值影響圖像增強(qiáng)的效果，因而采用改進(jìn)的樽海鞘群算法尋優(yōu)，實(shí)現(xiàn)參數(shù)的自適應(yīng)選值，提高圖像增強(qiáng)效果。

增強(qiáng)后的低頻子帶圖像和去噪后的高頻子帶圖像進(jìn)行NSST逆變換后代替原圖像的I分量，從HSI空間轉(zhuǎn)換到RGB空間得到增強(qiáng)后的圖像。

圖2 圖像增強(qiáng)過程

3 樽海鞘群算法及改進(jìn)

3.1 樽海鞘群算法(SSA)

樽海鞘鏈由兩種類型的樽海鞘組成：領(lǐng)導(dǎo)者和追隨者。領(lǐng)導(dǎo)者位于鏈的最前端，其它樽海鞘作為追隨者跟隨領(lǐng)導(dǎo)者依次移動。

SSA算法中樽海鞘的位置是用

N

×

D

的矩陣

X

表示，公式如式(8)所示，其中

N

為樽海鞘種群的數(shù)量，

D

為樽海鞘個體位置的維度。領(lǐng)導(dǎo)者的位置更新如式(9)所示。

(8)

(9)

c

=2

e

-(4max)

(10)

t

為當(dāng)前迭代數(shù)，

T

為最大迭代數(shù)。

追隨者跟隨領(lǐng)導(dǎo)者運(yùn)動，它的位置只與初始位置，速度，加速度有關(guān)。根據(jù)牛頓運(yùn)動定律更新追隨者的位置，如下所示(11)

(11)

(12)

3.2 改進(jìn)樽海鞘群優(yōu)化算法(ISSA)

SSA算法在解決具有大量局部最優(yōu)解的多峰、高維問題時，容易出現(xiàn)局部停滯現(xiàn)象，導(dǎo)致算法收斂精度低。針對上述問題，選取種群前一半個體為領(lǐng)導(dǎo)者，后一半個體為追隨者，并引入高斯變異、自適應(yīng)慣性權(quán)重和早熟判別機(jī)制對SSA算法進(jìn)行改進(jìn)，增強(qiáng)算法的全局搜索能力，提高算法的尋優(yōu)精度。

3.2.1 引入高斯變異

在算法中引入高斯變異對領(lǐng)導(dǎo)者位置更新公式進(jìn)行改進(jìn)，增強(qiáng)算法種群多樣性，提高算法的全局搜索能力，更新后的公式如式(13)所示

國慶節(jié)前，牛皮糖生日，照例請客。以前請客請的是三親四戚。這回他沒有請親戚，卻請了幾個八竿子打不著的人，但個個都是有點(diǎn)來頭的。席上一頓酒肉胡勸亂吃之后，一幫人就在地坪里就著暖洋洋的太陽扯白。

(13)

3

.

2

.

2 引入自適應(yīng)慣性權(quán)重根據(jù)追隨者位置更新公式可知，每個樽海鞘的位置只與其前一個樽海鞘的位置有關(guān)，為了提高追隨者搜索能力同時平衡算法的局部開發(fā)和全局探索，引入自適應(yīng)慣性權(quán)重系數(shù)

w

。算法初始，

w

較大，全局搜索性能較強(qiáng)，

w

隨著算法運(yùn)行次數(shù)增加而自適應(yīng)減小，算法局部開發(fā)能力提高。

(14)

α

為[15，30]之間取值；

w

=0

.

9、

w

=0

.

4分別為慣性權(quán)重系數(shù)的上限和下限。

改進(jìn)后的追隨者位置更新公式如下

(15)

3

.

2

.

3 早熟判別機(jī)制鑒于樽海鞘群算法在迭代過程中容易早熟，陷入局部最優(yōu)解，因此，利用當(dāng)前樽海鞘群體的方差

σ

是否小于某個閾值

ε

(

σ

反映樽海鞘群體中個體適應(yīng)度值的整體變化，值越小表示個體聚集度越高，種群多樣性低)和

f

/f

是否趨于1(趨于1時表明樽海鞘已經(jīng)陷入停滯狀態(tài))來判斷算法是否陷入局部最優(yōu)解。當(dāng)

σ

小于某個方差閾值

ε

和

f

/f

趨于1時，則判斷算法陷入局部最優(yōu)解。具體步驟如下：

將樽海鞘位置帶入適應(yīng)度函數(shù)，計(jì)算每個樽海鞘的適應(yīng)度值。

Fitness

=

f

(

x

)

(16)

找出最佳食物源。

(17)

記錄并保留最佳食物源

F

。

f

=max(

Fitness

)

(18)

f

=

F

(19)

δ

=

f

/f

(20)

(21)

(22)

若

σ

≤

ε

和

δ

趨于1，則對領(lǐng)導(dǎo)者進(jìn)行柯西變異，公式如下

(23)

如圖3所示，柯西分布函數(shù)兩端延伸比高斯分布函數(shù)長，該特點(diǎn)使柯西變異產(chǎn)生的個體比原個體有更大的差異，更易跳出局部最優(yōu)解，中心點(diǎn)處峰值較小的特點(diǎn)使其花費(fèi)更少的時間搜索領(lǐng)域空間，容易調(diào)節(jié)。

圖3 高斯分布和柯西分布曲線

4 實(shí)驗(yàn)結(jié)果分析

4.1 基于測試函數(shù)的仿真及分析

為了驗(yàn)證ISSA算法的適用性和優(yōu)越性，對PSO算法(P)、DE算法(D)、BA算法(B)、SSA算法(S)、ISSA算法(I)分別采取5個標(biāo)準(zhǔn)測試函數(shù)進(jìn)行仿真。各測試函數(shù)的具體參數(shù)設(shè)置如表1所示。

f

、

f

為單峰測試函數(shù)，

f

、

f

、

f

為多峰測試函數(shù)。各算法的參數(shù)設(shè)置如表2所示。實(shí)驗(yàn)環(huán)境：Intel(R) Core(TM) i5-8250U CPU @ 1.6GHZ 處理器；8.0GB運(yùn)行環(huán)境；64位操作系統(tǒng)，基于x64的處理器；Windows10下Matlab2014版本的電腦上測試。

1)Sphere函數(shù)

2)Schwelel2.22函數(shù)

3)Ackley函數(shù)

4) Rastrigin函數(shù)

5)Griewank函數(shù)

表1 測試函數(shù)參數(shù)設(shè)置

表2 實(shí)驗(yàn)參數(shù)設(shè)置

4.2 算法性能測試

圖4為不同算法在不同測試函數(shù)上迭代500次，適應(yīng)度值雖迭代次數(shù)增加的變化曲線，由圖可以看出，

ISSA

算法相比于其它四種算法下降趨勢更明顯，說明

ISSA

算法收斂性更好。

4.3 實(shí)驗(yàn)測試結(jié)果分析

為了進(jìn)一步驗(yàn)證算法的優(yōu)越性，對每個標(biāo)準(zhǔn)測試函數(shù)進(jìn)行30次測試和比較，其結(jié)果如表3所示。平均值為30次獨(dú)立運(yùn)行算法在測試函數(shù)中能夠達(dá)到的精度和的平均值。標(biāo)準(zhǔn)差可以反映算法的穩(wěn)定性效果

ISSA

算法標(biāo)準(zhǔn)差較低說明算法尋優(yōu)過程更穩(wěn)定。最優(yōu)值、最差值和平均值反映算法的收斂精度和尋優(yōu)能力，對比表3中數(shù)據(jù)可以看出，

ISSA

算法比

PSO

算法(

P

)、

DE

算法(

D

)、

BA

算法(

B

)、

SSA

算法(

S

)尋優(yōu)能力更強(qiáng)。綜上，

ISSA

算法(

I

)具有更好的全局搜索性，綜合性也更好，能精準(zhǔn)地找到最優(yōu)值。

圖4 收斂曲線圖

表3 測試實(shí)驗(yàn)結(jié)果對比

4.4 電力巡檢圖像增強(qiáng)分析

為了驗(yàn)證增強(qiáng)算法的有效性，獲取了四幅不同角度、不同距離的無人機(jī)電力巡檢圖片，分別用多尺度Retinex算法(MSRCR)、文獻(xiàn)[8]算法、直方圖均衡化算法(HE)、小波算法(WT)和ISSA-PCNN算法進(jìn)行圖像增強(qiáng)實(shí)驗(yàn)對比。

利用常用的兩個評價標(biāo)準(zhǔn)——平均梯度(AG)和空間頻率(SF)作為ISSA-PCNN算法的適應(yīng)度函數(shù)。平均梯度描述圖像的細(xì)節(jié)和紋理特征?？臻g頻率顯示圖像空間的活動程度，值越大，圖像增強(qiáng)的方法性能越好。適應(yīng)度函數(shù)如下：

min

f

=1

/

(

SF

+

AG

)

(24)

具體步驟如圖5所示。

圖5 改進(jìn)樽海鞘群算法優(yōu)化PCNN參數(shù)步驟

文獻(xiàn)[8]算法采用傳統(tǒng)PCNN參數(shù)設(shè)置

α

=1

.

000、

β

=0

.

100、

V

=20

.

000，

PCNN

其它參數(shù)設(shè)置與

ISSA

-

PCNN

算法相同：

V

=1、

T

=100、

W

=[0

.

5 1 0

.

5；1 0 1；0

.

5 1 0

.

5]。

圖6—圖9為不同圖像增強(qiáng)方法增強(qiáng)圖像結(jié)果。從視覺效果來看，MSRCR算法提高了圖像的對比度，但顏色較原圖像失真嚴(yán)重。文獻(xiàn)[8]算法、小波變換算法和直方圖均衡化算法增強(qiáng)效果不穩(wěn)定，對某些圖像會出現(xiàn)過度增強(qiáng)，導(dǎo)致圖像失真。ISSA-PCNN算法增強(qiáng)效果更穩(wěn)定，增強(qiáng)后的圖像更清晰自然，失真小。

圖6 第一組圖像增強(qiáng)方法對比

圖7 第二組圖像增強(qiáng)方法對比

圖8 第三組圖像增強(qiáng)方法對比

圖9 第四組圖像增強(qiáng)方法對比

運(yùn)用三個客觀評價指標(biāo)(圖像熵(IE)、空間頻率(SF)、平均梯度(AG))對其結(jié)果進(jìn)行客觀評價，數(shù)據(jù)如表4所示，分析發(fā)現(xiàn)ISSA-PCNN算法指標(biāo)總體高于另外四種算法的指標(biāo)，綜合以上分析，ISSA-PCNN算法具有良好的穩(wěn)定性和圖像增強(qiáng)能力，能更好地完成電力巡檢圖像的增強(qiáng)任務(wù)。

表4 圖像增強(qiáng)結(jié)果的客觀評價

5 總結(jié)

通過引入高斯變異算子、自適應(yīng)慣性權(quán)重和早熟判別機(jī)制，本文提出一種基于改進(jìn)樽海鞘群算法優(yōu)化脈沖耦合神經(jīng)網(wǎng)絡(luò)參數(shù)的圖像增強(qiáng)技術(shù)。從測試函數(shù)的仿真中可以看出，針對不同測試函數(shù)，改進(jìn)樽海鞘群算法均優(yōu)于標(biāo)準(zhǔn)樽海鞘群算法。將改進(jìn)樽海鞘群算法應(yīng)用于圖像增強(qiáng)模型中，仿真結(jié)果表明相比于其它圖像增強(qiáng)方法，該方法增強(qiáng)后的圖像從視覺角度看更清晰自然、對比度更高，從客觀角度分析圖像保真性更好。綜上，實(shí)驗(yàn)結(jié)果表明了改進(jìn)樽海鞘群算法具有更好的全局搜索和局部開發(fā)能力，算法的有效性和魯棒性也得到了驗(yàn)證。