成貴學(xué)，鄭曉楠

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

1 引言

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

常用的圖像增強(qiáng)方法有空域法，例如直方圖均衡化法、拉普拉斯算子法、中值濾波法等；還有頻域法，例如小波變換、傅里葉變換等。上述方法簡(jiǎn)單易行，但存在圖像細(xì)節(jié)損失、邊緣模糊，參數(shù)無(wú)法自適應(yīng)選擇等問(wèn)題。而脈沖耦合神經(jīng)網(wǎng)絡(luò)(Pulse coupled neural networks， PCNN)因其簡(jiǎn)單、處理速度快、性能穩(wěn)定等優(yōu)勢(shì)被廣泛應(yīng)用于圖像增強(qiáng)。文獻(xiàn)[6]提出基于PCNN內(nèi)部活動(dòng)項(xiàng)的彩色圖像增強(qiáng)；文獻(xiàn)[7]將改進(jìn)的PCNN 用于線陣CCD立靶測(cè)量系統(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ù)選擇問(wèn)題視為目標(biāo)函數(shù)優(yōu)化問(wèn)題，利用啟發(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)能力主要是個(gè)體間的相互作用和影響，個(gè)體本身沒(méi)有變異機(jī)制，易早熟收斂。本文提出一種改進(jìn)樽海鞘群優(yōu)化算法(Improved Salp Swarm Aalgorith， ISSA)并將其應(yīng)用于PCNN圖像增強(qiáng)領(lǐng)域。解決樽海鞘群算法收斂精度差且易陷入局部最優(yōu)解的缺陷，在保證算法運(yùn)行效率的同時(shí)提高圖像增強(qiáng)的效果，對(duì)算法進(jìn)行了三部分改進(jìn)：在領(lǐng)導(dǎo)者位置更新中引入高斯變異，提高種群多樣性；在追隨者位置更新中引入自適應(yīng)慣性權(quán)重，平衡算法局部開(kāi)發(fā)和全局搜索；設(shè)定早熟判別機(jī)制來(lái)判斷算法是否早熟收斂，并通過(guò)柯西變異跳出局部最優(yōu)解。最后，選取5個(gè)標(biāo)準(zhǔn)函數(shù)，對(duì)比ISSA、PSO、DE、BA、SSA算法，并在電力巡檢圖像上驗(yàn)證算法的有效性。

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

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

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

(1)

(2)

m

為

C

，，的均值。

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

簡(jiǎn)化的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

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

F

(

n

)；

L

(

n

)為第(

i

，

j

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

Y

(

n

)為輸出值。

V

為通道幅值；

W

為

L

(

n

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

β

為鏈接系數(shù)；

V

、

α

分別為

θ

(

n

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

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

利用PCNN增強(qiáng)低頻子帶圖像的過(guò)程中需確定參數(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)過(guò)程

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

3.1 樽海鞘群算法(SSA)

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

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

N

×

D

的矩陣

X

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

N

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

D

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

(8)

(9)

c

=2

e

-(4max)

(10)

t

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

T

為最大迭代數(shù)。

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

(11)

(12)

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

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

3.2.1 引入高斯變異

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

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

(13)

3

.

2

.

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

w

。算法初始，

w

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

w

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

(14)

α

為[15，30]之間取值；

w

=0

.

9、

w

=0

.

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

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

(15)

3

.

2

.

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

σ

是否小于某個(gè)閾值

ε

(

σ

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

f

/f

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

σ

小于某個(gè)方差閾值

ε

和

f

/f

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

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

Fitness

=

f

(

x

)

(16)

找出最佳食物源。

(17)

記錄并保留最佳食物源

F

。

f

=max(

Fitness

)

(18)

f

=

F

(19)

δ

=

f

/f

(20)

(21)

(22)

若

σ

≤

ε

和

δ

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

(23)

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

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

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

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

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

f

、

f

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

f

、

f

、

f

為多峰測(cè)試函數(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版本的電腦上測(cè)試。

1)Sphere函數(shù)

2)Schwelel2.22函數(shù)

3)Ackley函數(shù)

4) Rastrigin函數(shù)

5)Griewank函數(shù)

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

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

4.2 算法性能測(cè)試

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

ISSA

算法相比于其它四種算法下降趨勢(shì)更明顯，說(shuō)明

ISSA

算法收斂性更好。

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

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

ISSA

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

ISSA

算法比

PSO

算法(

P

)、

DE

算法(

D

)、

BA

算法(

B

)、

SSA

算法(

S

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

ISSA

算法(

I

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

圖4 收斂曲線圖

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

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

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

利用常用的兩個(gè)評(píng)價(jià)標(biāo)準(zhǔn)——平均梯度(AG)和空間頻率(SF)作為ISSA-PCNN算法的適應(yīng)度函數(shù)。平均梯度描述圖像的細(xì)節(jié)和紋理特征。空間頻率顯示圖像空間的活動(dòng)程度，值越大，圖像增強(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é)果。從視覺(jué)效果來(lái)看，MSRCR算法提高了圖像的對(duì)比度，但顏色較原圖像失真嚴(yán)重。文獻(xiàn)[8]算法、小波變換算法和直方圖均衡化算法增強(qiáng)效果不穩(wěn)定，對(duì)某些圖像會(huì)出現(xiàn)過(guò)度增強(qiáng)，導(dǎo)致圖像失真。ISSA-PCNN算法增強(qiáng)效果更穩(wěn)定，增強(qiáng)后的圖像更清晰自然，失真小。

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

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

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

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

運(yùn)用三個(gè)客觀評(píng)價(jià)指標(biāo)(圖像熵(IE)、空間頻率(SF)、平均梯度(AG))對(duì)其結(jié)果進(jìn)行客觀評(píng)價(jià)，數(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é)果的客觀評(píng)價(jià)

5 總結(jié)

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