基於近紅外成像高光譜技術的岩石光譜隨水滲入時間的變化研究

基於近紅外成像高光譜技術的岩石光譜隨水滲入時間的變化研究

四川雙利合譜科技有限公司-黃宇

　　一
、測試原理及方法：

　　高(gao)光(guang)譜(pu)成(cheng)像(xiang)技(ji)術(shu)是(shi)近(jin)二(er)十(shi)年(nian)來(lai)發(fa)展(zhan)起(qi)來(lai)的(de)基(ji)於(yu)非(fei)常(chang)多(duo)窄(zhai)波(bo)段(duan)的(de)影(ying)像(xiang)數(shu)據(ju)技(ji)術(shu)，其(qi)最(zui)突(tu)出(chu)的(de)應(ying)用(yong)是(shi)遙(yao)感(gan)探(tan)測(ce)領(ling)域(yu)
，並(bing)在(zai)越(yue)來(lai)越(yue)多(duo)的(de)民(min)用(yong)領(ling)域(yu)有(you)著(zhe)更(geng)大(da)的(de)應(ying)用(yong)前(qian)景(jing)。它(ta)集(ji)中(zhong)了(le)光(guang)學(xue)、光電子學
、電子學、信息處理
、計算機科學等領域的先進技術，是傳統的二維成像技術和光譜技術有機的結合在一起的一門新興技術。

　　高光譜成像技術的定義是在多光譜成像的基礎上
，在從紫外到近紅外（200-2500nm）deguangpufanweinei，liyongchengxiangguangpuyi，zaiguangpufugaifanweineideshushihuoshubaitiaoguangpuboduanduimubiaowutilianxuchengxiang。zaihuodewutikongjiantezhengchengxiangdetongshi
，yehuodelebeicewutideguangpuxinxi。

目標物體-成像物鏡-入射狹縫-準直透鏡-PGP-聚焦透鏡-CCD棱鏡-光柵-棱鏡：PGP

圖1 成像原理圖

　　光譜儀的光譜分辨率由狹縫的寬度和光學光譜儀產生的線性色散確定。最小光譜分辨率是由光學係統的成像性能確定的（點擴展大小）。

　　成像過程為：每次成一條線上的像後（X方向）
，在檢測係統輸送帶前進的過程中
，排列的探測器掃出一條帶狀軌跡從而完成縱向掃描（Y方向）。綜合橫縱掃描信息就可以得到樣品的三維高光譜圖像數據。

圖2 像立方體

　　二、材料與分析：

　　1
、實驗設備與材料

　　本文以中國礦業大學（北京）深部岩土力學與地下工程國家重點實驗室提供的砂岩、礫岩為研究對象，利用四川雙利合譜科技有限公司的近紅外高光譜成像儀GaiaField（光譜範圍900 nm - 1700 nm）采集測試對象的高光譜影像數據，從而分析水滲入岩石時
，砂岩、礫岩的光譜反射率隨水滲入時間的變化規律。表1為近紅外高光譜成像儀GaiaField係統參數。

表1  近紅外高光譜成像儀GaiaField係統參數

　　2、高光譜影像數據預處理

　　對成像光譜儀獲取的原始影像數據進行數據的預處理，預處理過程主要包括兩部分。第一部分是輻射定標；第二部分為噪聲去除。

首先進行輻射定標。輻射定標的計算公式如1所示。

　　其中

，Ref_target為目標物的反射率
，Ref_panel為標準參考板的反射率，DN_target為原始影像中目標物的的數值，DN_panel為原始影像中標準參考板的數值
，DN_dark為成像光譜儀係統誤差。

　　其次是噪聲去除，常用的方法有小波降噪、S-G降噪


