一種用于葉蠟石信息提取的高光譜影像處理方法
【技術(shù)領(lǐng)域】
[0001]本發(fā)明屬于高光譜影像處理方法，特別是涉及一種用于葉蠟石信息提取的高光譜影像處理方法。
【背景技術(shù)】
[0002]當(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)容】

[0003]本發(fā)明要解決的技術(shù)問題是提供一種降低其他地物或噪聲的影響，減少了處理的數(shù)據(jù)量的用于葉蠟石信息提取的高光譜影像處理方法。通過對高光譜影像重采樣，提取特定波段，進行一系列判斷和計算，計算出影像范圍內(nèi)的不同區(qū)域葉蠟石的豐度值。
[0004]為解決上述技術(shù)問題，本發(fā)明提供一種用于葉蠟石信息提取的高光譜影像處理方法，依次包括:
[0005]步驟一，預(yù)處理
[0006]獲取高光譜影像；對高光譜影像的進行預(yù)處理，進行大氣校正，獲取地面反射率的影像數(shù)據(jù)；
[0007]步驟二:采樣
[0008]對步驟一得到的影像數(shù)據(jù)進行采樣，對波段在2045nm，2075nm, 2090nm, 2105nm, 2165nm, 2195nm, 2240nm, 2315nm, 2330nm, 2375nm, 2390nm 的圖像采樣，并依次記錄為 bl ?bll，即bl為波段2045nm的采樣數(shù)據(jù)，b2為波段2075nm的采樣數(shù)據(jù)，以此類推；
[0009]步驟三:判斷
[0010]按進行下面一系列判斷，并記錄結(jié)果
[0011]al=(bl 大于 b2);
[0012]a2=(b3小于b4);
[0013]a3= (b4 大于 b5);
[0014]a4= (b5 小于 b6)；
[0015]a5= (b7 大于 b8)；
[0016]a6= (b8 小于 b9)；
[0017]a7= (blO 大于 bll);
[0018]上述判斷是針對每次判斷圖像的相應(yīng)像元進行判斷的；
[0019]本步驟結(jié)束后得到al?a7，共7個矩陣；
[0020]步驟四:計算
[0021]用下述公式進行計算bl2
[0022]bI2=bl+2*b4+b6+b7+b9+b10-b2-b3-2*b5-2*b8_bII
[0023]所述的*表示相乘；
[0024]用下面公式計算a0
[0025]a0=al*a2*a3*a4*a5*a6*a7*bl2
[0026]上述所有計算均為相應(yīng)像元計算。
[0027]本發(fā)明本方法的效果在于:只用了 11個波段，比如SASI全波段高光譜影像有101個波段，用這種方法處理數(shù)據(jù)量可以減少89%，并且由于是計算機自動一步提取，減少了主成分變換、端元波譜的選擇等操作步驟，運算速度可以提高了 9倍以上。由于去除了大部分對信息提取關(guān)系不大的波段，減少其他物質(zhì)或噪聲對其光譜的干擾，提高了信息提取的精度。對高光譜影像數(shù)據(jù)中葉蠟石信息的快速提取具有較好的作用和意義。
【具體實施方式】
[0028]本發(fā)明依次包括下述步驟:
[0029]步驟一，預(yù)處理
[0030]獲取高光譜影像；對高光譜影像的進行預(yù)處理，進行大氣校正，獲取地面反射率的影像數(shù)據(jù)；
[0031]本步驟所述的進行大氣校正是本領(lǐng)域公知的處理方法。
[0032]步驟二:采樣
[0033]對步驟一得到的影像數(shù)據(jù)進行采樣，對波段在2045nm，2075nm, 2090nm, 2105nm, 2165nm, 2195nm, 2240nm, 2315nm, 2330nm, 2375nm, 2390nm 的圖像采樣，并依次記錄為 bl ?bll,即bl為波段2045nm的采樣數(shù)據(jù)，b2為波段2075nm的采樣數(shù)據(jù)，以此類推。每次采樣得到的都是一幅灰度圖，圖中每一個像元的值都是其灰度值，即bl為一幅灰度圖，bl圖像的(1，O點的值為灰度值，其余點依此類推，其余采樣圖也依次類推。
[0034]步驟三:判斷
[0035]按進行下面一系列判斷，并記錄結(jié)果
[0036]al=(bl 大于 b2);
[0037]a2=(b3小于b4);
[0038]a3= (b4 大于 b5);
[0039]a4= (b5 小于 b6)；
[0040]a5= (b7 大于 b8)；
[0041]a6= (b8 小于 b9)；
[0042]a7= (blO 大于 bll);
[0043]上述判斷是針對每次判斷圖像的相應(yīng)像元進行判斷的，以al= (bl大于b2)為例，取bl圖像的某像元(例如(1，I)點)的灰度值，與b2圖像的相應(yīng)像元(當(dāng)bl圖像取(1，I)點，則b2圖像也必須取(1，I)點)的灰度值，然后根據(jù)判斷規(guī)則“bl大于b2”判斷，當(dāng)判斷結(jié)果為“是”時，記錄判斷結(jié)果為1，否則記錄結(jié)果為O。因此當(dāng)al= (bl大于b2)判斷完畢時，得到的al是與bl矩陣相同大小的矩陣，其中每個點的值是根據(jù)判斷規(guī)則得到的判斷結(jié)果(即每個點的值是O或I)。
[0044]其它判斷也按照類似的規(guī)則進行。本步驟結(jié)束后得到al?a7，共7個矩陣。
[0045]步驟四:計算
[0046]用下述公式進行計算bl2
[0047]bI2=bl+2*b4+b6+b7+b9+b10-b2-b3-2*b5-2*b8_bII
[0048]所述的*表示相乘。
[0049]用下面公式計算a0
[0050]a0=al*a2*a3*a4*a5*a6*a7*bl2
[0051]上述所有計算均為相應(yīng)像元計算，即使用不同圖像的相應(yīng)像元計算。以bl2=bl+2*b4+b6+b7+b9+bl0-b2-b3-2*b5-2*b8-bll 公式為例，當(dāng)計算點(x，y)時，取 bl、b4、b6、b7、b9、bl0、b2、b3、b5、b8、bll的點(x,y)的灰度值參與計算，得到的結(jié)果是bl7的點(x,y)的值。又例如 a0=al*a2*a3*a4*a5*a6*a7*bl2,當(dāng)計算點(x，y)時，取 al、a2、a3、a4、a5、a6、a7、bl2的點(x,y)的值參與計算，得到的結(jié)果是aO的點(x,y)的值。
[0052]計算得到的aO就是葉蠟石信息的豐度圖，即圖像中某區(qū)域的數(shù)值越大表示該區(qū)域葉蠟石的豐度越高。
【主權(quán)項】
1.一種用于葉蠟石信息提取的高光譜影像處理方法，其特征在于:包括下述步驟 步驟一，預(yù)處理 獲取高光譜影像；對高光譜影像的進行預(yù)處理，進行大氣校正，獲取地面反射率的影像數(shù)據(jù)； 步驟二:采樣 對步驟一得到的影像數(shù)據(jù)進行采樣，對波段在2045nm，2075nm, 2090nm, 2105nm, 2165nm，2195nm, 2240nm, 2315nm, 2330nm, 2375nm, 2390nm 的圖像采樣，并依次記錄為 bl ?bll，即bl為波段2045nm的采樣數(shù)據(jù)，b2為波段2075nm的采樣數(shù)據(jù)，以此類推； 步驟三:判斷 按進行下面一系列判斷，并記錄結(jié)果 al= (bl 大于 b2)； a2= (b3 小于 b4)； a3= (b4 大于 b5)； a4= (b5 小于 b6)； a5= (b7 大于 b8)； a6= (b8 小于 b9)； a7= (blO 大于 bll)； 上述判斷是針對每次判斷圖像的相應(yīng)像元進行判斷的； 本步驟結(jié)束后得到al?a7，共7個矩陣； 步驟四:計算 用下述公式進行計算bl2bI2=bl+2*b4+b6+b7+b9+bl0-b2-b3-2*b5-2*b8-bII所述的*表示相乘； 用下面公式計算aO a0=al*a2*a3*a4*a5*a6*a7*bl2 上述所有計算均為相應(yīng)像元計算。
【專利摘要】本發(fā)明屬于高光譜影像處理方法，特別是涉及一種用于葉蠟石信息提取的高光譜影像處理方法。它包括：步驟一，預(yù)處理；步驟二：采樣，對波段在2045nm，2075nm，2090nm，2105nm，2165nm，2195nm，2240nm，2315nm，2330nm，2375nm，2390nm的圖像采樣；步驟三：判斷，和步驟四：計算。本發(fā)明本方法的效果在于：只用了11個波段，比如SASI全波段高光譜影像有101個波段，用這種方法處理數(shù)據(jù)量可以減少89%，并且由于是計算機自動一步提取，減少了主成分變換、端元波譜的選擇等操作步驟，運算速度可以提高了9倍以上。由于去除了大部分對信息提取關(guān)系不大的波段，減少其他物質(zhì)或噪聲對其光譜的干擾，提高了信息提取的精度。對高光譜影像數(shù)據(jù)中葉蠟石信息的快速提取具有較好的作用和意義。
【IPC分類】G06K9-46, G06T5-00
【公開號】CN104574283
【申請?zhí)枴緾N201310503075
【發(fā)明人】楊燕杰, 趙英俊, 崔欣, 秦凱, 高仉生, 周覓, 孫雨, 宮寶昌, 錢坤, 楊國防, 崔建勇, 劉漢斌
【申請人】核工業(yè)北京地質(zhì)研究院
【公開日】2015年4月29日
【申請日】2013年10月23日