基于棱鏡掃描法的太陽光譜儀光譜定標

高震宇，方 偉，王玉鵬*，張 浩

1. 中國科學院長春光學精密機械與物理研究所, 吉林 長春 130033 2. 中國科學院大學, 北京 100049

基于棱鏡掃描法的太陽光譜儀光譜定標

高震宇1, 2，方 偉1，王玉鵬1*，張 浩1, 2

1. 中國科學院長春光學精密機械與物理研究所, 吉林 長春 130033 2. 中國科學院大學, 北京 100049

為了標定掃描式棱鏡太陽光譜儀的棱鏡不同轉動角度對應的中心波長和光譜帶寬，利用了一種棱鏡掃描方法對太陽光譜儀的光譜響應函數進行測量。該方法使用固定的單色光波長，控制棱鏡轉動實現單色光的像在探測器位置掃描，并通過坐標映射得到響應位置的光譜響應函數。文中根據光譜響應函數的定義，推導出棱鏡掃描法與單色儀波長掃描方法波長定標原理上的等效性。之后分別以532 nm固體激光器和632.8 nm氦氖激光器為光源，使用棱鏡掃描法測量太陽光譜儀對應波長位置的光譜響應函數，并以單色儀波長掃描法實驗作為對比。實驗結果表明，對于掃描式棱鏡太陽光譜儀，棱鏡掃描法測量的中心波長分別為531.86和632.67 nm，其準確度優(yōu)于單色儀波長掃描法測得的531.39和631.97 nm。由于不受單色儀性能的限制，前者測量的光譜帶寬值也優(yōu)于后者。最后以汞燈為光源使用棱鏡掃描法對太陽光譜儀進行了光譜定標實驗，實現了特征光譜定標法結合棱鏡掃描法對中心波長及光譜帶寬的標定。該方法同樣可以應用于掃描式光柵光譜儀以及單色儀的光譜定標。

光譜響應函數； 中心波長； 光譜帶寬； 坐標映射

引 言

作為地球能量的主要來源，太陽輻射對全球能量平衡的至關重要[1]。并且由于大氣、海洋以及不同地貌的陸地對太陽輻射的反射和吸收均具有不同的光譜選擇性[2]，因此精確監(jiān)測太陽光譜輻照度對地球氣候研究和大氣成分分析具有重要意義。太陽光譜儀是對太陽光譜輻照度進行連續(xù)定量化觀測的儀器，準確、高精度的定標是保證其獲得高質量監(jiān)測數據的關鍵。太陽光譜儀的定標主要包括光譜定標和輻射定標，其中光譜定標是輻射定標的前提，且前者的定標精度將對后者的定標精度產生直接影響[3]。

光譜定標可視為求取光譜儀器各光譜通道的光譜響應函數(SRF)的過程，根據光譜響應函數可得到對應的中心波長和光譜帶寬等特性。目前，光譜儀器常用的光譜定標方法主要有特征光譜定標法和單色儀波長掃描法[4]。其中，特征光譜定標法結構簡單操作方便，在HRIS[5]，SCIAMACHY[6]等儀器的在軌定標過程中都有所應用，但該方法一般用于具有線性色散、高光譜分辨率的儀器定標，且只能標定中心波長，不能實現光譜帶寬的標定[7]； 單色儀波長掃描法作為光譜定標常用方法，可同時實現寬光譜范圍的中心波長和光譜分辨率的標定，具有定標精度高，可全光譜定標等優(yōu)點[8]，在AVIRIS[9]，MODIS[10]等光譜成像儀器中得到應用。但由于實驗中光譜定標結果實際是定標光源光譜分布函數與光譜儀器光譜響應函數共同作用的結果[11]，其定標精度受限于所用單色儀的精度和波長準確度。另外，對于色散元件轉動式的光譜儀和單色儀，色散元件的轉動定位精度及定位準確度也會使該類儀器光譜定標產生不確定度及誤差項，這種不確定度與誤差無法通過提高定標用單色儀的精度來消除。

