[發明專利]一種用于陽起石信息提取的高光譜影像處理方法

有效

• 申請號: 201310712521.5
• 申請日: 2013-12-20
• 公開（公告）號: CN104732488B
• 公開（公告）日: 2018-07-06
• 發明（設計）人: 周覓;閔蘋;王鎖會;陳興旺;王俊虎;武鼎;黃艷菊
• 申請（專利權）人: 核工業北京地質研究院
• 主分類號: G06T5/00
• 分類號: G06T5/00;G06T7/00
• 代理公司: 核工業專利中心 11007
• 代理人: 高尚梅;劉昕宇

本發明屬于高光譜影像處理方法，特別是涉及一種用于陽起石信息提取的高光譜影像處理方法。它包括：步驟一，預處理；步驟二：采樣，對波段在1040nm，1700nm，2270nm，2300nm，2315nm，2345nm，2360nm，2390nm，2420nm的圖像采樣；步驟三：判斷，和步驟四：計算。本發明本方法的效果在于：只用了9個波段，如果用于高光譜影像SASI數據處理（有101個波段），需要處理的數據量減少了91%，并且由于是計算機自動一步提取，減少了主成分變換、端元波譜的選擇等操作步驟，運算速度可以提高了11倍以上。由于去除了大部分對信息提取關系不大的波段，減少其他物質或噪聲對其光譜的干擾，提高了信息提取的精度。對高光譜影像數據中陽起石信息的快速提取具有較好的作用和意義。

技術領域

本發明屬于高光譜影像處理方法，特別是涉及一種用于陽起石信息提取的高光譜影像處理方法。

背景技術

當前的高光譜遙感影像的陽起石信息提取方法主要是光譜全波段匹配或是部分連續波段的光譜匹配，具體算法有光譜角、混合解調濾波等，由于地表的物質組成很少是由單一礦物組成的，這些方法在信息提取的過程易受其他地物波譜或噪聲的影響，提取信息精度相對較低。其次現有的光譜提取方法人工操作步驟多，增加了人為判斷誤差。第三是高光譜數據波段多，數據量大，現有的方法處理時間長，降低了數據處理的速度和應用規模。因此，如何在陽起石信息提取的過程中減少其他地物或噪聲的影響、人工操作步驟和處理數據量，成為當前高光譜遙感影像處理的前沿之一。

發明內容

本發明要解決的技術問題是提供一種降低其他地物或噪聲的影響，減少了處理的數據量的用于陽起石信息提取的高光譜影像處理方法。通過對高光譜影像重采樣，提取特定波段，進行一系列判斷和計算，計算出影像范圍內的不同區域陽起石的豐度值。

為解決上述技術問題，本發明提供一種用于陽起石信息提取的高光譜影像處理方法，依次包括：

步驟一，預處理

獲取高光譜影像；對高光譜影像的進行預處理，進行大氣校正，獲取地面反射率的影像數據；

步驟二：采樣

對步驟一得到的影像數據進行采樣，對波段在1040nm，1700nm，2270nm，2300nm，2315nm，2345nm，2360nm，2390nm，2420nm的圖像采樣，并依次記錄為b1～b9，即b1為波段1040nm的采樣數據，b2為波段1700nm的采樣數據，以此類推；

步驟三：判斷

按進行下面一系列判斷，并記錄結果

a1=(b1小于b2)；

a2=(b3大于b4)；

a3=(b5小于b6)；

a4=(b7大于b8)；

a5=(b8小于b9)；

上述判斷是針對每次判斷圖像的相應像元進行判斷的；

本步驟結束后得到a1～a5，共5個矩陣；

步驟四：計算

用下述公式進行計算b10

b10=b2+b3+b6+b7+b9–b1-b4-b5-2*b8

所述的*表示相乘；

用下面公式計算a0

a0=a1*a2*a3*a4*a5*b10

上述所有計算均為相應像元計算。