[發明專利]一種基于孔雀石多波段邏輯運算模型的找礦方法及系統

在審

• 申請號: 202110506744.0
• 申請日: 2021-05-10
• 公開（公告）號: CN113239779A
• 公開（公告）日: 2021-08-10
• 發明（設計）人: 何海洋;李士杰;王博;唐衛東;高永寶;魏立勇;李云昊;劉蓓;朱一龍
• 申請（專利權）人: 中國人民武裝警察部隊黃金第五支隊
• 主分類號: G06K9/00
• 分類號: G06K9/00;G06K9/34;G06K9/62
• 代理公司: 北京維正專利代理有限公司 11508
• 代理人: 李傳亮

1.一種基于孔雀石多波段邏輯運算模型的找礦方法，其特征在于，所述方法包括：

獲取待識別區域的高分辨率光學衛星圖像數據及銅化探異常；

對所述高分辨率光學衛星圖像數據進行預處理，得到全波段高分辨率光學衛星圖像數據；

構建孔雀石多波段邏輯運算模型，利用所述全波段高分辨率光學衛星圖像數據提取得到孔雀石礦化異常；

根據所述銅化探異常及所述孔雀石礦化異常進行疊加，得到銅礦找礦靶區；

其中，所述構建孔雀石多波段邏輯運算模型，利用所述全波段高分辨率光學衛星圖像數據提取得到孔雀石礦化異常，包括：

根據ASTER標準波普數據庫確定孔雀石對應全波段WorldView-3衛星圖像數據中16個波段的波譜曲線，對比分析孔雀石波譜曲線特征，確定所述全波段WorldView-3衛星圖像數據中具有明顯吸收谷特征的可見光-近紅外第7波段VNIR7及短波紅外第7波段SWIR7；

從所述16個波段中選擇VNIR2及SWIR3，構建第一波段比值的表達式(VNIR2+SWIR3)*0.5/VNIR7，所述第一波段比值用于表征所述孔雀石在0.78μm的吸收谷特征；

從所述16個波段中選擇SWIR5，構建第二波段比值的表達式SWIR5/SWIR7，所述第二波段比值用于表征所述孔雀石在2.27μm的吸收谷特征；

對所述全波段WorldView-3衛星圖像數據進行所述第一波段比值和所述第二波段比值處理，得到圖像灰度值；

對所述圖像灰度值進行統計分析，分別得到所述第一波段比值的第一平均值和第一標準方差及所述第二波段比值的第二平均值和第二標準方差；

根據所述第一平均值和所述第一標準方差進行異常閾值分割，確定所述第一波段比值對應的第一最優測試閾值，根據所述第二平均值和所述第二標準方差進行異常閾值分割，確定所述第二波段比值對應的第二最優測試閾值；

根據所述第一波段比值的表達式、所述第一最優測試閾值、所述第二波段比值的表達式及所述第二最優測試閾值，建立孔雀石多波段邏輯運算模型，所述孔雀石多波段邏輯運算模型的表達式為：

((float(VNIR2+SWIR3)*0.5/(VNIR7))ge1.08）and?((float(SWIR5)/(SWIR7))ge1.10)

其中，所述ge用于表征大于或等于，所述第一最優測試閾值的取值為1.08，所述第二最優測試閾值的取值為1.10；

根據所述孔雀石多波段邏輯運算模型，對所述WorldView-3衛星圖像數據進行計算處理，得到所述待識別區域的孔雀石礦化異常。

2.根據權利要求1所述的基于孔雀石多波段邏輯運算模型的找礦方法，其特征在于，所述獲取待識別區域的高分辨率光學衛星圖像數據及銅化探異常，包括：

讀取高分辨率光學衛星WorldView-3對待識別區域進行掃描的衛星圖像數據，得到WorldView-3衛星圖像數據，所述WorldView-3衛星圖像數據包括8個可見光-近紅外波段和8個短波紅外波段；

根據已知地質資料獲取所述待識別區域的銅化探異常。

3.根據權利要求2所述的基于孔雀石多波段邏輯運算模型的找礦方法，其特征在于，所述對所述高分辨率光學衛星圖像數據進行預處理，得到全波段高分辨率光學衛星圖像數據，包括：

通過遙感圖像處理軟件的定標工具對所述WorldView-3衛星圖像數據進行輻射校正；

通過所述遙感圖像處理軟件的大氣校正功能對輻射校正后的所述WorldView-3衛星圖像數據進行大氣校正；

通過所述遙感圖像處理軟件的正射校正功能對大氣校正后的所述WorldView-3衛星圖像數據進行正射校正；

通過所述遙感圖像處理軟件的波段合成工具，對正射校正后的所述WorldView-3衛星圖像數據的8個可見光-近紅外波段和8個短波紅外波段進行波段合成，得到所述全波段WorldView-3衛星圖像數據。