快照式光譜成像技術(shù)綜述

高澤東，高洪興，朱院院，李 杰，郝 群，劉 宇，陳 超，成 剛，曹 杰，孟合民

(1.北京理工大學(xué) 光電學(xué)院，北京 100081；2.西安應(yīng)用光學(xué)研究所，陜西 西安 710065；3.中國人民解放軍32180部隊(duì)，北京 100012)

1 引 言

光譜成像儀通過采集場景的光譜輻射度收集三維(Three dimensional, 3D)數(shù)據(jù)集，稱作數(shù)據(jù)立方體。數(shù)據(jù)立方體包括二維(Two dimensional, 2D)空間信息和一維(One dimensional, 1D)光譜信息[1-2]。通過降維融合顯示[3]、地物分類[4]或目標(biāo)檢測識別[5-6]實(shí)現(xiàn)3D數(shù)據(jù)立方體的應(yīng)用。從數(shù)據(jù)獲取方式而言，光譜成像技術(shù)包括掃描式光譜成像與快照式光譜成像[1]。

掃描式光譜成像是指以時(shí)間序列測量數(shù)據(jù)立方體的1D或2D切片，其獲取3D數(shù)據(jù)立方體通常有3種方式：點(diǎn)掃描：一次獲取1D光譜曲線，通過2D空間掃描，獲得完整3D光譜數(shù)據(jù)立方體；線掃描：一次獲得1D空間信息+1D光譜信息，如帶狹縫的光譜成像方式，常見的有棱鏡或光柵色散型光譜成像方式、空間調(diào)制型光譜成像方式；譜掃描：一次獲得一個(gè)窄帶波段內(nèi)的2D空間信息，如濾光片輪切換方式，AOTF、LCTF等可調(diào)諧濾光成像方式[7]。掃描式光譜成像傳感器在采集有顯著空間變化或光譜變化的動(dòng)態(tài)場景時(shí)，由于需要掃描而無法獲取準(zhǔn)確的數(shù)據(jù)立方體[1]。然而，在實(shí)際應(yīng)用場合，人們希望光譜成像儀無需掃描，能在一次曝光周期內(nèi)獲取完整的數(shù)據(jù)立方體，因此快照式光譜成像技術(shù)被提出[8]。本文對快照式光譜成像技術(shù)進(jìn)行了系統(tǒng)的分析和研究，按照圖像分割、孔徑分割、光路分割和頻率分割法對該技術(shù)進(jìn)行了分類，研究了每種光譜成像技術(shù)的原理，分析了每種技術(shù)的優(yōu)缺點(diǎn)，并進(jìn)行對比及展望。

2 原理與分類

2.1 概念與原理

快照式光譜成像技術(shù)是指在探測器單個(gè)積分周期內(nèi)測量3D數(shù)據(jù)立方體。而成像探測器是2D探測器，單個(gè)積分周期只能獲取2D信息，因此快照式需要將數(shù)據(jù)立方體劃分為多個(gè)2D數(shù)據(jù)以同時(shí)測量，再通過后處理將2D數(shù)據(jù)重新組合為3D立方體。最簡易的快照式光譜成像儀為通道式光譜成像儀(Channel Spectral Imager, CSI)[9-10]。光譜數(shù)據(jù)立方體如圖1所示。

圖1 光譜數(shù)據(jù)立方體示意圖Fig.1 Schematic diagram of spectral data cube

CSI通過多通道并行采集，解決2D成像系統(tǒng)不能一次采集3D數(shù)據(jù)的問題，如圖2所示。每個(gè)通道由濾光片、光學(xué)系統(tǒng)和探測器構(gòu)成，每個(gè)通道配置不同中心波長的窄帶濾光片，一個(gè)通道獲取一個(gè)窄帶光譜圖像，形成3D數(shù)據(jù)立方體。CSI具有原理簡單、易于實(shí)現(xiàn)、相機(jī)可靠的優(yōu)點(diǎn)。但由于多通道并行采集，各通道存在視差，各通道采集的窄帶光譜圖像需要配準(zhǔn)；同時(shí)，一個(gè)成像通道對應(yīng)一個(gè)窄帶波段，光譜通道數(shù)取決于光學(xué)系統(tǒng)與探測器的數(shù)目，因此光譜通道數(shù)量較少，一般為3，4；采用多個(gè)探測器，體積、質(zhì)量均比較大。

圖2 通道式光譜成像原理Fig.2 Schematic diagram of channel spectral imager

由于通道式光譜成像技術(shù)存在光譜通道數(shù)較少、體積和質(zhì)量較大等問題，其他快照式光譜成像技術(shù)得以蓬勃發(fā)展。經(jīng)過近30年的發(fā)展，出現(xiàn)了多種快照式光譜成像技術(shù)，包括積分視場光譜儀(Integral Field Spectrometry，IFS)[8]、圖像映射光譜儀(Image Mapping Spectrometry, IMS)[11]、像素濾光的多光譜探測器(Multispectral Sensor with Per-pixel Filters, MS-PPF)[12]、基于光場結(jié)構(gòu)的光譜成像儀(Spectral Imager using Light Field Architecture, SI-LF)[13-14]、編碼孔徑快照式光譜成像儀(Coded Aperture Snapshot Spectral Imager, CASSI)[15]、基于復(fù)眼成像原理的多孔徑分塊濾光光譜成像儀(Multispectral Imager with Thin Observation module by Bound Optics, MSI-TOMBO)[16]、單片集成分塊濾光的多光譜探測器(Multispectral Sensor with Integrated Tilted Filters, MS-ITF)[17]、多孔徑漸變?yōu)V光的光譜成像儀(Spectral Imager with Multi-apertures and Linear Varying Filter, SI-MALVF)[18-19]、快照式高光譜成像傅里葉變換(Snapshot Hyperspectral Imaging Fourier Transform, SHIFT)光譜儀[20]、計(jì)算機(jī)層析光譜成像儀(Computed Tomographic Imaging Spectrometry, CTIS)[21]、基于分束器的光譜成像儀(Spectral Imager with Beam Splitting, SI-BS)[22-23]、基于濾光片堆棧的光譜成像儀(Spectral Imager with a Filter Stack, SI-FS)[24]、圖像復(fù)分光譜成像儀(Image-replicating Imaging Spectrometer, IRIS)[25]、多光譜Sagnac干涉儀(Multispectral Sagnac Interferometer, MSI)[26]等。

2.2 技術(shù)分類

快照式光譜成像分類如圖3所示，根據(jù)是否需要重構(gòu)，分為直接測量與間接測量[2]兩種方式。直接測量是指采集的2D圖像經(jīng)過簡單的數(shù)據(jù)重排就可以獲得3D數(shù)據(jù)立方體，而間接測量是指采集的2D圖像需要經(jīng)過算法重構(gòu)才能獲得數(shù)據(jù)立方體。根據(jù)是否需要迭代計(jì)算，將重構(gòu)方法分為直接重構(gòu)與迭代重構(gòu)。

掃描式光譜成像技術(shù)的研究重點(diǎn)在于如何獲取光譜信息，常從分光方式對它進(jìn)行系統(tǒng)研究，有色散型、濾光型及干涉型等[7]。而快照式光譜成像的關(guān)注點(diǎn)是如何將3D數(shù)據(jù)立方分割為多個(gè)2D圖像，根據(jù)分割方法歸納出圖像分割、孔徑分割、光路分割和頻域分割4類。

圖3 快照式光譜成像技術(shù)分類Fig.3 Classification of snapshot spectral imaging technologies

圖4 圖像分割原理Fig.4 Principle for image division

