基于分層稀疏濾波卷積神經(jīng)網(wǎng)絡(luò)的sar圖像分類方法
【技術(shù)領(lǐng)域】
[0001] 本發(fā)明屬于圖像處理技術(shù)領(lǐng)域，更進(jìn)一步涉及一種SAR圖像分類方法，可用于目 標(biāo)識別。
【背景技術(shù)】
[0002] 合成孔徑雷達(dá)SAR是一種微波成像雷達(dá)，具有良好的分辨率，不僅可以詳細(xì)地、準(zhǔn) 確地觀測地形、地貌，獲取地球表面信息，還可以透過地表和自然植被收集地表以下的信 息。SAR是從空間對地觀測的一種有效手段，能夠生成地面目標(biāo)區(qū)域或地域的高分辨率地 圖，提供類似于光學(xué)照片的雷達(dá)圖像，已廣泛應(yīng)用于軍事以及其它對地觀測領(lǐng)域。
[0003] 合成孔徑雷達(dá)的概念是1951年6月由美國Goodyear宇航公司的Carl Wiley首次 提出的。SAR是一種主動式微波成像傳感器，它利用脈沖壓縮技術(shù)提高距離分辨率，利用合 成孔徑原理提高方位分辨率，從而獲得大面積的高分辨率雷達(dá)圖像，具有全天時、全天候、 多波段、多極化、可變側(cè)視角及高分辨率等優(yōu)點(diǎn)，甚至在惡劣的環(huán)境下也能以較高的分辨率 提供詳細(xì)的地面測繪數(shù)據(jù)和圖像。我國從20世紀(jì)70年代中期開始SAR系統(tǒng)的研制工作， 先后取得了一定的研究成果，1979年9月，中科院電子所研制的機(jī)載SAR原理樣機(jī)試飛成 功，獲得我國第一批SAR圖像。我國第一顆SAR衛(wèi)星已躋身國際先進(jìn)行列，目前已進(jìn)入實(shí)際 應(yīng)用階段，并在國土測繪、資源普查、城市規(guī)劃、搶險(xiǎn)救災(zāi)等領(lǐng)域發(fā)揮了重要的作用。
[0004] SAR技術(shù)具有如下特有的優(yōu)勢：
[0005] I) SAR成像不依賴光照，而是靠自己發(fā)射的微波，能夠穿透云、雨、雪和煙霧，具有 全天時、全天候成像能力，這是SAR遙感最突出的優(yōu)勢。
[0006] 2)微波對地表有一定的穿透能力。
[0007] 3)對金屬目標(biāo)及地表紋理特征有較強(qiáng)的探測能力。
[0008] 現(xiàn)有的經(jīng)典的SAR圖像分類方法主要有以下兩類：
[0009] ( -）從特征入手。根據(jù)全極化SAR數(shù)據(jù)的特點(diǎn)，根據(jù)其數(shù)據(jù)分布特性或散射機(jī)理 提取包含極化信息的特征，設(shè)計(jì)分類方法以完成地物分類。該類算法大概可以細(xì)分為3種： 一種是基于極化SAR統(tǒng)計(jì)特性的分類方法；第二種是基于極化SAR散射機(jī)理的分類方法； 第三種是結(jié)合極化SAR統(tǒng)計(jì)分布和散射機(jī)理的分類方法。
[0010] (二）從處理方法入手。在已有的特征集上，引入更有效的處理方法，從而更充分 的利用現(xiàn)有的分類信息。SVM，Adab〇〇st以及神經(jīng)網(wǎng)絡(luò)等方法均屬于此類，目前它們在極化 SAR分類解譯方面都取得了大量優(yōu)秀的研究成果。
[0011] 但上述方法與光學(xué)圖像相比，由于SAR圖像視覺可讀性較差，使得SAR圖像信息處 理非常困難。另一方面，隨著SAR應(yīng)用的日趨廣泛以及技術(shù)的不斷成熟，其數(shù)據(jù)信息也在急 劇增長，SAR所收集的數(shù)據(jù)量之大已經(jīng)遠(yuǎn)遠(yuǎn)超出了人工作出迅速判斷的極限。這些因素都限 制了傳統(tǒng)的圖像分類技術(shù)如基于模板匹配、基于模型和基于核的分類技術(shù)在SAR圖像分類 中的應(yīng)用。目前SAR圖像識別技術(shù)主要有三個問題亟需解決：（1)由于SAR圖像中存在大 量的相干斑噪聲，采用常用的特征提取方法很難克服噪聲的影響，分類精度不高；(2)由于 SAR圖像中同類地物的場景復(fù)雜，傳統(tǒng)的特征提取方法在設(shè)計(jì)上費(fèi)時費(fèi)力，且具有較大的局 限性，不具備自適應(yīng)性。（3)由于對SAR圖像地物的標(biāo)注過程比較繁瑣費(fèi)力，因此需在標(biāo)記 樣本較少的情況下進(jìn)行分類，而傳統(tǒng)的分類方法在這種情況下，分類精度較低，分類結(jié)果不 穩(wěn)定。