、均值濾波
、最小噪聲分離等方法。本研究則運用國外較為常用的最小噪聲分離方法(Minimum Noise Fraction Rotation， MNF)進行噪聲去除。最小噪聲分離工具用於判定圖像數據內在的維數（即波段數），分離數據中的噪聲
，減少隨後處理中的計算需求量。MNF本質上是兩次層疊的主成分變換。第一次變換（基於估計的噪聲協方差矩陣）用於分離和重新調節數據中的噪聲，這步操作使變換後的噪聲數據隻有最小的方差且沒有波段間的相關。第二步是對噪聲白化數據（Noise-whitened）的標準主成分變換。為了進一步進行波譜處理
，通過檢查最終特征值和相關圖像來判定數據的內在維數。數據空間可被分為兩部分：一yi部bu分fen與yu較jiao大da特te征zheng值zhi和he相xiang對dui應ying的de特te征zheng圖tu像xiang相xiang關guan
，其qi餘yu部bu分fen與yu近jin似si相xiang同tong的de特te征zheng值zhi以yi及ji噪zao聲sheng占zhan主zhu導dao地di位wei的de圖tu像xiang相xiang關guan。由you於yu此ci次ci采cai集ji的de高gao光guang譜pu影ying像xiang沒mei有you白bai板ban校xiao正zheng
，因yin此ci數shu據ju預yu處chu理li的de第di一yi步bu輻fu射she定ding標biao沒mei有you進jin行xing分fen析xi處chu理li，直zhi接jie作zuoMNF降噪分析。圖3為MNF降噪前後的成像高光譜數據中DN值的變化。

圖3 MNF變換前(左)後(右)高光譜影像DN值的變化

　　3、特征波段的獲取

　　特te征zheng波bo段duan的de選xuan擇ze對dui於yu建jian立li穩wen定ding的de光guang譜pu模mo型xing是shi至zhi關guan重zhong要yao的de。采cai用yong全quan波bo段duan的de光guang譜pu數shu據ju建jian立li模mo型xing時shi，不bu僅jin計ji算suan工gong作zuo量liang大da，而er且qie校xiao正zheng模mo型xing的de預yu測ce精jing度du也ye很hen難nan達da到dao最zui優you值zhi。因yin此ci在zai建jian立li校xiao正zheng模mo型xing前qian進jin行xing特te征zheng波bo段duan選xuan擇ze是shi很hen有you必bi要yao的de。目mu前qian，常chang用yong的de特te征zheng波bo段duan選xuan擇ze方fang法fa有you相xiang關guan係xi數shu法fa


、方差分析、逐步多元線性回歸、粒子群算法、反向區間偏最小二乘法和連續投影算法等。其中連續投影算法是一種新型的波長選擇方法，得到了越來越廣泛的應用(成忠等, 2010; Moreira et al., 2009; Gomes et al., 2013; Talfan et al., 2015; Liu et al., 2009)
，本研究也采用該算法進行特征波段的提取。連續投影算法(Successive Projections Algorithm, SPA) (Araujo et al, 2001) shiyizhongnenggoucongguangpuxinxizhongchongfenxunzhaohanyouzuidixianduderongyuxinxidebianliangzu
，shidegebianliangzhijiandegongxianxingdadaozuixiao。tashiyizhongquedingxingdesousuofangfa，qibianliangxuanzedejieguoshikezhongxiande，zheduiyuyanzhengjidebianliangxuanzehuigengwenjian
，tanenggoucongyanzhongzhongdiedeguangpuxinxizhongtiquyouxiaoxinxi，jixiaohuaguangpubianliangzhijiandegongxianxingyingxiang，jianhuaxiaozhengmoxinghetigaojianmodesuduhexiaolv。

　　四、結果與分析

　　1、 砂岩光譜反射率隨著水滲入時間的變化

　　weilenengjiaoweikeguandelejieshayanguangpufanshelvsuishuishenrushijiandebianhuaguilv，benyanjiuzaigaoguangpuyingxiangshangrenyihuoqushayandelianggexiangyuandian，guanchaqisuishuishenrushijianguangpubianhuaguilv
，xiangyuandianxuanquweizhijiqiguangpufanshelvsuishuishenrushijiandebianhuarutu4所示。從圖4可知，從900-1700 nm範圍內，隨水滲入時間的推移，砂岩的光譜反射率總體趨勢為下降，在1420 nm附近較為顯著的水汽吸收穀，隨後反射率逐漸上升，在1650 nm後反射率呈指數上升趨勢。

圖4砂岩光譜反射率隨著水滲入時間的變化（左為紅色，右為綠色）

　　2、礫岩光譜反射率隨著水滲入時間的變化