圖像分割(Division of Image, DoI)方法將分割器件置于成像系統(tǒng)的像面上，并將像面上的圖像分成小塊，然后變形、收縮或定向到不同位置，在2D探測器上為光譜維色散留下足夠的空間。圖4為圖像分割的基本原理。值得注意的是：像面圖像分割是DoI方法的特征；DoI的圖像分塊、移動(dòng)留白及分光，相當(dāng)于在2D探測器上實(shí)現(xiàn)了點(diǎn)掃描或線掃描。常見的圖像分割方法有積分視場(IFS)、圖像映射(IMS)、像素濾光(MS-PPF)、基于光場結(jié)構(gòu)的SI-LF和壓縮編碼孔徑(CASSI)。其中IFS，IMS，MS-PPF，SI-LF屬于直接測量的圖像分割方法，而CASSI屬于迭代重構(gòu)的圖像分割方法。

孔徑分割(Division of Aperture, DoA)思想起源于通道式光譜成像。將多個(gè)孔徑置于成像系統(tǒng)的光瞳處，通過多孔徑場景復(fù)制，在2D探測器上得到多孔徑的復(fù)制圖像，如圖5所示。最常見的是在成像系統(tǒng)的光瞳處放置微透鏡陣列實(shí)現(xiàn)入射場景的多孔徑復(fù)制。值得注意的是：孔徑復(fù)制是DoA的特征，去除DoA的前端成像系統(tǒng)，DoA就蛻化成為CSI；DoA的孔徑復(fù)制與分光，相當(dāng)于在2D探測器上實(shí)現(xiàn)了譜掃描。常見的孔徑分割方法有復(fù)眼成像(MSI-TOMBO)、復(fù)制濾光(MS-ITF)、多孔徑漸變?yōu)V光(SI-MALVF)、傅里葉變換(SHIFT)和計(jì)算機(jī)層析(CTIS)。其中，MSI-TOMBO，MS-ITF，SI-MALVF屬于直接測量的孔徑分割方法，而SHIFT屬于直接重構(gòu)的孔徑分割方法，CTIS屬于迭代重構(gòu)的孔徑分割方法。

圖5 孔徑分割原理Fig.5 Principle for aperture division

光路分割(Division of Optical Path, DoOP)方法將入射光分成多個(gè)光路，如圖6所示，然后定位到不同方向，每個(gè)光路分別獲取不同的窄帶光譜圖像，如通過二色分光鏡將一個(gè)光路分成兩個(gè)光路，然后由兩個(gè)探測器進(jìn)行采集。DoOP的特征是將光路分成若干個(gè)子光路，相當(dāng)于在2D探測器上實(shí)現(xiàn)了譜掃描。常見的光路分割方法有基于分束器的SI-BS、基于濾光片堆棧的SI-FS和圖像復(fù)分IRIS，均屬于直接測量的光路分割方法。

圖6 光路分割原理Fig.6 Principle for optical path division

頻域分割(Division of Frequency Domain, DoFD)方法是指在空間域或譜域中復(fù)用具有不同特性的光子，然后在相應(yīng)的頻域中分割得到信號的捕獲方法。常見的頻域分割法是多光譜Sagnac干涉。

3 圖像分割

3.1 切片器積分視場光譜議

1938年，Bowen提出了首個(gè)快照式光譜成像儀：基于切片鏡的積分視場光譜儀(Integral Field Spectrometry with Faceted Mirrors, IFS-M)[8]，其原理如圖7所示。切片鏡置于前置光學(xué)系統(tǒng)的像面，切片鏡上有多個(gè)不同傾斜角度的長條鏡(切片元件)，將像面上的圖像切片成條帶，條帶圖像經(jīng)過折疊鏡反射到不同的光瞳位置。在每個(gè)光瞳后面有準(zhǔn)直鏡、色散器、二次成像鏡和探測器陣列，將每個(gè)瞳孔位置的條帶圖像色散，得到條帶圖像的色散圖像，經(jīng)過重排得到3D光譜數(shù)據(jù)立方體I((x,y,λ)。

圖7 基于切片鏡的積分視場光譜儀原理Fig.7 Schematic diagram of integral field spectrometry with faceted mirrors

由于切片鏡的切片元件是全反射方法，IFS-M技術(shù)的光能利用率高，非常適合于紅外測量。但卻包含如下缺點(diǎn)：切片元件的數(shù)量決定了一維空間像素?cái)?shù)，切片元件數(shù)量受到限制，使得空間像素?cái)?shù)較少；為了保持瞳孔之間的距離，切片元件的傾斜必須很大，易在切片元件末端產(chǎn)生散焦；多個(gè)探測器導(dǎo)致其體積和質(zhì)量很大。由于存在空間像素少、體積質(zhì)量大等缺點(diǎn)，IFS-M沒有得到廣泛應(yīng)用。

圖8 微切片器的積分視場光譜儀原理Fig.8 Schematic diagram of integral field spectrometry with micro-slicer

基于微切片器的積分視場光譜儀(Integral Field Spectrometry with Micro-slicer，IFS-μ)[27]能夠較好地克服上述問題，以減少光譜采樣為代價(jià)，可采集更多的空間元素。IFS-μ使用與IFS-M相同的切片鏡，但切片元件的平面更大。IFS-μ首先通過變形前置光學(xué)系統(tǒng)將輸入圖像的一個(gè)軸拉伸變形、放大，然后通過微切片器的切片、變形、縮小，使得圖像的各個(gè)空間元素以窄條的形式呈現(xiàn)。這些窄條充分分離，可以在沒有重疊的情況下進(jìn)行光譜色散，如圖8所示。通過微切片技術(shù)可實(shí)現(xiàn)100×100×120的光譜數(shù)據(jù)立方體。盡管如此，IFS-μ卻只適用于空間像素?cái)?shù)要求低、光譜通道數(shù)要求高的場合。

3.2 可變形光纖束積分視場光譜儀

基于可變形光纖束的積分視場光譜儀(Integral Field Spectrometer with Deformable fiber bundle, IFS-F)如圖9所示。采用可變形光纖束，將成像物鏡像面上的2D空間圖像“拉直”為1D線圖像，再經(jīng)過光柵、棱鏡等色散元件進(jìn)行色散，獲得1D線圖像的色散圖像，最后經(jīng)過重排獲得3D光譜數(shù)據(jù)立方體[28]。

圖9 基于可變形光纖束的積分視場光譜儀原理Fig.9 Schematic diagram of integral field spectrometer with deformable fiber bundle

IFS-F的特點(diǎn)是：空間像素?cái)?shù)取決于光纖的根數(shù)，空間像素?cái)?shù)越多，光纖束越多，IFS-F的體積質(zhì)量越大。因此，IFS-F具有較低的空間像素?cái)?shù)，較高的光譜通道數(shù)，主要應(yīng)用于高光譜分辨率、低空間分辨率需求的場合，如實(shí)時(shí)自動(dòng)目標(biāo)識別、IED/爆炸物/化學(xué)探測、生物識別生物醫(yī)學(xué)、導(dǎo)彈攔截、爆炸閃電分析和燃燒診斷等。

2006年，美國Opto-knowledge公司研制了0.4～1.0 μm的4-DIS光譜儀，空間分辨率較低，僅有21×19，光譜通道數(shù)為300，時(shí)間分辨率達(dá)10 000 Hz，如圖10(a)所示。2012年，Opto-knowledge公司采用4個(gè)可變形光纖束提升了空間分辨率，空間分辨率也僅44×40，如圖10(b)所示[29]。

圖10 Opto-knowledge光譜成像儀Fig.10 Opto-knowledge spectral imager

3.3 微透鏡陣列積分視場光譜儀

圖11 基于微透鏡陣列的積分視場光譜儀系統(tǒng)布置Fig.11 Schematic diagram of integral field spectrometer with lenslet arrays

