一種高光譜成像光譜儀光譜定標方法

撖芃芃

（中國科學院長春光學精密機械與物理研究所，吉林 長春 130033）

0 引 言

高光譜成像光譜儀是一種“圖譜合一”的光學遙感儀器，其作用是獲取地球目標的詳細光譜景像，在陸地、海洋的輻射信息及大氣等方面的監(jiān)測有較多的應用［1－5］。高光譜成像光譜儀的特點是成像技術(shù)和光譜技術(shù)有機地結(jié)合，能夠獲取大量的窄波段連續(xù)光譜圖像數(shù)據(jù)。成像光譜儀的工作原理如圖1所示。

圖1 成像光譜儀工作原理

前端成像鏡頭將來自目標的光信息收集并成像于像面處，像面的任意一點均包含目標對應空間點的光譜信息。通過一定的方式對第二維空間圖像進行掃描，即可采集到由目標二維空間圖像信息及其各點光譜信息所構(gòu)成的三維數(shù)據(jù)立方［6－9］。

隨著成像光譜儀技術(shù)的發(fā)展，對成像光譜儀定標技術(shù)提出了較為嚴格的要求，尤其是對光譜定標的準確度，光譜定標結(jié)果的準確性直接影響到成像光譜儀光譜分析的結(jié)果。文中正是在這樣一個背景下提出了一種高光譜成像光譜儀的光譜定標方法，并搭建了一套光譜定標系統(tǒng)對某種高光譜成像光譜儀進行光譜定標，該系統(tǒng)采用高斯曲線擬合的方法確定中心波長的位置，提取光譜維方向的半高寬作為光譜帶寬，同時采用最小二乘法進行全波段光譜定標，經(jīng)過擬合計算得到該光譜定標方法的標準差為0.23nm，滿足應用要求。

1 高光譜成像光譜儀光譜定標原理

根據(jù)高光譜成像光譜儀的工作原理可知，條帶目標經(jīng)前端成像鏡頭成像高光譜成像光譜儀入射狹縫處，后經(jīng)成像光譜儀色散分光并成像，從而得到按波長排列的條帶目標的單色像所組成的二維圖譜，其中與波長排列方向相垂直的方向稱為空間維，而波長排列方向稱為光譜維。采用面陣探測器接收高光譜成像光譜儀所輸出的二維圖譜時，須將探測器的行像元和列像元與二維圖譜的空間維和光譜維進行校準，以保證二維圖譜信息的準確提取。因此，探測器每一行像元即可采集到入射狹縫不同波長的單色像，則稱采集同一波長單色像的每一行像元為一個光譜通道；而探測器每一列像元所采集到的不同波長的光信息均來自于同一空間點，則稱采集同一空間點不同波長光信息的每一列像元為一個空間通道。綜上所述，探測器的像元P（i，j）所輸出的亮度值可表示為［10－12］：

式中：Eλ——光源在波長λ處的輻射照度；

ρ（λ）——被探測目標的光譜反射率；

τa（λ）——大氣的光譜透過率；

wIFOV——高光譜成像光譜儀系統(tǒng)的瞬時視場；

A——高光譜成像光譜儀系統(tǒng)的有效光學面積；

τ0（λ）——高光譜成像光譜儀系統(tǒng)的光譜透過率；

Ss（λ）——高光譜成像光譜儀系統(tǒng)的傳遞函數(shù)；

Rd（λ）——高光譜成像光譜儀探測器的光譜響應率；

Re（λ）——接收系統(tǒng)的電子學光譜響應率。

對于式（1）而言，取決于高光譜成像光譜儀自身的物理量僅有τ0（λ）、Ss（λ）、Rd（λ）、Re（λ），由它們所構(gòu)建的物理量則高光譜成像光譜儀像元P（i，j）的相對光譜響應函數(shù)Si，j（λ），即：

在理想情況下，系統(tǒng)各部分的光譜響應是均勻的，即τ0（λ）、Ss（λ）、Rd（λ）、Re（λ）為不隨波長變化的常數(shù)。當入射狹縫無限窄時，則像元P（i，j）光譜響應函數(shù)的歸一化結(jié)果可表示為矩形函數(shù)，如圖2所示。

圖2 理想光譜響應函數(shù)曲線

即：

然而，在實際工程應用中，入射狹縫具有一定寬度，且τ0（λ）、Ss（λ）、Rd（λ）、Re（λ）并非常數(shù)，這使得式（3）矩形函數(shù)與實際光譜響應函數(shù)不符。經(jīng)實踐證明，成像光譜儀的實際光譜響應函數(shù)與高斯函數(shù)較為近似［13－15］，即：

此外，采用高斯函數(shù)的半高寬BFWHM表征成像光譜儀的光譜帶寬δλ與實際應用相符，即：

