一種用于水鋁礦信息提取的高光譜影像處理方法
【技術(shù)領(lǐng)域】
[0001]本發(fā)明屬于高光譜影像處理方法，特別是涉及一種用于水鋁礦信息提取的高光譜影像處理方法。
【背景技術(shù)】
[0002]當前的高光譜遙感影像的水鋁礦信息提取方法主要是光譜全波段匹配或是部分連續(xù)波段的光譜匹配，具體算法有光譜角、混合解調(diào)濾波等，由于地表的物質(zhì)組成很少是由單一礦物組成的，這些方法在信息提取的過程易受其他地物波譜或噪聲的影響，提取信息精度相對較低。其次現(xiàn)有的光譜提取方法人工操作步驟多，增加了人為判斷誤差。第三是高光譜數(shù)據(jù)波段多，數(shù)據(jù)量大，現(xiàn)有的方法處理時間長，降低了數(shù)據(jù)處理的速度和應用規(guī)模。因此，如何在水鋁礦信息提取的過程中減少其他地物或噪聲的影響、人工操作步驟和處理數(shù)據(jù)量，成為當前高光譜遙感影像處理的前沿之一。

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

[0003]本發(fā)明要解決的技術(shù)問題是提供一種降低其他地物或噪聲的影響，減少了處理的數(shù)據(jù)量的用于水鋁礦信息提取的高光譜影像處理方法。通過對高光譜影像重采樣，提取特定波段，進行一系列判斷和計算，計算出影像范圍內(nèi)的不同區(qū)域水鋁礦的豐度值。
[0004]為解決上述技術(shù)問題，本發(fā)明提供一種用于水鋁礦信息提取的高光譜影像處理方法，依次包括:
[0005]步驟一，預處理
[0006]獲取高光譜影像；對高光譜影像的進行預處理，進行大氣校正，獲取地面反射率的影像數(shù)據(jù)；
[0007]步驟二:采樣
[0008]對步驟一得到的影像數(shù)據(jù)進行采樣，對波段在1355nm，1445nm, 1490nm, 1520nm, 1550nm, 1625nm, 1775nm, 1940nm, 2150nm, 2255nm 的圖像采樣，并依次記錄為 bl ?blO,即 bl為波段1355nm的采樣數(shù)據(jù)，b2為波段1445nm的采樣數(shù)據(jù)，以此類推；
[0009]步驟三:判斷
[0010]按進行下面一系列判斷，并記錄結(jié)果
[0011]al = (bl 大于 b2)；
[0012]a2 = (b2 小于 b3)；
[0013]a3 = (b3 大于 b4)；
[0014]a4 = (b5 小于 b6)；
[0015]a5 = (b7 大于 b8)；
[0016]a6 = (b9 大于 blO)；
[0017]上述判斷是針對每次判斷圖像的相應像元進行判斷的；
[0018]本步驟結(jié)束后得到al?a6，共6個矩陣；
[0019]步驟四:計算
[0020]用下述公式進行計算bll
[0021]bll = bl+2*b3+b6+b7+b9-2*b2-b4-b5-b8_bl0
[0022]所述的*表示相乘；
[0023]用下面公式計算aO
[0024]aO = al*a2*a3*a4*a5*a6*bll
[0025]上述所有計算均為相應像元計算。
[0026]如上所述的一種用于水鋁礦信息提取的高光譜影像處理方法，其中，所述的步驟三中的判斷是針對每次判斷圖像的相應像元進行判斷的。
[0027]如上所述的一種用于水鋁礦信息提取的高光譜影像處理方法，其中，所述的對圖像的相應像元進行判斷是按照下述過程進行的:
[0028]取需要判斷圖像的對應像元，判斷像元灰度是否符合判斷標準，如果符合，判定該像元位置的結(jié)果為“ 1 ”，否則判定該像元位置的結(jié)果為“0”，然后判斷下一個對應像元，直到完成對整個圖像的判斷。
[0029]本發(fā)明本方法的效果在于:只用了 10個波段，相對于高光譜影像SASI全波段101個波段，需要處理的數(shù)據(jù)量減少了 90 %，并且由于是計算機自動一步提取，減少了主成分變換、端元波譜的選擇等操作步驟，運算速度可以提高了 10倍以上。由于去除了大部分對信息提取關(guān)系不大的波段，減少其他物質(zhì)或噪聲對其光譜的干擾，提高了信息提取的精度。對高光譜影像數(shù)據(jù)中水鋁礦信息的快速提取具有較好的作用和意義。
【具體實施方式】
[0030]本發(fā)明依次包括下述步驟:
[0031]步驟一，預處理
[0032]獲取高光譜影像；對高光譜影像的進行預處理，進行大氣校正，獲取地面反射率的影像數(shù)據(jù)；
[0033]本步驟所述的進行大氣校正是本領(lǐng)域公知的處理方法。
[0034]步驟二:采樣
[0035]對步驟一得到的影像數(shù)據(jù)進行采樣，對波段在1355nm，1445nm, 1490nm, 1520nm, 1550nm, 1625nm, 1775nm, 1940nm, 2150nm, 2255nm 的圖像采樣，并依次記錄為 bl ?blO,即 bl為波段1355nm的采樣數(shù)據(jù)，b2為波段1445nm的采樣數(shù)據(jù)，以此類推。每次采樣得到的都是一幅灰度圖，圖中每一個像元的值都是其灰度值，即bl為一幅灰度圖，bl圖像的(1，1)點的值為灰度值，其余點依此類推，其余采樣圖也依次類推。
[0036]步驟三:判斷
[0037]按進行下面一系列判斷，并記錄結(jié)果
[0038]al = (bl 大于 b2)；
[0039]a2 = (b2 小于 b3)；
[0040]a3 = ((b3 大于 b4)；
[0041]a4 = (b5 小于 b6)；
[0042]a5 = (b7 大于 b8)；
[0043]a6 = (b9 大于 blO)；
[0044]上述判斷是針對每次判斷圖像的相應像元進行判斷的，以al = (bl大于b2)為例，取bl圖像的某像元(例如(1，1)點)的灰度值，與b2圖像的相應像元(當bl圖像取(1，1)點，則b2圖像也必須取(1，1)點)的灰度值，然后根據(jù)判斷規(guī)則“bl大于b2”判斷，當判斷結(jié)果為“是”時，記錄判斷結(jié)果為1，否則記錄結(jié)果為0。因此當al = (bl大于b2)判斷完畢時，得到的al是與bl矩陣相同大小的矩陣，其中每個點的值是根據(jù)判斷規(guī)則得到的判斷結(jié)果(即每個點的值是0或1)。
[0045]其它判斷也按照類似的規(guī)則進行。本步驟結(jié)束后得到al?a6，共6個矩陣。
[0046]步驟四:計算
[0047]用下述公式進行計算bll
[0048]bll = bl+2*b3+b6+b7+b9-2*b2-b4-b5-b8_bl0 所述的 * 表不相乘。
[0049]用下面公式計算aO
[0050]aO = al*a2*a3*a4*a5*a6*bll
[0051]上述所有計算均為相應像元計算，即使用不同圖像的相應像元計算。以bll = bl+2*b3+b6+b7+b9-2*b2-b4-b5-b8-bl0 公式為例，當計算點(x，y)時，取 bl、b3、b6、b7、b9、b2、b4、b5、b8、bl0的點(x, y)的灰度值參與計算,得到的結(jié)果是bll的點(x, y)的值。又例如 aO = al*a2*a3*a4*a5*a6*bll,當計算點(x, y)時，取 al、a2、a3、a4、a5、a6、bll 的點(x，y)的值參與計算，得到的結(jié)果是aO的點(X, y)的值。
[0052]計算得到的aO就是水鋁礦信息的豐度圖，即圖像中某區(qū)域的數(shù)值越大表示該區(qū)域水鋁礦的豐度越高。
【主權(quán)項】
1.一種用于水鋁礦信息提取的高光譜影像處理方法，其特征在于:包括下述步驟 步驟一，預處理 獲取高光譜影像；對高光譜影像的進行預處理，進行大氣校正，獲取地面反射率的影像數(shù)據(jù)； 步驟二:采樣 對步驟一得到的影像數(shù)據(jù)進行采樣，對波段在1355nm，1445nm, 1490nm, 1520nm, 1550nm，1625nm, 1775nm, 1940nm, 2150nm, 2255nm的圖像采樣，并依次記錄為bl?blO,即bl為波段2045nm的采樣數(shù)據(jù)，b2為波段2075nm的采樣數(shù)據(jù)，以此類推； 步驟三:判斷 按進行下面一系列判斷，并記錄結(jié)果 al = (bl 大于 b2)； a2 = (b2 小于 b3)； a3 = ((b3 大于 b4)； a4 = (b5 小于 b6)； a5 = (b7 大于 b8)； a6 = (b9 大于 blO)； 上述判斷是針對每次判斷圖像的相應像元進行判斷的； 本步驟結(jié)束后得到al?a6，共6個矩陣； 步驟四:計算 用下述公式進行計算bllbll = bl+2*b3+b6+b7+b9-2*b2-b4-b5-b8-bl0 所述的*表示相乘； 用下面公式計算aOaO = al*a2*a3*a4*a5*a6*bll 上述所有計算均為相應像元計算。2.如權(quán)利要求1所述的一種用于水鋁礦信息提取的高光譜影像處理方法，其特征在于:所述的步驟三中的判斷是針對每次判斷圖像的相應像元進行判斷的。3.如權(quán)利要求2所述的一種用于水鋁礦信息提取的高光譜影像處理方法，其特征在于:所述的對圖像的相應像元進行判斷是按照下述過程進行的: 取需要判斷圖像的對應像元，判斷像元灰度是否符合判斷標準，如果符合，判定該像元位置的結(jié)果為“1”，否則判定該像元位置的結(jié)果為“0”，然后判斷下一個對應像元，直到完成對整個圖像的判斷。
【專利摘要】本發(fā)明屬于高光譜影像處理方法，特別是涉及一種用于水鋁礦信息提取的高光譜影像處理方法。它包括：步驟一，預處理；步驟二：采樣，對波段在1355nm，1445nm，1490nm，1520nm，1550nm，1625nm，1775nm，1940nm，2150nm，2255nm的圖像采樣；步驟三：判斷，和步驟四：計算。本發(fā)明本方法的效果在于：只用了10個波段，相對于高光譜影像SASI全波段101個波段，需要處理的數(shù)據(jù)量減少了90％，并且由于是計算機自動一步提取，減少了主成分變換、端元波譜的選擇等操作步驟，運算速度可以提高了10倍以上。由于去除了大部分對信息提取關(guān)系不大的波段，減少其他物質(zhì)或噪聲對其光譜的干擾，提高了信息提取的精度。對高光譜影像數(shù)據(jù)中水鋁礦信息的快速提取具有較好的作用和意義。
【IPC分類】G06T5/00, G06T7/00
【公開號】CN105405102
【申請?zhí)枴緾N201410431253
【發(fā)明人】楊燕杰, 趙英俊, 武鼎, 宗自華, 張川, 陸冬華, 錢坤, 車永飛, 葉發(fā)旺, 李瀚波, 高仉生, 吳文歡
【申請人】核工業(yè)北京地質(zhì)研究院
【公開日】2016年3月16日
【申請日】2014年8月28日