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

有效

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

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

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

背景技術(shù)

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

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

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

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

步驟一：波段采樣

對(duì)已有的地面反射率的影像數(shù)據(jù)進(jìn)行采樣，提取在1415nm、1430nm、1940nm、1970nm、2165nm、2195nm、2225nm、2240nm、2285nm、2345nm的影像，并依次記錄為b1～b10，每次采樣得到的都是一幅灰度圖，圖中每一個(gè)像元的值都是其灰度值，

步驟二：判斷

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

a1＝(b1小于b2)；

a2＝(b3小于b4)；

a3＝(b5大于b6)；

a4＝(b7小于b8)；

a5＝(b9大于b10)；

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

步驟三：計(jì)算

用下述公式進(jìn)行計(jì)算b11

b11＝b2+b4+b5+b8+b9-b1-b3-b6-b7-b10

所述的*表示相乘，

用下面公式計(jì)算a0，

a0＝a1*a2*a3*a4*a5*b11

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

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

具體實(shí)施方式

下面結(jié)合實(shí)施例對(duì)本發(fā)明作進(jìn)一步說(shuō)明。一種用于伊利石信息提取的高光譜影像處理方法，包括

步驟一：波段采樣

對(duì)已有的地面反射率的影像數(shù)據(jù)進(jìn)行采樣，提取在1415nm、1430nm、1940nm、1970nm、2165nm、2195nm、2225nm、2240nm、2285nm、2345nm的影像，并依次記錄為b1～b10，即b1為波段1415nm的采樣數(shù)據(jù)，b2為波段1430nm的采樣數(shù)據(jù)，以此類推。每次采樣得到的都是一幅灰度圖，圖中每一個(gè)像元的值都是其灰度值，即b1為一幅灰度圖，b1圖像的(1,1)點(diǎn)的值為灰度值，其余點(diǎn)依此類推，其余采樣圖也依次類推。

步驟二：判斷

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

a1＝(b1小于b2)；

a2＝(b3小于b4)；

a3＝(b5大于b6)；

a4＝(b7小于b8)；

a5＝(b9大于b10)；

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

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