1. 程式人生 > >高光譜遙感影象相關知識梳理大全

高光譜遙感影象相關知識梳理大全

前言

​ 本資料整理了高光譜遙感影象概念定義、分析處理與分類識別的基本知識。第一部分介紹高光譜影象的一般性原理和知識,第二部分介紹了高光譜影象的噪聲問題;第三部分介紹高光譜影象資料冗餘問題以及資料降維解決冗餘的方法;第四部分介紹高光譜影象的混合像元問題,對光譜解混做了一定介紹;第五部分和第六部分分別介紹了高光譜影象的監督分類和分監督分類的特點、流程和常用演算法。

1.基本介紹

高光譜遙感(Hyperspectral remote sensing) 是將成像技術和光譜技術相結合的多維資訊獲取技術,同時探測目標的二維集合空間與一維光譜資訊,獲取高光譜解析度的連續、窄波段影象資料。

高光譜影象高解析度影象多光譜影象不同。

高光譜識別優勢:

  • 光譜解析度高、波段眾多,能夠獲取地物幾乎連續的光譜特徵曲線,並可以根據需要選擇或提取特定波段來突出目標特徵;

  • 同一空間解析度下,光譜覆蓋範圍更寬,能夠探測到地物更多對電磁波的響應特徵;

  • 波段多,為波段之間的相互校正提供了便利;

  • 定量化的連續光譜曲線資料為地物光譜機理模型引入影象分類提供了條件;

  • 包含豐富的輻射、空間和光譜資訊,是多種資訊的綜合載體。

高光譜在識別方面的困難:

  1. 資料量大,影象包含幾十個到上百個波段,資料量是單波段遙感影象的幾百倍;資料存在大量冗餘,處理不當,反而會影響分類精度;
  2. 對高光譜影象的分類一方面要求更高的光譜定標和反射率轉換的精度,另一方面又因為成像機理複雜,資料量巨大而導致對影象資料預處理困難,包括大氣矯正、幾何校正、光譜定標和反射率轉換等;
  3. 波段多、波段間的相關性高,因此分類需要的訓練樣本數目大大增多,往往因訓練樣本不足導致得到的訓練引數不可靠(維數災難);
  4. 針對常規遙感的處理模型和方法不能滿足高光譜影象分類的需要。主要問題之一是統計學分類模型的引數估計問題,其對光譜特徵的選擇要求很高。

高光譜影象分類中的Hughes 現象:

Hughes現象:在機器學習問題中,需要在高維特徵空間(每個特徵都能夠取一系列可能值)的有限資料樣本中學習一種“自然狀態”(可能是無窮分佈),要求有相當數量的訓練資料含有一些樣本組合。給定固定數量的訓練樣本,其預測能力隨著維度的增加而減小。

​ 在高光譜遙感影象中,當訓練樣本數目有限時,分類精度隨著影象波段數目的增加先增加,在到達一定極值後,分類精度隨這波段數目的增加而下降。

​ 傳統遙感影象分析是利用影象空間資訊,高光譜影象分析的核心是光譜分析。高光譜的遙感資料是一個光譜影象立方體,其最主要的特點將影象空間維與光譜維資訊合為一體,與單波段相比,多出了一維光譜資訊。在獲取地表空間影象同時,會得到每個像元對應的地物光譜資訊。

單波段遙感影像圖

像元光譜特徵圖

上圖分別是某波段的遙感影像圖和某像元的光譜特徵圖。

2.高光譜資料噪聲處理

​ 高光譜影象精細光譜資料能夠充分反映地物光譜的細微特徵。依據不同地物在光譜特徵上的相對差異就可實現地物類別區分,實現目標探測和精細分類。

​ 但是,由於成像光譜儀波段通道很密而造成光成像能量不足,故相對全色圖而言,高光譜影象的信噪比(SNR)提高比較困難。在影象資料的獲取過程中,地物光譜特徵在噪聲的影響下容易產生“失真”。另外由於高光譜資料量大,在精細分類過程中,往往需要對資料進行降維處理,而在降維過程中需要最大限度保留訊號和壓縮噪聲,所以精確的噪聲評估很有必要。另外噪聲對精細分類的結果也有最直接影響。所以需要對高光譜資料進行噪聲評估。

​ 主要的噪聲認為是高斯白噪聲,分為加性噪聲和乘性噪聲。

​ 噪聲評估主要有三種方法:實驗室法、暗電流法和影象法。由於前兩種方法在實驗中難以實現,廣泛應用的是第三種——影象法

