基于已知球面波前的Hartmann-Shack傳感器結構參量標定

周曉斌，欒亞東，史雷蕾，雷增強

(西安應用光學研究所，陜西 西安 710065)

?

基于已知球面波前的Hartmann-Shack傳感器結構參量標定

周曉斌，欒亞東，史雷蕾，雷增強

(西安應用光學研究所，陜西 西安 710065)

微透鏡陣列到CCD的距離是影響Hartmann-Shack波前探測器精度的主要裝配誤差之一。對該平移裝配誤差的修正，能夠有效減小波前探測誤差。理論求解了球面波前通過微透鏡引起的子孔徑光斑質(zhì)心移動量與波前探測器結構參數(shù)之間的關系，借助該關系能夠求出微透鏡到CCD之間的實際距離，以其改進波前斜率的計算。實驗驗證了理論推導的合理性，并對實際裝配參數(shù)進行標定，得實際距離為24.2mm。利用標定后的參數(shù)重建波前，其相對誤差減小20.4%。實驗表明該標定方法能夠有效提高波前傳感器測量準確性。

自適應光學；Hartmann-Shack傳感器；波前重構；平移誤差

引言

Hartmann-Shack波前探測器(H-S)普遍應用于自適應光學系統(tǒng)、激光光束質(zhì)量判斷、光學元件及光學系統(tǒng)檢測、大氣擾動探測和人眼像差的測量[1-2]。影響H-S探測精度的因素除光斑質(zhì)心位置和波前重構算法外，還有其自身裝配誤差[3]。微透鏡陣列到CCD的距離l是限制H-S測量準確性裝配誤差之一。一般波前斜率計算都是利用微透鏡焦距fn代替該參數(shù)，但是并不能保證CCD靶面與微透鏡焦面完全重合，如果簡單認為l=fn會引入波前斜率計算誤差，文獻[4-5]對此有詳細的分析。因此，對于裝調(diào)完成的H-S，需要進行標定以消除系統(tǒng)誤差。

傳統(tǒng)標定方法是利用修正系數(shù)對波前測量值進行校正。該方法使用fn計算斜率，用該斜率對標準像差板進行波前重建，對比測量結果與實際值差異，計算通用修正系數(shù)。此方法對實驗要求較高，需要多次測量統(tǒng)計后才能計算出系數(shù)。梁春等人提出利用點光源提供離焦波前[6]，再根據(jù)Zernike模式法重建波前中離焦項系數(shù)推算出距離l。該方法依賴于波前重建算法精度。文獻[7]提出利用已知斜率的平面波校準l的方法，其依靠點光源的位移量保證校準精度。

本文結合球面波前特點，理論推導了已知半徑球面波前通過微透鏡陣列的光斑位移量，從而反推出H-S裝配結構參量l，并以實際距離代替微透鏡焦距計算波前子孔徑斜率，重構波前。

1 Hartmann-Shack工作原理

平面波垂直入射到微透鏡陣列時，光斑分別位于子孔徑中心。當畸變波前入射后，子孔徑光斑位置會發(fā)生變化(Δx,Δy)。根據(jù)光斑位移和微透鏡陣列到CCD距離l，可以計算波面子區(qū)域內(nèi)的平均波前斜率g：

(1)

再由波前復原算法即可重構畸變波前。為了匹配微透鏡陣列尺寸，通常在系統(tǒng)前端設置縮束系統(tǒng)，并與微透鏡陣列和CCD探測器構成典型的Hartmann-Shack波前探測器，其工作原理如圖1所示。

圖1 Hartmann-Shack工作原理Fig.1 Principle diagram of Hartmann-Shack

利用模式法重建波前[8-9]，一個完整的波前φ(x,y)可以用Zernike多項式展開為

(2)

式中：Zk(x,y)為各項Zernike多項式；ak為各項系數(shù)。要重構波前只需求出系數(shù)ak。對上式微分有：

(3)

可以得到矩陣式：

G=FA

(4)

式中：G為波前斜率向量；A為Zernike多項式的系數(shù)向量；F稱為波前重構矩陣，表示Zernike多項式在x,y方向的平均斜率。由于F通常非方陣，所以需求出它的廣義逆F+，從而計算出系數(shù)矩陣A。

A=F+G

(5)

2 基于球面波前的Hartmann-Shack結構參量標定

文獻[10]推導出半徑為R的球面波前φ(x,y)=φo-k[(x2+y2)/2R]，其中波矢k=2π/λ，通過焦距fn和子孔徑間隔d的微透鏡陣列后，光斑質(zhì)心相對于平面波前入射時的位移量。微透鏡子孔徑(i，j)內(nèi)X方向相位差為