　　與(yu)砂(sha)岩(yan)分(fen)析(xi)相(xiang)似(si)
，為(wei)了(le)能(neng)較(jiao)為(wei)客(ke)觀(guan)的(de)了(le)解(jie)礫(li)岩(yan)光(guang)譜(pu)反(fan)射(she)率(lv)隨(sui)水(shui)滲(shen)入(ru)時(shi)間(jian)的(de)變(bian)化(hua)規(gui)律(lv)
，本(ben)研(yan)究(jiu)在(zai)礫(li)岩(yan)的(de)高(gao)光(guang)譜(pu)影(ying)像(xiang)上(shang)任(ren)意(yi)獲(huo)取(qu)兩(liang)個(ge)像(xiang)元(yuan)點(dian)，觀(guan)察(cha)其(qi)隨(sui)水(shui)滲(shen)入(ru)時(shi)間(jian)光(guang)譜(pu)變(bian)化(hua)規(gui)律(lv)，像(xiang)元(yuan)點(dian)選(xuan)取(qu)位(wei)置(zhi)及(ji)其(qi)光(guang)譜(pu)反(fan)射(she)率(lv)隨(sui)水(shui)滲(shen)入(ru)時(shi)間(jian)的(de)變(bian)化(hua)如(ru)圖(tu)5所示。從圖5可知，從總體上
，兩個像元的光譜反射率隨水滲入時間
，礫岩的光譜反射率變化趨勢相一致，不同的是反射率高低不一致（這是因為礫岩每個像元點之間也是存在差異的）。從900-1700 nm範圍內，隨水滲入時間的推移
，礫岩的光譜反射率總體趨勢為下降，在1430 nm附近較為顯著的水汽吸收穀，隨後反射率回升。

圖5礫岩光譜反射率隨著水滲入時間的變化（左為紅色，右為綠色）

　　3、砂岩光譜反射率與水滲入岩石時間的相關性分析

　　圖4分別列舉了砂岩高光譜影像中兩像元點隨水滲入時間的光譜反射率變化規律，從圖4可知，隨著水滲入時間的推移
，砂岩的光譜反射率呈一定的變化規律。圖6則(ze)分(fen)別(bie)分(fen)析(xi)砂(sha)岩(yan)高(gao)光(guang)譜(pu)影(ying)像(xiang)中(zhong)兩(liang)個(ge)像(xiang)元(yuan)光(guang)譜(pu)反(fan)射(she)率(lv)與(yu)水(shui)滲(shen)入(ru)時(shi)間(jian)的(de)相(xiang)關(guan)性(xing)，其(qi)中(zhong)左(zuo)邊(bian)為(wei)紅(hong)色(se)像(xiang)元(yuan)點(dian)，右(you)邊(bian)為(wei)綠(lv)色(se)像(xiang)元(yuan)點(dian)的(de)相(xiang)關(guan)性(xing)分(fen)析(xi)圖(tu)。從(cong)圖(tu)6可(ke)知(zhi)，砂(sha)岩(yan)高(gao)光(guang)譜(pu)上(shang)的(de)不(bu)同(tong)像(xiang)元(yuan)點(dian)
，其(qi)波(bo)段(duan)相(xiang)關(guan)性(xing)與(yu)水(shui)滲(shen)入(ru)有(you)著(zhe)顯(xian)著(zhu)差(cha)異(yi)
，但(dan)是(shi)相(xiang)關(guan)性(xing)的(de)變(bian)化(hua)趨(qu)勢(shi)及(ji)最(zui)相(xiang)關(guan)波(bo)段(duan)的(de)範(fan)圍(wei)與(yu)相(xiang)關(guan)係(xi)數(shu)等(deng)相(xiang)似(si)。

圖6砂岩光譜反射率與水滲入岩石時間的相關性分析（左為紅色像元，右為綠色像元）

　　4、礫岩光譜反射率與水滲入岩石時間的相關性分析

　　圖5分別列舉了礫岩高光譜影像中兩像元點隨水滲入時間的光譜反射率變化規律，從圖5可知
，隨著水滲入時間的推移，雖然礫岩的像元點光譜反射率高低不太相似
，但是呈現變化規律一致。圖7則(ze)分(fen)別(bie)分(fen)析(xi)礫(li)岩(yan)高(gao)光(guang)譜(pu)影(ying)像(xiang)中(zhong)兩(liang)個(ge)像(xiang)元(yuan)光(guang)譜(pu)反(fan)射(she)率(lv)與(yu)水(shui)滲(shen)入(ru)時(shi)間(jian)的(de)相(xiang)關(guan)性(xing)
，其(qi)中(zhong)左(zuo)邊(bian)為(wei)紅(hong)色(se)像(xiang)元(yuan)點(dian)，右(you)邊(bian)為(wei)綠(lv)色(se)像(xiang)元(yuan)點(dian)的(de)相(xiang)關(guan)性(xing)分(fen)析(xi)圖(tu)。從(cong)圖(tu)7kezhi，liyangaoguangpushangdebutongxiangyuandian
，qiboduanxiangguanxingyushuishenruyouzhexianzhuchayi，danshixiangguanxingdebianhuaqushijizuixiangguanboduandefanweiyuxiangguanxishudengxiangsi。

