[發(fā)明專利]一種用于綠泥石信息提取的高光譜影像處理方法

在審

• 申請?zhí)枺?/td>: 201210580124.2
• 申請日: 2012-12-27
• 公開（公告）號: CN103902998A
• 公開（公告）日: 2014-07-02
• 發(fā)明（設計）人: 楊燕杰;趙英俊
• 申請（專利權）人: 核工業(yè)北京地質(zhì)研究院
• 主分類號: G06K9/46
• 分類號: G06K9/46
• 代理公司: 核工業(yè)專利中心 11007
• 代理人: 高尚梅;劉昕宇

技術領域

本發(fā)明屬于信息提取方法，具體涉及一種用于綠泥石信息提取的高光譜影像處理方法。?

背景技術

通過對圖像信息處理來提取圖像區(qū)域內(nèi)不同類型礦物的豐度，是目前勘探領域常用的方法。這種方法不需要工作人員直接到野外勘探、采樣，可以提取不同類型礦物的豐度。當前的高光譜遙感影像的綠泥石信息提取方法主要是光譜全波段匹配或是部分連續(xù)波段的光譜匹配，具體算法有光譜角、混合解調(diào)濾波等，由于地表的物質(zhì)組成很少是由單一礦物組成的，這些方法在信息提取的過程易受其他地物波譜或噪聲的影響，提取信息精度相對較低。其次現(xiàn)有的光譜提取方法人工操作步驟多，增加了人為判斷誤差。第三是高光譜數(shù)據(jù)波段多，數(shù)據(jù)量大，現(xiàn)有的方法處理時間長，降低了數(shù)據(jù)處理的速度和應用規(guī)模。因此，如何在綠泥石信息提取的過程中減少其他地物或噪聲的影響、人工操作步驟和處理數(shù)據(jù)量，成為當前高光譜遙感影像處理的前沿之一。?

發(fā)明內(nèi)容

本發(fā)明針對現(xiàn)有技術的缺陷，提供一種用于綠泥石信息提取的高光譜影像處理方法。?

本發(fā)明是這樣實現(xiàn)的：一種用于綠泥石信息提取的高光譜影像處?理方法，包括，?

步驟一：波段采樣?

對已有的地面反射率的高光譜影像數(shù)據(jù)進行采樣，提取波段在1130nm,1790nm,2165nm,2255nm,2270nm,2330nm,2360nm的影像，并依次記錄為b1～b7，每次采樣得到的都是一幅灰度圖，圖中每一個像元的值都是其灰度值，?

步驟二：判斷?

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

a1=（b1小于b2）；?

a2=（b3大于b4）；?

a3=（b5大于b6）；?

a4=（b6小于b7）；?

上述判斷是針對每次判斷圖像的相應像元進行判斷的，當判斷結果為“是”時，記錄判斷結果為1，否則記錄結果為0，?

步驟三：計算?

用下述公式進行計算b8?

b8=b2+b3+b5+b7-b1-b4-2*b6?

所述的*表示相乘，?

用下面公式計算a0?

a0=a1*a2*a3*a4*b8?

上述所有計算均為相應像元計算，即使用不同圖像的相應像元計算。?

本發(fā)明的效果是：本發(fā)明提取特定波段，進行一系列判斷和計算，可以準確的計算出影像范圍內(nèi)的不同區(qū)域綠泥石的豐度值。保留綠泥石的光譜特征明顯的波段，去除其他特征不明顯的波段，從而在信息提取的過程突出綠泥石的光譜特征，降低其他地物或噪聲的影響，減少了處理的數(shù)據(jù)量，并可以用IDL程序達到用較少的人工操作實現(xiàn)最終結果信息提取的目的，提高了綠泥石信息提取的精度和速度。?

具體實施方式

下面結合實施例對本發(fā)明作進一步說明。一種用于綠泥石信息提取的高光譜影像處理方法，包括，?

步驟一：波段采樣?

對已有的地面反射率的高光譜影像數(shù)據(jù)進行采樣，提取波段在1130nm,1790nm,2165nm,2255nm,2270nm,2330nm,2360nm的影像，并依次記錄為b1～b7，即b1為波段1130nm的采樣數(shù)據(jù)，b2為波段1790nm的采樣數(shù)據(jù)，以此類推。每次采樣得到的都是一幅灰度圖，圖中每一個像元的值都是其灰度值，即b1為一幅灰度圖，b1圖像的（1、1）點的值為灰度值，其余點依此類推，其余采樣圖也依次類推。?

步驟二：判斷?

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

a1=（b1小于b2）；?

a2=（b3大于b4）；?

a3=（b5大于b6）；?

a4=（b6小于b7）；?

上述判斷是針對每次判斷圖像的相應像元進行判斷的，以a1=（b1小于b2）為例，取b1圖像的某像元（例如（1、1）點）的灰度值，與b2圖像的相應像元（當b1圖像取（1、1）點，則b2圖像也必須?。?、1）點）的灰度值，然后根據(jù)判斷規(guī)則“b1小于b2”判斷，當判斷結果為“是”時，記錄判斷結果為1，否則記錄結果為0。因此當a1=（b1小于b2）判斷完畢時，得到的a1是與b1矩陣相同大小的矩陣，其中每個點的值是根據(jù)判斷規(guī)則得到的判斷結果（即每個點的值是0或1）。?

其它判斷也按照類似的規(guī)則進行。本步驟結束后得到a1～a4，共4個矩陣（影像）。?

步驟三：計算?

用下述公式進行計算b8?

b8=b2+b3+b5+b7-b1-b4-2*b6?