任夢(mèng)如，劉洪成，葉發(fā)旺，張 川，武 鼎，李瀚波

（核工業(yè)北京地質(zhì)研究院 遙感信息與圖像分析技術(shù)國(guó)家級(jí)重點(diǎn)實(shí)驗(yàn)室，北京 100029）

高光譜成像儀（hyperspectral imager）又稱成像光譜儀，具有 “圖譜合一”的特點(diǎn)，即在獲取目標(biāo)地物的窄波段連續(xù)光譜圖像的同時(shí)，也獲取每個(gè)像元幾乎連續(xù)的光譜曲線。經(jīng)過精確的輻射定標(biāo)后，將獲取的目標(biāo)地物光譜輻射值轉(zhuǎn)換成反射率或輻亮度信息?；谀繕?biāo)的光譜曲線，可以實(shí)現(xiàn)對(duì)地物組成成分的識(shí)別。經(jīng)過幾十年的發(fā)展，高光譜成像技術(shù)愈發(fā)成熟與先進(jìn)，光譜分辨率可以達(dá)到納米級(jí)，被廣泛應(yīng)用于社會(huì)各個(gè)領(lǐng)域，如農(nóng)業(yè)生產(chǎn)管理、生態(tài)環(huán)境監(jiān)測(cè)、地質(zhì)礦產(chǎn)勘查等［1-4］。按照巖礦光譜形成基本原理和不同蝕變礦物對(duì)光譜吸收機(jī)制的不同，目前在可見光-短波紅外（400～2 500 nm）譜段，成像光譜技術(shù)可識(shí)別30～40種礦物的種類、礦物組成成分以及礦物相對(duì)含量。比如與成礦作用密切相關(guān)的蝕變礦物的識(shí)別，對(duì)于尋找地質(zhì)信息所蘊(yùn)涵的礦物共生組合規(guī)律，對(duì)圈定礦化蝕變帶，減少和降低傳統(tǒng)地質(zhì)工作中可能因地質(zhì)路線、采樣部署不合理以及肉眼的局限性而遺漏的一些主要地質(zhì)現(xiàn)象，以及追索礦化熱液蝕變中心和圈定找礦靶區(qū)等都有重要作用［4-5］。

目前，國(guó)外高光譜成像儀產(chǎn)品已逐漸標(biāo)準(zhǔn)化，機(jī)載光譜成像技術(shù)也發(fā)展到商業(yè)運(yùn)行階段，具有代表性的高光譜成像儀有加拿大的 CASI/SASI、 美 國(guó) 的 AVIRIS、 德 國(guó) 的ROSIS、澳大利亞的HyMap等，但體積和重量偏大，不適于無人機(jī)搭載，無法進(jìn)行近距離高光譜成像測(cè)量。而國(guó)內(nèi)起步相對(duì)較晚，光譜成像技術(shù)不夠成熟，且受元件制造技術(shù)水平的限制， 成型的產(chǎn)品還 比 較 少［6-7］。 而Fe2＋、Fe3＋、Mn3＋等陽(yáng)離子礦物與含羥基 O—H、硫基S—H等陰離子基團(tuán)的礦物在短波紅外（1 000～2 500 nm）譜段具有明顯的特征譜帶與吸收特征，可以有效識(shí)別綠泥石、綠簾石、蛇紋石等礦物［8-9］。因此開展輕小型短波紅外高光譜成像儀研制對(duì)野外地質(zhì)礦產(chǎn)勘查具有重要意義，可以有效促進(jìn)蝕變礦物精細(xì)化識(shí)別水平的進(jìn)步。

鑒于此，筆者采用模塊化光譜儀制造方式研制了一款短波紅外高光譜成像儀，具有光譜分辨率高、結(jié)構(gòu)緊湊、輕小型化、便攜型性強(qiáng)等優(yōu)點(diǎn)，適用于無人機(jī)載成像，可以較低的成本獲得熱液蝕變的詳細(xì)礦物信息，在鈾礦地質(zhì)勘查領(lǐng)域有廣闊的應(yīng)用前景。