影象法分為以下幾種:

  • 均勻區域法(HA)

    主要思想:從影象中選取四個以上均勻區域,通過計算這些均勻區域標準差的平均值獲取影象噪聲的估計值。

    不足:需要人工進行均勻區域選擇,無法自動化;滿足條件的均勻區域在大部分遙感圖中並不存在;子區域噪聲估計並不能代表整幅影象的噪聲。

  • 地學統計法(GS)

    主要思想:從影象中選擇幾條均勻的窄條帶,通過對這些窄條帶的半方差函式的計算實現對影象噪聲的估算。一定程度上利用了成像光譜資料的空間相關性。

    不足:與均勻區域法相似。

  • 區域性均值與區域性標準差法(LMLSD)

    主要思想:假定影象由大量均勻的小塊構成,且噪聲以加性噪聲為主。使用區域性均值和區域性標準差的概念,將影象分割為很多小塊,然後計算這些子塊的標準差作為區域性噪聲大小,並選擇包含子塊數最多區間的區域性標準差的平均值作為整個影象的最佳噪聲估計。

    不足:只對高斯白噪聲有效,對高斯隨機噪聲的影象,訊號被噪聲干擾。

  • 空間/光譜維去相關法(SSDC)

    主要思想:是一種專門針對高光譜影象的噪聲評估方法,利用高光譜影象空間維和光譜維存在高相關性的特點,通過多元線性迴歸去除具有高相關性的訊號,利用得到的殘差影象對噪聲進行估算。

    評價:該方法受地物覆蓋型別影響小,並且可以自動執行,是目前較為穩定的高光譜影象噪聲評價方法。

3.高光譜影象資料降維的常用方法

​ 高光譜遙感影象所具有的大量光譜波段為地物資訊提取提供了極其豐富的資訊,有利於更精細的地物分類,然而波段的增多也會導致資訊的冗餘和資料處理複雜性的提高。

資料降維滿足下面條件:儘可能保留資料的特徵資訊;去除資料冗餘與相關性。

主要從特徵選擇特徵提取兩方面進行降維處理。

3.1 特徵選擇

特徵選擇是光譜組合,即從原光譜波段數為N的波段中選擇M個波段(N>M),可能的光譜組合數為

N!M!(NM)!
在這裡可以採用Band Index 方法進行光譜波段選擇降維。

*Band Index 方法:*Hyperspectral 遙感影象根據相關性分為K組(如shortwave light、visible light、near-infrared),設每個組中的波段數為 nl(l=1,2,,k)。用pij 表示波段i 與波段j的相關係數 ,σi 表示波段i的方差,Ra表示波段i與不同組其他波段相關係數的絕對值的和,Rw 表示波段i與同組其他波段的相關係數的絕對值的均值,則有

Rw=1nl|pij|
從而波段i的Band Index 可表示為:
Pi=σiRa+Rw
通過觀察,我們知道隨著方差增大,波段包含的資訊越多;隨著相關係數降低,波段的獨立性越高。

Band Index 是一個重要引數,其反映了波段總體包含特徵資訊和相關性。

​ 通過結合Band Index目標物體有效的光譜範圍(effiective spectral scope of object),我們可以進行波段選擇,進而下一步識別分類。

3.2 特徵提取

​ 高光譜的資料降維技術是以影象特徵提取為目的,利用低維資料來有效地表達高維資料的特徵,同時也壓縮了資料量,更有利於資訊的快速提取。資料降維包含的內容非常廣泛,高光譜遙感影象主要以降低光譜維度提取光譜維度特徵為主。
​ 上一節的特徵選擇就講到了在原始特徵空間進行特徵選擇形成原始空間的一個子空間的特徵選擇方法,接下來介紹線性變換方法Y=BX ,從高維資料空間中,產生一個合適的低維子空間(不是簡單的特徵選擇組合),使資料在這個空間的分佈可以在某種最優意義上描述原來的資料。

3.2.1 主成分分析

主成分分析(PCA)是最基本的高光譜資料降維方法,在高光譜資料壓縮、去相關、消噪和特徵提取中發揮了巨大的作用。PCA 變換又稱為霍特林變換(hotelling transform)K-L (karhunen-loeve)變換。變換後的各主成分分量彼此不相關,且隨主成分編號的增加該分量包含的資訊量減少。

在高光譜遙感資料的PCA變換中,一般將每個波段當作一個向量來處理,設高光譜遙感資料有p個波段,影象空間維度為m×n,則具體處理流程:

  1. 影象向量化:輸入影象資料可以表示成X=(x1,x2,,xp)T,其中 xi 表示為一個 N×1 列向量,這裡有N=m×n。即將影象按行或按列展開有規則連線起來,稱為一個向量。
  2. 向量中心化:將向量組中的所有向量減去向量組的均值向量,即 Y=XE(X)
  3. 計算向量組 Y 的協方差矩陣 Σ
  4. 求協方差矩陣 Σ的特徵值矩陣 Λ 和特徵向量矩陣 A
  5. 進行主成分變換, Z=ATY

