[發(fā)明專(zhuān)利]一種復(fù)種指數(shù)提取方法及裝置

有效

• 申請(qǐng)?zhí)枺?/td>: 201810096305.5
• 申請(qǐng)日: 2018-01-31
• 公開(kāi)（公告）號(hào): CN108345992B
• 公開(kāi)（公告）日: 2021-07-09
• 發(fā)明（設(shè)計(jì)）人: 王昊宇;趙祥;高濤;杜曉錚
• 申請(qǐng)（專(zhuān)利權(quán)）人: 北京師范大學(xué)
• 主分類(lèi)號(hào): G06Q10/06
• 分類(lèi)號(hào): G06Q10/06;G06Q50/02
• 代理公司: 北京路浩知識(shí)產(chǎn)權(quán)代理有限公司 11002
• 代理人: 王瑩;吳歡燕

1.一種復(fù)種指數(shù)提取方法，其特征在于，包括：

S1，根據(jù)目標(biāo)區(qū)域中農(nóng)作物耕種的起始時(shí)間，獲取從所述起始時(shí)間開(kāi)始預(yù)設(shè)時(shí)長(zhǎng)內(nèi)所述目標(biāo)區(qū)域的第一植被指數(shù)時(shí)間序列曲線(xiàn)；

S2，使用二階差分方法，獲取所述第一植被指數(shù)時(shí)間序列曲線(xiàn)中的波峰，若各所述波峰為單干擾峰，則對(duì)所述單干擾峰進(jìn)行條件平滑；

S3，再次使用所述二階差分方法獲取條件平滑后的所述第一植被指數(shù)時(shí)間序列曲線(xiàn)中的峰值，并統(tǒng)計(jì)條件平滑后的所述第一植被指數(shù)時(shí)間序列曲線(xiàn)中的峰值總數(shù)，將所述峰值總數(shù)作為所述目標(biāo)區(qū)域的復(fù)種指數(shù)；

其中，所述步驟S2中判斷各所述波峰是否為單干擾峰的步驟具體包括：

對(duì)于所述第一植被指數(shù)時(shí)間序列曲線(xiàn)中的每個(gè)波峰，若該波峰不滿(mǎn)足第一預(yù)設(shè)條件或第二預(yù)設(shè)條件，則該波峰為單干擾波峰；

其中，所述第一預(yù)設(shè)條件為該波峰的后一個(gè)波谷與該波峰的前一個(gè)波谷之間的間隔時(shí)長(zhǎng)大于所述目標(biāo)區(qū)域中農(nóng)作物的生長(zhǎng)周期；

所述第二預(yù)設(shè)條件為該波峰對(duì)應(yīng)的前獨(dú)立長(zhǎng)勢(shì)與該波峰對(duì)應(yīng)的植被指數(shù)年際變化之間的比值大于或等于第一預(yù)設(shè)閾值，且該波峰對(duì)應(yīng)的后獨(dú)立長(zhǎng)勢(shì)與該波峰對(duì)應(yīng)的植被指數(shù)年際變化之間的比值大于或等于所述第一預(yù)設(shè)閾值；

該波峰對(duì)應(yīng)的前獨(dú)立長(zhǎng)勢(shì)為該波峰與該波峰的前一個(gè)波谷之間的差值；

該波峰對(duì)應(yīng)的后獨(dú)立長(zhǎng)勢(shì)為該波峰與該波峰的后一個(gè)波谷之間的差值；

該波峰對(duì)應(yīng)的植被指數(shù)年際變化為該波峰所在的年份中最大植被指數(shù)與最小植被指數(shù)之間的差值。

2.根據(jù)權(quán)利要求1所述的方法，其特征在于，所述步驟S1之前還包括：

根據(jù)目標(biāo)區(qū)域的坐標(biāo)，獲取所述坐標(biāo)所屬的預(yù)設(shè)坐標(biāo)范圍和所述預(yù)設(shè)坐標(biāo)范圍對(duì)應(yīng)的農(nóng)作物耕種的起始時(shí)間和農(nóng)作物的生長(zhǎng)周期；其中，所述預(yù)設(shè)坐標(biāo)范圍與所述起始時(shí)間預(yù)先關(guān)聯(lián)存儲(chǔ)，所述預(yù)設(shè)坐標(biāo)范圍與所述生長(zhǎng)周期預(yù)先關(guān)聯(lián)存儲(chǔ)。

3.根據(jù)權(quán)利要求1所述的方法，其特征在于，所述步驟S1之前還包括：

獲取目標(biāo)區(qū)域至少兩年的第二植被指數(shù)時(shí)間序列曲線(xiàn)；其中，所述第二植被指數(shù)時(shí)間序列曲線(xiàn)以小于或等于半個(gè)月為一個(gè)波段；

使用SG濾波算法對(duì)所述第二植被指數(shù)時(shí)間序列曲線(xiàn)進(jìn)行平滑處理；其中，所述SG濾波算法的半波寬小于或等于所述第二植被指數(shù)時(shí)間序列曲線(xiàn)中波段總數(shù)的二分之一；

相應(yīng)地，所述步驟S1具體包括：

根據(jù)目標(biāo)區(qū)域中農(nóng)作物耕種的起始時(shí)間，從平滑處理后的所述第二植被指數(shù)時(shí)間序列曲線(xiàn)中選擇從所述起始時(shí)間開(kāi)始預(yù)設(shè)時(shí)長(zhǎng)內(nèi)的第一植被指數(shù)時(shí)間序列曲線(xiàn)。

4.根據(jù)權(quán)利要求1-3任一所述的方法，其特征在于，所述步驟S2中使用二階差分方法獲取所述第一植被指數(shù)時(shí)間序列曲線(xiàn)中的波峰的計(jì)算公式為：

d_i¹＝NDVI_i-NDVI_i-1；

若d_i¹0，d_i²＝1；

若d_i¹0，d_i²＝-1；

d_i³＝d_i²-d_i-1²；

其中，NDVI_i為任一所述第一植被指數(shù)時(shí)間序列曲線(xiàn)中第i時(shí)刻的植被指數(shù)，NDVI_i-1為第i時(shí)刻前一時(shí)刻的植被指數(shù)，d_i¹為第i時(shí)刻的一階差分，d_i²為第i時(shí)刻的中間屬性，d_i-1²為第i時(shí)刻前一時(shí)刻的中間屬性，d_i³為第i時(shí)刻的二階差分；

若d_i³＝-2，則獲知第i時(shí)刻的植被指數(shù)為波峰；

若d_i³＝2，則獲知第i時(shí)刻的植被指數(shù)為波谷。