1 短波紅外高光譜成像儀光學(xué)設(shè)計(jì)

1.1 高光譜成像儀工作原理

高光譜成像儀主要由物鏡、狹縫、準(zhǔn)直鏡、透射式光柵、棱鏡、紅外面陣探測(cè)器和數(shù)據(jù)存儲(chǔ)系統(tǒng)組成（圖1）。光線經(jīng)過成像物鏡射入并會(huì)聚到入射狹縫，經(jīng)準(zhǔn)直鏡后入射到與透射式光柵和棱鏡組成的分光模塊，復(fù)色光線經(jīng)過色散部件后變?yōu)閱紊饩€，此時(shí)地物目標(biāo)成為一系列線陣列像元，之后通過成像鏡的會(huì)聚作用，成像在短波紅外焦平面探測(cè)器上成像，最后經(jīng)過光電轉(zhuǎn)換由數(shù)據(jù)存儲(chǔ)系統(tǒng)接收，從而獲得目標(biāo)的二維空間信息。

圖1 推掃式高光譜成像儀原理圖Fig.1 The schematic diagram of push broom hyperspectral imager

空間維在記錄目標(biāo)的二維空間信息的同時(shí)由光譜維獲取地物目標(biāo)的光譜信息，最終得到完整的數(shù)據(jù)立方體，如圖2所示。

1.2 光學(xué)系統(tǒng)設(shè)計(jì)

物鏡的設(shè)計(jì)遵循儀器結(jié)構(gòu)緊湊、輕小型化、大視場(chǎng)、低F數(shù)、與光譜儀光瞳匹配的原則，兼顧成本和裝調(diào)難度。由于地物在光譜范圍為1 000～2 500 nm的短波紅外區(qū)間反射能量較低，因此物鏡光譜透過率要高，采用透射式鏡組結(jié)構(gòu)，像方遠(yuǎn)心設(shè)計(jì)。為能夠?qū)δ繕?biāo)精細(xì)化成像，鏡頭焦距盡量擴(kuò)大且能夠在短距離情況下成像［10］。

分光模塊采用棱鏡-光柵-棱鏡（PGP）型透射式色散系統(tǒng)，利用光柵和棱鏡構(gòu)成復(fù)合分光器件，解決單一分光器件存在的譜線彎曲和色散非線性問題，達(dá)到簡(jiǎn)化系統(tǒng)結(jié)構(gòu)，提高光譜分辨率，降低數(shù)據(jù)處理復(fù)雜度的目的［11-12］。設(shè)計(jì)結(jié)合了棱鏡和光柵色散分光技術(shù)的優(yōu)點(diǎn)，具有光譜分辨率高、光學(xué)效率高、光譜線性度好、重量輕便等優(yōu)點(diǎn)。PGP型透射式色散系統(tǒng)為同軸系統(tǒng)，類似管狀結(jié)構(gòu)，有利于光校操作，系統(tǒng)的結(jié)構(gòu)設(shè)計(jì)簡(jiǎn)單，鏡頭接口和探測(cè)器接口易搭配。光譜儀光學(xué)系統(tǒng)結(jié)構(gòu)圖如圖3所示。

輕小型短波紅外光譜儀技術(shù)參數(shù)如表1所示。

1.3 光學(xué)系統(tǒng)評(píng)價(jià)

圖2 推掃式高光譜成像儀數(shù)據(jù)立方體Fig.2 Data cube of push broom hyperspectral imager

圖3 輕小型短波紅外光譜儀光學(xué)系統(tǒng)結(jié)構(gòu)圖Fig.3 Optical system structure of compact SWIR hyperspectral imager

表1 輕小型短波紅外光譜儀技術(shù)參數(shù)Table 1 Technical parameters of compact SWIR hyperspectral imager

