[發明專利]基于MODIS遙感影像的高分辨率遙感影像植被指數時間序列校正方法

有效

• 申請號: 201310690726.8
• 申請日: 2013-12-13
• 公開（公告）號: CN103617629A
• 公開（公告）日: 2014-03-05
• 發明（設計）人: 劉永學;汪業成;張群;李滿春;王加勝;程亮;張荷霞;孫超;趙賽帥;成王玉;鄒偉
• 申請（專利權）人: 南京大學
• 主分類號: G06T7/00
• 分類號: G06T7/00;G01S7/497;G01S17/88
• 代理公司: 南京同澤專利事務所(特殊普通合伙) 32245
• 代理人: 石敏

1.一種基于MODIS遙感影像的高分辨率遙感影像植被指數時間序列校正方法，包括以下步驟:

1）獲取預定區域、預定時間段的HJ-1A/1B遙感影像數據以及該區域同時間段的MODIS遙感數據；

2）讀取MODIS遙感數據，獲得MODIS遙感影像中每個像元的NDVI值、實際成像日和質量評價值，其中所述質量評價值分為以下幾類：“-1”表示沒有評價結果，“0”表示數據質量非常好，“1”表示數據可用但是質量一般，“2”表示有雪冰覆蓋，“3”表示有云層干擾；

3）分別將所有MODIS遙感影像中具有相同坐標值的像元的NDVI值按實際成像日的時間順序組合成第一MODIS?NDVI時間序列，將所述第一MODIS?NDVI時間序列中質量評價值不等于0和1的像元的NDVI值用該序列中其相鄰的NDVI值的插值進行代替，得到第二MODIS?NDVI時間序列；

4）對每幅HJ-1A/1B遙感影像像元的植被指數進行歸一化處理，得到為HJ?NDVI影像；

5）分別將所有HJ?NDVI影像中具有相同坐標值的像元的NDVI值按實際成像日的時間順序組合成與坐標值對應的第一HJ?NDVI時間序列，并依照第二MODIS?NDVI時間序列的實際成像日對第一HJ?NDVI時間序列按最相鄰的NDVI值進行插值，得到與第二MODIS?NDVI時間序列的實際成像日相一致的第二HJ?NDVI時間序列；

6）對所述預定區域中MODIS遙感影像的每一個像元，找到在空間上被該像元覆蓋的HJ-1A/1B遙感影像中的像元，兩者存在覆蓋關系的第二MODIS?NDVI時間序列和第二HJ?NDVI時間序列構成時間序列對，對所述時間序列對以第二MODIS?NDVI時間序列為標準序列做協整檢驗；

7）通過了協整檢驗的時間序列對的第二MODIS?NDVI時間序列減去相對應的第二HJ?NDVI時間序列得到波動序列；對該波動序列進行異常區段檢驗得到正常波動區段和異常波動區段，并將第二HJ?NDVI時間序列按正常波動區段和異常波動區段的實際成像日劃分為正常HJ?NDVI時間區段和異常HJ?NDVI時間區段；

8）對第二HJ?NDVI時間序列中的異常HJ?NDVI時間區段進行校正，得到所述預定區域、預定時間段內具有較高時間分辨率的高分辨率遙感影像時間序列。

2.根據權利要求1所述的基于MODIS遙感影像的高分辨率遙感影像植被指數時間序列校正方法，其特征在于：在執行第6）步之前，首先將第二MODIS?NDVI時間序列和和第二HJ?NDVI時間序列進行Savitzky-Golay濾波。

3.根據權利要求1所述的基于MODIS遙感影像的高分辨率遙感影像植被指數時間序列校正方法，其特征在于：第8）步中，利用正常HJ?NDVI時間區段與相應時間段的第二MODIS?NDVI時間序列對所述異常HJ?NDVI時間區段進行校正。

4.根據權利要求1所述的基于MODIS遙感影像的高分辨率遙感影像植被指數時間序列校正方法，其特征在于：第4）步中對每幅HJ-1A/1B遙感影像進行歸一化前，先對每幅HJ-1A/1B遙感影像進行輻射校正和幾何精校正。

5.根據權利要求1所述的基于MODIS遙感影像的高分辨率遙感影像植被指數時間序列校正方法，其特征在于，第6）步中對所述時間序列對進行協整檢驗的方法是：

1A）使用ADF檢驗方法判斷對應的所述時間序列對中的第二MODIS?NDVI時間序列和和第二HJ?NDVI時間序列是否為同階單整，若同階單整，則轉到第1B）步；若非同階單整，則所述時間序列對不具備協整關系；

1B）對同階單整的時間序列對進行協整回歸，使用ADF方法檢驗回歸殘差的平穩性，若殘差平穩，則所述時間序列對是協整的，兩者具有相依關系；否則兩者不具備協整關系。