PCA變換是基於資訊量的一種正交線性變換,變換後的影象資訊主要集中在前幾個主成分分量中,在變換域中丟棄資訊量小的成分分量,經過反變換後仍能得到復原影象的近似影象。

在PCA變換的基礎上提出了分塊主成分分析方法定向主成分分析方法(DPCA)選擇主成分方法(selective PCA)

但PCA 變換存在兩個明顯的缺陷:一是影象數值變換影響明顯;二是變換後的信噪比並不一定隨著主成分編號的增加而降低。針對這兩個問題分別發展了標準化的PCA(standardized PCA, SPCA)殘差調整的PCA(residual-scaled PCA, RPCA)

3.2.2 最大噪聲分數變換

​ 當噪聲方差或噪聲在影象各波段分佈不均勻時,基於方差最大化的PCA方法並不能保證影象質量隨著主成分的增大而降低。所以這裡引入最大噪聲分數(maximum noise fraction,MNF) 變換 ,該變換根據影象質量排列成分。MNF 方法主要採用SNR和噪聲比例來描述影象質量引數。

假設X=[x1,x2,,xp]Tp×N 矩陣,行向量組的均值向量 E(X)=0,協方差矩陣 D(X)=Σ ,假設

X=S+N
其中SN分別指影象中的訊號和噪聲,且兩者不相關。

ΣSΣN 分別為SN的協方差矩陣。這裡假設噪聲為加性噪聲,則噪聲比例可表示為:

Var{N}/Var{X}
MNF 變換是一種線性變換,則有
Zi=aTiX,i=1,,p

相關推薦

光譜遙感影象相關知識梳理大全

前言 ​ 本資料整理了高光譜遙感影象概念定義、分析處理與分類識別的基本知識。第一部分介紹高光譜影象的一般性原理和知識,第二部分介紹了高光譜影象的噪聲問題;第三部分介紹高光譜影象資料冗餘問題以及資料降維解決冗餘的方法;第四部分介紹高光譜影象的混合像元問題,對

光譜遙感影象處理(14)-----ENVI使用教程之監督分類與非監督分類

     監督分類:又稱訓練分類法,用被確認類別的樣本像元去識別其他未知類別像元的過程。它是在分類之前通過目視判讀和野外調查,對遙感影象上某些樣區中影像地物的類別屬性有了先驗知識,對每一種類別選取一定

基於深度卷積神經網路的光譜遙感影象分類---PCA+2D-CNN(偽空譜特徵)

論文地址基於深度卷積神經網路的高光譜遙感影象分類西華大學學報同時利用高光譜影象的光譜資訊和空間資訊的深度卷積神經網路分類模型。基於深度學習到的深度特徵,用邏輯迴歸分類器進行分類訓練。高光譜影象面臨的挑戰

光譜遙感影象空譜聯合分類方法研究--基於光譜,空間,空譜的DBN+LR

論文地址高光譜遙感影象空譜聯合分類方法研究遼寧工程技術大學基於DBN的特徵提取方法和影象分類架構用於高光譜資料分析,將譜域-空域特徵提取和分類器相結合提高了分類精度。1.引言針對高光譜資料量大,特徵維數

光譜遙感影象處理(12)----ENVI使用教程之波段比與波段運算

一、波段比          波段比又稱比值增強(ratio enhancement)。是計算同一幅影象不同波段間相應像元的像元值之比的影象增強方法。          比值增強的主要作用有:可增強地

zip相關知識梳理(一)

