根據Landset 8計算植被覆蓋度

    之前處理了山東省的Modis合成產品,這次用landset8的原始數據算下植被覆蓋度練練手。

簡要說下步驟:

  1. 數據集介紹
  2. 輻射定標
  3. 大氣校正
  4. 計算NDVI
  5. 導出Excel
  6. 計算植被覆蓋度

1.數據集介紹

先看一下landset的參數和各波段的信息,本想用landset7的數據,衆所周知03年以後壞掉了,看了下各位大神開發的修復程序,還是決定先不搞它了,8和7的區別如下
8和7的對比圖

2.輻射定標

然後預處理唄,輻射定標,大氣校正,輻射定標還是很簡單的
簡單設置個參數

3.大氣校正

但是但是!到了大氣校正,flaash的參數設置還挺多的,而且導入文件的時候還蹦出來個選項卡,粗略查了一下資料問了下老師結論如下:在這裏插入圖片描述
在這裏插入圖片描述
這兩個分別代表:1.從文件中讀取,不同的波段設置不同的參數。
2.對所有波段用一個單一的參數。
老師告訴我根據數據的類型不同,選項不同,比如地理位置和植被類型
在這裏插入圖片描述
一開始還忘了設置K-T變換參數報錯,在Defaults裏面選擇OLI(600…),其餘K-T設置如圖
K-T變換
另外網上說要修改波譜庫的指向,我修改完發現兩次大氣校正結果一樣
修改波譜庫
處理吧,時間還是蠻久的!
完成
第一次校正,發現結果的DN值被擴大了,又是各種找資料,結論是:都變大,算植被覆蓋度的話並不會受影響。
被擴大的DN值

4.計算NDVI

接下來band math,歸一化算NDVI唄,常規操作,另外大家要是有一步計算出來的公式歡迎分享給我,band math是最浪費我時間的操作了,想公式有點困難。。NDVI

5.導出到Excel

根據Excel計算植被覆蓋度,我取5%置信區間
Excel

6.計算植被覆蓋度並查看結果

VFC計算公式

輸入公式計算,基本完成,這一景圖是選的青藏高原一角,時間是18年2月份,所以肯定冰雪覆蓋,水體我改成了白色。
over

發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章