通過對(duì)光學(xué)系統(tǒng)的性能評(píng)價(jià)可以判斷光學(xué)系統(tǒng)的各個(gè)指標(biāo)是否符合設(shè)計(jì)要求。通過點(diǎn)列圖不僅可以衡量光學(xué)系統(tǒng)的成像質(zhì)量，還能夠判斷光譜儀的光譜分辨率［6］。從望遠(yuǎn)鏡系統(tǒng)點(diǎn)列圖（圖4）可以看出，在一個(gè)像元內(nèi)，點(diǎn)列圖形狀比較規(guī)整，且相對(duì)集中，滿足設(shè)計(jì)要求。適配焦平面探測(cè)器像元為30 μm，而全視場(chǎng)內(nèi)物鏡的RMS半徑為18 μm，可以保證入射光線內(nèi)的所有景物均可以在探測(cè)器平面上成像。

望遠(yuǎn)鏡系統(tǒng)的光學(xué)調(diào)制傳遞函數(shù)（MTF）曲線如圖5所示，系統(tǒng)在不同中心波長(zhǎng)、不同視場(chǎng)條件下，奈奎斯特頻率16.7 line/mm處的MTF值均高于0.85，系統(tǒng)的實(shí)際光譜分辨率優(yōu)于8 nm，滿足設(shè)計(jì)要求。

從分光模塊在1.00、1.75和2.50 μm處的單色點(diǎn)列圖（圖6）看出，不同波長(zhǎng)不同視場(chǎng)下點(diǎn)列圖比較集中，分布在一個(gè)像元之內(nèi)。彌散斑RMS半徑均小于16 μm，小于短波紅外焦平面探測(cè)器的像元尺寸，滿足光譜分辨率要求。

圖4 望遠(yuǎn)鏡系統(tǒng)點(diǎn)列圖Fig.4 Beam path of the optical system

圖5 望遠(yuǎn)鏡系統(tǒng)的光學(xué)調(diào)制傳遞函數(shù)曲線Fig.5 MTF curve of the telescopic system

圖6 單色點(diǎn)列圖Fig.6 The monochromatic dot plot

2 近景高光譜成像數(shù)據(jù)獲取

2.1 近景高光譜成像空間技術(shù)指標(biāo)

近景高光譜成像儀設(shè)備通過推掃式移動(dòng)，空間成像幾何關(guān)系如圖7所示，瞬時(shí)視場(chǎng)與探測(cè)器單個(gè)像元尺寸相關(guān)，總視場(chǎng)與掃描行的全部像元數(shù)相關(guān)，根據(jù)探測(cè)器的焦平面尺寸，光譜維設(shè)置為256維，空間維設(shè)置為320維。

圖7 推掃式高光譜成像儀空間成像幾何關(guān)系Fig.7 The geometry signature of push broom hyperspectral imager

假定焦平面探測(cè)器的像元為正方形，在不考慮光學(xué)系統(tǒng)像差的情況下，瞬時(shí)視場(chǎng)是指探測(cè)器單個(gè)像元對(duì)應(yīng)的視場(chǎng)，與像元大小和光學(xué)系統(tǒng)焦距相關(guān)。

一般情況下，由于瞬時(shí)視場(chǎng)較小，正切值可用弧度代替：

空間分辨率GR可表示為：

總視場(chǎng)由探測(cè)器空間維像元總數(shù)以及光學(xué)系統(tǒng)焦距決定，即：

研制的近景高光譜成像儀光學(xué)系統(tǒng)焦距f為 30 mm， 像 元 尺 寸 為 30 μm，IFOV＝0.058°，F(xiàn)OV＝18.2°。在近景高光譜成像條件下，當(dāng)相對(duì)距離為1.0 m時(shí)，根據(jù)式（3），計(jì)算巖石樣品的空間分辨率為1 mm。同理，相對(duì)距離為2.0 m時(shí)，巖石樣品的空間分辨率為2 mm。

2.2 數(shù)據(jù)存儲(chǔ)