zip相關知識梳理(一) 經過對zip檔案的長時間研究,對zip檔案進行相關知識進行梳理,雖然網上很多牛人對其做了相關基礎解析,但是對於特殊情況沒有進行說明,比如超過4G的zip檔案該以什麼格式進行編譯,現本人對其進行詳細講解。如有不足,歡迎各位大神指點! 一、zip檔案格式 zip檔案格式: [本地檔

Web應用伺服器 相關知識梳理(二)Tomcat的系統架構

    Tomcat非常模組化,其總體架構本人歸納為以下:              話說 有這樣的一家人 被稱為Tomcat:       (一)Server           父親被稱為Server,其權利可謂至高無上,控制家中每對孩子;作為每對夫妻的生存環境

MySQL資料庫相關知識梳理

MySQL基本知識梳理 int STDCALLmysql_query(MYSQL*mysql,constchar*q):連線上資料庫後,將連線上資料庫額控制代碼和執行的SQL語句傳入到函式中即可執行該語句。 MYSQL_RES*res; res = mysql_stor

Web應用伺服器 相關知識梳理(四)Tomcat的其他元件——Session

         由於HTTP是一種無狀態協議,當用戶的一次訪問請求結束後,後端伺服器就無法知道下一次來訪問的還是不是上次訪問的使用者;提到Session與Cookie,二者的作用都是為了保持訪問使用者與後端伺服器的互動狀態,解決了HTTP是一種無狀態協議的弊端。 Cook

Python 處理遙感影象光譜輻射定標、大氣校正和計算反射率

唔,最近在做作業的時候,一些實驗內容涉及到了用ENVI處理遙感影象,然後自己手動操作軟體一遍遍的輸入各種引數神馬的感覺挺無聊。。。。然後決定自己用python裡面的opencv庫寫個指令碼批處理影象反射率的計算試試~ 核心步驟就是 遙感影像光譜輻射定標 →大氣校正→計算反射率

資料庫及相關知識詳解大全

友情提示以下內容是對資料庫以及相關知識所做的一個闡述,內容較多需多花費些各位看官寶貴時間,不過全是乾貨以及一些小的例項,我相信各位點進來的看官都會有所收穫。 資料庫介紹 什麼是資料庫? 資料庫(Database)是按照資料結構來組織、儲存和管理資料的倉庫, 每個資料庫都有一個或多個不

影象超解析度及相關知識

                                  影象超解析度及相關知識   初學者可能往往會

數字影象處理筆記——影象處理綜述與相關知識介紹

影象獲取方式 伽馬光 伽馬光成像原理是給人體注入同位素,然後用檢測器來檢測同位素的放射量來成像 X光 X光成像原理是利用人體不同部位密度不同,X光透過性不同,使得成像的灰度值不同 CT CT成像是一組X光不停旋轉以後得到不同角度的成像,從而可以實現

電商之梳理rabbitmq相關知識---通訊--訊息佇列

rabbitmq 編輯 MQ全稱為Message Queue, 訊息佇列(MQ)是一種應用程式對應用程式的通訊方法。應用程式通過讀寫出入佇列的訊息(針對應用程式的資料)來通訊,而無需專用連線來連結它們。訊息傳遞指的是程式之間通過在訊息中傳送資料進行通訊,而不是

MySQL 基礎知識梳理學習(三)----InnoDB日誌相關的幾個要點

完整 就是 相等 獲取 通過 over 插入 redo 完整性 1.REDO LOG 是用來做數據庫crash recovery的,是保證數據庫數據正確性、完整性及持久性的重要機制之一。在數據庫操作中,它保存了對InnoDB表中數據的修改記錄,所以也叫日誌文件。 2.日誌文

js中與元素寬位置相關知識彙總

常見clientWidth、clientHeight、offsetWidth、offsetLeft,clientX、scrollTop等詞語,比較混亂,現在總結下他們的區別。 1. clientWidth:元素的寬度(width+padding),offsetWidth相比clientWidth,

電商之梳理mongodb相關知識---分散式資料庫

mongodb 《mongodb》是2013年中國鐵道出版社出版的圖書,作者是鄒貴金。該書介紹了MongoDB基礎、應用、管理、效能優化、資料庫的架構,環境搭建例項,程式設計例項等內容。 書 名 mongodb 作 者 鄒貴金 出版社中國鐵道出

電商之梳理Hbase相關知識---分散式資料庫

HBase 本詞條由“科普中國”百科科學詞條編寫與應用工作專案 稽核 。 HBase是一個分散式的、面向列的開源資料庫,該技術來源於 Fay Chang 所撰寫的Google論文“Bigtable:一個結構化資料的分散式儲存系統”。就像Bigtable利用了

光譜影象基於MATLAB的PCA降維

工具:matlab,本人使用的是2016a 使用資料集用公開的Salinas資料集為例,Salinas資料集為512*217*204,即它有204個波段,我們要把它從204維降至3維。 matlab中內含了進行PCA降維的函式,但這個函式輸進去的資料要是二維

電商之梳理UML相關知識-------建模使用

統一建模語言 編輯 同義詞 UML(統一建模語言)一般指統一建模語言 本詞條由“科普中國”百科科學詞條編寫與應用工作專案 稽核 。 Unified Modeling Language (UML)又稱統一建模語言或標準建模語言,是始於1997年一個OMG標準