【發(fā)明內(nèi)容】

[0012] 本發(fā)明的目的在于針對上述已有技術(shù)的不足，提出一種基于分層稀疏濾波卷積神 經(jīng)網(wǎng)絡(luò)的SAR圖像分類方法，利用深度神經(jīng)網(wǎng)絡(luò)提取SAR圖像局部和全局特征，提高SAR圖 像的分類精度。
[0013] 本發(fā)明的技術(shù)方案是：通過逐層訓(xùn)練適應(yīng)于SAR圖像的稀疏濾波器，構(gòu)建多層稀 疏濾波卷積神經(jīng)網(wǎng)絡(luò)，用于提取SAR圖像局部和全局的特征，進(jìn)而訓(xùn)練分類器，達(dá)到對SAR 圖像分類的目的。其實(shí)現(xiàn)步驟包括如下：
[0014] (1)將SAR圖像數(shù)據(jù)庫樣本集劃分為訓(xùn)練數(shù)據(jù)集X和測試樣本集y ;
[0015] ⑵訓(xùn)練SVM分類器：
[0016] 2a)從訓(xùn)練數(shù)據(jù)集X中隨機(jī)抽取m塊尺寸dX d的訓(xùn)練圖像塊，并進(jìn)行全局對比度 歸一化，構(gòu)成訓(xùn)練圖像塊集JTe ?
[0017] 2b)利用訓(xùn)練圖像塊集X訓(xùn)練第一層稀疏字典巧，.其中N表示X中每個圖 像塊的特征數(shù)目；
[0018] 2c)利用第一層的稀疏字典01求訓(xùn)練集X的第一層稀疏特征圖： Z e #><(^+1)><('^+1)，其中 11^分別表示圖片的高度和寬度；
[0019] 2d)對第一層稀疏特征圖Z進(jìn)行非線性變換，得到特征圖=C1G RNX(U d+1)/wX(v d+1)/w， 其中w表示池化的比例；
[0020] 2e)從訓(xùn)練集X的特征圖C1上隨機(jī)抽取m2塊尺寸NXd 2Xcy^訓(xùn)練圖像塊，構(gòu)成 訓(xùn)練集
[0021] 2f)利用訓(xùn)練集X2采用與2b)相同的方法，訓(xùn)練第二層稀疏字典：其 中N 2表示X 2中每個圖像塊的特征數(shù)量；
[0022] 2g)利用第二層的稀疏字典D2采用與2c)相同的方法，求訓(xùn)練集X的第二層稀疏 特征圖
[0023] 2h)對第二層稀疏特征圖乏進(jìn)行與2d)相同的非線性變換，得到非線性變換特征 圖C 2;
[0024] 2i)級聯(lián)(^和C 2構(gòu)成一維向量c，訓(xùn)練線性核SVM分類器；
[0025] (3)抽取測試集y的特征并進(jìn)行分類，得到分類結(jié)果：
[0026] 3a)對測試集y利用訓(xùn)練階段獲得的第一層稀疏字典D1和第二層稀疏字典D 2，采 用與訓(xùn)練集X相同的非線性變換方法抽取測試集第一層和第二層的非線性變換特征 < 和 ζ，級聯(lián)石和6構(gòu)成一維向量L
[0027] 3b)將一維向量?輸入到SVM分類器進(jìn)行分類，得到最終分類結(jié)果。
[0028] 與現(xiàn)有技術(shù)相比，本發(fā)明具有以下優(yōu)點(diǎn)：
[0029] 本發(fā)明通過逐層無監(jiān)督的訓(xùn)練獲得適應(yīng)SAR圖像特征分布的稀疏濾波器，相比費(fèi) 時費(fèi)力通過手工設(shè)計(jì)的特征提取方法如，SIFT，HOG等更具有普適性，能夠很好的克服相干 斑噪聲的影響，同時通過提取SAR圖像深層的特征，在標(biāo)記樣本很少的情況下，仍能達(dá)到很 高分類精度和非常穩(wěn)定的分類結(jié)果。
【附圖說明】
[0030] 圖1是本發(fā)明的實(shí)現(xiàn)流程圖；
[0031] 圖2是本發(fā)明仿真使用的SAR圖像。
【具體實(shí)施方式】
[0032] 參照圖1，本發(fā)明的實(shí)現(xiàn)步驟如下。
[0033] 步驟1 :將SAR圖像數(shù)據(jù)庫樣本集劃分為訓(xùn)練數(shù)據(jù)集X和測試樣本集y。
[0034] 首先，在包含6類SAR圖像數(shù)據(jù)庫樣本集的每類樣本集中各取大小為256X256的 1000張圖片，然后，再從每一類圖片中隨機(jī)抽取200張構(gòu)成訓(xùn)練集X，剩余作為測試集y。
[0035] 步驟2 :從訓(xùn)練數(shù)據(jù)集X中隨機(jī)抽取m塊尺寸dX d的訓(xùn)練圖像塊，并進(jìn)行全局對 比度歸一化，構(gòu)成訓(xùn)練圖像塊集JTe JTrf。
[0036] 步驟3 :利用訓(xùn)練圖像塊集X訓(xùn)練第一層稀疏字典。
[0037] 3a)將訓(xùn)練圖像塊集X的特征矩陣表示為：
[0038]
[0039] 其中表示字典，N表示每個圖像塊的特征數(shù)量，ε是極小的常數(shù)，F(xiàn) e R?xn 表示特征矩陣。矩陣F第i行的值對應(yīng)第i個圖像塊的特征值，第j列的值表示不同圖像 塊的第j類特征；
[0040] 3b)根據(jù)特征矩陣F求稀疏字典D1:
[0041 ] 常用的字典學(xué)習(xí)方法有稀疏編碼算法，稀疏自編碼算法，稀疏RBM算法，OMP正交 匹配追蹤算法，ICA獨(dú)立成分分析算法，稀疏濾波算法等，本實(shí)例中采用但不局限于稀疏濾 波算法求稀疏字典。即：
[0042] 首先，按照公式/ = //14對特征矩陣F的每一列進(jìn)行歸一化處理，再對每一行 進(jìn)行歸一化處理，得到歸一化后的特征矩陣F2;
[0043]
然后，對歸一化后的特征矩陣F2進(jìn)行稀疏約束，求得第一層稀疏字典：
[0044] 步驟4 :利用第一層的稀疏字典01求訓(xùn)練集X的第一層稀疏特征圖Z。
[0045] 利用稀疏字典求解整幅輸入圖片的稀疏特征圖的常用方法有：隨機(jī)抽取處理合成 法，重疊卷積算法和分片卷積算法，本實(shí)例中采用但不局限于重疊卷積算法，其步驟如下：
[0046] 4a)求解輸入圖片的第i張稀疏特征圖Z1:
[0047] Zt =I Φ：Κ? ：,
[0048] 其中K1G R dxd表示第i個卷積核，i = 〇~N，?表示卷積操作，I e Ruxv表示訓(xùn) 練集X的一張圖片，UX V為圖片尺寸，卷積核K1由稀疏字典D i的第i列/變換得到， rj\ ((u d+l) X (v d+l)
[0049] 4b)利用N個不同的卷積核!^對輸入圖片進(jìn)行卷積操作，得到第一層稀疏特征圖： ^j^NX (u d+l) X (v d+l)
[0050] 步驟5 :對第一層稀疏特征圖進(jìn)行非線性變換。
[0051] 非線性變換包括稀疏特征圖歸一化和池化操作，常用的歸一化方法有局部響應(yīng)歸 一化法和局部對比度歸一化法，常用的池化操作有平均池化，最大池化和隨機(jī)池化，本實(shí)例 中采用但不局限于局部響應(yīng)歸一化法和最大池化。其步驟如下：
[0052] 5a)求解第i張局部響應(yīng)歸一化特征圖B1在（X，y)位置上的值
[0053]
[0054] 其中 表示第i張稀疏特征圖Z1在（X，y)位置上的值，α，β，c分別表示不同數(shù) 值的常量，η表示與第i張稀疏特征圖相鄰的稀疏特征圖數(shù)目；
[0055] 5b)對第i張稀疏特征圖Z1中所有坐標(biāo)上的值進(jìn)行局部響應(yīng)歸一化操作，得到第 i張稀疏特征圖Z1的局部響應(yīng)歸一化特征圖B S
[0056]
[0057] 5c)對N張稀疏特征圖采用5a) _5b)的操作，得到第一層局部響應(yīng)歸一化特征圖： B |~gl gNJ ^ rNX (u d+l) X (V d+l