基于微透鏡陣列的積分視場光譜儀(Integral Field Spectrometer with Lenslet Arrays，IFS-L)將微透鏡陣列放置在物鏡的像面上，如圖11所示，微透鏡空間采樣像面上的圖像得到中間像，并在光瞳平面將單個(gè)微透鏡空間范圍內(nèi)的所有光束會聚成一點(diǎn)；光束會聚點(diǎn)之間的空白區(qū)域，為光柵色散留下了足夠的空間；最后，在2D探測器陣列上形成會聚點(diǎn)的色散條帶[30]。

系統(tǒng)的性能主要取決于微透鏡陣列的設(shè)計(jì)。IFS-L的空間像素?cái)?shù)取決于微透鏡陣列的個(gè)數(shù)，光譜通道數(shù)取決于微透鏡之間的間隙。其優(yōu)點(diǎn)在于：微透鏡產(chǎn)生的會聚點(diǎn)光束比入射光束更快(F數(shù)更小)；相對于IFS-F，IFS-L的微透鏡陣列體積要小于可變形光纖束，有效減小了體積質(zhì)量。折衷了空間像素?cái)?shù)和光譜通道數(shù)，數(shù)據(jù)立方體尺寸可達(dá)200×200×27。但是，微透鏡的使用降低了后面光學(xué)系統(tǒng)的F#數(shù)，對后端的光學(xué)器件提出了更高的要求。此外，為了有效利用探測器，會聚點(diǎn)色散留下更大的區(qū)域，需要微透鏡陣列繞光軸傾斜，增加了系統(tǒng)的裝調(diào)難度。

3.4 針孔陣列積分視場光譜儀

圖12 超像素陣列相機(jī)原理Fig.12 Schematic diagram of hyper-pixel array camera

2009年，Bodkin提出了一種基于針孔陣列采樣的積分視場光譜成像儀(Integral Field Spectrometer with Pinholes, IFS-P)，又稱作超像素陣列(Hyper-pixel Array，HPA)相機(jī)。如圖12所示，針孔陣列置于物鏡的像面，替代IFS-L的微透鏡陣列進(jìn)行離散采樣，針孔之間的間隔使每個(gè)采樣點(diǎn)之間有足夠的間隔，保證色散后的光譜不重疊[31]。

2009年，Du等提出了類似于HPA方案的基于棱鏡的多光譜視頻成像系統(tǒng)(Prism-based Multispectral Video Imaging System, PMVIS)[32]。2011年，Cao 等針對PMVIS空間分辨率低的問題，提出了一種基于混合相機(jī)的高光譜采集系統(tǒng)(Hybird Spectral Video Imaging System，HVIS)[33]，其基本原理如圖13所示。一路利用灰度相機(jī)捕獲光譜視頻流，一路利用RGB相機(jī)捕獲高空間分辨率的彩色視頻，兩路視頻同步校準(zhǔn)，通過實(shí)時(shí)光譜傳播算法進(jìn)行融合，最終獲得既具有高空間像素?cái)?shù)又具有高光譜通道數(shù)的光譜數(shù)據(jù)立方體。

圖13 基于混合相機(jī)的高光譜采集系統(tǒng)原理Fig.13 Schematic diagram of hybird spectral video imaging system

IFS-P的優(yōu)點(diǎn)是原理簡單、計(jì)算量極小、實(shí)時(shí)性好，適用于探測大目標(biāo)，常用于農(nóng)業(yè)調(diào)查、環(huán)境監(jiān)測。其缺點(diǎn)是：空間采樣不連續(xù)，難以適用于小目標(biāo)探測識別；針孔采樣的光能利用率較低。

3.5 圖像映射光譜儀

IFS-M技術(shù)適合于低空間像素?cái)?shù)(10數(shù)量級)和高光譜通道數(shù)(100數(shù)量級)的應(yīng)用場合。然而，對于許多應(yīng)用，如顯微鏡，要求較高的空間像素?cái)?shù)(100數(shù)量級)，光譜通道數(shù)只需10～40個(gè)，因此IFS-M難以大量應(yīng)用。雖然微切片器(IFS-μ)能夠?qū)崿F(xiàn)更高的空間采樣，但嚴(yán)重增加了系統(tǒng)設(shè)計(jì)的復(fù)雜性。

圖像映射光譜儀(Image Mapping Spectrometry, IMS)類似于IFS-M，也使用切片鏡，切片鏡置于顯微鏡像面，將顯微圖像切片。IMS與IFS-M的顯著區(qū)別是：IFS-M的切片鏡中，一個(gè)切片元件一個(gè)傾斜角，而IMS的多個(gè)切片元件共享相同的傾斜角度，圖像的多個(gè)切片被映射到同一個(gè)光瞳[11]，提高了空間分辨率。

進(jìn)一步，切片元件長度N決定了x向空間像素?cái)?shù)，切片元件的數(shù)量M決定了y向的空間像素?cái)?shù)，傾斜角度種類L決定了光譜通道數(shù)λ，最終3D(x,y,λ)數(shù)據(jù)立方體的大小為N×M×L。圖14所示的切片元件長度未知，8個(gè)切片元件將圖像切割成8個(gè)條帶，4個(gè)傾斜角使得光譜通道數(shù)為4，故最終的光譜數(shù)據(jù)立方體的大小為x×8×4。

圖14 圖像映射光譜儀原理Fig.14 Schematic diagram of image mapping spectrometry

IMS的優(yōu)點(diǎn)是：每個(gè)光瞳由多個(gè)切片元件共享，允許更高的空間分辨率，通過高分辨率探測器可采集100×100×25或350×350×46的光譜數(shù)據(jù)立方體；微切片鏡將整個(gè)系統(tǒng)的體積縮小到一個(gè)緊湊的尺寸，體積重量較小。其缺點(diǎn)是：切片元件的加工難度大，空間像元數(shù)y越多，切片元件數(shù)越多，切片元件寬度越小(μm級)，加工越難；小切片容易造成衍射，降低成像質(zhì)量；在顯微圖像的衍射點(diǎn)進(jìn)行切片時(shí)，一些反射光線會偏離幾何光路，導(dǎo)致相鄰切片圖像間的串?dāng)_，串?dāng)_成為主要的背景噪聲，嚴(yán)重影響光譜測量精度。

IMS類似于IFS-M，將圖像進(jìn)行空間分割，只需要簡單的圖像重排就可以得到3D數(shù)據(jù)立方體，因此歸類為直接測量的圖像分割技術(shù)。

3.6 像素濾光探測器陣列

圖15 4×4濾鏡單元拼接的馬賽克快照式光譜成像原理Fig.15 Principle for mosaic snapshot spectral imaging composed of 4 × 4 pixel filter array

1976年，Bayer采用相鄰2×2像素級濾鏡單元獲取RGB彩色圖像，其中每個(gè)單元包括1個(gè)紅色、2個(gè)綠色和1個(gè)藍(lán)色的像素級濾鏡。2014年，比利時(shí)微電子研究中心IMEC將此概念推廣，提出了一種基于像素濾光的探測器陣列(Multispectral Sensor with Per-pixel Filters, MS-PPF)，將Fabry-Pérot像素級濾光片陣列整體集成在傳感器表面。將相鄰n×n像素濾鏡陣列視作一個(gè)光譜體素單元，每個(gè)濾鏡的中心波長不同時(shí)，就可同時(shí)獲得n2個(gè)通道的光譜數(shù)據(jù)。圖15是在CMOS傳感器上方制作了4×4個(gè)濾鏡單元，可同時(shí)獲得16個(gè)通道的光譜數(shù)據(jù)[12]。不難知道，這種濾鏡單元拼接的多光譜快照式成像方案通過犧牲空間分辨率來獲取具有一定空間分辨率和光譜分辨率的光譜數(shù)據(jù)立方體。

