[發(fā)明專利]一種森林動(dòng)態(tài)變化模式自動(dòng)提取方法

在審

• 申請(qǐng)?zhí)枺?/td>: 201610069963.6
• 申請(qǐng)日: 2016-02-02
• 公開（公告）號(hào): CN105718936A
• 公開（公告）日: 2016-06-29
• 發(fā)明（設(shè)計(jì)）人: 邱炳文;劉哲
• 申請(qǐng)（專利權(quán)）人: 福州大學(xué)
• 主分類號(hào): G06K9/46
• 分類號(hào): G06K9/46;G06K9/62
• 代理公司: 福州元?jiǎng)?chuàng)專利商標(biāo)代理有限公司 35100
• 代理人: 蔡學(xué)俊

1.一種森林動(dòng)態(tài)變化模式自動(dòng)提取方法，其特征在于：包括如下步驟，

S01：逐像元建立研究區(qū)多年逐日增強(qiáng)型植被指數(shù)以及積雪指數(shù)時(shí)序數(shù)據(jù)集；

S02：逐像元逐年分別提取植被覆蓋強(qiáng)度指標(biāo)、離散度指標(biāo)以及持續(xù)度指標(biāo)；

S03：逐像元逐年提取亮度指標(biāo)；

S04：逐像元檢測(cè)植被覆蓋強(qiáng)度指標(biāo)、離散度指標(biāo)、持續(xù)度指標(biāo)以及亮度指標(biāo)的多年變化趨勢(shì)；

S05：建立森林動(dòng)態(tài)變化模式識(shí)別流程圖；

S06：實(shí)現(xiàn)森林變化模式自動(dòng)提取，獲得研究區(qū)森林動(dòng)態(tài)變化分布圖。

2.根據(jù)權(quán)利要求1所述的一種森林動(dòng)態(tài)變化模式自動(dòng)提取方法，其特征在于：在所述步驟S01中，利用MODIS波段反射率數(shù)據(jù)，通過如下方式建立研究區(qū)每個(gè)柵格單元的多年逐日增強(qiáng)型植被指數(shù)MODISEVI以及歸一化差分積雪指數(shù)MODISNDSI時(shí)序數(shù)據(jù)，計(jì)算MODISEVI2以及NDSI時(shí)序數(shù)據(jù)，分別為：

EVI2=2.5*(Nir-Red)(Nir+Red+1),]]>

NDSI=Green-SWIR6Green+SWIR6,]]>

其中，Nir，Red，Green，SWIR6分別為MODIS的近紅外、紅光、綠光以及中紅外6波段的反射率；基于無云日計(jì)算得到MODISEVI以及NDSI時(shí)序數(shù)據(jù)，利用線性插值方法獲得研究區(qū)多年的逐日MODISEVI以及NDSI時(shí)序數(shù)據(jù)，采用Whittakersmoother數(shù)據(jù)平滑方法，逐像元構(gòu)建多年時(shí)空連續(xù)的逐日MODISEVI和NDSI時(shí)序數(shù)據(jù)集。

3.根據(jù)權(quán)利要求1所述的一種森林動(dòng)態(tài)變化模式自動(dòng)提取方法，其特征在于：在所述步驟S02中，逐像元逐年基于植被指數(shù)時(shí)序數(shù)據(jù)的高值區(qū)域的均值，構(gòu)建該像元每年的植被覆蓋強(qiáng)度指標(biāo)；逐像元逐年基于植被指數(shù)時(shí)序數(shù)據(jù)的中高值區(qū)域的極差以及標(biāo)準(zhǔn)差，構(gòu)建該像元每年的植被覆蓋離散度指標(biāo)，其中，DM＝(Max-Quan₂)(SD_m)，DM為植被指數(shù)中高值區(qū)域內(nèi)數(shù)據(jù)的極差與標(biāo)準(zhǔn)差的乘積，Max為對(duì)應(yīng)年份增強(qiáng)型植被指數(shù)時(shí)序數(shù)據(jù)的最大值，Quan₂為對(duì)應(yīng)年份增強(qiáng)型植被指數(shù)時(shí)序數(shù)據(jù)的二分位數(shù)；逐像元逐年基于植被指數(shù)時(shí)序曲線上大于等于預(yù)設(shè)閾值的持續(xù)時(shí)間，構(gòu)建該像元每年的植被覆蓋持續(xù)度指標(biāo)。

4.根據(jù)權(quán)利要求1所述的一種森林動(dòng)態(tài)變化模式自動(dòng)提取方法，其特征在于：在所述步驟S03中，逐像元逐年基于歸一化差分積雪指數(shù)的年均值，構(gòu)建該像元每年的亮度指標(biāo)。