針對上述情況，本文分析并使用了一種棱鏡掃描方法對一架棱鏡色散型太陽光譜儀進行光譜響應函數測量，該方法與單色儀波長掃描法過程不同但原理類似[12]，以固定波長的單色光為光源，利用轉動棱鏡的方法實現單色光的像在探測器位置掃描，并通過坐標映射測得相應位置的光譜響應函數。文中利用一組單色儀波長掃描法測量光譜響應函數實驗作為對比，結果表明，對于文中所述太陽光譜儀，棱鏡掃描法測量的中心波長與光譜帶寬均與單色儀波長掃描法的結果相符且具有更好的準確度與精度。本文最后基于棱鏡掃描法以汞燈為光源對太陽光譜儀進行了光譜定標，驗證了棱鏡掃描法光譜定標可行性。

1 棱鏡掃描法定標原理

圖1為自行研發(fā)并制造的棱鏡型太陽光譜儀的結構原理圖，由入射狹縫、Fery棱鏡、轉動電機、出射狹縫和光電二極管探測器組成。該裝置狹縫和探測器固定，通過電機控制棱鏡轉動獲得工作光譜范圍內的全光譜信息。棱鏡的轉動角度信息通過讀取參考光路中的線陣CCD響應獲得，并將線陣探測器像元編號作為角度代碼。該太陽光譜儀分為三個通道對300～2 500 nm光譜范圍太陽光進行監(jiān)測，光譜分辨率為1～42 nm。本文主要針對300～1 000 nm通道進行分析討論。

圖1 棱鏡太陽光譜儀光路結構

上述系統(tǒng)中，探測器的響應是探測器在對應位置接收到的所有波長能量的疊加[13]，光譜響應函數即為光電二極管探測器對這些不同波長單色光能量響應隨波長分布的函數。針對該系統(tǒng)，單色儀波長掃描法需要固定太陽光譜儀的分光棱鏡在被測波長位置，通過連續(xù)改變單色儀輸出光波長，測量并記錄探測器輸出則獲得該波長位置處的光譜響應函數。由于不同波長對應的狹縫像沿著色散方向分布，此過程即為不同單色波長的像函數對固定的探測器的平移掃描。在所測的中心波長附近一定程度范圍內，像差對單色像能量分布影響有限，可以考慮不同波長單色光的像能量分布相同，僅在色散方向上有位移。則像函數可以表示為

(1)

式中，y為像面坐標，λ為波長，λ0為棱鏡固定位置對應的中心波長，Δλ為波長與中心波長λ0的波長差，Δy為波長為λ的像中心距探測器中心的平移量，Δλ與Δy一一對應。光譜響應函數可表示為

(2)

式中，s(y)為出射狹縫函數。

現在固定輸入的單色光波長，通過電機控制分光棱鏡旋轉，使單色光的狹縫像在探測器附近移動完成掃描。此過程單色像能量分布變化可以忽略，僅在色散方向產生位移，可以等效為固定棱鏡，改變單色光波長的掃描過程。該過程可以看作單色光的狹縫像固定，控制出射狹縫平移完成掃描。此時光譜響應函數可表示為

(3)

由式(2)和式(3)可見，針對本文所述的太陽光譜儀，在假設中心波長附近的小光譜范圍內單色像函數無明顯變化時，固定輸入單色光波長的棱鏡掃描法具有和單色儀波長掃描法相同的推導結果，因此可以考慮上述兩種方法等效。事實上，在中心波長的光譜帶寬的范圍內單色像函數的變化很小，允許上述假設成立。如圖2所示為根據圖1中太陽光譜儀點擴散函數(PSF)模擬[14]的太陽光譜儀632.8 nm的單色光像函數與620 nm以及645 nm單色光像函數的相對偏差，可知最大相對偏差小于1%。

圖2 632.8 nm與附近波長像函數相對偏差

2 棱鏡掃描法測量光譜響應函數

根據棱鏡掃描法的定標原理，針對棱鏡太陽光譜儀，搭建的定標裝置示意圖如圖3所示，由定標光源、石英會聚透鏡、光纖、反射式拋物面準直鏡構成。定標光源根據需要可以是單色儀出射的單色光、譜線燈或者激光，其中激光額外加入一片毛玻璃用來去偏振和勻光，將經過毛玻璃作用后出射的光作為定標光源。石英會聚透鏡和反射式拋物面準直鏡與光纖之間通過SMA905接口連接，并且可以通過微調接口處的距離來調節(jié)出射光的發(fā)散角。