MS-PPF的優(yōu)點(diǎn)為：將單片濾鏡陣列集成到標(biāo)準(zhǔn)CMOS探測器上，大大減少了系統(tǒng)雜散光、提高了靈敏度和成像速度；采用CMOS工藝技術(shù)，大大降低了成本，具有大規(guī)模制造的潛力；每個(gè)像素濾鏡的晶片級設(shè)計(jì)可實(shí)現(xiàn)緊湊的快照式光譜采集，極大地增加了攝像機(jī)的便攜性、靈活性和易用性，特別適用于尺寸和質(zhì)量受限的應(yīng)用領(lǐng)域，如小型無人機(jī)偽裝目標(biāo)偵察識別應(yīng)用場合。其缺點(diǎn)在于：濾鏡單元內(nèi)不同濾鏡對應(yīng)的位置近似相同，當(dāng)光譜通道數(shù)增大時(shí)，空間誤差增大；光譜通道數(shù)較低，目前的技術(shù)水平是5×5，共25個(gè)光譜通道。

3.7 光場結(jié)構(gòu)多光譜成像技術(shù)

2009年，Levoy在光場相機(jī)的基礎(chǔ)上，提出了基于光場結(jié)構(gòu)的多光譜成像技術(shù)(Spectral Imaging with Light Field Architecture, SI-LF)[13]。其原理如圖16所示，SI-LF包括物鏡、濾光片陣列、針孔陣列、探測器。針孔陣列放置在物鏡的像面上，探測器放置在微透鏡陣列的后方，濾光片陣列布置在物鏡的孔徑處。該結(jié)構(gòu)可分解為兩部分：第一部分是由目標(biāo)、物鏡、針孔陣列和探測器構(gòu)成的光場相機(jī)結(jié)構(gòu)；第二部分是由濾光片陣列、針孔陣列和探測器構(gòu)成的小孔成像結(jié)構(gòu)。

圖16 光場結(jié)構(gòu)的光譜成像原理Fig.16 Spectral imaging principle of light field structure

第一部分中，物鏡將目標(biāo)成像到針孔陣列，因探測器位于針孔陣列后方，會聚到針孔上的光線會繼續(xù)傳播，根據(jù)入射光的角度，重新在探測器表面分配入射光，在探測器上形成一個(gè)體素，即一個(gè)針孔對應(yīng)多個(gè)探測器像元。第二部分中，針孔將濾光片陣列成像到探測器上，使得第一部分形成的體素內(nèi)的小像元調(diào)制有濾光片陣列的光譜信息。將一個(gè)體素視為一個(gè)子單元，提取子單元相同位置的像素構(gòu)成一幅窄帶光譜圖像。

SI-LF設(shè)計(jì)的變體包括用微透鏡陣列代替針孔陣列和用線性可變?yōu)V波器代替濾波器陣列[14]。SI-LF的特點(diǎn)是：針孔數(shù)量等于3D光譜數(shù)據(jù)立方體的空間像素?cái)?shù)；單個(gè)體素包含的像元個(gè)數(shù)等于光譜通道數(shù)；單個(gè)體素內(nèi)各個(gè)像元的光譜配置與濾光片陣列的光譜配置一致；光場相機(jī)的結(jié)構(gòu)實(shí)現(xiàn)了空間像素的復(fù)制功能。

SI-LF的優(yōu)點(diǎn)是濾光片陣列相對于在平鋪濾光的MS-IFT要大得多，制造難度大大降低，更容易實(shí)現(xiàn)。缺點(diǎn)在于：采用針孔，光能量損失較大；體素內(nèi)各個(gè)像素的光來自于目標(biāo)的不同方向，為使光譜調(diào)制效果較好，需滿足朗伯體反射假設(shè)，即各個(gè)方向均勻反射。

Hagen與Gao把SI-LF歸類為基于孔徑分割的光譜成像技術(shù)，認(rèn)為SI-LF類似于Tanida提出的MSI-TOMBO和IMEC提出的MS-ITF，均含有針孔或微透鏡陣列[1-2]。本文則將SI-LF歸類為基于圖像分割的光譜成像技術(shù)，原因是：針孔位于物鏡的像面上，功能上實(shí)現(xiàn)了對圖像的空間采樣，類似于IFS-L，IFS-P，而不是MSI-TOMBO和MS-ITF的場景復(fù)制；SI-LF在探測器上的采集圖像類似于基于像素濾光的光譜成像技術(shù)MS-PPF。

3.8 壓縮編碼孔徑快照式光譜成像儀

2007年，美國杜克大學(xué)Brady研究小組基于壓縮感知理論，提出了雙色散元件的編碼孔徑快照式光譜成像儀(Dual Disperser Coded Aperture Snapshot Spectral Imager，DD-CASSI)[15]。2008年，針對DD-CASSI雙臂4f系統(tǒng)架構(gòu)元件多、長度長、裝調(diào)難等難點(diǎn)，Brady提出了單色散元件的編碼孔徑快照式光譜成像儀(Single Disperser Coded Aperture Snapshot Spectral Imager,SD-CASSI)[34]。

2009年，Brady針對SD-CASSI光路折轉(zhuǎn)不易裝調(diào)的問題，提出了直視型編碼孔徑快照式光譜成像儀(Direct View Coded Aperture Snapshot Spectral Imager,DV-CASSI)[35]。DV-CASSI原理如圖17所示，包括物鏡、編碼孔徑、帶通濾光片、F/8中繼透鏡、Amici棱鏡和CCD探測器。CASSI基于壓縮成像理論，使用編碼孔徑在物鏡的像面上對所成圖像進(jìn)行隨機(jī)二值化編碼，然后通過Amic棱鏡對編碼圖像進(jìn)行色散，由CCD探測器陣列采集空間與光譜混疊的復(fù)合圖像，最后通過圖像重構(gòu)算法將采集的復(fù)合圖像重構(gòu)成3D數(shù)據(jù)立方體。圖像重構(gòu)是圖像形成過程反問題的求解，常采用稀疏重建的梯度投影[36]或兩步迭代收縮/閾值算法[37]。

圖17 雙色散元件的編碼孔徑快照式光譜成像儀原理Fig.17 Schematic diagram of dual disperser coded aperture snapshot spectral imager

CASSI不需要掃描，就可由一幀2D圖像反演出場景的3D光譜數(shù)據(jù)立方體，成像過程具有強(qiáng)實(shí)時(shí)的顯著優(yōu)點(diǎn)；并且基于壓縮感知理論，不犧牲空間分辨率換取光譜分辨率，因此空間分辨率和光譜分辨率均較高。然而，CASSI建立在壓縮感知假設(shè)的基礎(chǔ)上，它要求輸入場景在梯度域中是稀疏的；而且，基于壓縮感知理論的重構(gòu)算法為迭代逼近算法，算法復(fù)雜度高、收斂慢，難以重構(gòu)實(shí)時(shí)；此外，壓縮感知概念采用信號的空間分辨率與圖像質(zhì)量的權(quán)衡，替代空間分辨率與光譜分辨率的折衷，使得結(jié)果數(shù)據(jù)質(zhì)量不可預(yù)測，重構(gòu)常會引入空間和光譜偽影，降低了空間分辨率較高的優(yōu)勢。

近年來，研究人員使用多次編碼快照[38]、高階圖像重建模型[39]、優(yōu)化編碼孔徑[40]和使用兩個(gè)攝像機(jī)的混合設(shè)計(jì)[41]等方式提高CASSI的重建質(zhì)量。

Hagen將CASSI歸類為迭代重構(gòu)的光路分割[2]，本文則將CASSI歸類為迭代重構(gòu)的圖像分割技術(shù)，理由是：(1)編碼孔徑位于物鏡的像面上，在像面對圖像進(jìn)行編碼空間采樣，類似于IFS-L、IFS-P；(2)如果場景完全稀疏，此時(shí)編碼透光孔之間就有了足夠的色散空間，類似于IFS-P。

3.9 小 結(jié)