某待測高光譜成像光譜儀像元P（82，241）的實際光譜響應數(shù)據(jù)與高斯擬合曲線對比，實際光譜響應函數(shù)曲線如圖3所示。

圖3 實際光譜響應函數(shù)曲線

從圖中可以看出，高斯函數(shù)能夠準確地表征實際光譜響應函數(shù)。

綜上所述，采用高斯函數(shù)對高光譜成像光譜儀各像元的實際光譜響應函數(shù)進行擬合，即可得到描述系統(tǒng)光譜性能的參量，即各像元的中心波長λc、光譜帶寬δλ（即光譜分辨率）、光譜通道間的光譜取樣間隔Δλ以及系統(tǒng)的工作波長范圍。其中，高斯擬合曲線的峰值波長和半高寬即為相應像元的中心波長和光譜帶寬；相鄰光譜通道的中心波長差即為光譜取樣間隔；最大光譜通道的中心波長與最小光譜通道的中心波長差即為系統(tǒng)的工作波長范圍。

2 高光譜成像光譜儀光譜定標系統(tǒng)設計

光譜定標的任務就是確定各個光譜通道的光譜響應函數(shù)，從而得到描述系統(tǒng)光譜性能的參量。根據(jù)高光譜成像光譜儀的工作原理和成像系統(tǒng)物像共軛原理可知，若目標為單色點物，則探測器某一像元必然與之共軛，共軛像元所在的行位置取決于單色點物的波長，列位置取決于單色點物所在的視場。由此可知，若要獲取像元P（i，j）的光譜響應函數(shù)，首先調(diào)整單色點物所在視場，令其共軛像元落于探測器第j列；然后，連續(xù)改變點物的波長，同時讀取像元P（i，j）的亮度值，即可得到該像元的實際光譜響應數(shù)據(jù)；最后，通過高斯擬合即可得到該像元的光譜響應函數(shù)，進而求得其中心波長λc、光譜帶寬δλ。重復上述方法，即可獲得各像元的光譜響應函數(shù)。光譜定標過程如圖4所示。

以上述光譜定標原理為依據(jù)設計高光譜成像光譜儀的光譜定標系統(tǒng)，令可調(diào)的連續(xù)光源發(fā)出的光聚焦于單色儀入射狹縫上，經(jīng)單色儀色散分光后，于出射狹縫輸出單色光，出射狹縫即為單色點物，出射狹縫經(jīng)平行光管變?yōu)闊o窮遠單色點物后被高光譜成像光譜儀探測接收。調(diào)整視場角令像點位于探測器第j列，控制單色儀逐漸增大輸出波長，令像點較為接近像元P（i，j），且保證像元未接收到像點信號。然后，以步進方式控制單色儀進行波長掃描，波長每遞進一次，探測器像元P（i，j）就讀取一次亮度值，直至像點掃過像元P（i，j），即得到了像元P（i，j）的實際光譜響應數(shù)據(jù)。通過高斯擬合，從而獲得該像元的光譜響應函數(shù)。因此，高光譜成像光譜儀的光譜定標系統(tǒng)由定標光源、單色儀、平行光管、待測高光譜成像光譜儀組成，光譜定標系統(tǒng)的實物圖如圖5所示。

圖4 光譜定標過程示意圖

圖5 高光譜成像光譜儀光譜定標系統(tǒng)實物圖

3 高光譜定標數(shù)據(jù)分析

測試的某型高光譜成像光譜儀使用時進行了像素合并以提高系統(tǒng)信噪比，即4個像素合并為一個像元，因此，該成像光譜儀具有480個空間通道，270個光譜通道，波長范圍為400～800nm。圖5中的光譜定標系統(tǒng)開發(fā)了適用于高光譜成像光譜儀的光譜定標軟件，軟件界面如圖6所示。

圖6 高光譜成像光譜儀光譜定標軟件

圖中左上部分區(qū)域顯示CCD像面的實時圖像，通過實時圖像可以獲取單色光斑的大小、形狀、亮度及坐標位置等信息。中上部分區(qū)域為檢測參數(shù)設置區(qū)域，首先設置CCD探測器取樣次數(shù)；然后設置所要監(jiān)測的像元的序號；最后設置波長參數(shù)，使記錄的波長信息與單色儀輸出的波長信息相匹配，包括定標起始波長、定標步進波長及當前定標波長。右上部分區(qū)域為CCD參數(shù)設置區(qū)域，通過改變幀頻和積分時間，控制采集到的單色光斑的亮度，確保定標數(shù)據(jù)具有良好的信噪比且避免飽和。下方區(qū)域為光譜曲線監(jiān)視區(qū)域，隨著波長的掃描動態(tài)地繪制出被監(jiān)測像元的光譜響應曲線。

