[發(fā)明專利]一種區(qū)域生態(tài)建設見效時間的空間分異檢測方法

有效

• 申請?zhí)枺?/td>: 202010250696.9
• 申請日: 2020-04-01
• 公開（公告）號: CN111460386B
• 公開（公告）日: 2022-02-11
• 發(fā)明（設計）人: 劉彥隨;曹智;李裕瑞;陳玉福;王永生
• 申請（專利權）人: 中國科學院地理科學與資源研究所
• 主分類號: G06F17/18
• 分類號: G06F17/18
• 代理公司: 北京律譜知識產(chǎn)權代理有限公司 11457
• 代理人: 黃云鐸

1.一種區(qū)域生態(tài)建設見效時間的空間分異檢測方法，其特征在于，所述見效時間的空間分異檢測方法包括以下步驟：

步驟S1，獲取待檢測區(qū)域生態(tài)建設工程執(zhí)行前后數(shù)年的植被指數(shù)年度最大值合成數(shù)據(jù)集；

步驟S2，將所述待檢測區(qū)域劃分為有限數(shù)量的宏觀單元，生成各宏觀單元的植被指數(shù)/年份圖，并通過各所述宏觀單元中生成的植被指數(shù)擬合線變化率變化確定待檢測區(qū)域擬合線變化率變化閾值；

步驟S3，依據(jù)各所述宏觀單元的植被指數(shù)/年份圖設定時間段閾值；

步驟S4，依據(jù)設定的植被指數(shù)擬合線變化率變化閾值和時間段閾值，判斷在所述宏觀單元內(nèi)基本單元的植被指數(shù)/年份圖中是否具有明顯的轉折時間點；

步驟S5，依據(jù)所述轉折時間點及其兩側所述擬合線變化率變化，判斷該時間點是否為生態(tài)建設的明顯見效時間，并分別從明顯見效時間計算生態(tài)恢復速率；

在步驟S2中，還包括以下步驟：

步驟S21，將待檢測區(qū)域依據(jù)行政區(qū)或自然區(qū)劃分為宏觀單元，計算合成數(shù)據(jù)集中各宏觀單元范圍內(nèi)每一年份像素單元植被指數(shù)值的平均值；

步驟S22，依據(jù)待檢測區(qū)域植被指數(shù)的年份，生成各宏觀單元的所述植被指數(shù)/年份圖；

步驟S23，依據(jù)生態(tài)建設工程執(zhí)行期設置轉折點，將所述植被指數(shù)/年份圖中曲線分為2段，生成每段所述曲線的線性回歸擬合線；

步驟S24，依據(jù)兩段所述擬合線的變化率，求取所述擬合線變化率變化，并設定變化率變化閾值。

2.根據(jù)權利要求1所述見效時間的空間分異檢測方法，其特征在于，在步驟S24中，所述擬合線變化率變化為兩段所述擬合線的變化率求差；選取所述擬合線變化率變化中最小值，并設定該最小值的50％-70％的數(shù)值為所述變化率變化閾值。

3.根據(jù)權利要求1所述見效時間的空間分異檢測方法，其特征在于，在步驟S4中，還包括以下步驟：

步驟S41，依次設定各基本單元的植被指數(shù)/年份圖中具有任意位置的2個轉折時間點，生成所有情況的3段擬合線，并從中選取估計誤差平方之和最小的一組數(shù)據(jù)；判斷該組數(shù)據(jù)中，由2個轉折時間點構成的3個時間段中每個時間段是否都大于所述時間段閾值，如果是，執(zhí)行步驟S42，如果否，執(zhí)行步驟S43；

步驟S42，判斷起始年份至第一轉折時間點和第一轉折時間點至第二轉折時間點、以及第一轉折時間點至第二轉折時間點和第二轉折時間點至結束年份形成的植被指數(shù)擬合線變化率變化是否都大于所述變化率變化閾值，如果是，確定第一轉折時間點和第二轉折時間點是該基本單元的植被指數(shù)/年份圖的轉折時間點，并執(zhí)行步驟S5，如果否，執(zhí)行步驟S43；

步驟S43，設定該基本單元的植被指數(shù)/年份圖中具有任意位置的2個轉折時間點，生成所有情況的3段擬合線，并從中選取估計誤差平方與擬合優(yōu)度比之和最小的一組數(shù)據(jù)；判斷該組數(shù)據(jù)中，由2個轉折時間點構成的3個時間段中每個時間段是否都大于所述時間段閾值，如果否，執(zhí)行步驟S44，如果是，執(zhí)行步驟S42，步驟S42執(zhí)行為否時，執(zhí)行步驟S44；

步驟S44，排除該基本單元的植被指數(shù)/年份圖中具有2個轉折時間點，設定該基本單元的植被指數(shù)/年份圖中具有任意位置的1個轉折時間點，生成所有情況的2段擬合線，并從中選取估計誤差平方之和最小的一組數(shù)據(jù)；判斷該組數(shù)據(jù)中，由1個轉折時間點構成的2個時間段中每個時間段是否都大于所述時間段閾值，如果是，執(zhí)行步驟S45，如果否，執(zhí)行步驟S46；

步驟S45，判斷該轉折時間點兩側線性函數(shù)形成的擬合線變化率變化是否大于所述變化率變化閾值，如果是，確定該轉折時間點是該基本單元的植被指數(shù)/年份圖的轉折時間點，并執(zhí)行步驟S5，如果否，執(zhí)行步驟S46；

步驟S46，該基本單元的植被指數(shù)/年份圖中沒有明顯的轉折時間。

4.根據(jù)權利要求1所述見效時間的空間分異檢測方法，其特征在于，在步驟S4中，所述基本單元是所述宏觀單元內(nèi)一像素點所表示區(qū)域。