圖像分割方法目前已形成了3個(gè)技術(shù)方向：(1)基于視場分割器件的IFS；(2)基于像素濾光的MS-PPF；(3)基于壓縮感知理論的CASSI。

IFS通過分割器件，如切片鏡、可變形光纖束、微透鏡陣列、針孔為分割器件等將視場分割。IFS技術(shù)成熟度較高，商業(yè)化程度高，在天文、工業(yè)、醫(yī)療、農(nóng)業(yè)、安防等領(lǐng)域都得到了廣泛應(yīng)用。視場分割器件是關(guān)鍵核心部件，IFS的研究熱點(diǎn)當(dāng)前仍然聚焦于此。

MS-PPF依賴于先進(jìn)的微電子工藝技術(shù)，它集成了薄膜技術(shù)與光刻技術(shù)，研究門檻高，主要受限于研究機(jī)構(gòu)的基礎(chǔ)設(shè)備、基礎(chǔ)工藝和設(shè)計(jì)理念[42]，目前主要集中在具有高精尖微電子設(shè)計(jì)制作能力的研究機(jī)構(gòu)，如比利時(shí)IMEC、芬蘭SPECIM、法國SILIOS、美國PixelIteq公司等。

CASSI建立在壓縮感知理論上，采用簡單的實(shí)驗(yàn)裝置便可實(shí)現(xiàn)快照式成像，受到許多研究人員的偏愛。改進(jìn)編碼方式、重建成像光路以及優(yōu)化重構(gòu)算法均可以提升CASSI性能，是近年來的研究熱點(diǎn)。

4 孔徑分割

4.1 多孔徑分塊濾光

2000年，日本學(xué)者Tanida提出了一種基于復(fù)眼成像原理的光電觀測模塊(Thin Observation Module by Bound Optics，TOMBO)架構(gòu)[43]。2003年，Tanida將TOMBO架構(gòu)應(yīng)用到彩色成像，提出color-TOMBO架構(gòu)，該方法是在微透鏡上方附著R，G，B寬波段濾光片陣列，實(shí)現(xiàn)顏色濾波[44]。2004年，Tanida又將color-TOMBO中的濾光分割進(jìn)一步推廣到多光譜成像，提出MSI-TOMBO，如圖18所示，即在TOMBO系統(tǒng)上配置空間分割的濾光片陣列[16]。

MSI-TOMBO的微透鏡陣列實(shí)現(xiàn)了場景圖像的光學(xué)復(fù)制功能，故歸類為直接測量的孔徑分割技術(shù)。

圖18 基于復(fù)眼成像原理的光電觀測模塊架構(gòu)原理Fig.18 Schematic diagram of thin observation module by bound optics

4.2 單片集成分塊濾光

2013年，比利時(shí)IMEC發(fā)布了復(fù)制濾光的多光譜探測器(Multispectral Sensor with Integrated Tilted Filters, MS-ITF)[17]。MS-ITF是單片集成式器件，包括光學(xué)濾波器與CMOS探測器。如圖19所示，光學(xué)濾波器含有4行8列，共32個(gè)濾鏡單元，并直接制作在CMOS表面。32個(gè)的濾鏡單元窄帶寬度小于10 nm，覆蓋了CMOS探測器敏感的600～1 000 nm波段。當(dāng)采用2 048×1 088像素?cái)?shù)的CMOS探測器時(shí)，4×8個(gè)濾鏡單元實(shí)現(xiàn)了快照式獲取256×256像素、32個(gè)通道的光譜數(shù)據(jù)立方體。

IMEC同時(shí)還提出了基于MS-ITF的快照式多光譜攝像機(jī)，如圖19所示，它由一個(gè)物鏡、一個(gè)光闌、一個(gè)微透鏡陣列和一個(gè)多光譜探測器組成。物鏡將遠(yuǎn)處的場景成像在光闌處，在物鏡的光瞳處放置4行8列的微透鏡陣列，微透鏡陣列將光闌處的成像復(fù)制并中繼到多光譜探測器，多光譜探測器上的4×8個(gè)濾鏡單元分別濾光，并由CMOS傳感器實(shí)時(shí)采集。原理樣機(jī)在日光條件下的采集幀頻大約為30 frame/s，在更高照度水平的條件下，如機(jī)器視覺應(yīng)用中，可達(dá)到340 frame/s。

圖19 基于復(fù)制濾光多光譜探測器的多光譜成像原理Fig.19 Schematic diagram of multispectral imager using multispectral sensor with integrated tilted filters

單片式MS-ITF探測器集成度高，使得光譜攝像機(jī)體積小、質(zhì)量輕，適合大規(guī)模制造，大大降低了成本；此外，立方體重建簡單，允許以視頻速率實(shí)時(shí)、低延遲操作。但是采用的單片集成式光學(xué)濾波器，加工制作難度大。

值得注意的是：IFS-L的微透鏡陣列實(shí)現(xiàn)的是光學(xué)采樣(圖像分割)功能，而MS-ITF多光譜像機(jī)的微透鏡陣列實(shí)現(xiàn)場景圖像的光學(xué)復(fù)制功能。因此，歸類為直接測量的孔徑分割技術(shù)。

4.3 多孔徑漸變?yōu)V光

2018年，Hubold提出了多孔徑漸變?yōu)V光的光譜成像技術(shù)(Spectral Imaging with multi-apertures and Linear Varying Filter, SI-MALVF)，其原理如圖20所示，包括線性漸變?yōu)V光片(Linear Varying Filter, LVF)、微透鏡陣列(Micro-lens Array, MLA)、隔板陣列、玻璃蓋片與FPA成像器件。其中，MAL用于成像場景的多孔徑光學(xué)復(fù)制，LVF實(shí)現(xiàn)光譜濾光，隔板陣列用于孔徑之間的雜散光隔離[18]。

首先，假設(shè)LVF的變譜方向在較短的距離內(nèi)，光譜近似不變。設(shè)計(jì)的關(guān)鍵是線性漸變?yōu)V光片LVF、微透鏡陣列MLA、FPA探測器的位置關(guān)系。在Hubold的設(shè)計(jì)中，將MLA與FPA對正放置，將LVF放置在MLA的上方，如圖20所示。如圖21所示，在LVF與MAL對正放置時(shí)，MLA的高度方向(LVF的恒譜方向)對應(yīng)區(qū)域所成的子圖像是相同的，只能獲得變譜方向微透鏡個(gè)數(shù)的光譜圖像，造成FPA探測器極大浪費(fèi)。Hubold創(chuàng)新地提出將LVF相對MLA傾斜一定角度θ，如圖22所示，此時(shí)MLA的水平x方向和垂直y方向?qū)?yīng)區(qū)域的光譜均不相同，將光譜通道數(shù)提升至微透鏡的個(gè)數(shù)。

圖20 多孔徑漸變?yōu)V光原理Fig.20 Schematic diagram of spectral imaging with multi-apertures and linear varying filter

圖21 漸變?yōu)V光片與微透鏡陣列對正放置Fig.21 Aligned position between LVF, MLA and detector

圖22 漸變?yōu)V光片相對微透鏡陣列旋轉(zhuǎn)Fig.22 LVF rotation relative to MLA

Hubold采用3 243×4 864的CCD探測器，6×11微透鏡陣列，可在450～850 nm波段內(nèi)獲取400×400×66大小的3D光譜數(shù)據(jù)立方體，體積緊湊，達(dá)到60 mm×60 mm×28 mm。

Hubold方法的優(yōu)點(diǎn)是：通過傾斜LVF，大幅提升了光譜通道數(shù)；LVF的制作難度遠(yuǎn)低于分塊濾光片，采用商用LVF將大幅降低光譜儀的制造難度；體積小、質(zhì)量輕。其缺點(diǎn)是：雖然假設(shè)LVF變譜方向的光譜在較短的距離內(nèi)近似不變，但是當(dāng)單個(gè)微透鏡覆蓋的LVF區(qū)域較大時(shí)，光譜擾動(dòng)明顯，覆蓋的區(qū)域越大，光譜擾動(dòng)越大，不利于大孔徑的光譜成像。