高光譜成像儀在推掃獲取數(shù)據(jù)過程中，影像空間維（spatial dimensions）依次記錄每一個(gè)波段的掃描行的，直到覆蓋全部的256維。圖像的長(zhǎng)度（lines）等于探測(cè)器掃描單張影像的數(shù)量，總數(shù)據(jù)量大小除以每張影像的數(shù)據(jù)量大小來計(jì)算得到。數(shù)據(jù)存儲(chǔ)采用BIL格式，并配置*.hdr頭文件。

2.3 輻射校正和光譜重建

輻射校正的主要目的是將高光譜影像的像素值轉(zhuǎn)換為地物的反射率值；而輻射定標(biāo)則是確定焦平面探測(cè)單元獲取的影像灰度值與其接收到的電磁波信號(hào)之間的數(shù)學(xué)關(guān)系，并以此來確定不同像元的光譜特征曲線。

由于無人機(jī)等搭載的高光譜成像儀，探測(cè)距離通常較近，大氣對(duì)高光譜成像儀的輻照度影響較小，可忽略不計(jì)。在進(jìn)行數(shù)據(jù)采集時(shí)，可以在同一視場(chǎng)內(nèi)放置反射率已知的標(biāo)準(zhǔn)反射率板白板作為參照，對(duì)獲得的數(shù)據(jù)進(jìn)行相對(duì)輻射校正和光譜重建，利用下式可以直接獲得被測(cè)目標(biāo)的反射率：

式（5）、 （6）中： ρobj（m）—待測(cè)目標(biāo)的第 m 波段反射率；ρstd（m）—標(biāo)準(zhǔn)反射率板的第m波段反射率；DNstd（m）—測(cè)量標(biāo)準(zhǔn)反射率板得到的第m波段DN值；DNobj（m）—測(cè)量待測(cè)目標(biāo)得到的第m波段DN值。

這里，我們可以反射率為0.99的標(biāo)準(zhǔn)白板作為標(biāo)準(zhǔn)反射率板，根據(jù)式（6）計(jì)算每個(gè)像元的反射率，便可獲取近景高光譜成像數(shù)據(jù)和光譜曲線，然后利用標(biāo)準(zhǔn)反射板進(jìn)行校正處理，最后將校正后的發(fā)射率曲線與ENVI光譜庫(kù)中的對(duì)應(yīng)標(biāo)準(zhǔn)曲線對(duì)比，就可以較準(zhǔn)確地對(duì)目標(biāo)物進(jìn)行識(shí)別。

3 近景高光譜成像礦物信息識(shí)別

3.1 礦物識(shí)別算法

基于地物光譜反射率曲線，利用標(biāo)準(zhǔn)光譜庫(kù)或地面光譜儀同步測(cè)量的光譜曲線作為已知光譜數(shù)據(jù)，采用光譜匹配的方法，通過巖石光譜與參考光譜的相似性度量或巖石光譜與光譜庫(kù)中標(biāo)準(zhǔn)光譜曲線的比較，計(jì)算成像像元光譜與參考光譜在光譜空間中的差異性來識(shí)別近景高光譜影像中的蝕變礦物類型［8，13］。

巖礦的光譜曲線包含一系列特征吸收譜帶，且每一個(gè)特征吸收譜帶或譜帶組合與巖石內(nèi)部微粒的位置屬性存在一定的對(duì)應(yīng)關(guān)系，如：Fe2＋在1 000 nm附近產(chǎn)生特征譜帶，CO32-在2 300 nm處產(chǎn)生特征吸收，Al-OH在2 210 nm處產(chǎn)生特征吸收等。根據(jù)不同礦物在短波紅外波段區(qū)間的吸收峰位（表2），即可以區(qū)分不同含羥基或陰性離子團(tuán)的礦物。利用陰離子的特征譜帶或其譜帶組合可以對(duì)礦物進(jìn)行識(shí)別。

表2 部分礦物的特征吸收峰位置Table 2 Infrared spectral absorption position of source minerals