圖7礫岩光譜反射率與水滲入岩石時間的相關性分析（左為紅色像元
，右為綠色像元）

　　5
、砂岩隨水滲入岩石時間的光譜特征波長選擇及模型構建

　　本研究利用Matlab編寫的SPA來提取砂岩基於水滲入時間變化的特征波長，以減少變量的輸入個數簡化後期建模過程。SPA運算過程所設定的特征個數為3-10，特征波個數根據RMSE確定，如圖8的Figure1所示，RMSE越小，水滲入時間的估算越精確，Figure2為確定特征波段的波段位置。表2分別列出了基於連續投影算法的砂岩隨水侵時間的特征波長的波段位置

，越靠前，其特征波段重要性越高
，特征波段分別為1552.65、1658.47、1698.22
、903.28 nm。

圖8砂岩隨水滲入岩石時間的光譜特征波長選擇參數設置及特征波長確定

表2 基於連續投影算法的砂岩隨水侵時間的特征波長

　　根據選取的砂岩特征波段運用逐步多元回歸分析SMR構建模型並運用獨立的數據對模型進行檢驗
，實測值與預測值之間的PRESS為207.49，RMSEP為4.16 min，SDV為4.13
，BIAS為-1.30，相關係數r為0.997.

圖9 基於特征波段構建模型及檢驗效果圖

　　6、礫岩隨水滲入岩石時間的光譜特征波長選擇及模型構建

　　與砂岩研究相似，SPA運算過程所設定的特征個數為3-10，特征波個數根據RMSE確定，如圖10的Figure1所示，RMSE越小，水滲入時間的估算越精確，Figure2為確定特征波段的波段位置。表3分別列出了基於連續投影算法的礫岩隨水侵時間的特征波長的波段位置，越靠前，其特征波段重要性越高，特征波段分別為1381.20、1430.59、909.80
、1698.22
、903.28 nm。

圖10礫岩隨水滲入岩石時間的光譜特征波長選擇參數設置及特征波長確定

表3 基於連續投影算法的礫岩隨水侵時間的特征波長

　　根據選取的礫岩特征波段運用逐步多元回歸分析SMR構建模型並運用獨立的數據對模型進行檢驗
，實測值與預測值之間的PRESS為1672.61
，RMSEP為13.63 min

，SDV為13.67，BIAS為-4.46
，相關係數r為0.984。

圖11 基於特征波段構建模型及檢驗效果圖

　　五、討論

　　本研究運用單個像元點來分析砂岩
、礫li岩yan受shou水shui滲shen入ru後hou，其qi光guang譜pu反fan射she率lv的de變bian化hua趨qu勢shi
，並bing運yun用yong連lian續xu投tou影ying算suan法fa提ti取qu了le特te征zheng波bo長chang，同tong時shi利li用yong逐zhu步bu多duo元yuan回hui歸gui分fen析xi構gou建jian光guang譜pu與yu水shui滲shen入ru時shi間jian的de模mo型xing，結jie果guo表biao明ming，無wu論lun是shi砂sha岩yan或huo者zhe礫li岩yan，所suo構gou建jian的de模mo型xing精jing度du都dou較jiao高gao
，其qi中zhong砂sha岩yan的deRMSP為RMSEP為4.16 min
，相關係數r為0.997
；礫岩的RMSEP為13.63 min，相關係數r為0.984。

　　研究者可根據自己的研究方向，針對局部（整體）砂岩或者礫岩，研究其隨水滲入時間光譜反射率的變化規律。除了逐步多元回歸分析算法
，研究者還可以運用光譜指數、主成分分析、偏最小二乘法、人工神經網絡、支持向量機
、方法有相關係數法、方差分析
、逐步多元線性回歸
、粒子群算法
、反向區間偏最小二乘法和連續投影算法等，研究者均可嚐試分析。常用的數據處理軟件有ENVI、ERDAS、MATLAB

、excel等
，研究者往往不可能隻通過一個軟件就把數據分析的完全徹底
，需要綜合幾個軟件進行分析

，而且要求研究者有較好的編程能力。