[發明專利]一種適合于熱紅外高光譜遙感的石英含量定量計算方法

在審

• 申請號: 201610958025.1
• 申請日: 2016-10-27
• 公開（公告）號: CN107991244A
• 公開（公告）日: 2018-05-04
• 發明（設計）人: 郭幫杰;張杰林;武鼎;王俊虎;周覓
• 申請（專利權）人: 核工業北京地質研究院
• 主分類號: G01N21/25
• 分類號: G01N21/25;G01J5/00
• 代理公司: 核工業專利中心11007
• 代理人: 高尚梅

1.一種適合于熱紅外高光譜遙感的石英礦物含量定量計算方法，其特征在于：具體包括以下步驟：

步驟一、對研究區的熱紅外高光譜遙感數據預處理；

步驟二、熱紅外高光譜遙感數據大氣校正，具體包括以下子步驟：

步驟2.1熱紅外高光譜遙感數據進行大氣透過率曲線、上行輻射和下行輻射模擬，得到研究區大氣透過率模擬曲線及上下行輻射模擬曲線；

步驟2.2按照步驟一所采用的熱紅外高光譜遙感數據波段進行重采樣，得到重采樣后大氣透過率曲線，進而得到熱紅外高光譜遙感數據大氣校正圖像；

步驟三、熱紅外高光譜數據溫度/發射率分離，利用歸一化法進行溫度/發射率分離，得到發射率圖像；

步驟四、石英礦物含量定量計算經驗公式確定，具體包括以下子步驟：

步驟4.1熱紅外光譜曲線的觀察分析，找出與石英礦物含量相關的波段范圍；

步驟4.2采集不同巖性的巖漿巖標本進行熱紅外光譜數據采集；

步驟4.3將步驟4.2采集到的熱紅外光譜數據進行重采樣，使其與步驟一采用的熱紅外高光譜遙感數據相匹配，并結合步驟4.1找出石英礦物含量密切相關的波段組合；

步驟4.4對步驟4.2所述巖石標本進行地球化學測試，獲取全量分析數據；

步驟4.5根據各樣品全量分析數據，使用CIPW標準礦物計算方法計算出各樣品的石英礦物含量；

步驟4.6將步驟4.3得到的與石英含量密切相關的波段組合在不同組合方式下，進行四則運算、對數運算及指數運算，然后結合步驟4.4得到的石英含量進行回歸分析，獲取相關性最高的波段組合形式；

步驟4.7根據石英礦物含量及其對應的發射率光譜曲線中相關性最高的波段組合，得出石英礦物含量定量計算經驗公式；

步驟五、石英礦物含量圖獲取，將步驟4.7所述經驗公式代入發射率圖像，進行波段運算，獲取石英礦物含量圖。

2.根據權利要求1所述的一種適合于熱紅外高光譜遙感的石英礦物含量定量計算方法，其特征在于：所述的步驟三中，石英含量進行回歸分析，相關性最高的波段組合形式為ln[B2^2/(B1*B3)]；其中，B1＝B_8.49，B2＝B_8.60μm，B3＝B_8.71μm，B_8.82可作為B_8.49的備用波段，下標表示該波段的波長。

3.根據權利要求2所述的一種適合于熱紅外高光譜遙感的石英礦物含量定量計算方法，其特征在于：得出石英礦物含量定量計算經驗公式，y＝8208.9x³-5525x²+1228.2x+0.3139，x＝ln[B2^2/(B1*B3]。

4.根據權利要求3所述的一種適合于熱紅外高光譜遙感的石英礦物含量定量計算方法，其特征在于：結合了CIPW標準礦物計算方法，以獲取SiO₂、TiO₂、Al₂O₃、TFe₂O₃、MnO、MgO、CaO、Na₂O、K₂O及P₂O₅的全量分析數據計算出各樣品石英礦物含量。

5.根據權利要求4所述的一種適合于熱紅外高光譜遙感的石英礦物含量定量計算方法，其特征在于：與石英礦物含量相關的波段范圍具體為8.60μm的發射峰至兩側發射谷范圍的光譜特征。