ENVI5對landsat8 OLI數(shù)據(jù)輻射定標(biāo)與大氣校正

1、ENVI5.1 landsat 8OLI數(shù)據(jù)輻射定標(biāo)與大氣校正ENVI 5.1 FLAASH大氣校正工具中新增了一些多光譜傳感器類型，如：Landsat 8 OLI、Pleiades-1A、Pleiades-1B、Rasat、Spot6。支持多光譜傳感器類型明細(xì)見下圖1：圖1.1            ENVI 5.0與ENVI 5.1 FLAASH中支持傳感器類型對比1   打開Landsat 8數(shù)據(jù)在這里以Landsat 8數(shù)據(jù)“LC8123

2、0322013132LGN02”為例進行介紹。File -> Open，選擇LC81230322013132LGN02_MTL.txt文件，點擊OK打開。分為五個數(shù)據(jù)集：多光譜數(shù)據(jù)（1-7波段），全色波段數(shù)據(jù)（8波段），卷云波段數(shù)據(jù)（9波段），熱紅外數(shù)據(jù)（10，11波段）和質(zhì)量波段數(shù)據(jù)（12波段）。 圖2.1            Data Manager面板2   輻射定標(biāo)ENVI 5.1 工具箱中查找工具：/Radi

3、ometric Correction/Radiometric Calibration；雙擊此工具，選擇要校正的多光譜數(shù)據(jù)“LC81230322013132LGN02_MTL_MultiSpectral”進行輻射定標(biāo)。1)   選擇定標(biāo)的文件；圖2.2            打開需要定標(biāo)的多光譜數(shù)據(jù)2)   定標(biāo)參數(shù)設(shè)置FLAASH大氣校正對于radiance數(shù)據(jù)的要求是：BIL存儲格式，單位轉(zhuǎn)換“Scale Fact

4、or”的設(shè)置，單擊Apply FLAASH Settings得到相應(yīng)的參數(shù)。具體的參數(shù)設(shè)置如下圖2.3；（定標(biāo)類型：輻射亮度值或大氣表觀反射率，F(xiàn)LAASH大氣校正需要輸入BIL格式的輻射亮度值，故此出選擇radiance，BIL）圖2.3            輻射定標(biāo)參數(shù)設(shè)置3)   查看定標(biāo)結(jié)果通過定標(biāo)以后，查看影像的DN值是不是輻射定標(biāo)后的結(jié)果，可以采用 查看定標(biāo)后的DN值，查詢結(jié)果見下圖2.4； 圖2.4 

5、           輻射定標(biāo)結(jié)果查看3   FLAASH大氣校正ENVI 5.1 工具箱中查找工具：/Radiometric Correction/Atmospheric Correction Module/FLAASH Atmospheric Correction，雙擊此工具，打開輻射定標(biāo)的數(shù)據(jù)，進行相關(guān)的參數(shù)設(shè)置進行大氣校正。1)   Input Radiance Image：打開輻射定標(biāo)結(jié)果數(shù)據(jù)；2)

6、0; 設(shè)置輸出反射率的路徑；3)  設(shè)置輸出FLAASH校正文件的路徑，最優(yōu)狀態(tài)：路徑所在磁盤空間足夠大；4)  中心點經(jīng)緯度Scene Center Location：自動獲取；5)   選擇傳感器類型：Landsat-8 OLI；其對應(yīng)的傳感器高度以及影像數(shù)據(jù)的分辨率自動讀??；6)  設(shè)置研究區(qū)域的地面高程數(shù)據(jù)，網(wǎng)上查詢研究區(qū)域平均海拔；7)   影像生成時的飛行過境時間：在layer manager中的Lc8數(shù)據(jù)圖層右鍵選擇View Metadata，瀏覽ti

7、me字段獲取成像時間；注：也可以從元文件“LC81230322013132LGN02_MTL.txt”中找到，具體名稱：DATE_ACQUIRED = 2013-05-12；SCENE_CENTER_TIME = 02:55:26.6336980Z；8)  大氣模型參數(shù)選擇：Sub-Arctic Summer（根據(jù)成像時間和緯度信息選擇）；9)   氣溶膠模型Aerosol Model：Urban，氣溶膠反演方法Aerosol Retrieval：2-band（K-T）；10)  其他參數(shù)按照默認(rèn)設(shè)置即可。 圖2.

8、5            FLAASH校正主界面參數(shù)設(shè)置11)  多光譜參數(shù)設(shè)置中，K-T反演選擇默認(rèn)模式：Defaults->Over-Land Retrieval standard（600:2100）波譜響應(yīng)函數(shù)：默認(rèn)指向.Program Files-Exelis-ENVI51-classic-filt_func-landsat8_oli.sli該波譜相應(yīng)函數(shù)包含了卷云波段和全色波段，且卷云波段位于全色7個波段的中間位置，如果利用該波譜相應(yīng)函數(shù)，

9、則會導(dǎo)致用錯波段，SWIR1波段大氣校正后的結(jié)果全為0。因此，需要重新制作與大氣校正輸入的多光譜數(shù)據(jù)的7波段對應(yīng)的波譜相應(yīng)函數(shù)。步驟如下：（1）           啟動ENVI classic，選擇Window->Start New Plot Window。（2）           在ENVI Plot Window窗口中，選擇File->Input Data

10、>Spectral Library，打開ldcm_oli.sli波譜響應(yīng)文件。（3）           如下圖所示選擇7個波段的波譜響應(yīng)函數(shù)，選擇File->Save plot as-> Spectral Library，按照默認(rèn)參數(shù)保存為.sli文件（文件名：landsat8_mul7band_spectrfun.sli）。 圖1：選擇7個波段的波譜響應(yīng)文件將波譜相應(yīng)函數(shù)filter function file重新指向剛才導(dǎo)出的波譜相應(yīng)函數(shù)。見下圖2.6

11、；，點擊OK； 圖2.6            多光譜參數(shù)設(shè)置12)     高級參數(shù)設(shè)置：根據(jù)內(nèi)存大小設(shè)置Tile Size（Mb）：100(8g物理內(nèi)存)，其他參數(shù)默認(rèn)即可，詳細(xì)見下圖2.7； 圖2.7            高級參數(shù)設(shè)置13)  點擊“Apply”運行FLAASH校正。14)  運行完后能得到反演的能見度和水汽柱圖2.8            能見度和水汽柱反演結(jié)果4   查看校正結(jié)果