基于人工神經(jīng)網(wǎng)絡(luò)的遙感圖像去條帶方法

引言

紅外焦平面探測(cè)是一種兼具輻射敏感和信號(hào)處理功能的新一代紅外探測(cè)技術(shù)，但是由于制造過(guò)程和工作環(huán)境的影響， 使得焦平面陣列(FPA ) 各個(gè)陣列元即使在相同的輻射通量照射下，也會(huì)輸出不相同的響應(yīng)電壓。這種紅外響應(yīng)引起的遙感圖像的失真被稱作紅外圖像的非均勻性。為了提高觀測(cè)頻率、掃描范圍和空間分辨率，航天遙感一般采用推掃式的多元敏感線陣列對(duì)地物成像，通過(guò)觀察發(fā)現(xiàn)，推掃得到的遙感圖像出現(xiàn)有規(guī)律的條帶失真，條帶寬度與遙感器多元敏感元個(gè)數(shù)的掃描線寬度一致，而且隨著時(shí)間的推移，條帶現(xiàn)象日趨嚴(yán)重，與單敏感元掃描圖像中的噪聲相比有明顯差異，這種失真其實(shí)是焦平面非均勻性的一種表現(xiàn)形式。條帶失真是影響線陣列紅外遙感圖像質(zhì)量的主要因素，必須要用諸如定標(biāo)的方法去除，但是由于探測(cè)器單元響應(yīng)會(huì)隨著時(shí)間和工作環(huán)境的變化改變，所以僅僅用定標(biāo)的方法來(lái)校正條帶失真有很大的局限性。

本文對(duì)焦平面非均勻性校正的神經(jīng)網(wǎng)絡(luò)法進(jìn)行改進(jìn)，介紹了一種基于人工神經(jīng)網(wǎng)絡(luò)的遙感圖像條帶消除的方法。這種方法可以完全不對(duì)FPA進(jìn)行標(biāo)定(或自動(dòng)標(biāo)定)，并且可以通過(guò)線性和非線性模型校正，是紅外成像系統(tǒng)非均勻性校正的理想方法。

傳統(tǒng)校正方法

線性校正假設(shè)探測(cè)器單元的響應(yīng)呈線性：

y = ax + b

式中，x 為某一探測(cè)器單元的輸入信號(hào)，y 為可測(cè)的輸出信號(hào)。如果能求出增益因子a 和偏移因子b，就可求得無(wú)畸變的輸入信號(hào)x 。

傳統(tǒng)的非均勻性校正方法是在紅外焦平面成像系統(tǒng)使用前，用標(biāo)準(zhǔn)的兩個(gè)或多個(gè)參考溫度源，對(duì)每一個(gè)陣列單元響應(yīng)進(jìn)行定標(biāo)，以保證每個(gè)陣列單元在兩個(gè)或多個(gè)參考溫度之間有相同的響應(yīng)，其校正值被存儲(chǔ)起來(lái)，在進(jìn)行數(shù)字處理時(shí)固定地將其疊加上去。如果每個(gè)陣列單元的輸出特性隨時(shí)間是完全線性和穩(wěn)定的，那么，在上述定標(biāo)溫度范圍內(nèi)，這種校正是有效的，不過(guò)隨著陣列數(shù)的增加，存儲(chǔ)校正系數(shù)所需要的存儲(chǔ)容量就大為增加。再加上系統(tǒng)的不穩(wěn)定性、陣列單元的非線性和1/ f 噪聲等因素的影響，使得經(jīng)過(guò)一段時(shí)間后，陣列單元特性會(huì)發(fā)生漂移或溫度背景范圍出現(xiàn)變化，必須對(duì)紅外焦平面陣列進(jìn)行再定標(biāo)。顯然，這類(lèi)校正方法不但麻煩，而且可能并不符合實(shí)際使用情況，從而導(dǎo)致校正效果不佳，因此，必須研究自適應(yīng)的非均勻性校正方法。

神經(jīng)網(wǎng)絡(luò)法