圖3 棱鏡掃描法定標裝置示意圖

使用棱鏡掃描法對太陽光譜儀532和632.8 nm處光譜響應函數進行測量，定標光源分別選擇532 nm的固體激光器和632.8 nm的He-Ne激光器，二者的光譜帶寬均遠小于0.1 nm，與單色儀出射的單色光相比，可以忽略定標光源本身光譜帶寬對光譜響應函數測量的影響，且波長準確度和功率穩(wěn)定性好。分別在532和632.8 nm對應的棱鏡位置處，控制棱鏡精密轉動完成掃描并記錄棱鏡轉角對應的響應值，得到探測器響應隨太陽光譜儀棱鏡轉角變化的曲線，其中棱鏡轉角以一塊線陣CCD像元編號表示。對上述測量結果進行一次坐標映射，將探測器像元坐標映射為波長坐標，則可以獲得儀器的光譜響應函數。此過程可以通過光譜儀光學系統(tǒng)的相關參數模擬獲得，但加工裝調過程中引入的像差無法包含在模擬結果之中。因此本文采用實驗方法，測量多組被測中心波長附近的波長的光譜響應曲線，分別擬合并記錄波長對應的探測器像元，最后通過多項式擬合獲得波長與探測器像元的對應關系。對532和632.8 nm波長附近分別進行上述測量，并對數據進行三次多項式擬合，可得對應532和632.8 nm的擬合結果分別為

y=-28 949.27+20.27x-4.67×10-3x2+3.61×10-7x3

y=-77 827.24+47.54x-9.64×10-3x2+6.55×10-7x3

擬合結果對應的擬合殘差如圖4所示，其中代表橫坐標的棱鏡轉動角和波長的對應關系基于三次多項式擬合，因此并非是線性的。由圖可知，在532和632.8 nm附近位置擬合殘差分別小于0.06和0.09 nm。

圖4 多項式擬合對應的擬合殘差

根據擬合結果對棱鏡掃描法的直接測量的探測器響應隨棱鏡轉角變化曲線進行坐標映射，將角度坐標映射為波長坐標后獲得光譜響應函數曲線，532和632.8 nm波長位置的光譜響應函數如圖5所示，其中波長坐標與探測器像元坐標的對應關系由上文所述的非線性擬合獲得。對兩波長的光譜響應函數分別進行高斯函數擬合，可得實驗測量的中心波長分別為531.86和632.67 nm，對應的函數FWHM值即光譜帶寬分別為9.39和15.14 nm。

圖5 棱鏡掃描法測量的太陽光譜儀光譜響應函數

與棱鏡掃描法作對比，使用單色儀波長掃描法對太陽光譜儀532和632.8 nm波長位置的光譜響應函數進行測量，實驗中使用高功率鹵鎢燈作為光源，實驗所用單色儀最優(yōu)光譜分辨率為0.1 nm，準確度為0.2 nm，重復性為0.1 nm。測量結果如圖6所示，其擬合后的高斯曲線中心波長分別為531.39和631.97 nm，對應的光譜帶寬為9.53和15.36 nm。

對比上述兩種方法的實驗結果，可知對于文中所述的棱鏡轉動式光譜儀，棱鏡掃描法的結果與單色儀波長掃描法的結果相符，且對中心波長的測量更具優(yōu)勢。因為單色儀波長掃描法測量的中心波長會受到擬合過程以及所用單色儀的準確度和重復性影響，除此之外，對于本文中的色散元件掃描式的光譜儀，使用單色儀波長掃描法測量光譜響應函數時需要固定色散元件在定標波長位置處，色散元件實際固定位置與標稱位置之間的誤差也對中心波長的測量產生直接影響。而棱鏡掃描法是使用激光作為光源，光源本身造成的偏移可以忽略，并且該方法不存在色散棱鏡固定位置的誤差，因而避免了相應的對不確定度及準確度的影響，其不確定度主要來自擬合過程和坐標映射過程。在光譜帶寬測量方面，由于激光光源本身帶寬引入的影響幾乎可以忽略，而單色儀出射光的光譜帶寬會對標定的光譜帶寬產生影響[15]，因此單色儀波長掃描法測得光譜帶寬數值大于棱鏡掃描法測的數值。