2019年，Mu等改變了系統(tǒng)結(jié)構(gòu)，將MLA與CVF的位置顛倒，MLA放置在外側(cè)，LVF放置在FPA探測器前面，不再旋轉(zhuǎn)LVF，而是旋轉(zhuǎn)MLA，大幅抑制了光譜擾動(dòng)。Mu將該結(jié)構(gòu)命名為光學(xué)復(fù)制與重繪成像光譜儀(Optically Replicating and Remapping Imaging Spectrometer, ORRIS)[19]。

MLA旋轉(zhuǎn)一定角度后，每個(gè)微透鏡沿著變譜方向移動(dòng)Δx，如圖23所示。平行于恒譜方向設(shè)計(jì)長條帶，用于從子圖像中取樣數(shù)據(jù)進(jìn)行重建，該長條帶的寬度等于微透鏡變譜方向的移動(dòng)距離Δx，此時(shí)長條帶區(qū)域內(nèi)子圖像將形成無縫平移。當(dāng)長條帶寬度與微透鏡縱向個(gè)數(shù)的乘積等于光譜圖像寬度W時(shí)，拼接后會形成完整的場景圖像。因長條帶的寬度很窄，滿足LVF光譜不變假設(shè)，拼接重建后的圖像光譜帶寬很窄，光譜擾動(dòng)微弱。

圖23 光學(xué)復(fù)制與重繪成像光譜儀的透鏡陣列旋轉(zhuǎn)與圖像重建Fig.23 MLA rotation and image reconstruction of optically replicating and remapping imaging spectrometer

性能上ORRIS原理樣機(jī)能夠在單個(gè)積分時(shí)間內(nèi)捕獲380～850 nm波段內(nèi)80個(gè)連續(xù)光譜通道，空間分辨率為400×400像素，體積大約為50 mm (長)×70 mm (寬) ×70 mm(高)，質(zhì)量為0.5 kg。微透鏡陣列實(shí)現(xiàn)了場景圖像的多孔徑光學(xué)復(fù)制功能，因此將SI_MALVF歸類為直接測量的孔徑分割方法。

4.4 多孔徑傅里葉變換

快照式高光譜成像傅里葉變換(Snapshot Hyperspectral Imaging Fourier Transform, SHIFT)光譜儀的想法是對Hirai于1994年開發(fā)的多圖像傅里葉變換光譜儀(Multiple-image Fourier Transform Spectrometer, MIFTS)的擴(kuò)展[20]。最初的MIFTS是基于邁克爾遜干涉儀和透鏡陣列，而SHIFT光譜儀是基于透鏡陣列后面的一對雙折射Nomarski棱鏡。

如圖24所示，目標(biāo)物經(jīng)過N×M透鏡陣列，復(fù)制為N×M個(gè)子圖像陣列，再經(jīng)過由兩個(gè)Nomarski棱鏡組成的雙折射偏振干涉儀(Birefringent Polarization Interferometer,BPI)，最終探測器上的子圖像產(chǎn)生光程差(Optical Path Difference, OPD)。由于BPI相對于探測器陣列旋轉(zhuǎn)一個(gè)小角度δ，每個(gè)子圖像具有不同的OPD，N×M個(gè)干涉子圖像具有不同的光程差。通過依次提取每個(gè)干涉子圖像來組裝3D干涉圖立方體。沿著干涉圖立方體的OPD軸進(jìn)行傅里葉變換，重建3D光譜數(shù)據(jù)立方體。

圖24 快照式高光譜成像傅里葉變換成像原理Fig.24 Schematic diagram of Snapshot Hyperspectral Imaging Fourier Transform

SHIFT的優(yōu)點(diǎn)是：基于Nomarski棱鏡的設(shè)計(jì)能夠有效減少相機(jī)的體積；雙折射棱鏡便于制造；共光路一體化設(shè)計(jì)，提高了對振動(dòng)的魯棒性。缺點(diǎn)是：采用微透鏡陣列，使得干涉圖子圖像之間存在視差；由于依賴雙折射效應(yīng)，在成像非偏振場景時(shí)，SHIFT的光通量至少降低50%。

因SHIFT的微透鏡陣列的作用是復(fù)制場景，而最終的3D光譜數(shù)據(jù)立方體需要傅里葉變換重構(gòu)，為此將SHIFT歸類為直接重構(gòu)的孔徑分割技術(shù)。

4.5 計(jì)算層析

1997年，Descour和Dereniak提出計(jì)算機(jī)層析成像光譜儀(Computed Tomographic Imaging Spectrometry, CTIS)，它由物鏡、視場光闌、準(zhǔn)直透鏡、2D光柵、二次成像透鏡和CCD等元件構(gòu)成[21]，如圖25所示。全息2D光柵色散元件將準(zhǔn)直后的光束投影色散到3×3方向，通過最大期望迭代反演算法，可由獲取的投影圖像反演到3D光譜數(shù)據(jù)立方體。

圖25 計(jì)算機(jī)層析成像光譜儀原理Fig.25 Schematic diagram of computed tomographic imaging spectrometry

Vandervlugt等人通過設(shè)計(jì)更高階的2D光柵[45]，如5×5和7×7圖案，實(shí)現(xiàn)更多方向的投影，產(chǎn)生更高的光譜分辨率。對于一個(gè)2 048×2 048的FPA，典型的CTIS設(shè)計(jì)將光闌成像到100×100的區(qū)域，沿著626像素長的外對角線給出投影，可生成134通道的光譜圖像。

CTIS實(shí)現(xiàn)了瞬時(shí)快照式的光譜采集，且不需要任何濾光設(shè)備，具有很高的光轉(zhuǎn)換效率。但是，CTIS只能將3D光譜數(shù)據(jù)投影到有限的若干個(gè)方向上，導(dǎo)致采集數(shù)據(jù)對應(yīng)的傅里葉變換域存在一個(gè)圓錐形的采樣丟失，這一圓錐丟失問題嚴(yán)重影響了重構(gòu)的準(zhǔn)確性；并且在實(shí)際應(yīng)用中該系統(tǒng)具有校準(zhǔn)困難，不穩(wěn)定等問題，只適用于精密制造的測量儀器，實(shí)用性不強(qiáng)。

CTIS的2D光柵將場景投影到不同方向，類似于微透鏡陣列的場景復(fù)制，又因需要迭代反演重構(gòu)，因此將CTIS歸類為迭代重構(gòu)的孔徑分割技術(shù)。

4.6 小 結(jié)

孔徑分割方法主要包括場景復(fù)制與光譜分光兩個(gè)過程。場景復(fù)制一般通過微透鏡陣列實(shí)現(xiàn)。按照光譜分光方式可將現(xiàn)有的孔徑分割技術(shù)分為濾光式、干涉式和色散式。其中，干涉式基于傅里葉變換原理，如SHIFT，通過干涉條紋反演數(shù)據(jù)立方體，要求儀器能獲得高質(zhì)量的干涉條紋，裝調(diào)難度大。而濾光式原理簡單，如采用微濾光片陣列的MSI-TOMBO、MS-ITF，采用漸變?yōu)V光片的SI-MALVF。微濾光片陣列方案的設(shè)計(jì)制作難度較高，新型漸變?yōu)V光片方案極大地降低了快照式光譜儀的研制難度，具有很好的研究前景。

綜上不難看出，新型光譜分光的孔徑分割快照式光譜儀是孔徑分割技術(shù)的發(fā)展趨勢。

5 光路分割

5.1 分束器