3.2 典型礦物特征信息提取

在進(jìn)行數(shù)據(jù)采集前，首先用ASD對(duì)實(shí)驗(yàn)樣品進(jìn)行光譜測(cè)量，并對(duì)礦物含量進(jìn)行判斷。利用研制的成像儀對(duì)典型礦物含量高的巖石面進(jìn)行掃描，采用光譜角礦物識(shí)別算法進(jìn)行礦物識(shí)別，并獲得了對(duì)應(yīng)礦物的光譜曲線。之后對(duì)獲得結(jié)進(jìn)行處理，并與標(biāo)準(zhǔn)光譜庫(kù)對(duì)應(yīng)曲線進(jìn)行對(duì)比（圖 8）。

圖8 典型礦物精細(xì)化識(shí)別影像及對(duì)應(yīng)光譜曲線Fig.8 Refined identification image and corresponding spectral curve of typical minerals

從實(shí)測(cè)光譜曲線上看出，石膏在短波紅外波段有4個(gè)明顯的吸收峰（1 448、1 748、1 940、2 216 nm），方解石在2 335 nm處有明顯吸收特征，與表2中吸收峰位置一致；光譜形態(tài)與標(biāo)準(zhǔn)光譜庫(kù)中對(duì)應(yīng)光譜曲線的吻合性好。表明研制的成像光譜儀可以通過礦物的反射率曲線來進(jìn)行不同礦物的精細(xì)化識(shí)別。

3.3 混合礦物信息識(shí)別

為進(jìn)一步驗(yàn)證光譜儀對(duì)同一巖石內(nèi)部不同礦物識(shí)別效果，在實(shí)驗(yàn)室內(nèi)利用輕小型短波紅外高光譜成像儀近景測(cè)量獲取了某巖石樣品（圖9 a）的高光譜圖像，采用蝕變礦物提取算法進(jìn)行混合礦物識(shí)別，成功對(duì)巖石樣品中蝕變礦物進(jìn)行提取，并可見明顯的孔雀石化和絹云母化的分界面，其分布情況見圖9 b。

輕小型短波紅外高光譜成像儀的成功研制為礦物精細(xì)化識(shí)別工作提供了技術(shù)支持，但部分礦物在短波紅外波段光譜特征不明顯，如赤鐵礦僅在可見光-近紅外波段范圍內(nèi)有明顯吸收特征，因此在使用過程中需要結(jié)合實(shí)際情況；實(shí)驗(yàn)室測(cè)量時(shí)，外界影響較小，儀器姿態(tài)基本無變化，在進(jìn)行無人機(jī)載成像時(shí)外界擾動(dòng)大，對(duì)實(shí)際結(jié)果有一定的影響，需設(shè)計(jì)專用穩(wěn)定平臺(tái)，并進(jìn)一步完善礦物提取方法，以提高礦物識(shí)別的準(zhǔn)確性。

圖9 孔雀石化巖石樣品（a）和精細(xì)化高光譜成像識(shí)別結(jié)果（b）Fig.9 The malachite specimen（a）and refined mineral identification results（b）

4 結(jié)語(yǔ)

采用透射式光柵研制的輕小型短波紅外高光譜成像儀可以獲取較高光譜分辨率和空間分辨率的近景高光譜成像數(shù)據(jù)。通過性能測(cè)試表明本儀器不僅能夠在更高的光譜分辨率和空間分辨率的尺度上對(duì)礦物進(jìn)行精細(xì)化識(shí)別，還可以在一個(gè)影像平面上進(jìn)行混合礦物識(shí)別。此外，由于本光譜儀具有結(jié)構(gòu)緊湊、輕小型化、便攜型性強(qiáng)等優(yōu)點(diǎn)，適用于無人機(jī)載成像，可以提高野外地質(zhì)勘探的工作效率，在地質(zhì)資源勘查領(lǐng)域有廣闊的應(yīng)用前景。