圖6 單色儀波長掃描法測量太陽光譜儀光譜響應函數

3 棱鏡掃描法的光譜定標

傳統(tǒng)的特征光譜定標法只能實現中心波長的標定，無法標定光譜帶寬。但是應用汞燈等具有特征譜線的光源，結合棱鏡掃描法可以同時實現中心波長和光譜帶寬的標定，并且由于特征譜線具有較高的穩(wěn)定性和幾乎可以忽略的光譜帶寬，與單色儀波長掃描法相比可以減少波長偏移和光源本身光譜帶寬對光譜定標的影響，提高定標的精度與準確度。

圖7 太陽光譜儀對汞燈光源精細掃描響應曲線

太陽光譜儀對汞燈光源光譜進行精細掃描的響應曲線如圖7所示，圖中響應峰值對應波長由太陽光譜儀測量波長與色散元件角度對應關系結合汞燈特征譜線推斷獲得。由于太陽光譜儀光譜分辨率有限，圖中312.57，313.15 nm以及576.96，579.07 nm無法有效分開，因此不作為有效數據。對圖中每一個譜線的響應數據進行高斯函數擬合，可以獲得與該譜線波長對應的色散元件轉動角，將這些定標數據進行多項式擬合，則實現太陽光譜儀在365.02～1 014 nm波段的波長定標。另外對擬合好的高斯函數進行坐標映射，求對應波長位置的光譜響應函數，進而獲得對應波長位置的光譜分辨率。以汞燈為光源的光譜定標數據如表1所示。對表中光譜分辨率數據進行多項式擬合，可得太陽光譜儀365.02～1 014 nm波段范圍內光譜分辨率變化曲線如圖8所示，圖中數據點為實驗測試獲得，曲線是各數據點經三次多項式擬合而成。

表1 汞燈光源光譜定標有效數據

圖8 太陽光譜儀365.02～1 014 nm光譜分辨率曲線

由上述過程可見棱鏡掃描法可以實現文中太陽光譜儀的光譜定標，但由于汞燈在波長范圍內特征譜線數量較少，因此上述光譜定標過程數據量較小，對此可以將氖、氬、氙等光譜校準燈組合使用，則特征譜線從紫外到近紅外均有分布，結合棱鏡掃描法可應用于寬光譜范圍的色散元件掃描式光譜儀器的光譜定標。

4 結 論

針對棱鏡掃描式太陽光譜儀，分析并使用了棱鏡掃描法測量其光譜響應函數，并設置了單色儀波長掃描法測量光譜響應函數作為對比實驗。對比結果表明，對于棱鏡掃描式光譜儀，棱鏡掃描法的測量結果與單色儀波長掃描法測量結果相符，并且由于不受單色儀性能限制，在測定中心波長和光譜帶寬方面相比單色儀波長掃描法具有更好的準確度。最后，本文以汞燈為光源利用棱鏡掃描法對太陽光譜儀進行光譜定標，驗證了棱鏡掃描法光譜定標的可行性。對于色散元件轉動式的光柵光譜儀和單色儀，將本文中的棱鏡掃描改為色散元件掃描，上述方法同樣有效。

[1] LIU En-chao, ZHENG Xiao-bing, LI Xin, et al(劉恩超，鄭小兵，李 新，等). Opt. Precision Eng.(光學精密工程), 2013, 21(3): 608.

[2] ZHANG Hao, FANG Wei, YE Xin, et al(張 浩，方 偉，葉 新，等). Acta Optica Sinica(光學學報), 2013, 33(2): 178.

[3] ZHENG Yu-quan(鄭玉權). Opt. Precision Eng. (光學精密工程), 2010, 18(11): 2347.

[4] ZHANG Chun-lei, XIANG Yang(張春雷，向 陽). Spectroscopy and Spectral Analysis(光譜學與光譜分析), 2012, 32(1): 268.

