? ? 前面介紹了Landsat8數(shù)據(jù)下載护赊,那么今天介紹一下给赞,在數(shù)據(jù)下載后如何進(jìn)行處理。
一蓬坡、處理流程
? ? 根據(jù)高永剛等研究結(jié)果為更好地保持原始光譜的信息量猿棉,建議先進(jìn)行圖像融合磅叛,后輻射校正處理。同時萨赁,也解決了大氣校正(FLAASH)不能對全色波段進(jìn)行輻射校正的問題弊琴。
? ? 同時,地物分類或者做地圖使用那完全沒必要做大氣校正杖爽;而要通過反射情況提取水體信息那就要做大氣校正敲董。
? ? 因?yàn)長1T文件已經(jīng)經(jīng)過了帶DEM的地形校正,所以坐標(biāo)精度基本能滿足中小比例尺的要求掂林,所以加下來要進(jìn)行圖像如何臣缀、輻射校正和大氣校正。
二泻帮、圖像融合
? ? ENVI 5.3新增了GS融合法精置,將8波段30 m分辨率的多光譜數(shù)據(jù)和15 m的全色數(shù)據(jù)進(jìn)行融合,從而得到15 m分辨率的多光譜數(shù)據(jù)锣杂。
? ? 1脂倦、啟動ENVI 5.3,選擇File——Open元莫,選擇MTL.txt文件打開赖阻;選擇“B8.TIF”文件打開。
? ? 2踱蠢、在工具箱(Toolbox)中火欧,點(diǎn)擊Image Sharpening,雙擊Gram-Schmidt Pan Sharpening打開對話框茎截。
? ? 3苇侵、在對話框中選擇低分辨率的多光譜數(shù)據(jù)文件,點(diǎn)擊OK企锌,再選擇高分辨全色波段數(shù)據(jù)文件榆浓,點(diǎn)擊OK。
? ? 4撕攒、在Pan Sharpening Parameters參數(shù)面板如圖下圖陡鹃,選擇傳感器類型為:Landsat8_oli,重采樣方法選擇Cubic Convolution抖坪,設(shè)置輸出路徑和文件名萍鲸。
? ? 5、輸出結(jié)果: 具體結(jié)果如下圖6擦俐,右側(cè)為原始多光譜猿推,左側(cè)為融合影像,可見分辨率顯著提高。
三蹬叭、輻射校正
? ? 1、在ENVI 5.3的工具箱中查找工具:Radiometric Correction——Radiometric Calibration雙擊状知,選擇要校正的多譜數(shù)據(jù)(而本次試驗(yàn)選擇以上融合的多光譜)秽五,“LC81200332020291LGN00_rh.dat”進(jìn)行輻射定標(biāo)。
? ? 2饥悴、選擇參數(shù)對話框中點(diǎn)擊:Apply Flaash Setting后會自動讀取元數(shù)據(jù)中的正射參數(shù)坦喘,也可以自己輸入?yún)?shù),主要參數(shù)如下圖:
? ? 3西设、輸出結(jié)果: 右側(cè)為多光譜融合影像瓣铣,而左側(cè)為輻射校正后影像,可以看出輻射校正后影像更加明鮮亮贷揽,植被更加突出棠笑。
四、大氣校正
? ? 目前通常采用輸模型MODTRA對Landsat8數(shù)據(jù)進(jìn)行大氣校正禽绪,具體采用ENVI5.3中FLAASH模塊蓖救。
1、選擇Toolbox—Radiometric Correction——Atmospheric Correction Module——FLAASH Atmospheric Correction印屁,打開FLAASH大氣校正工具循捺。參數(shù)設(shè)置如下:
? 2、參數(shù)說明
? 2.1基本參數(shù)設(shè)置(見圖10)
? 2.1.1文件輸入和輸出:單擊Input Radiance Image雄人,選擇以上輻射亮度值數(shù)據(jù)LC81200332020291LGN00_fs.dat从橘。打開Radiance Scale Factors 對話框,選擇Use single scale factor for all bands(Single scale factor:1.000000)础钠。
? ? 單擊 Output Reflectance File 按鈕選擇輸出文件名和路徑恰力。
? ? 2.1.2傳感器與圖像信息:緯度(Lat)、經(jīng)度(Lon)珍坊,F(xiàn)LAASH自動獲任础;傳感器類型(Sensor Type)選擇LandSat-8 OLI阵漏;海拔(Ground Elevation/km)驻民,可以從百度查詢區(qū)域平均海拔;影像日期履怯、時間(Flight Date回还、Flight Time(格林尼治時間)),一般自動獲取叹洲。
? ? 2.1.3大氣模型(Atmospheric Model):Mid-Latitude Summer(5月份緯度:40-50)柠硕。
? ? 2.1.4氣溶膠模型(Aerosol Model):以城市(Urban)、郊外(Rural)為主。
? ? 2.1.5氣溶膠反演(Aerosol Retrieval):2-Band(K-T)蝗柔。
? ? 2.1.6初始能見度(Initial Visibility/km):40闻葵。
? ? 2.2多光譜設(shè)置(Multispectral Settings)
? ? 2.2.1 Defaults 下拉框: Over-Land Retrieval Standard (660:2100)。 Filter Function File :選擇landsat8_oli.sli波譜響應(yīng)文件癣丧。
? ? 2.3高級設(shè)置(Advanced Settings)
? ? 2.3.1Tile Size設(shè)置為100M槽畔,其余按照默認(rèn)設(shè)置。
? ? 3胁编、結(jié)果輸出: 如圖13顯示厢钧,有圖為大氣校正結(jié)果,左圖為輻射校正結(jié)果嬉橙,可見大氣校正圖像使水體更易辨識早直,而對植被分類無影響,所以一般情況下做輻射校正就可以滿足要求市框。
? ? 以上就是遙感數(shù)據(jù)的預(yù)處理過程霞扬,如果有什么問題和建議歡迎留言。
參考文獻(xiàn):
[1] 高永剛, 周文靜, 徐涵秋. 遙感影像融合與輻射校正順序關(guān)系探討[J]. 福州大學(xué)學(xué)報:自然科學(xué)版, 2017, 045(006):833-839.
[2] 初慶偉 ,張洪群 ,吳業(yè)煒 . Landsat-8 衛(wèi)星數(shù)據(jù)應(yīng)用探討. 遙感信息, 2013, 4(28): 110-114.
[3] 袁金國, 牛錚, 王錫平. 基于FLAASH的Hyperion高光譜影像大氣校正. 光譜學(xué)與光譜分析, 2009, 5(29): 1181-1185.
[4] 王磊, 陳一江. 武漢中心 LandSat-8 數(shù)據(jù)處理流程 v0.11, 2015年1月20日.