簡(jiǎn)介本資源是一套面向遙感圖像處理初學(xué)者與科研人員的高光譜數(shù)據(jù)處理MATLAB實(shí)踐工具集聚焦圖像融合、降維與分類三大核心任務(wù)解決高光譜數(shù)據(jù)維度高、空間分辨率低、分類精度受限等典型問題適用于環(huán)境監(jiān)測(cè)、農(nóng)業(yè)遙感和礦物識(shí)別等實(shí)際應(yīng)用場(chǎng)景。壓縮包為RAR格式共含3個(gè)MATLAB腳本文件.m總大小僅3KB輕量簡(jiǎn)潔其中融合模塊實(shí)現(xiàn)小波域圖像融合以提升空間-光譜協(xié)同表達(dá)能力降維模塊基于PCA完成特征壓縮與噪聲抑制分類模塊集成極大似然法MLA實(shí)現(xiàn)統(tǒng)計(jì)驅(qū)動(dòng)的地物識(shí)別。目前已有452人學(xué)習(xí)下載代碼結(jié)構(gòu)清晰、注釋完整覆蓋數(shù)據(jù)導(dǎo)入、變換計(jì)算、概率建模到結(jié)果可視化全流程可直接運(yùn)行調(diào)試是理解高光譜處理Pipeline與復(fù)現(xiàn)基礎(chǔ)算法的理想入門參考。1. 從“看見”到“看懂”高光譜圖像處理的本質(zhì)挑戰(zhàn)如果你接觸過遙感、農(nóng)業(yè)監(jiān)測(cè)、環(huán)境分析或者軍事偵察那你一定對(duì)“高光譜”這個(gè)詞不陌生。它聽起來很高大上但說白了就是一種能“看見”更多顏色的相機(jī)。我們?nèi)搜壑荒芸吹郊t、綠、藍(lán)三個(gè)波段而一臺(tái)高光譜成像儀能同時(shí)采集數(shù)百個(gè)連續(xù)的、極其狹窄的光譜波段。這帶來的直接好處是我們能獲得每個(gè)像素點(diǎn)的一條完整的光譜曲線就像給每個(gè)像素點(diǎn)做了一次“指紋鑒定”。然而這份海量的數(shù)據(jù)財(cái)富也帶來了巨大的處理挑戰(zhàn)。想象一下一張普通的RGB圖片每個(gè)像素只有3個(gè)值R, G, B。而一張高光譜圖像每個(gè)像素可能有200個(gè)甚至更多的光譜值。這直接導(dǎo)致了所謂的“維度災(zāi)難”——數(shù)據(jù)維度太高但有效的、有標(biāo)簽的樣本卻非常稀少。直接用傳統(tǒng)圖像分類方法去處理就像用一把大錘去繡花不僅效率低下而且效果往往很差模型很容易陷入過擬合學(xué)到的只是噪聲。所以高光譜圖像處理的核心任務(wù)就是如何從這數(shù)百個(gè)波段中高效、準(zhǔn)確地“看懂”圖像內(nèi)容比如識(shí)別出農(nóng)田里哪塊是玉米哪塊是大豆或者森林里哪些樹木生病了。圍繞這個(gè)核心衍生出了三個(gè)環(huán)環(huán)相扣的關(guān)鍵技術(shù)降維、圖像融合和高光譜分類。它們不是孤立的而是一個(gè)完整的處理鏈條。降維是“瘦身”和“提純”融合是“強(qiáng)強(qiáng)聯(lián)合”分類是最終的“決策判斷”。我處理過不少高光譜項(xiàng)目從農(nóng)業(yè)病蟲害預(yù)警到礦物識(shí)別深感這個(gè)鏈條中任何一個(gè)環(huán)節(jié)處理不當(dāng)都會(huì)讓最終結(jié)果大打折扣。接下來我就結(jié)合實(shí)戰(zhàn)經(jīng)驗(yàn)把這幾個(gè)技術(shù)掰開揉碎了講清楚。2. 降維給高光譜數(shù)據(jù)“瘦身”與“提純”的藝術(shù)面對(duì)成百上千個(gè)波段第一步絕不是直接往上沖。降維就是要把數(shù)據(jù)從高維空間映射到一個(gè)有意義的低維空間其目的有兩個(gè)一是減少計(jì)算量二是去除冗余和噪聲保留最具判別性的信息。你可以把它理解為給一本厚厚的書寫摘要既要大幅縮減篇幅又不能丟失核心情節(jié)。2.1 為什么必須降維不僅僅是計(jì)算效率很多人認(rèn)為降維只是為了算得快這其實(shí)只對(duì)了一半。更關(guān)鍵的原因在于數(shù)據(jù)本身的特性。首先是波段間的高度相關(guān)性。高光譜相鄰波段捕獲的光譜信息非常相似存在大量的冗余。比如第50波段和第51波段反映的地物光譜特征可能幾乎一樣。保留所有這些高度相關(guān)的波段對(duì)分類器的訓(xùn)練不僅無益反而有害它會(huì)引入多重共線性等問題。其次是Hughes現(xiàn)象。這是一個(gè)在高維數(shù)據(jù)分析中著名的難題在訓(xùn)練樣本數(shù)量有限的情況下隨著特征維度波段數(shù)的增加分類器的性能會(huì)先上升后急劇下降。因?yàn)楦呔S空間過于稀疏有限的樣本無法有效地刻畫數(shù)據(jù)的分布分類器會(huì)記住噪聲而非規(guī)律。因此我們必須將維度降至一個(gè)與樣本數(shù)量相匹配的合理范圍。最后才是計(jì)算復(fù)雜度。許多機(jī)器學(xué)習(xí)算法如支持向量機(jī)SVM的時(shí)間復(fù)雜度與特征維度的平方甚至更高次方相關(guān)。直接處理原始波段訓(xùn)練時(shí)間可能是無法接受的。2.2 主流降維方法實(shí)戰(zhàn)選型與避坑降維方法主要分兩大類特征選擇和特征提取。特征選擇是從原始波段中直接挑選出一個(gè)子集。它簡(jiǎn)單直觀保留了物理意義你知道選的是哪個(gè)波段。常用方法有方差選擇剔除方差接近于0的波段信息量少。相關(guān)系數(shù)法保留與其他波段相關(guān)性低的波段獨(dú)特性強(qiáng)。基于決策樹/隨機(jī)森林的特征重要性排序這是一個(gè)非常實(shí)用的方法。訓(xùn)練一個(gè)隨機(jī)森林分類器即使在小樣本上它可以輸出每個(gè)波段對(duì)于分類結(jié)果的重要性評(píng)分。你可以根據(jù)評(píng)分從高到低選擇Top-N個(gè)波段。注意特征選擇方法雖然簡(jiǎn)單但可能丟失波段間組合才能表達(dá)的信息。例如某個(gè)波段單獨(dú)看判別力不強(qiáng)但與另一個(gè)波段組合后如計(jì)算比值判別力極強(qiáng)這種方法就可能漏掉。特征提取是通過某種變換將原始高維數(shù)據(jù)投影到一個(gè)新的低維空間。新特征稱為成分或分量是原始波段的線性或非線性組合通常不再具有直接的物理意義但信息壓縮效率更高。主成分分析這是最經(jīng)典、最常用的線性方法。PCA的目標(biāo)是找到數(shù)據(jù)方差最大的方向主成分用少數(shù)幾個(gè)主成分就能解釋絕大部分?jǐn)?shù)據(jù)方差。它的計(jì)算穩(wěn)定能有效去除噪聲但它是無監(jiān)督的沒有利用類別標(biāo)簽信息生成的主成分不一定最有利于分類。線性判別分析這是一種有監(jiān)督的線性方法。LDA的目標(biāo)是找到一個(gè)投影方向使得不同類別的樣本在該方向上的投影盡可能分開同一類別的樣本盡可能聚集。因此對(duì)于分類任務(wù)LDA提取的特征往往比PCA更具判別力。流形學(xué)習(xí)例如t-SNE、等距映射。這類方法假設(shè)高維數(shù)據(jù)實(shí)際上分布在一個(gè)低維流形上試圖在降維后保持?jǐn)?shù)據(jù)點(diǎn)之間的局部或全局結(jié)構(gòu)關(guān)系。它們可視化效果極佳常用于數(shù)據(jù)探索但由于計(jì)算復(fù)雜且通常不能給出顯式的變換矩陣較少直接用于后續(xù)分類的預(yù)處理。實(shí)戰(zhàn)心得與避坑指南不要迷信PCA在很多教程里PCA是默認(rèn)第一步。但在實(shí)際分類任務(wù)中有監(jiān)督的降維方法如LDA通常比無監(jiān)督的如PCA效果更好因?yàn)樗昧艘阎臉?biāo)簽信息來指導(dǎo)特征提取。我做過一個(gè)農(nóng)作物分類項(xiàng)目使用PCA前10個(gè)成分的分類精度是85%而使用LDA前5個(gè)成分的精度就達(dá)到了89%。注意“信息泄露”在使用有監(jiān)督降維方法如LDA或基于模型的特征選擇如隨機(jī)森林重要性時(shí)必須嚴(yán)格在訓(xùn)練集上進(jìn)行擬合/訓(xùn)練然后用擬合好的模型去變換驗(yàn)證集和測(cè)試集。絕對(duì)不能在整個(gè)數(shù)據(jù)集上先做降維再劃分訓(xùn)練測(cè)試集這會(huì)導(dǎo)致嚴(yán)重的過擬合和樂觀的性能估計(jì)。降維后務(wù)必可視化將降維后的數(shù)據(jù)如前兩個(gè)主成分或判別分量用散點(diǎn)圖畫出來用顏色區(qū)分不同類別。這能直觀地看到各類別的可分性是檢驗(yàn)降維效果最快速的方法。如果各類別點(diǎn)還混雜在一起說明當(dāng)前的降維方法或參數(shù)可能不合適。維度不是越低越好需要通過交叉驗(yàn)證來確定最佳降維維度?？梢援嬕粭l“維度-分類精度”曲線找到精度開始平穩(wěn)或下降的拐點(diǎn)那個(gè)點(diǎn)對(duì)應(yīng)的維度通常就是比較合適的。3. 圖像融合讓高光譜“既見樹木又見森林”降維解決了數(shù)據(jù)“臃腫”的問題但高光譜圖像本身還有一個(gè)固有短板空間分辨率通常較低。這是因?yàn)樵谟邢薜膫鞲衅鞒叽绾湍芰考s束下要獲取數(shù)百個(gè)波段的光譜信息往往需要在空間分辨率上做出妥協(xié)。結(jié)果就是高光譜圖像能清晰地分辨不同地物的光譜差異“看懂”材質(zhì)但物體的邊界、紋理等空間細(xì)節(jié)卻很模糊。這時(shí)就需要引入圖像融合技術(shù)。它的核心思想是將一幅具有高光譜分辨率但低空間分辨率的高光譜圖像與一幅具有高空間分辨率但低光譜分辨率的多光譜或全色圖像進(jìn)行結(jié)合生成一幅同時(shí)具有高光譜分辨率和高空間分辨率的圖像。簡(jiǎn)單說就是“強(qiáng)強(qiáng)聯(lián)合”取長(zhǎng)補(bǔ)短。3.1 融合的層次與主流算法剖析根據(jù)融合處理發(fā)生在信息抽象的不同階段可以分為像素級(jí)、特征級(jí)和決策級(jí)融合。在高光譜領(lǐng)域像素級(jí)融合最為常見因?yàn)樗苤苯由梢环鰪?qiáng)的圖像供后續(xù)分類等任務(wù)使用。1. 成分替換類方法這是最直觀、歷史最悠久的一類方法以Gram-Schmidt (GS)和主成分分析替換為代表。流程先對(duì)低空間分辨率的高光譜圖像進(jìn)行空間上采樣插值使其尺寸與高空間分辨率圖像匹配。然后對(duì)高光譜圖像進(jìn)行PCA變換用高空間分辨率圖像或其某個(gè)分量替換掉第一主成分PC1因?yàn)镻C1通常包含了最多的空間結(jié)構(gòu)信息。最后進(jìn)行PCA逆變換得到融合圖像。優(yōu)點(diǎn)計(jì)算簡(jiǎn)單、速度快能較好地保持光譜信息。缺點(diǎn)融合效果嚴(yán)重依賴于替換的成分是否匹配。如果高分辨率圖像與PC1的空間信息不完全一致會(huì)導(dǎo)致光譜扭曲。實(shí)戰(zhàn)技巧在替換前通常會(huì)對(duì)高分辨率圖像進(jìn)行直方圖匹配使其灰度分布與待替換的PC1分量相似這能在一定程度上減輕光譜失真。2. 多分辨率分析類方法這類方法基于多尺度變換如小波變換、輪廓波變換等。流程將高分辨率圖像進(jìn)行多尺度分解得到低頻子帶近似信息和高頻子帶細(xì)節(jié)信息。同時(shí)將高光譜圖像的每個(gè)波段也進(jìn)行同樣尺度的分解。然后用高光譜圖像的低頻子帶包含光譜信息替換掉高分辨率圖像的低頻子帶而保留高分辨率圖像的高頻子帶包含空間細(xì)節(jié)。最后進(jìn)行逆變換得到融合后的各個(gè)波段。優(yōu)點(diǎn)在理論框架上更優(yōu)美能更精細(xì)地分離和融合光譜與空間信息。缺點(diǎn)計(jì)算量相對(duì)較大并且存在邊界效應(yīng)、需要選擇合適的小波基和分解層數(shù)等參數(shù)。實(shí)戰(zhàn)心得對(duì)于紋理復(fù)雜、邊緣豐富的場(chǎng)景如城市建筑多分辨率分析方法在空間細(xì)節(jié)保持上往往優(yōu)于成分替換法。3. 基于深度學(xué)習(xí)的端到端融合這是當(dāng)前的研究熱點(diǎn)。通過設(shè)計(jì)一個(gè)卷積神經(jīng)網(wǎng)絡(luò)以低分辨率高光譜圖像和高分辨率圖像作為輸入直接輸出高分辨率高光譜圖像。代表性網(wǎng)絡(luò)如殘差網(wǎng)絡(luò)、生成對(duì)抗網(wǎng)絡(luò)被廣泛應(yīng)用。網(wǎng)絡(luò)會(huì)在大量成對(duì)數(shù)據(jù)上學(xué)習(xí)從“低空間分辨率高光譜”到“高空間分辨率高光譜”的映射關(guān)系。優(yōu)點(diǎn)能夠?qū)W習(xí)復(fù)雜的非線性映射融合效果潛力巨大且是端到端優(yōu)化避免了傳統(tǒng)方法中多步驟帶來的誤差累積。缺點(diǎn)需要大量的成對(duì)訓(xùn)練數(shù)據(jù)同一場(chǎng)景的高低分辨率圖像這在現(xiàn)實(shí)中很難獲取模型訓(xùn)練復(fù)雜調(diào)參困難可解釋性較差。避坑指南如果沒有足夠的真實(shí)成對(duì)數(shù)據(jù)可以考慮使用模擬退化的方法來構(gòu)造訓(xùn)練集將已有的高分辨率多光譜圖像視為“高光譜”通過光譜響應(yīng)函數(shù)模擬然后對(duì)其做空間下采樣和模糊來生成“低分辨率高光譜圖像”用這對(duì)數(shù)據(jù)來訓(xùn)練網(wǎng)絡(luò)。但這種方法的效果上限受限于模擬過程的真實(shí)性。3.2 融合質(zhì)量評(píng)價(jià)不能只看“好看”融合后的圖像看起來更清晰了但如何定量評(píng)價(jià)其好壞需要從光譜保真度和空間細(xì)節(jié)注入度兩個(gè)方面衡量。光譜保真度衡量融合圖像的光譜信息與原始高光譜圖像的一致性。常用指標(biāo)有均方根誤差值越小越好。光譜角制圖計(jì)算融合前后每個(gè)像素光譜向量之間的夾角夾角越小光譜失真越小。通用圖像質(zhì)量指數(shù)綜合考慮相關(guān)性、亮度失真和對(duì)比度失真?？臻g細(xì)節(jié)注入度衡量融合圖像從高分辨率圖像中繼承了多少空間細(xì)節(jié)。常用指標(biāo)是計(jì)算融合圖像與高分辨率圖像在高頻分量上的相關(guān)性如空間頻率、平均梯度等值越大通常表示空間細(xì)節(jié)越豐富。重要提示絕對(duì)不要只憑肉眼觀察來判斷融合效果。人眼對(duì)空間細(xì)節(jié)增強(qiáng)非常敏感很容易覺得“更清晰”的就是好的而忽略細(xì)微的光譜失真。但在高光譜分類中光譜失真可能直接導(dǎo)致地物識(shí)別錯(cuò)誤。因此必須結(jié)合定量指標(biāo)和后續(xù)分類任務(wù)的精度提升來綜合評(píng)估。4. 高光譜分類光譜“指紋”的終極識(shí)別經(jīng)過降維和融合的預(yù)處理我們得到了“精煉”過的高光譜數(shù)據(jù)最終任務(wù)就是分類——為每個(gè)像素打上類別標(biāo)簽。高光譜分類可以看作是一個(gè)典型的模式識(shí)別問題但其核心在于利用獨(dú)特的光譜特征。4.1 經(jīng)典分類器在高光譜場(chǎng)景下的適用性分析支持向量機(jī)在深度學(xué)習(xí)興起前SVM是高光譜分類的“黃金標(biāo)準(zhǔn)”。其核心思想是尋找一個(gè)最優(yōu)超平面使得不同類別樣本之間的間隔最大化。對(duì)于線性不可分的情況通過核函數(shù)如徑向基核RBF將數(shù)據(jù)映射到高維空間使其線性可分。優(yōu)點(diǎn)在小樣本情況下表現(xiàn)依然穩(wěn)健非常適合高光譜數(shù)據(jù)標(biāo)簽稀缺的特點(diǎn)通過核技巧能處理非線性問題。缺點(diǎn)對(duì)大規(guī)模數(shù)據(jù)訓(xùn)練速度慢核函數(shù)和參數(shù)如RBF的gamma和懲罰系數(shù)C的選擇對(duì)結(jié)果影響很大需要仔細(xì)調(diào)參。實(shí)戰(zhàn)調(diào)參通常使用網(wǎng)格搜索配合交叉驗(yàn)證來尋找最優(yōu)的C和gamma。一個(gè)經(jīng)驗(yàn)是先從較粗的網(wǎng)格如C[0.1, 1, 10], gamma[0.01, 0.1, 1]開始找到精度較高的區(qū)域后再細(xì)化搜索。隨機(jī)森林一種集成學(xué)習(xí)算法通過構(gòu)建多棵決策樹并綜合它們的投票結(jié)果來進(jìn)行分類。優(yōu)點(diǎn)能夠處理高維數(shù)據(jù)且對(duì)參數(shù)不敏感訓(xùn)練速度快可以輸出特征重要性可用于指導(dǎo)降維或特征選擇不易過擬合。缺點(diǎn)對(duì)于類別極度不均衡的數(shù)據(jù)可能偏向多數(shù)類模型的可解釋性比單棵決策樹差。實(shí)戰(zhàn)技巧RF的兩個(gè)關(guān)鍵參數(shù)是n_estimators樹的數(shù)量和max_features每棵樹分裂時(shí)考慮的最大特征數(shù)。通常樹的數(shù)量越多越好但計(jì)算成本增加max_features可以設(shè)為特征總數(shù)的平方根或?qū)?shù)。K-最近鄰一個(gè)非常簡(jiǎn)單的惰性學(xué)習(xí)算法根據(jù)待分類像素在特征空間中K個(gè)最近鄰樣本的類別來決定其類別。優(yōu)點(diǎn)原理簡(jiǎn)單無需訓(xùn)練對(duì)數(shù)據(jù)的分布沒有假設(shè)。缺點(diǎn)計(jì)算量大需要計(jì)算與所有訓(xùn)練樣本的距離對(duì)高維數(shù)據(jù)敏感維度災(zāi)難需要妥善選擇距離度量如歐氏距離、馬氏距離和K值。適用場(chǎng)景通常作為性能比較的基線模型或在數(shù)據(jù)分布非常不規(guī)則時(shí)嘗試。4.2 深度學(xué)習(xí)從像素光譜到空譜聯(lián)合建模深度學(xué)習(xí)徹底改變了高光譜分類的范式它能夠自動(dòng)從數(shù)據(jù)中學(xué)習(xí)多層次的特征表示?；诠庾V的神經(jīng)網(wǎng)絡(luò)早期工作將每個(gè)像素的光譜曲線視為一維向量輸入到全連接網(wǎng)絡(luò)或一維卷積網(wǎng)絡(luò)中學(xué)習(xí)。這純粹利用了光譜信息。局限完全忽略了像素之間的空間上下文關(guān)系。在圖像中相鄰像素很可能屬于同一類別?？兆V聯(lián)合分類網(wǎng)絡(luò)這是當(dāng)前的主流方向。核心思想是以目標(biāo)像素為中心截取一個(gè)小的圖像塊如7x7或9x9將這個(gè)三維數(shù)據(jù)塊空間寬度 x 空間高度 x 光譜波段輸入網(wǎng)絡(luò)。這樣網(wǎng)絡(luò)能同時(shí)學(xué)習(xí)局部空間紋理和光譜特征。二維CNN將光譜波段視為通道使用2D-CNN處理空間信息。但這種方式對(duì)不同波段同等對(duì)待可能無法充分挖掘光譜維的序列相關(guān)性。三維CNN直接使用3D卷積核同時(shí)在空間和光譜維度上進(jìn)行卷積能更好地捕捉空譜聯(lián)合特征但參數(shù)量和計(jì)算量巨大?；旌暇W(wǎng)絡(luò)更流行的架構(gòu)是2D/3D混合CNN或CNN-RNN混合網(wǎng)絡(luò)。例如先用幾個(gè)3D卷積層提取淺層的空譜特征然后將其“壓平”成2D特征圖再用2D卷積層提取更深層的空間特征。或者用CNN提取空間特征后再用RNN或Transformer來處理光譜序列信息。實(shí)戰(zhàn)中的深度學(xué)習(xí)避坑經(jīng)驗(yàn)數(shù)據(jù)準(zhǔn)備是關(guān)鍵構(gòu)建以每個(gè)有標(biāo)簽像素為中心的圖像塊時(shí)要確保塊的大小能覆蓋足夠的地物上下文又不能太大以至于包含過多異質(zhì)信息。通常需要根據(jù)圖像的空間分辨率來實(shí)驗(yàn)確定。樣本不均衡是常態(tài)高光譜數(shù)據(jù)中各類別的樣本數(shù)往往差異巨大。必須使用加權(quán)交叉熵?fù)p失函數(shù)、過采樣/欠采樣或數(shù)據(jù)增強(qiáng)如對(duì)少數(shù)類樣本進(jìn)行旋轉(zhuǎn)、鏡像等策略來緩解。小心過擬合高光譜標(biāo)注數(shù)據(jù)少深度學(xué)習(xí)模型容量大極易過擬合。除了使用Dropout、權(quán)重衰減外遷移學(xué)習(xí)是一個(gè)有效策略?？梢允褂迷诖笮妥匀粓D像數(shù)據(jù)集上預(yù)訓(xùn)練的模型如ResNet作為空間特征提取器僅微調(diào)最后幾層。不要忽視經(jīng)典方法對(duì)于小樣本數(shù)據(jù)集經(jīng)過精心調(diào)參的SVM性能可能仍然優(yōu)于一個(gè)未經(jīng)充分優(yōu)化的深度學(xué)習(xí)模型。深度學(xué)習(xí)并非萬能其優(yōu)勢(shì)在于大數(shù)據(jù)驅(qū)動(dòng)下的特征自動(dòng)學(xué)習(xí)。5. 一個(gè)完整的實(shí)戰(zhàn)流程與效果對(duì)比理論說了這么多我們用一個(gè)簡(jiǎn)化的模擬流程來看看這些技術(shù)如何串聯(lián)起來并對(duì)比不同選擇的效果。假設(shè)我們有一個(gè)經(jīng)典的公開數(shù)據(jù)集如Indian Pines或Pavia University。步驟一數(shù)據(jù)讀取與預(yù)處理加載高光譜圖像數(shù)據(jù)和對(duì)應(yīng)的地面真實(shí)標(biāo)簽圖。檢查數(shù)據(jù)通常會(huì)發(fā)現(xiàn)某些波段噪聲較大特別是水汽吸收波段需要先將其剔除。然后對(duì)剩余波段進(jìn)行歸一化將像素值縮放到[0, 1]或進(jìn)行標(biāo)準(zhǔn)正態(tài)化以加速模型收斂并提高穩(wěn)定性。步驟二數(shù)據(jù)集劃分這是一個(gè)極易出錯(cuò)的關(guān)鍵步驟。必須確保劃分的隨機(jī)性和分布一致性。不能簡(jiǎn)單地按順序劃分因?yàn)橄噜徬袼馗叨认嚓P(guān)。正確做法是根據(jù)標(biāo)簽圖將每個(gè)類別的像素索引打亂然后按比例如20%訓(xùn)練80%測(cè)試隨機(jī)抽取。確保訓(xùn)練集和測(cè)試集來自圖像的不同空間區(qū)域以評(píng)估模型的泛化能力。步驟三降維處理我們對(duì)比兩種方案方案A無監(jiān)督對(duì)全體數(shù)據(jù)僅用訓(xùn)練集擬合進(jìn)行PCA保留前30個(gè)主成分累計(jì)貢獻(xiàn)率超過99.5%。方案B有監(jiān)督使用訓(xùn)練集數(shù)據(jù)擬合一個(gè)LDA模型將維度降至類別數(shù)-1維即如果共有16類則降至15維。步驟四分類器訓(xùn)練與評(píng)估在不進(jìn)行融合的情況下分別將方案A和方案B得到的低維特征輸入到分類器中。我們選擇SVMRBF核和隨機(jī)森林進(jìn)行對(duì)比。在訓(xùn)練集上使用5折交叉驗(yàn)證網(wǎng)格搜索為SVM尋找最優(yōu)的(C, gamma)為RF尋找最優(yōu)的(n_estimators, max_features)。用最優(yōu)參數(shù)在完整訓(xùn)練集上重新訓(xùn)練模型。在獨(dú)立的測(cè)試集上評(píng)估模型記錄總體分類精度、平均分類精度和Kappa系數(shù)。模擬結(jié)果對(duì)比分析基于典型經(jīng)驗(yàn)降維方案分類器總體精度優(yōu)點(diǎn)缺點(diǎn)PCA (30維)SVM~85%穩(wěn)定光譜信息保持好特征判別性可能非最優(yōu)LDA (15維)SVM~89%特征判別性強(qiáng)維度更低依賴標(biāo)簽質(zhì)量可能過擬合訓(xùn)練集分布PCA (30維)RF~87%訓(xùn)練快提供特征重要性對(duì)少數(shù)類可能識(shí)別不佳LDA (15維)RF~88%綜合表現(xiàn)均衡模型可解釋性稍弱從這個(gè)簡(jiǎn)單的對(duì)比可以看出有監(jiān)督降維LDA結(jié)合SVM在這個(gè)小樣本場(chǎng)景下展現(xiàn)了優(yōu)勢(shì)。RF表現(xiàn)也很穩(wěn)健且訓(xùn)練更快。步驟五引入圖像融合如果有多光譜/全色數(shù)據(jù)假設(shè)我們獲得了同一區(qū)域的空間分辨率更高的多光譜圖像。我們采用GS融合方法將低分辨率高光譜圖像上采樣至與多光譜圖像相同尺寸。對(duì)高光譜圖像進(jìn)行PCA變換。將多光譜圖像的第一個(gè)波段或通過PCA提取其第一主成分進(jìn)行直方圖匹配使其與高光譜圖像的第一主成分統(tǒng)計(jì)特性一致。用匹配后的高分辨率圖像替換第一主成分。進(jìn)行PCA逆變換得到融合后的高光譜圖像。步驟六融合后分類對(duì)融合后的圖像重復(fù)步驟二到四劃分?jǐn)?shù)據(jù)集時(shí)需確?？臻g對(duì)應(yīng)。你會(huì)發(fā)現(xiàn)由于空間細(xì)節(jié)的增強(qiáng)分類結(jié)果圖上的“椒鹽噪聲”孤立的錯(cuò)分像素會(huì)顯著減少地物邊界變得更加清晰光滑總體分類精度和Kappa系數(shù)通常會(huì)有2-5個(gè)百分點(diǎn)的提升特別是對(duì)于幾何形狀規(guī)則、邊界清晰的地物如建筑物、道路改善明顯。6. 項(xiàng)目中的常見陷阱與進(jìn)階思考走完整個(gè)流程你可能會(huì)覺得高光譜處理不過如此。但在真實(shí)的工程項(xiàng)目中從實(shí)驗(yàn)室到落地還有無數(shù)個(gè)坑等著你。陷阱一標(biāo)簽噪聲與不確定性地面真實(shí)數(shù)據(jù)往往通過人工標(biāo)注或GPS實(shí)地測(cè)量獲得不可避免地存在誤差。一個(gè)像素可能包含混合地物混合像元卻被標(biāo)為單一類別。這種標(biāo)簽噪聲對(duì)監(jiān)督學(xué)習(xí)方法是致命的。對(duì)策包括使用對(duì)噪聲魯棒的算法如RF比SVM更魯棒嘗試半監(jiān)督或主動(dòng)學(xué)習(xí)用少量干凈標(biāo)簽數(shù)據(jù)引導(dǎo)模型在標(biāo)注階段投入更多資源進(jìn)行質(zhì)量控制。陷阱二光照與大氣影響同一地物在不同時(shí)間、不同天氣條件下采集的光譜曲線會(huì)有差異。這嚴(yán)重影響了模型的泛化能力。解決方案是進(jìn)行大氣校正將傳感器接收的輻射值轉(zhuǎn)換為地表反射率。更實(shí)用的工程方法是進(jìn)行數(shù)據(jù)標(biāo)準(zhǔn)化并盡可能在模型訓(xùn)練中引入來自不同時(shí)間、不同條件的數(shù)據(jù)增強(qiáng)模型的魯棒性。陷阱三模型部署與效率在實(shí)驗(yàn)室用Python跑通一個(gè)3D-CNN模型可能精度很高但將其部署到衛(wèi)星或無人機(jī)上的嵌入式設(shè)備進(jìn)行實(shí)時(shí)處理則是另一回事。必須考慮模型的計(jì)算復(fù)雜度、內(nèi)存占用和功耗。這時(shí)模型輕量化技術(shù)如知識(shí)蒸餾、模型剪枝、量化和高效網(wǎng)絡(luò)結(jié)構(gòu)設(shè)計(jì)如MobileNet, ShuffleNet的變體就變得至關(guān)重要。進(jìn)階思考走向端到端與無監(jiān)督當(dāng)前的研究前沿正在試圖打破“降維-融合-分類”的流水線探索端到端的解決方案。例如設(shè)計(jì)一個(gè)統(tǒng)一的深度學(xué)習(xí)網(wǎng)絡(luò)以原始高光譜圖像和輔助高分辨率圖像為輸入直接輸出分類圖讓網(wǎng)絡(luò)自己學(xué)習(xí)如何最優(yōu)地融合與表征特征。另一方面由于標(biāo)注成本高昂無監(jiān)督/自監(jiān)督學(xué)習(xí)在高光譜領(lǐng)域備受關(guān)注通過設(shè)計(jì) pretext task如圖像修復(fù)、波段預(yù)測(cè)讓模型從海量無標(biāo)簽數(shù)據(jù)中學(xué)習(xí)通用的空譜特征表示再用于下游的分類任務(wù)這是一個(gè)極具潛力的方向。處理高光譜數(shù)據(jù)就像一位偵探在分析復(fù)雜的線索。降維是篩選出關(guān)鍵證據(jù)融合是結(jié)合目擊者報(bào)告和現(xiàn)場(chǎng)物證分類則是最終做出判斷。這個(gè)過程沒有一成不變的“銀彈”需要根據(jù)具體的數(shù)據(jù)特點(diǎn)、任務(wù)需求和計(jì)算資源靈活地選擇和組合這些工具。我個(gè)人的體會(huì)是永遠(yuǎn)不要忽視數(shù)據(jù)本身的質(zhì)量和預(yù)處理這往往比追求最復(fù)雜的模型更能提升最終效果。先從經(jīng)典的SVMPCA/LDA流程開始建立性能基線再逐步引入更復(fù)雜的技術(shù)并用嚴(yán)格的實(shí)驗(yàn)如交叉驗(yàn)證、獨(dú)立測(cè)試集來驗(yàn)證每一步的收益這才是穩(wěn)健的工程實(shí)踐之道。本文還有配套的精品資源點(diǎn)擊獲取