ENVI下不變目標法相對大氣校正詳解

2022-10-11 01:15:08 字數 3369 閱讀 9812

相對大氣校正按照數學基礎可以分為2種,非線性校正法和線性校正法。非線性校正法最典型的是直方圖匹配,影象的直方圖是影象中所有灰度值的概率分布。即將校正影象的直方圖與參考影象的直方圖進行匹配,使兩幅影象具有相同或相近的灰度值概率分布,達到兩幅影象上同名地物具有相同灰度值的目的。

線性校正法有個前提假設:不同時相的影象灰度值之間滿足線性關係,這種假設在近似情況下是成立的。這樣就可以通過線性等式來描述不同時相間的灰度關係,用x表示參考影象,y表示待校正影象,他們之間的線性關係可描述為:

y=a*x+b                              (13.5)

其中:a、b為線性等式中的引數,即為增益和偏移量。線性校正的關鍵是確定公式(13.5)中的兩個引數,通常是根據兩時相光譜性質相對穩定地物樣本點的dn值,利用線性回歸的方法求得公式(13.5)中的引數a、b。之後,用公式(13.5)對影象的各波段dn值做線性變換,完成相對大氣校正。

根據前述原理,完成線性相對校正需要以下3個步驟:

第一步,在兩幅影象中搜尋相對固定目標即光譜穩定的地物樣本點,即偽不變特徵要素(pif:pseudo-invariant features);

第二步,運用這些偽不變特徵點的dn值,利用線性回歸的方法求解公式(13.5)中的引數,得到影象間的線性關係;

第三步,根據該關係式,通過波段運算,得到與參考影象具有相同或相近輻射值的結果影象,完成相對大氣校正。整個過程的關鍵是pif的選擇。

下面以兩個不同大氣環境下成像、已經經過精確配準、envi標準格式的landsat  tm5資料為例(2023年和2023年),介紹線性校正法的操作步驟(參見實驗資料光碟\第13章 輻射定標與大氣校正\2-不變目標法相對大氣校正):

注:1、以下操作在envi5.3中完成,envi5.4及以上版本類似

2、可以使用envi擴充套件工具:

基於不變目標法的相對大氣校正完成,本文介紹envi中根據基本原理和流程操作完成,能更好的理解整個處理過程。

第一步、pif選擇

選擇一幅目視質量較好的影象作為基準影象(2023年),另外一幅作為待校正影象(2023年)。在兩個影象上選擇相同區域的瀝青房頂、礫石面、混凝土停機坪、潔淨水體、混凝土、沙地等地物作為pif,這些地物不會隨時間的變化而變化。

(1)           在主介面中,選擇file→open,開啟兩幅影象「2000.dat」和「2001.dat」。

(2)           在主介面中,選擇views→two vertical views,開啟兩個垂直顯示視窗。

(3)           在圖層管理(layer manger)中,選中其中乙個view,在工具欄中單擊 開啟資料管理(data manger),在「2000.dat」右鍵選擇rgb:543,單擊load data按鈕。同樣的方法為另外乙個view載入顯示「2001.dat」。

(4)           在主介面中,選擇views→link views,在link views面板中,單擊new link按鈕,在右邊兩個檢視中單擊滑鼠左鍵,單擊ok按鈕,將兩個影象進行地理鏈結顯示。

(5)           在圖層管理(layer manger)中,在「2000.dat」年圖層中右鍵選擇new region of interest,開啟roi tool面板。

(6)           通過目視方式,從兩幅影象找到光譜穩定、相同地物作為樣本,用polygon或者point型別繪製感興趣區。

(7)           在繪製一定數量感興趣後(不宜太多,太多後面的回歸運算量會很大),在roi tool面板中,選擇file→export→export to csv,在export rois to csv面板中,選擇輸出的roi以及輸出的檔案路徑及檔名(圖13.13),單擊ok按鈕將roi內對應位置和像元值輸出為csv檔案。

圖13.13 export rois to csv面板

(8)           回到主介面中,在圖層管理(layer manger)中選中顯示「2001.dat」的view,在工具欄中單擊 開啟資料管理(data manger),找到前面繪製的roi資料右鍵選擇load按鈕,將roi載入到「2001.dat」圖層上。

(9)           回到主介面中,在圖層管理(layer manger)中的「2001.dat」上的roi圖層上雙擊滑鼠開啟roi tool面板。重複(7)步驟匯出為csv檔案。

分別用microsoft office excel開啟上面步驟得到兩個文字檔案,這樣我們得到了2023年的影象和2023年的影象相對應偽不變特徵要素(pif)的畫素值。從excel中可以看到,兩個時相影象中每乙個波段的畫素值是一一對應關係,剛好對應公式(13.5)中的x和y。第二步中就是利用這些畫素值,根據最小二乘回歸分析法獲得公式(13.5)中的a和b兩個引數。

第二步、線性關係式求解

使用最小二乘回歸的方法來求解線性回歸式,根據最小二乘回歸分析法從公式(13.5)可得到

其中xi、yi分別表示第i個待校正影象(2023年)和參考影象(2023年)對應的pif的像元值。

這裡直接使用microsoft office excel求解線性回歸式,將第一步獲得pif的像元值分別匯入同乙個excel電子**中(同乙個sheet),利用excel電子**計算散點圖功能很容易計算公式(13.5)中的a和b。這樣就得到待校正影象每個波段的線性變換關係式(如表13.7)。

表13.7 回歸解算的a和b值

波段增益(a)

偏移(b)

band1

0.39

90band2

0.489

32.2

band3

0.54

37.4

band4

0.496

46.2

band5

0.63

60.6

band7

0.51

51第三步、線性變換

(2)           如圖13.14所示,根據表13.7中的a和b值分別輸入gain和offset。選擇output data type:byte。

圖13.14 gain and offset values面板

(3)    選擇輸出路徑和檔名。單擊ok按鈕執行運算。

不變目標法相對大氣校正操作過程已經完成。

摘自《envi遙感影象處理方法 第二版》高等教育出版社

ENVI擴充套件工具 PS操作技巧 ENVI設定標頭檔案

envi擴充套件工具 柵格影象批處理工具包 工具在網盤 2 將extensions和custom code兩個資料夾拷貝到envi安裝路徑,覆蓋同名資料夾即可。l envi53 c program files exelis envi53 l envi54 c program files harris...

ENVI紋理計算

假彩色合成 是指多光譜遙感影象彩色合成處理時,如果參與合成的三個波段的波長與對應的紅 綠 藍三種原色的波長不同,那麼合成影象的顏色就不可能是地面景物的真是顏色,這種合成就是假彩色合成。在彩色合成中,若將 tm4,3,2 對應 r,g,b,則稱為標準假彩色成。即將近紅外波段對應紅色,紅光波段對應於綠色...

ENVI 非監督分類

非監督分類的概念 非監督分類,又稱 聚類分析或者點群分析 在多光譜影象中搜尋 定義其自然相似光譜集群的過程。它不必對影象地物獲取先驗知識,僅依靠影象上不同地物光譜資訊進行特徵提取,在統計特徵的差別來達到分類的目的,最後對已分出的各個類別的實際屬性進行確認。在envi中isodata和k means兩...