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

有效

• 申請?zhí)枺?/td>: 201210579987.8
• 申請日: 2012-12-27
• 公開（公告）號: CN103900966B
• 公開（公告）日: 2017-09-29
• 發(fā)明（設(shè)計）人: 楊燕杰;趙英俊
• 申請（專利權(quán)）人: 核工業(yè)北京地質(zhì)研究院
• 主分類號: G01N21/25
• 分類號: G01N21/25
• 代理公司: 核工業(yè)專利中心11007
• 代理人: 高尚梅,劉昕宇

技術(shù)領(lǐng)域

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

背景技術(shù)

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

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

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

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

步驟一：波段采樣

對已有的地面反射率的影像數(shù)據(jù)進行采樣，提取波段在1085nm、1190nm、2225nm、2255nm、2270nm、2285nm、2330nm、2345nm、2375nm的影像，并依次記錄為b1～b9，每次采樣得到的都是一幅灰度圖，圖中每一個像元的值都是其灰度值，

步驟二：判斷

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

a1=（b1小于b2）；

a2=（b3大于b4）；

a3=（b4小于b5）；

a4=（b6大于b7）；

a5=（b8小于b9）；

上述判斷是針對每次判斷圖像的相應(yīng)像元進行判斷的，當(dāng)判斷結(jié)果為“是”時，記錄判斷結(jié)果為1，否則記錄結(jié)果為0，

步驟三：計算

用下述公式進行計算b10

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

所述的*表示相乘，

用下面公式計算a0，

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

上述所有計算均為相應(yīng)像元計算，即使用不同圖像的相應(yīng)像元計算。

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

具體實施方式

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

步驟一：波段采樣

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

步驟二：判斷

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

a1=（b1小于b2）；

a2=（b3大于b4）；

a3=（b4小于b5）；

a4=（b6大于b7）；

a5=（b8小于b9）；

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

其它判斷也按照類似的規(guī)則進行。本步驟結(jié)束后得到a1～a5，共5個矩陣。

步驟三：計算