[5] Barry P S, Shepanski J, Segal C. SPIE, 2002, 4480: 231.

[6] Mager R, Fricke W, Burrows J P. SPIE, 1997, 3106: 84.

[7] QI Xiang-dong, HAN Peng-peng, PAN Ming-zhong, et al(齊向東，撖芃芃，潘明忠，等). Opt. Precision Eng. (光學精密工程), 2011, 19(12): 2870.

[8] JIN Hui, JIANG Hui-lin, ZHENG Yu-quan, et al(金 輝，姜會林，鄭玉權，等). Chinese Journal of Luminescence(發(fā)光學報), 2013, 34(2): 235.

[9] Vane G, Chrien T G, Miller E A, et al. SPIE, 1987, 834: 91.

[10] Amold G T, Fitzgerald M F, Grant P S, et al. SPIE, 1996, 2820: 56.

[11] WANG Ming-zhi, YAN Lei, YANG Bin, et al(王明志，晏 磊，楊 彬，等). Spectroscopy and Spectral Analysis(光譜學與光譜分析), 2013, 33(8): 2280.

[12] Harder J W, Thuillier G, Richard E C, et al. Solar Physics, 2010, 263(1): 3.

[13] Mouroulis P, Sellar R G, Wilson D W, et al. Optical Engineering, 2007, 46(6): 063001.

[14] Mouroulis P, Green R O. SPIE, 2007, 6667: 66670.

[15] SUN Ci, Bayanheshig, CUI Ji-cheng, et al(孫 慈，巴音賀希格，崔繼承，等). Spectroscopy and Spectral Analysis(光譜學與光譜分析), 2014, 34(8): 2289.

*Corresponding author

(Received May 4, 2015; accepted Sep. 20, 2015)

Spectral Calibration of a Solar Spectrometer Based on the Prism-Rotating Method

GAO Zhen-yu1, 2, FANG Wei1, WANG Yu-peng1*, ZHANG Hao1, 2

1. Changchun Institute of Optics, Fine Mechanics and Physics, Chinese Academy of Sciences, Changchun 130033, China 2. University of Chinese Academy of Sciences, Beijing 100049, China

In order to calibrate the corresponding center wavelength and bandwidth to the prism’s rotating angle of a scanning prism-dispersive solar spectrometer, a prism-rotating method to measure the spectral response function (SRF) of the solar spectrometer is suggested. The measuring process is as follows. With the wavelength of monochromatic light invariant, the prism is rotated for scanning the monochromatic image at the location of detector. Then the spectral response function is obtained by mapping the coordinate in position to the coordinate in wavelength. At first in this paper, by analyzing the definition of SRF, the conclusion is deduced that the prism-rotating method is equivalent to the wavelength-scanning method in practice. Then the 532 nm solid-state laser and 632.8 nm He-Ne laser are used as light source. The measurement of SRF of the solar spectrometer based on prism-rotating method is performed. A measurement of SRF with wavelength-scanning method is also performed to be used as a comparison. Experimental results indicate that the center wavelength of the scanning prism-dispersive solar spectrometer measured with prism-rotating method is 531.86 and 632.67 nm respectively. On the other hand, the result is 53139 and 631.97 nm with wavelength-scanning method, which is less precise than the result of prism-rotating method. The values of bandwidth measured with prism-rotating method are also more precise than the latter owing to avoiding the performance deficiency of monochromator. At last, using mercury vapor lamp as light source, an experiment for spectral calibration of the solar spectrometer with prism-rotating method is performed. The values of center wavelength and bandwidth are both achieved with the method combining the prism-rotating method and chracteristic spectrum. This method can also be used to calibrate the monochromator and the grating-dispersive spectrometer whose grating is rotatable.

Spectral response function (SRF); Center wavelength; Spectral bandwidth; Coordinate mapping

2015-05-04，

2015-09-20

國家自然科學基金項目(41474161)資助

高震宇，1988年生，中國科學院長春光學精密機械與物理研究所博士研究生 e-mail: lnsygzy@126.com *通訊聯(lián)系人 e-mail: wangyp@ciomp.ac.cn

TH774

A

10.3964/j.issn.1000-0593(2016)06-1930-06