基于分束器的光譜成像(Spectral Imaging with Beam Splitting, SI-BS)源于通道式光譜成像，其思路是采用分束器將不同中心波長的光束轉(zhuǎn)移到不同的成像通道，由多個(gè)探測器成像獲取[9]。圖26將一組二色濾光片作為分束器，每個(gè)二色濾光片反射預(yù)定義顏色的光，透射其他顏色的光，圖中采用4個(gè)濾光片獲取5通道光譜圖像。這種方式的優(yōu)點(diǎn)是：共光路，不存在視差，可以通過裝調(diào)實(shí)現(xiàn)窄帶光譜圖像的精確配準(zhǔn)，不需要通過后期的圖像處理方式實(shí)現(xiàn)。但是，相機(jī)體積較大，質(zhì)量較重。

圖26 基于多個(gè)光譜分束器的光譜成像原理Fig.26 Spectral imaging principle using multiple spectral-beam splitters

1978年，Bleukens提出了一種用于彩色攝像機(jī)的分色棱鏡組件，如圖27所示，3個(gè)膠結(jié)的分束器立方體將入射光分成3個(gè)色帶，每個(gè)色帶由獨(dú)立的探測器觀察[22]。該組件使彩色攝像機(jī)結(jié)構(gòu)更為緊湊，廣泛用于RGB彩色攝像機(jī)。通過改變分束棱鏡設(shè)計(jì)可以調(diào)整分束后的中心波長，實(shí)現(xiàn)3通道快照式光譜成像。2012年，Murakami 提出了四分色的分色棱鏡，如圖28所示，可實(shí)現(xiàn)4通道的快照式光譜成像，但是通過棱鏡將入射光分成4束以上，在現(xiàn)階段仍沒有更好的辦法實(shí)現(xiàn)[23,46]。SI-BS將光路通過分束器進(jìn)行分割，不需要計(jì)算重構(gòu)，歸類為直接測量的光路分割技術(shù)。

圖27 三通道分色棱鏡組件Fig.27 Dichroic prism assembly

圖28 基于四通道分色棱鏡組件的光譜成像原理Fig.28 Spectral imaging principle using quarter prism assembly

5.2 濾光片堆棧

基于濾光片堆棧的光譜成像(Spectral Imaging with a Filter Stack, SI-FS)如圖29所示。首先用準(zhǔn)直鏡對入射光進(jìn)行準(zhǔn)直，然后濾光片堆棧將不同光譜的光束反射到不同位置，并由聚焦鏡收集，最后在探測器的不同位置形成不同波段的光譜圖像[22]。其中，濾光片堆棧類似于SI-BS的分束器，由一組傾斜的雙色濾光片構(gòu)成，每個(gè)雙色濾光片反射設(shè)計(jì)波段的光束，透射其他波段的光束。

與SI-BS不同的是，SI-FS濾波器堆棧集成在一起，縮短了光路，大幅減小了體積和質(zhì)量；并且由于相鄰之間的角度很小，因此不同波長之間的光路差異可以忽略不計(jì)。然而，濾波器堆棧中可以容納傾斜角度的范圍有限，很難進(jìn)一步增加SI-FS的光譜通道數(shù)，最先進(jìn)的SI-FS可以快照捕獲12個(gè)光譜通道[24]。

由于SI-FS采用濾光片堆棧將入射復(fù)色光反射到探測器的不同位置，因此將該技術(shù)歸類為直接測量的光路分割方法。

圖29 濾光片堆棧原理示意圖Fig.29 Schematic diagram of spectral imaging with a filter stack

5.3 圖像復(fù)分

2003年，美國Heriot-Watt大學(xué)的 Harvey提出了圖像復(fù)分成像光譜儀(Image-replicating Imaging Spectrometer, IRIS)，利用雙折射介質(zhì)的雙光束偏振干涉法實(shí)現(xiàn)光譜濾波，利用Wollaston棱鏡實(shí)現(xiàn)雙光束的分離[25]。

圖30 級聯(lián)雙折射干涉儀的光譜濾波和光線分離原理Fig.30 Spectral filtering and light separation principle of cascaded birefringent interferometer

最簡單的雙光束偏振干涉儀原理如圖30所示，由兩級雙折射干涉儀級聯(lián)而成，每個(gè)雙折射儀由一個(gè)波片延遲器和Wollaston棱鏡組成。輸入光由一個(gè)偏振片線性極化，并傳輸?shù)降谝徊ㄆ舆t器上，波片的快軸與偏振片的光軸呈45°對準(zhǔn)。在尋常光和非常光兩個(gè)偏振分量之間引入光程差b1t1，其中b1是第一波片的折射率，t1是第一波片的厚度。這兩個(gè)偏振分量與波長相關(guān)，可描述為：

T‖=cos2(b1t1k/2),

T⊥=sin2(b1t1k/2),

(1)

其中k=2π/λ是波數(shù)。這兩個(gè)偏振分量被Wollaston棱鏡分離，并指向不同方向。然后第一束偏振光T‖進(jìn)入第二波片延遲器，又產(chǎn)生兩束正交的偏振分量，與原有的偏振分量產(chǎn)生干涉，經(jīng)過第二Wollaston棱鏡后指向不同的方向，描述為：

T‖,‖=cos2(b1t1k/2)cos2(b2t2k/2),

T‖,⊥=cos2(b1t1k/2)sin2(b2t2k/2),

(2)

其中：b2是第二波片的折射率，t2是第二波片的厚度。同樣第二束偏振光T⊥進(jìn)入第二波片延遲器、第二Wollaston棱鏡后，也產(chǎn)生兩束指向不同方向的偏振光：

T⊥,‖=sin2(b1t1k/2)cos2(b2t2k/2),

T⊥,⊥=sin2(b1t1k/2)sin2(b2t2k/2).

(3)

至此，產(chǎn)生了4束指向不同方向，具有不同光譜帶的偏振光，如圖30所示。

光束被分光為2n(n表示W(wǎng)ollaston棱鏡的個(gè)數(shù))束，同時(shí)濾波成2n光譜帶，最后2n束光經(jīng)成像鏡成像在探測器上。

光譜圖像直接記錄，不需要反演計(jì)算，也不存在伴隨反演的噪聲放大。但由于需使用線偏振片，當(dāng)成像非偏振場景時(shí)，會損失一半的光通量；而且需要具有充分折射率的大規(guī)格Wollaston棱鏡，因此難以測量很多的光譜帶，目前最多可測量16個(gè)通道的光譜帶；此外，難以校正與偏振相關(guān)的色差。

IRIS通過級聯(lián)多個(gè)Wollaston棱鏡，將光路拆分為2n個(gè)，可歸類為直接測量的光路分割技術(shù)。

5.4 小 結(jié)

光路分割方法基于通道式光譜成像方法，主要通過分束器或堆棧將光路一分為多。該類技術(shù)光路復(fù)雜、光譜通道數(shù)較少、體積質(zhì)量較大，因此在快照式光譜成像領(lǐng)域的研究不活躍?？梢灶A(yù)測研究新型緊湊、巧妙的分束器，增加光譜通道，減少體積和質(zhì)量仍然是該方法的研究方向。

6 頻域分割

多光譜Sagnac干涉儀(Multispectral Sagnac Interferometer, MSI)是色散補(bǔ)償?shù)腟agnac偏振干涉技術(shù)(Dispersion-Compensated Polarization Sagnac Interferometer, DCPSI)向光譜成像技術(shù)的延伸，其原理如圖31所示[26]。

圖31 多光譜Sagnac干涉儀原理Fig.31 Schematic diagram of multispectral Sagnac interferometer

入射光被線柵分束器WGBS分成兩束。其中，反射光束(偏振方向垂直于紙面)通過閃耀光柵G2以與光軸呈正衍射角衍射，經(jīng)過反射鏡M2和M1的連續(xù)反射，放大橫向剪切量，至閃耀光柵G1，經(jīng)過G1負(fù)衍射角衍射，平行于光軸到達(dá)線性偏振片。透射光束(偏振光譜垂直于光軸)以反向光路最后也以平行于光軸方向到達(dá)線性偏振片。最終兩束光與光軸平行、橫向偏移地離開干涉儀，橫向剪切量為：