(6)

則由相位差引起光斑位置，在子孔徑(i,j)焦面上X方向光斑中心偏移量為

(7)

所以，當匯聚或發(fā)散波前入射H-S時，光斑間距會相應縮小或增大，并且光斑仍均勻分布。

相鄰兩子孔徑的光斑位移量之差與球面波前半徑之間關系為

(8)

為了符號統(tǒng)一,δx引入負號即匯聚R為正，光斑間隔變小，利用(8)式便可以計算出實際裝調(diào)完成后的距離l。

可由平行光管提供已知半徑的球面波前。當平行光管移焦x，由牛頓公式得焦像距x′=-f2/x。根據(jù)H-S前段望遠(縮束)系統(tǒng)光學結構參數(shù)，物鏡fo、目鏡fe、H-S距平行光管出瞳L等，計算得最終入射到微透鏡的球面波前半徑為

(9)

式中，M為H-S縮束系統(tǒng)倍率。

為了減小系統(tǒng)誤差，采用差值法處理。結合(8)式、(9)式可知，當平行光管移焦量由x1變化為x2時，

(10)

從公式(10)可以看出，l與具體光路參數(shù)無關，只要已知f、d、M，測量平行光管2次移焦后相鄰光斑間距變化量δx1和δx2，就可以計算出l。而l的標定誤差主要由光斑間距變化量所確定。根據(jù)球面波的對稱性可知Y方向關系與X方向相同。

3 實驗驗證

利用以上算法對一個已經(jīng)裝調(diào)完成的H-S進行標定。平行光管焦距f=1 600mm；H-S縮束倍率M=4；微透鏡陣列焦距fn=32mm，d=1mm。光斑質(zhì)心提取軟件是基于CCD控制軟件包自行開發(fā)的，利用Zernike模式法重建波前。

圖2 實驗裝置示意圖Fig.2 Schematic diagram of experiment setup

從圖3和圖4可以看出，光斑移動規(guī)律與理論推導一致，其中圖4中OPD(optional path difference)表示波面光程差,發(fā)現(xiàn)平行光管光軸并未與H-S光軸重合。H-S測量的為球面波非對稱部分，但這并不影響l的計算。因為標準球面波引起光斑移動是均勻的，只要H-S垂直于平行光管光軸即可。而平行光管和H-S視場都很小，只要保證平面波入射子區(qū)域光斑分布均勻就能滿足光路要求。

圖3 平行光管移焦100 mm后光斑變化Fig.3 Spots changed before and after collimator defocus 100 mm

圖4 移焦100 mm后波前復原圖Fig.4 OPD distribution after defocus 100 mm

選取平行光管移焦40 mm、60 mm、80 mm 3組數(shù)據(jù)兩兩差值，由公式(10)計算得平均值l=24.2 mm。為了檢驗標定后的H-S測量準確性，分別用fn=32.0 mm和l=24.2 mm計算斜率，對同一波前進行重建，對比重構波前PV值和理論計算值。球面波矢高為其理論PV值。

圖5 H-S截取球面波矢高計算Fig.5 Rise calculation for spherical wave measured by Hartmann-Shack

表1 重構誤差對比

從上表可以看出，利用l計算斜率重構波前，相比于直接利用微透鏡焦距重構，波前擬合PV值更接近理論值，相對誤差減小20.4%，證明標定后波前測量準確度提高。然而，無論利用fn還是l重構，測量值相比于理論值都存在較大誤差，一方面是由于球面波球心點選擇造成理論計算PV值存在誤差。另一方面，在重構波前系數(shù)中X方向有一定的傾斜量，會影響測量PV值。但是，在相同球心選擇誤差和波前傾斜量下，對l進行標定能夠有效提高H-S測量準確性。

4 總結

本文研究利用已知球面波前對Hartmann-Shack波前探測器平移裝配誤差的標定方法。對此參數(shù)的標定可以有效提高H-S測量的可靠性和準確性。相比于傳統(tǒng)的像差板確定修正系數(shù)，此方法原理簡單，與波前重構算法無關，只需要簡單的實驗儀器即可實現(xiàn)，能夠有效指導實際裝調(diào)過程，降低裝調(diào)精度要求。

[1] Kou Songfeng，Liu Genrong，Niu Dongsheng，et al. Design of SONG Shack-Hartmann optical system[J]. Journal of Applied Optics，2014，35(1)：1-6. 寇松峰，劉根榮，牛冬生，等.“宋”望遠鏡夏克哈特曼光學系統(tǒng)設計[J]. 應用光學，2014，35(1)：1-6.