神經(jīng)網(wǎng)絡(luò)法的主要特征是通過(guò)自學(xué)習(xí)模擬信息內(nèi)部關(guān)系，進(jìn)而獲得系統(tǒng)特征參數(shù)。假設(shè)輸入x和輸出y之間有一種復(fù)雜的關(guān)系f，神經(jīng)網(wǎng)絡(luò)通過(guò)不斷調(diào)整結(jié)構(gòu)的權(quán)重系數(shù)和閾值得到逼近的關(guān)系f?，蕼厦y=f’(x)； Rumelhart 和Mcclalland提出的多層前饋網(wǎng)絡(luò)的反向傳播算法(BP算法) ，由于解決了感知器不能解決的多層網(wǎng)絡(luò)學(xué)習(xí)算法的問(wèn)題，可以很好地對(duì)復(fù)雜函數(shù)進(jìn)行逼近，在工程中得到了廣泛的應(yīng)用。一般使用的BP網(wǎng)絡(luò)是一個(gè)三層前向網(wǎng)絡(luò)，結(jié)構(gòu)如圖1所示。

假定校正輸出為Y(n)，輸入為X(n)，則：

Y(n)＝WT(n)X(n)+VT(n)

其中W和V是增益矢量和截距矢量，神經(jīng)網(wǎng)絡(luò)法就是不斷依據(jù)實(shí)際景像調(diào)整W和V，來(lái)去除條帶失真。根據(jù)三層BP網(wǎng)絡(luò)結(jié)構(gòu)，在中間層根據(jù)一定的法則計(jì)算某像素輸出，并以此作為該像素的輸出，反饋給線性校正神經(jīng)元來(lái)調(diào)整W和V。調(diào)整以誤差信號(hào)均方值最小為準(zhǔn)則。

對(duì)每一次迭代，令期望響應(yīng)與輸出響應(yīng)之差為誤差，用e(n)表示，則：

e(n)=f(n)-Y(n)=f(n)-WT(n)X(n)+VT(n)

其中，f(n)表示期望的校正后輸出，則誤差函數(shù)為：F(W,V)=(Wx+V-f)2，利用最陡下降法，可以得到計(jì)算W和V的迭代公式：

Wn+1=Wn-2ax(y-f)
Vn+1=Vn-2a(y-f)

式中：n為幀數(shù),a為步長(zhǎng)。

神經(jīng)網(wǎng)絡(luò)算法的改進(jìn)

從上面的分析可以得到，神經(jīng)網(wǎng)絡(luò)法對(duì)非均勻性的校正的關(guān)鍵是如何建立期望的校正輸出模型，在非線陣列的焦平面非均勻性校正中，一般將校正元相鄰元的輸出平均值作為本元的期望輸出帶入網(wǎng)絡(luò)進(jìn)行網(wǎng)絡(luò)訓(xùn)練，發(fā)展出了4領(lǐng)域和8領(lǐng)域等方法。對(duì)于線陣列，由于探測(cè)單元只有兩個(gè)相鄰的探測(cè)元，直接應(yīng)用上面的方法進(jìn)行網(wǎng)絡(luò)訓(xùn)練，校正效果不是太好。鑒于此，對(duì)算法進(jìn)行如下改進(jìn)。

假設(shè)一幅圖像有n條掃描線組成，對(duì)于每一條掃描線響應(yīng)，可以用Yk(i)表示，其中k表示第k條掃描線 ，i表示線陣列的第i個(gè)探測(cè)元?？梢詫⒕€陣列擴(kuò)展為有三條線陣列的焦平面，在第k次成像時(shí)，焦平面成像為[Yk-1(i):Yk(i):Yk+1(i)]，這樣就可以假定Yk(i)的期望校正輸出為：

Y’k(i)＝1/8(Yk(i-1)+Yk(i+1)+Yk-1(i-1)+Yk-1(i)+Yk-1(i+1)+Yk+1(i-1)+Yk+1(i)+Yk+1(i+1))

算法過(guò)程如下：

1.計(jì)算鄰域平均值:

Y’k(i)＝1/8(Yk(i-1)+Yk(i+1)+Yk-1(i)+Yk+1(i)+Yk-1(i-1)+Yk+1(i+1)+Yk-1(i-1)+Yk+1(i+1))

2. 令y = Wx + V ,其中W為增益校正因子,V為偏移量校正因子。誤差函數(shù):

F(W,V)=(Wx+V-f)2

利用此函數(shù)的梯度函數(shù)和最陡下降法,可以得到計(jì)算和的迭代公式：

Wn+1=Wn-2ax(y-f)
Vn+1=Vn-2a(y-f)

式中 n為幀數(shù)，a為步長(zhǎng)。

3.利用線性校正算法得到：

Yn+1=Wn+1