對像元P（82，241）進行光譜定標時，光譜定標軟件參數(shù)設置見圖6，所列像元P（82，241）采集到的光譜定標數(shù)據(jù)見表1。

表1 像元P（82，241）的光譜定標數(shù)據(jù)

采用最小二乘法對表1所列數(shù)據(jù)進行高斯擬合，得到像元P（82，241）的光譜響應曲線（見圖3）。

高斯函數(shù)能夠準確地表征高光譜成像光譜儀的實際光譜響應曲線，進而求得該像元的中心波長為503nm，光譜帶寬為1.9nm。

對中心視場所對應的241空間通道的各光譜通道進行光譜定標，部分光譜通道的光譜定標結(jié)果見表2。

表2 部分光譜通道的光譜定標結(jié)果 nm

采用最小二乘法將表2所列的各光譜通道與其中心波長的關(guān)系進行擬合，得到線性、二次、三次、四次擬合多項式。經(jīng)計算驗證，二次擬合多項式的標準誤差最小，僅為0.23nm，二次擬合多項式為：

式中：λc——各光譜通道的中心波長；

N——光譜通道系數(shù)。

擬合曲線如圖7所示。

4 結(jié) 語

提出的高光譜成像光譜儀的光譜定標方法主要是針對高光譜成像光譜儀較為嚴格的光譜定標，根據(jù)該方法搭建了一套光譜定標系統(tǒng)，并通過對某型高光譜成像光譜儀進行光譜定標，驗證了該方法的準確性，經(jīng)過擬合計算得到該光譜定標方法的標準差為0.23nm，滿足應用要求。該方法的提出為我國成像光譜技術(shù)及標定方法的研究提供了理論基礎，通過實際的測量取得了理想的結(jié)果，為下一步成像光譜定標技術(shù)的深入研究奠定了基礎。

圖7 241空間通道的各光譜通道的中心波長及其擬合曲線

［1］ Goetz A F H，Vane G，Solmon J E，et al.Imaging spectrometry for earth remote sensing［J］.Science，1985，228：1147－1153.

［2］ 王建宇.機載成像光譜儀性能分析和發(fā)展現(xiàn)狀：遙感新進展與發(fā)展戰(zhàn)略［M］.北京：中國科學技術(shù)出版社，1996：99－107.

［3］ 鄭玉權(quán).小型Offner光譜成像系統(tǒng)的設計［J］.光學精密工程，2005，13（6）：650－657.

［4］ Green R O，Chrlen T G，Niei Son I J，et al.Airborne visIb1e／infrared imaging spectrometer（AVIRIS）：recent improvements to the sensor and data facility［J］.SPIE，1993，1937：180－190.

［5］ Rickard L J，Basedow R W，Zai Ewski E F，et al.HYDICE：an airborne system for hyperspectral imaging［J］.SPIE，1993，1937：173－179.

［6］ 齊向東，撖芃芃，潘明忠，等.凸面光柵成像光譜儀的光譜定標［J］.光學精密工程，2011，19（12）：2870－2876.

［7］ 劉倩倩，鄭玉權(quán).超高分辨率光譜定標技術(shù)發(fā)展概況［J］.中國光學，2012，5（6）：566－577.

［8］ 金輝，姜會林，鄭玉權(quán)，等.高光譜遙感器的光譜定標［J］.發(fā)光學報，2013，34（2）：235－239.

［9］ 楊增鵬.棱鏡－光柵型高光譜成像光譜儀優(yōu)化設計及集成技術(shù)研究［D］.北京：中國科學院大學，2015.

［10］ Babey S K，Anger C D.Compact airborne spectrographic imager（CASI）：aprogress review［J］.SPIE，1993，1937：152－163.

［11］ Davis C O，Kappus M E，Bowles J H，et al.Calibration，characterization，and first results with the Ocean PHILIS hyperspectral imager［J］.SPIE，1999，3753：160－168.

［12］ Folkman M A，Peari M J，Liao I B，et al.EO－1／Hyperion hyperspectral imager design，development，characterization and calibration［J］.SPIE，2001，4151：40－51.

［13］ Wiison T H，Davis C O.Hyperspectral remote sensing technology（HRST）program and the Naval Earth Map Observer（NEM）satellite［J］.SPIE，1998，3437：2－11.

［14］ Fietcher P A.Image acquisition planning for the CHRIS sensor onboard PROBA［J］.SPIE，2004，5546：141－150.

［15］ Thomas G C，Robert O G，Michael L.Accuracy of the spectral and radiometric laboratory calibration of the airborne visible／infrared imaging spectrometer（AVIRIS）［J］.SPIE，1990，1298：37－49.