之前處理了山東省的Modis合成產品,這次用landset8的原始數據算下植被覆蓋度練練手。
簡要說下步驟:
- 數據集介紹
- 輻射定標
- 大氣校正
- 計算NDVI
- 導出Excel
- 計算植被覆蓋度
1.數據集介紹
先看一下landset的參數和各波段的信息,本想用landset7的數據,衆所周知03年以後壞掉了,看了下各位大神開發的修復程序,還是決定先不搞它了,8和7的區別如下
2.輻射定標
然後預處理唄,輻射定標,大氣校正,輻射定標還是很簡單的
3.大氣校正
但是但是!到了大氣校正,flaash的參數設置還挺多的,而且導入文件的時候還蹦出來個選項卡,粗略查了一下資料問了下老師結論如下:
這兩個分別代表:1.從文件中讀取,不同的波段設置不同的參數。
2.對所有波段用一個單一的參數。
老師告訴我根據數據的類型不同,選項不同,比如地理位置和植被類型
一開始還忘了設置K-T變換參數報錯,在Defaults裏面選擇OLI(600…),其餘K-T設置如圖
另外網上說要修改波譜庫的指向,我修改完發現兩次大氣校正結果一樣
處理吧,時間還是蠻久的!
第一次校正,發現結果的DN值被擴大了,又是各種找資料,結論是:都變大,算植被覆蓋度的話並不會受影響。
4.計算NDVI
接下來band math,歸一化算NDVI唄,常規操作,另外大家要是有一步計算出來的公式歡迎分享給我,band math是最浪費我時間的操作了,想公式有點困難。。
5.導出到Excel
根據Excel計算植被覆蓋度,我取5%置信區間
6.計算植被覆蓋度並查看結果
輸入公式計算,基本完成,這一景圖是選的青藏高原一角,時間是18年2月份,所以肯定冰雪覆蓋,水體我改成了白色。