[發明專利]基于遙感數據的礦山和固廢填埋場自動提取方法和裝置

有效

• 申請號: 202110307785.7
• 申請日: 2021-03-23
• 公開（公告）號: CN113033381B
• 公開（公告）日: 2021-09-10
• 發明（設計）人: 余嘉琦;李營;劉玉平;饒瑞璇;殷守敬;王晨;陳輝;王雪蕾;肖如林;賈興;靳川平;孫浩;徐丹;滕佳華
• 申請（專利權）人: 生態環境部衛星環境應用中心
• 主分類號: G06K9/00
• 分類號: G06K9/00;G06K9/20;G06K9/62
• 代理公司: 北京恩赫律師事務所 11469
• 代理人: 劉守憲;李善學

1.一種基于遙感數據的礦山和固廢填埋場自動提取方法，其特征在于，所述方法包括：

S1：獲取待檢測區的亞米級遙感影像和熱紅外遙感影像；

S2：對所述亞米級遙感影像和熱紅外遙感影像進行預處理；

S3：根據歸一化植被指數方法，從預處理后的亞米級遙感影像上提取非植被覆蓋區；

S4：根據K-T變換方法，從預處理后的亞米級遙感影像上提取非建筑用地區；

S5：根據預處理后的熱紅外遙感影像計算地表溫度，并從預處理后的熱紅外遙感影像上提取高溫異常區；

S6：將高溫異常區與非植被覆蓋區、非建筑用地區進行空間相交分析，其交集為礦山和固廢填埋場圖斑；

S7：根據礦山和固廢填埋場圖斑的亞米級遙感影像的光譜特征對所述礦山和固廢填埋場圖斑進行分類，得到礦山區域和固廢填埋場區域。

2.根據權利要求1所述的基于遙感數據的礦山和固廢填埋場自動提取方法，其特征在于，所述S3包括：

S31：基于預處理后的亞米級遙感影像，計算每個像元的歸一化植被指數NDVI；

其中，NDVI為歸一化植被指數，ρ_NIR、ρ_RED分別為預處理后的亞米級遙感影像的近紅外波段、紅波段的遙感反射率；

S32：將NDVI＜0的像元提取為非植被覆蓋區。

3.根據權利要求2所述的基于遙感數據的礦山和固廢填埋場自動提取方法，其特征在于，所述S4包括：

S41：基于預處理后的亞米級遙感影像，根據K-T變換模型Y＝R^Tx+r進行主成分分析，提取每個像元的第一主成分分量Y1和第二主成分分量Y2；

其中，Y為變換后的主成分分量，Y＝{Y1、Y2......Yn}；R為K-T變換系數，R＝{R1、R2......Rn}，R^T為R的轉置，x為預處理后的亞米級遙感影像各波段的灰度值，x＝{x1、x2......xn}，r為常數偏移量，r＝{r1、r2......rn}；

S42：根據每個像元的第一主成分分量Y1和第二主成分分量Y2計算每個像元的比值建筑指數RBI；

RBI＝Y₁/Y₂

S43：將RBI＜-0.05的像元提取為非建筑用地區。

4.根據權利要求3所述的基于遙感數據的礦山和固廢填埋場自動提取方法，其特征在于，所述S5包括：

S51：根據預處理后的熱紅外遙感影像計算每個像元的地表溫度T_s；

其中，C₁、C₂為Plank函數參數，λ為有效波長，B(T_s)＝[L_λ-L_up-τ*(1-ε)*L_down]/(ε*τ)，L_λ為預處理后的熱紅外遙感影像λ波段的輻射亮度，L_up、L_down、τ分別為大氣上行輻射、下行輻射和大氣透過率，ε為比輻射率；

S52：對每個像元的地表溫度T_s通過如下公式進行正態分布擬合，并計算擬合優度p-value，得到p-value∈(0.75，1)時對應的標準差σ；

其中，σ為標準差，μ為T_s的均值；

S53：將T_s＞1.96σ的像元提取為高溫異常區。