[2] Guo Guangyan，F(xiàn)an Zhongwei，Yu Jin，et al. Wavefront detection technology based on Shark-Hartmann theory[J]. Journal of Applied Optics，2014，35(5)：823-829. 郭廣妍，樊仲維，余錦，等.基于Shark-Hartmann理論的波前探測技術研究[J]. 應用光學，2014，35(5)：823-829.

[3] Xia Mingliang，Li Chao，Liu Zhaonan，et al. Adaptive threshold selection method for Shack-Hartmann wavefront sensors[J]. Optics and Precision Engineering，2010，18(2)：334-340. 夏明亮，李抄，劉肇南，等.Shack-Hartmann波前傳感器圖像自適應閾值的選取[J]. 光學精密工程，2010，18(2)：334-340.

[4] Ye Hongwei，Xian Hao，Zhang Yudong. Study on translation error of Hartmann-Shack wavefront sensor[J]. Opto-Eletronnic Engineering，2003，30(2)：1-10. 葉紅衛(wèi)，鮮浩，張雨東. 對Hartmann-Shack波前傳感器平移誤差的研究[J]. 光電工程，2003，30(2)：1-10.

[5] Dai Yang，Li Faquan，Cheng Xuewu，et al. Assembly error of Hartmann-Shack sensor[J]. High Power Laser and Particle Beams，2006，18(9)：1469-1474. 戴陽，李發(fā)泉，程學武，等. Hartmann-Shack傳感器組裝誤差分析[J]. 強激光與粒子束，2006，18(9)：1469-1474.

[6] Liang Chun，Sheng Jianxin，Tong Gui，et al. A self-reference method for measuring Hartmann-Shack warefront sensor parameter[J]. Acta photonica sinica，2009，38(4)：780-784. 梁春，沈建新，童桂，等. Hartmann-Shack傳感器結構參量的自基準標定[J]. 光子學報，2009，38(4)：780-784.

[7] Wang Yanping，Wang Qianqian，Ma Chong. Study on hartmann wavefront analyzer calibration method[J]. Chinese Journal of Lasers，2015，42(1)：0108003-1-0108003-7. 王艷萍，王茜蒨，馬沖. 哈特曼波前分析儀校準方法研究[J]. 中國激光，2015，42(1)：0108003-1-0108003-7.

[8] Zhou Renzhong. Adaptive optics[M]. Beijing: National Defence Industry Press，1996. 周仁忠. 自適應光學[M]. 北京：國防工業(yè)出版社，1996.

[9] Xia Aili. The study of adaptive optical system based on Shack-Hartmann wavefront sensor[D]. Xi’an:Xi’an Institute of Optics and Precision Mechanics of CAS，2011. 夏愛麗. 基于Shack-Hartmann波前傳感器的自適應光學系統(tǒng)研究[D]. 西安：中國科學院西安光學精密機械研究所，2011.

[10]Wan Min，Zhang Jiaru，Zhou Ke’en，et al. Calibrate plane beam using Hartmann-Shack method[J]. High Power Laser and Particle Beams，1997，9(3)：443-446. 萬敏，張家如，周克恩，等. 用哈特曼-夏克法標定平面光束[J]. 強激光與粒子束，1997，9(8)：443-446.

Structural parameters calibration of Hartmann-Shack sensor based on known spherical wavefront

Zhou Xiaobin，Luan Yadong，Shi Leilei，Lei Zengqiang

(Xi’an Institute of Applied Optics，Xi’an 710065，China)

The distance between the charge couple-device (CCD) and the micro lens array is one of the assembling errors which affects the accuracy of Hartmann-Shack wavefront detector. The correction of this motion assembly parameter can effectively reduce the wavefront detector error. The relationship between the subaperture spot centroid shift amount,caused by the spherical wavefront passing through the micro lens, and the structural parameters of wavefront sensor was derived theoretically. By utilizing the relationship, it is able to work out the actual distance from the CCD to the micro lens array and improve the calculation of the slope of wavefront. Experiments verified the rationality of theoretical derivation, and the actual assembly parameter calibrated is 24.2 mm. The wavefront was reconstructed by the calibrated parameter, the relative error can reduce by 20.4%. Experiment results show that the calibration method can effectively improve the measurement accuracy of wavefront sensor.

adaptive optics; Hartmann-Shack sensor; wavefront reconstruction; motion error

1002-2082(2015)06-0909-04

2015-06-11；

2015-08-30

周曉斌(1989-)，男，陜西岐山人，碩士研究生，主要從事光學系統(tǒng)設計工作。E-mail:xbin205@163.com

TP

A

10.5768/JAO201536.0603002