[發明專利]一種基于EVI2季節變化曲線的年GPP估算方法及系統

有效

• 申請號: 201810076612.7
• 申請日: 2018-01-26
• 公開（公告）號: CN108387525B
• 公開（公告）日: 2020-01-21
• 發明（設計）人: 彭代亮
• 申請（專利權）人: 中國科學院遙感與數字地球研究所
• 主分類號: G01N21/17
• 分類號: G01N21/17;G01N21/55
• 代理公司: 11227 北京集佳知識產權代理有限公司
• 代理人: 古利蘭;王寶筠

1.一種基于EVI2季節變化曲線的年GPP估算方法，其特征在于，包括：

獲取全年完整時間序列遙感影像，其中，所述遙感影像中包含近紅外波段和紅光波段；

對所述遙感影像進行定量化處理，獲取時間序列地表反射率數據；

對所述時間序列地表反射率數據進行裁剪，提取出研究區范圍內的地表反射率數據；

對所述研究區范圍內的地表反射率數據中的近紅外波段和紅光波段進行波段運算，逐像素計算出每一期所述研究區遙感影像的EVI2時間序列數據；

對所述EVI2時間序列數據進行濾波處理；

對經過濾波處理后的EVI2時間序列數據逐像元計算曲線積分，得到∑EVI2；

基于植被類型，提取出預先建立的∑EVI2與GPP樣本之間的估算模型；

基于所述得到的∑EVI2和所述估算模型估算出所述研究區范圍內的年GPP。

2.根據權利要求1所述的方法，其特征在于，對所述研究區范圍內的地表反射率數據中的近紅外波段和紅光波段進行波段運算，逐像素計算出每一期所述研究區遙感影像的EVI2時間序列數據包括：

基于公式計算出所述EVI2時間序列數據，其中，ρ_NIR為近紅外波段的反射率，ρ_RED為紅光波段的反射率。

3.根據權利要求1所述的方法，其特征在于，所述對經過濾波處理后的EVI2時間序列數據逐像元計算曲線積分，得到∑EVI2包括：

基于公式計算出∑EVI2，其中，L為整個時間序列的長度，f(x_i)為一年中的第x_i天，N為將濾波后的時間序列曲線下的陰影面積分解成的份數。

4.根據權利要求1所述的方法，其特征在于，所述對所述EVI2時間序列數據進行濾波處理包括：

基于S-G濾波器對所述EVI2時間序列數據進行濾波處理。

5.根據權利要求1所述的方法，其特征在于，所述對所述遙感影像進行定量化處理，獲取時間序列地表反射率數據包括：

對所述遙感影像進行正射校正、幾何校正、輻射定標以及大氣校正，獲取時間序列地表反射率數據。

6.一種基于EVI2季節變化曲線的年GPP估算系統，其特征在于，包括：

獲取模塊，用于獲取全年完整時間序列遙感影像，其中，所述遙感影像中包含近紅外波段和紅光波段；

處理模塊，用于對所述遙感影像進行定量化處理，獲取時間序列地表反射率數據；

裁剪模塊，用于對所述時間序列地表反射率數據進行裁剪，提取出研究區范圍內的地表反射率數據；

計算模塊，用于對所述研究區范圍內的地表反射率數據中的近紅外波段和紅光波段進行波段運算，逐像素計算出每一期所述研究區遙感影像的EVI2時間序列數據；

濾波模塊，用于對所述EVI2時間序列數據進行濾波處理；

積分模塊，用于對經過濾波處理后的EVI2時間序列數據逐像元計算曲線積分，得到∑EVI2；

提取模塊，用于基于植被類型，提取出預先建立的∑EVI2與GPP樣本之間的估算模型；

估算模塊，用于基于所述得到的∑EVI2和所述估算模型估算出所述研究區范圍內的年GPP。

7.根據權利要求6所述的系統，其特征在于，所述計算模塊具體用于：

基于公式計算出所述EVI2時間序列數據，其中，ρ_NIR為近紅外波段的反射率，ρ_RED為紅光波段的反射率。

8.根據權利要求6所述的系統，其特征在于，所述積分模塊具體用于：

基于公式計算出∑EVI2，其中，L為整個時間序列的長度，f(x_i)為一年中的第x_i天，N為將濾波后的時間序列曲線下的陰影面積分解成的份數。

9.根據權利要求6所述的系統，其特征在于，所述濾波模塊具體用于：

基于S-G濾波器對所述EVI2時間序列數據進行濾波處理。

10.根據權利要求9所述的系統，其特征在于，所述處理模塊具體用于：

對所述遙感影像進行正射校正、幾何校正、輻射定標以及大氣校正，獲取時間序列地表反射率數據。