[發明專利]一種火后森林砍伐分布識別的制圖方法及裝置

有效

• 申請號: 202211296896.3
• 申請日: 2022-10-21
• 公開（公告）號: CN115601596B
• 公開（公告）日: 2023-04-07
• 發明（設計）人: 許爾琪;李科為
• 申請（專利權）人: 中國科學院地理科學與資源研究所
• 主分類號: G06V10/764
• 分類號: G06V10/764;G06V20/13;G06V10/28;G06V10/36
• 代理公司: 北京云科知識產權代理事務所(特殊普通合伙) 11483
• 代理人: 張飆

1.一種火后森林砍伐分布識別的制圖方法，其特征在于，包括：

獲取并基于目標區域的遙感影像數據，計算所述遙感影像數據對應的燃燒指數、穗帽變換亮度指數和穗帽變換濕度指數；

設置影像合成的起始日期范圍，獲取多個年份中所述起始日期范圍內的所有目標區域遙感影像，并對所述所有目標區域遙感影像進行影像合成，得到長時間序列年度無云最優影像；

基于預設的燃燒指數閾值和差分燃燒指數閾值，從所述長時間序列年度無云最優影像中提取出未過火森林區域和過火森林區域，對所述未過火森林區域進行像元篩選，得到所述長時間序列年度無云最優影像中的成熟森林像元；

逐年獲取并基于所述成熟森林像元的穗帽變換亮度指數和穗帽變換濕度指數對應的均值及標準差，計算不同年份中所述過火森林區域中每個像元對應的標準穗帽變換亮度指數和標準穗帽變換濕度指數，并根據所述標準穗帽變換亮度指數和所述標準穗帽變換濕度指數，得到所述過火森林區域中每個像元對應的森林擾動指數；

獲取火后未砍伐樣區和火后砍伐樣區，對所述火后未砍伐樣區和所述火后砍伐樣區進行分層隨機抽樣處理，得到火后未砍伐樣區樣本像元和火后砍伐樣區樣本像元，計算并基于每個火后未砍伐樣區樣本像元和每個火后砍伐樣區樣本像元對應的樣本森林擾動指數，設置發生砍伐森林擾動指數閾值；

根據所述森林擾動指數和所述發生砍伐森林擾動指數閾值，逐年對所述過火森林區域中每個像元進行分類，得到不同年份對應的初始火后砍伐森林區域，并生成初始火后森林砍伐時空分布圖像。

2.如權利要求1所述的一種火后森林砍伐分布識別的制圖方法，其特征在于，生成初始火后森林砍伐時空分布圖像后，還包括：

基于圖像差分法計算相鄰年份所有火后未砍伐樣區樣本像元和所述火后砍伐樣區樣本像元的差分森林擾動指數，根據所述差分森林擾動指數，設置差分擾動指數閾值；

基于所述過火森林區域中每個像元對應的森林擾動指數，計算相鄰年份中所述過火森林區域中每個像元對應的差分森林擾動指數；

根據所述發生砍伐森林擾動指數閾值和所述差分森林擾動指數，逐年對所述過火森林區域中每個像元進行分類，得到火后森林初次砍伐區域及火后森林二次砍伐區域，并對應生成火后森林初次砍伐時空分布圖像和火后森林二次砍伐時空分布圖像。

3.如權利要求1所述的一種火后森林砍伐分布識別的制圖方法，其特征在于，得到不同年份對應的初始火后砍伐森林區域后，還包括：

將所述初始火后砍伐森林區域中的每個像元設置為火后砍伐森林像元，獲取當前年份中所有火后砍伐森林像元對應的第一森林擾動指數，同時獲取預設年份中所有火后砍伐森林像元對應的第二森林擾動指數；

判斷所述第一森林擾動指數和所述第二森林擾動指數的擾動差值是否在預設差值范圍內，若否，則認為所述火后砍伐森林像元分類錯誤，并將所述火后砍伐森林像元修正為火后未砍伐森林像元，若是，則認為所述火后砍伐森林像元分類正確。

4.如權利要求1所述的一種火后森林砍伐分布識別的制圖方法，其特征在于，基于預設的燃燒指數閾值和差分燃燒指數閾值，從所述長時間序列年度無云最優影像中提取出未過火森林區域和過火森林區域，具體包括：

獲取所述長時間序列年度無云最優影像中每個像元對應的燃燒指數，將所述燃燒指數與所述燃燒指數閾值進行對比；

若所述燃燒指數大于所述燃燒指數閾值，則將所述像元設置為森林像元，整合所有森林像元，得到森林區域，若所述燃燒指數小于或等于所述燃燒指數閾值，則將所述像元設置為非森林像元，整合所有非森林像元，得到非森林區域；

計算相鄰年份中所述森林區域的每個像元對應的差分燃燒指數，將所述差分燃燒指數與所述差分燃燒指數閾值進行對比；

若所述差分燃燒指數小于或等于所述差分燃燒指數閾值，則將所述森林像元設置為未過火森林像元，整合所有未過火森林像元，得到未過火森林區域，若所述燃燒指數大于或等于所述燃燒指數閾值，則將所述森林像元設置為過火森林像元，整合所有過火森林像元，得到過火森林區域。