(4)

其中：m是衍射級次，d是閃耀光柵周期，a+b+c是G1與G2的距離，λ為波長。再經(jīng)過成像鏡，在探測器上形成與波長無關(guān)的干涉條紋，其空間條紋載波頻率為：

(5)

在DCPSI中通過引入單階閃耀光柵(m=1)，形成了與波長無關(guān)的載波頻率干涉條紋，實(shí)現(xiàn)了色散補(bǔ)償。通過空間濾波和反演計(jì)算，可實(shí)時(shí)獲取3個(gè)Stokes參數(shù)。MSI通過引入多階閃耀光柵(m=2,3,4,5)，形成了與波長無關(guān)，但與衍射級次相關(guān)的干涉條紋；并且通過傅里葉逆變換，可以獲得與衍射級相關(guān)的反演子圖像。因多階閃耀光柵的衍射效率的極大值與衍射級次相關(guān)，衍射級越高，衍射效率的主峰(閃耀波長)越向短波偏移，如圖32所示。最終，每級衍射的反演子圖像與光譜帶唯一對應(yīng)，實(shí)現(xiàn)了快照式四通道光譜成像。

圖32 多階閃耀光柵的衍射效率Fig.32 Diffraction efficiency of multi-order blazed grating

多光譜Sagnac干涉技術(shù)本質(zhì)上是一種在重合的干涉場上同時(shí)測量多個(gè)不同光譜切片的多光譜方法。其優(yōu)點(diǎn)是多個(gè)譜帶之間本質(zhì)是空間配準(zhǔn)的，不需要進(jìn)行配準(zhǔn)處理，降低了后續(xù)處理的難度。其缺點(diǎn)是：光譜通道數(shù)較少，光譜帶必須對應(yīng)于光柵衍射級的閃耀波長；由于線性偏振輸入，光吞吐率減半。因此，該技術(shù)目前的實(shí)用化程度不高。MSI在頻域上進(jìn)行分光，需要通過傅里葉變換重構(gòu)，故將該技術(shù)歸類為直接重構(gòu)的頻域分割技術(shù)。

7 對比分析

表1對比了17種現(xiàn)有光譜成像儀，列出了不同光譜成像儀的分割方法、重構(gòu)方法、空間分辨率、光譜分辨率的受限因素、光通量利用率以及目前各種方案能達(dá)到的數(shù)據(jù)立方體規(guī)格。

表1 快照式光譜成像技術(shù)對比

7.1 空間像素?cái)?shù)與光譜通道數(shù)的權(quán)衡

大多數(shù)快照式光譜成像技術(shù)采用空間分辨率換取光譜分辨率的方法，由2D空間圖像重排或重構(gòu)出3D立方體，為此快照式光譜成像難以同時(shí)獲得較多的空間像素?cái)?shù)和光譜通道數(shù)。

在實(shí)際應(yīng)用中，應(yīng)根據(jù)實(shí)際需求，對空間像素?cái)?shù)與光譜通道數(shù)進(jìn)行權(quán)衡?？勺冃喂饫wIFS-F由于低空間像素?cái)?shù)、高光譜通道數(shù)，適用于動(dòng)態(tài)目標(biāo)跟蹤識別需求。IFS-L，IFS-P，IMS，MS-PPF，SI-LF和MS-ITF由于像素?cái)?shù)和光譜通道數(shù)適中，適用于光譜通道數(shù)與空間像素?cái)?shù)需求均不高的探測識別場合。SI-BS，SI-FS及MSI一般應(yīng)用于光譜通道數(shù)需求較少的場合。

7.2 圖譜匹配

在實(shí)際應(yīng)用中，希望獲得的每個(gè)像素與光譜匹配，如果圖譜之間失配，將為實(shí)際應(yīng)用帶來較大困難。如MS-PPF技術(shù)每個(gè)體素中的光譜實(shí)際上不能正確匹配到同一空間位置，光譜通道數(shù)越多，圖譜失配越嚴(yán)重?；诳讖綇?fù)制的MSI-TMOBO，MS-ITF，SI-MALVF和SHIFT，由于采用了微透鏡陣列復(fù)制場景，場景之間存在視角場景，極易造成圖譜失配，為此要進(jìn)一步研究光譜通道圖像的配準(zhǔn)方法，提高圖譜匹配性能。

7.3 空間采樣的連續(xù)性

IFS-P和SI-LF采用針孔進(jìn)行空間離散采樣，難以用于小目標(biāo)的偵察。雖然IFS-L采用微透鏡會聚采樣空間，但易造成空間光譜的混疊，為后續(xù)應(yīng)用帶來不便。因此，這些方案不適用于小目標(biāo)偵察和光譜精度要求較高的場合。

7.4 光能利用率

光通量一直是制約光譜成像應(yīng)用到遠(yuǎn)距離探測識別的關(guān)鍵因素。SI-LF通過針孔成像，將一個(gè)采樣空間復(fù)制到多個(gè)像素通道。SHIFT，IRIS和MSI由于采用了偏振片，將光通量降至50%以下，大幅降低了光譜成像儀的探測識別能力。在一些應(yīng)用場合，可通過高靈敏度探測器提高光譜圖像的成像質(zhì)量。

7.5 動(dòng)態(tài)范圍

在光譜成像中，成像系統(tǒng)、色散器和探測器的光譜響應(yīng)率不同，導(dǎo)致不同光譜通道的圖像亮度和對比度有很大差異。同一個(gè)成像探測器獲取的不同光譜通道的圖像，會出現(xiàn)有的像素過飽和，有的像素光照不足，難以兼顧不同光譜通道的成像質(zhì)量，這給探測器的增益調(diào)整帶來了難度。因此，快照式光譜成像需要采用具有大動(dòng)態(tài)范圍的探測器。

8 總 結(jié)

快照式光譜成像能夠在單個(gè)積分時(shí)間內(nèi)獲取完整的3D光譜數(shù)據(jù)立方體，不需要通過時(shí)間序列點(diǎn)掃描、線掃描或譜掃描的方式獲得立方體，滿足動(dòng)態(tài)場景的探測識別需求。根據(jù)分割策略，將快照式光譜成像分為圖像分割、孔徑分割、光路分割及頻域分割4類。圖像分割的特征是將分割元件置于成像系統(tǒng)像面；孔徑分割的特征是在成像系統(tǒng)的出瞳處復(fù)制場景；光路分割的特征是在通過分光元件將入射光分成不同的光路?；诜指畈呗缘姆诸惙椒芨采w現(xiàn)有的各種快照式光譜成像技術(shù)方案。

本文綜述了17種不同的快照式光譜成像技術(shù)方案的原理、優(yōu)點(diǎn)及缺點(diǎn)，最后對比了它們的性能。在應(yīng)用快照式光譜成像儀時(shí)，應(yīng)根據(jù)實(shí)際需要權(quán)衡空間像素?cái)?shù)與光譜通道數(shù)。對于像素與光譜匹配精度較高的場合，應(yīng)避免使用像素濾光的MS-PPF，如果采用孔徑分割的MSI-TMOBO，MS-ITF，SI-MALVF，SHIFT需進(jìn)行高精度的光譜圖像配準(zhǔn)。當(dāng)探測識別小目標(biāo)時(shí)，應(yīng)避免選擇空間采樣不連續(xù)的方案IFS-P，SI-LF。對于遠(yuǎn)距離探測識別，應(yīng)采用高靈敏度探測器。為獲得較高質(zhì)量的光譜圖像，應(yīng)采用具有大動(dòng)態(tài)范圍的探測器。