采樣點(diǎn)分布對(duì)基于面形斜率徑向基模型的自由曲面擬合精度的影響

武鵬飛，張　贊，鄭　義，張　娟，王堯堯，趙　星*，王靈杰，伍雁雄

(1.南開(kāi)大學(xué) 現(xiàn)代光學(xué)研究所 光學(xué)信息技術(shù)科學(xué)教育部重點(diǎn)實(shí)驗(yàn)室，天津 300071；2.中國(guó)科學(xué)院 長(zhǎng)春光學(xué)精密機(jī)械與物理研究所光學(xué)系統(tǒng)先進(jìn)制造技術(shù)中國(guó)科學(xué)院重點(diǎn)實(shí)驗(yàn)室， 吉林 長(zhǎng)春 130033)

?

采樣點(diǎn)分布對(duì)基于面形斜率徑向基模型的自由曲面擬合精度的影響

武鵬飛1，張贊1，鄭義1，張娟1，王堯堯1，趙星1*，王靈杰2，伍雁雄2

(1.南開(kāi)大學(xué) 現(xiàn)代光學(xué)研究所 光學(xué)信息技術(shù)科學(xué)教育部重點(diǎn)實(shí)驗(yàn)室，天津 300071；2.中國(guó)科學(xué)院 長(zhǎng)春光學(xué)精密機(jī)械與物理研究所光學(xué)系統(tǒng)先進(jìn)制造技術(shù)中國(guó)科學(xué)院重點(diǎn)實(shí)驗(yàn)室， 吉林 長(zhǎng)春 130033)

鑒于自由曲面模型的面形擬合精度在自由曲面表征以及面形初始結(jié)構(gòu)選取等研究中的重要性，本文針對(duì)基于面形斜率的高斯徑向基表征模型，研究了不同的采樣點(diǎn)分布類型對(duì)該模型面形擬合精度的影響。采用不同采樣點(diǎn)分布擬合離軸二次曲面和帶凸起的拋物面，結(jié)果表明采用均勻隨機(jī)分布的采樣點(diǎn)有利于實(shí)現(xiàn)高精度的面形擬合，且達(dá)到一定的擬合精度后，采樣點(diǎn)的數(shù)目對(duì)擬合精度的影響有限。以離軸三反系統(tǒng)為設(shè)計(jì)實(shí)例，對(duì)比了由不同采樣方式生成初始面形后系統(tǒng)的像質(zhì)優(yōu)化結(jié)果。結(jié)果顯示，采用均勻隨機(jī)型采樣方式得到的初始面形進(jìn)行系統(tǒng)優(yōu)化，最終全視場(chǎng)平均調(diào)制傳遞函數(shù)(MTF)可以達(dá)到0.72以上,遠(yuǎn)高于由邊緣集中采樣方式生成初始面形后系統(tǒng)像質(zhì)的優(yōu)化結(jié)果，從而印證了理論研究結(jié)果。

光學(xué)設(shè)計(jì)；自由曲面；擬合精度；采樣點(diǎn)分布

1　引　言

隨著以超精密金剛石切削技術(shù)為代表的現(xiàn)代加工與檢測(cè)技術(shù)的發(fā)展，光學(xué)自由曲面[1]的研究逐漸走向成熟。光學(xué)自由曲面可以在一個(gè)表面同時(shí)實(shí)現(xiàn)多種功能，在提高光學(xué)性能、簡(jiǎn)化結(jié)構(gòu)方面有著非球面無(wú)法比擬的優(yōu)勢(shì)[2]。然而，自由曲面的面形自由度高，表面矢高變化復(fù)雜，對(duì)面形表征精度的要求極高，所以面形擬合精度成為自由曲面表征模型的重要研究?jī)?nèi)容之一[3-4]。自由曲面的表征模型分為兩大類：一類是以泰勒函數(shù)為代表的有理函數(shù)；另一類是以高斯徑向基為代表的無(wú)理函數(shù)[5]。作為典型的有理函數(shù)自由曲面表征模型，Zernike多項(xiàng)式的發(fā)展較為成熟，但它對(duì)大離軸量非對(duì)稱表面的像差補(bǔ)償能力有限，很難達(dá)到滿意的成像質(zhì)量[6-7]。因此，具有局部面形表征能力的高斯徑向基模型成為近些年的研究熱點(diǎn)[6-9]。

根據(jù)高斯徑向基模型的局域性特點(diǎn)，本課題組提出了基于面形斜率的高斯徑向基模型(Gaussian Radial Basis Function Based on Slope， Gaussian RBF-BS)[10]。不同于傳統(tǒng)高斯徑向基模型中所有基函數(shù)形狀因子數(shù)值相同的特點(diǎn)，Gaussian RBF-BS模型中決定模型局域性的形狀因子與面形斜率相關(guān)聯(lián)，因此該模型對(duì)非對(duì)稱性表面的擬合精度高出傳統(tǒng)高斯徑向基模型1-2個(gè)數(shù)量級(jí)，具有更高的面形表征能力[10]。

在光學(xué)設(shè)計(jì)過(guò)程中，面形初始結(jié)構(gòu)在一定程度上決定了系統(tǒng)的優(yōu)化結(jié)果和成像質(zhì)量[11-13]，一般利用自由曲面數(shù)理模型通過(guò)面形擬合方法獲得光學(xué)自由曲面系統(tǒng)的初始面形。所以，作為光學(xué)系統(tǒng)高效優(yōu)化、快速收斂的一個(gè)重要前提，面形擬合精度不僅反映了自由曲面表征模型的面形表征能力，而且影響了系統(tǒng)優(yōu)化設(shè)計(jì)中自由曲面的面形初始結(jié)構(gòu)，進(jìn)而影響最終的優(yōu)化設(shè)計(jì)結(jié)果。大量研究表明，光學(xué)表面的面形擬合精度受到基函數(shù)數(shù)目和采樣點(diǎn)特性等多種因素的影響[14-18]。而高精度面形擬合中，采樣點(diǎn)特性對(duì)面形擬合精度的影響更加顯著，因此國(guó)內(nèi)許多研究組對(duì)采樣點(diǎn)特性與面形擬合精度的關(guān)系進(jìn)行了研究。馮婕等研究了Zernike多項(xiàng)式項(xiàng)數(shù)及采樣點(diǎn)數(shù)目對(duì)擬合精度的影響[14]，孫學(xué)真研究了抽樣點(diǎn)數(shù)目對(duì)Zernike曲面面形擬合的影響[15]，而謝蘇瀧等則從樣項(xiàng)比方面對(duì)Zernike多項(xiàng)式曲面擬合進(jìn)行了研究[16]。這些研究表明，采樣點(diǎn)數(shù)目是影響擬合精度的重要因素。國(guó)外研究組也研究了采樣點(diǎn)數(shù)目以及分布方式對(duì)擬合精度的影響。G W Forbes研究了基于正交基底的自由曲面擬合方法，并分析了面形擬合誤差的影響因素[17-18]。Rolland則從擬合精度與采樣點(diǎn)分布關(guān)系的角度對(duì)Zernike多項(xiàng)式以及Q多項(xiàng)式的自由曲面擬合進(jìn)行了研究。她通過(guò)分析比較不同采樣點(diǎn)分布類型及采樣點(diǎn)數(shù)目對(duì)Zernike多項(xiàng)式模型面形擬合精度的影響，指出在基函數(shù)數(shù)量足夠的條件下，曲面采樣點(diǎn)的分布和數(shù)量對(duì)擬合過(guò)程的收斂速度、計(jì)算效率以及面形擬合精度有很大的影響[3-4]。

作者曾對(duì)比研究了Gaussian RBF-BS模型和Zernike模型的面形表征能力，并且利用三反系統(tǒng)實(shí)例研究了兩模型在光學(xué)系統(tǒng)設(shè)計(jì)中的像差平衡能力。研究表明，Zernike模型與Gaussian RBF-BS模型不僅數(shù)學(xué)表達(dá)式截然不同，而且兩模型的面形表征能力和像差平衡能力也不同[10]。因此要使具有高自由度的Gaussian RBF-BS模型在光學(xué)系統(tǒng)設(shè)計(jì)中得到更好的應(yīng)用，研究采樣點(diǎn)分布對(duì)該模型面形擬合精度的影響具有重要意義。

本文研究了兩種采樣點(diǎn)分布類型以及采樣點(diǎn)數(shù)目對(duì)Gaussian RBF-BS模型面形擬合精度的影響，通過(guò)對(duì)比和分析擬合結(jié)果確定了適用于Gaussian RBF-BS模型的采樣點(diǎn)分布類型，并以離軸三反系統(tǒng)設(shè)計(jì)為例，分析比較了由不同采樣方式生成初始面形后系統(tǒng)的像質(zhì)優(yōu)化結(jié)果，驗(yàn)證了理論研究結(jié)果的正確性。

2　Gaussian RBF-BS模型

光學(xué)自由曲面的Gaussian RBF-BS模型表達(dá)式為[5]：

(1)

針對(duì)傳統(tǒng)高斯徑向基形狀因子取值相同導(dǎo)致的面形擬合能力受限問(wèn)題[19]，本文根據(jù)基函數(shù)所處位置的面形矢高平均變化率來(lái)調(diào)整基函數(shù)形狀因子εi，使εi正比于自由曲面的局部面形斜率。面形矢高變化較大的區(qū)域，εi值較大，從而基函數(shù)的形狀會(huì)更為細(xì)銳，εi對(duì)基函數(shù)中心及鄰近區(qū)域內(nèi)的矢高擬合貢獻(xiàn)較大，而對(duì)口徑內(nèi)其他區(qū)域矢高擬合的影響較小，擬合行為表現(xiàn)出明顯的局域性，因此容易獲得較高的面形擬合精度。

3　采樣點(diǎn)分布類型及待擬合表面

(a)邊緣集中型采樣點(diǎn)分布

(b)均勻隨機(jī)型采樣點(diǎn)分布

為了研究不同采樣點(diǎn)分布類型和數(shù)量對(duì)Gaussian RBF-BS模型面形擬合精度的影響，選擇兩種面形進(jìn)行擬合實(shí)驗(yàn)研究。

隨著IP技術(shù)的發(fā)展，大容量計(jì)算能力將會(huì)帶來(lái)信號(hào)處理方面革新。傳統(tǒng)的播出系統(tǒng)在發(fā)展到云播出之前，中間的一個(gè)重要步驟是一體化播出服務(wù)器，其核心技術(shù)是將傳統(tǒng)播出系統(tǒng)的播出信號(hào)、信號(hào)切換、鍵控等功能軟件化，運(yùn)行在通用計(jì)算平臺(tái)上，接口雖然是SDI，但其內(nèi)部處理完全I(xiàn)P化，配合云計(jì)算平臺(tái)和云處理軟件完成整個(gè)工作流程。

第一種面形是離軸二次曲面，其面形表達(dá)式如下[5]：

(2)

其中：x方向的頂點(diǎn)曲率和非球面系數(shù)分別為cx=-1/455.4 mm-1，kx=-1.306；y方向的頂點(diǎn)曲率和非球面系數(shù)分別為cy=-1/455.4 mm-1，ky=-1.306；有效口徑為110 mm，離軸量為50 mm。該自由曲面的矢高如圖2所示。

圖2　離軸二次曲面的面形矢高

第二種面形是帶凸起的拋物面。該面形具有一定的普適性[15]，不僅包含二次曲面項(xiàng)，還包含不同函數(shù)中心的高斯項(xiàng)，位置分別位于0.92口徑處、0.36口徑處(凸起程度最大)和0.86口徑處。因此，這個(gè)表面既有位于口徑邊緣的凸起，又有位于口徑內(nèi)部的凸起。其面形表達(dá)式如下[5]：

0.6e-0.49[(x+3)2+(y-2)2]+0.03e-0.81[(x-5)2+(y-7)2].

(3)

單位口徑內(nèi)，該自由曲面的矢高如圖3所示。與離軸二次曲面相比，該表面在口徑內(nèi)部具有較大的矢高變化，即面形斜率變化較大，因此其高精度的面形擬合的難度也更大。

圖3　帶凸起的拋物面面形矢高

4　面形擬合結(jié)果及分析

實(shí)驗(yàn)中，均勻隨機(jī)分布和邊緣集中分布兩種采樣點(diǎn)分布類型生成采樣點(diǎn)坐標(biāo)，然后利用Gaussian RBF-BS模型對(duì)上述兩種面形進(jìn)行面形擬合，最后通過(guò)計(jì)算所有采樣點(diǎn)擬合誤差的均方根來(lái)評(píng)估模型的面形擬合精度。

(a)均勻隨機(jī)分布采樣點(diǎn)擬合后的面形矢高差

(b)邊緣集中分布采樣點(diǎn)擬合后的面形矢高差

(c)不同采樣點(diǎn)分布類型的面形擬合精度隨采樣點(diǎn)數(shù)目的變化

(c)Fitting precision versus number of sampling points for different distribution of sampling points

圖4使用不同采樣點(diǎn)分布的離軸二次曲面面形擬合結(jié)果

Fig.4Fitting results of off-axis quatric surface by using sampling points in different distributions

圖4(c)是不同采樣點(diǎn)分布類型的面形擬合精度隨采樣點(diǎn)數(shù)目的變化曲線。從圖4(c) 可以看出，在Gaussian RBF-BS模型擬合離軸二次曲面的過(guò)程中，上述兩種采樣點(diǎn)分布類型最終均可得到比較滿意的擬合結(jié)果，擬合精度達(dá)到了10-11m 量級(jí)。在采樣點(diǎn)數(shù)目相同的條件下，與均勻隨機(jī)型采樣點(diǎn)分布相比，邊緣集中型采樣點(diǎn)分布面形擬合后的RMS值更高，即采用均勻隨機(jī)采樣點(diǎn)分布進(jìn)行面形擬合能夠獲得更高的面形擬合精度。這一結(jié)果從圖4(a) 和4(b)的矢高差分布情況也可明顯看出。同時(shí)隨著采樣點(diǎn)數(shù)目的增多，使用兩種采樣點(diǎn)分布類型進(jìn)行面形擬合后的RMS數(shù)值變化幅度較小，這說(shuō)明達(dá)到一定的擬合精度后，采樣點(diǎn)數(shù)量對(duì)模型面形擬合精度的影響并不明顯。

圖5是Gaussian RBF-BS模型擬合帶凸起的拋物面的結(jié)果。由于面形矢高變化較大，擬合難度較高，需要使用更多的基函數(shù)才能獲得較高的擬合精度。圖5(a)、5(b)分別為采樣點(diǎn)數(shù)目為3 000 個(gè)，基函數(shù)為729項(xiàng)時(shí)，以均勻隨機(jī)型分布和邊緣集中型分布方式在全口徑內(nèi)擬合得到的帶凸起拋物面的面形矢高差。圖5(c)是不同采樣點(diǎn)分布類型的面形擬合精度隨采樣點(diǎn)數(shù)目的變化曲線。從圖5(c)可以看出，在Gaussian RBF-BS模型擬合帶凸起的拋物面時(shí)，隨著采樣點(diǎn)數(shù)目的增多，擬合精度呈現(xiàn)先提高后不變的趨勢(shì)。采樣點(diǎn)數(shù)目不足會(huì)導(dǎo)致凸起拋物面面形矢高變化較大的凸起處欠采樣，從而影響擬合精度。當(dāng)采樣點(diǎn)數(shù)目超過(guò)2 000個(gè)時(shí)，最終的擬合精度能達(dá)到nm量級(jí)，擬合精度不再受采樣點(diǎn)數(shù)目的影響，采樣點(diǎn)分布成為決定擬合精度的主要因素之一。因此，在此擬合精度下，相同的采樣點(diǎn)數(shù)目，邊緣集中型采樣點(diǎn)分布進(jìn)行面形擬合后的RMS數(shù)值同樣高于均勻隨機(jī)型采樣點(diǎn)分布，即邊緣集中型采樣點(diǎn)分布的擬合精度低于均勻隨機(jī)型采樣點(diǎn)分布的擬合精度。

綜合兩類面形的擬合分析結(jié)果可知，Gaussian RBF-BS模型要達(dá)到一定的面形擬合精度，需要足夠的采樣點(diǎn)數(shù)量，并且擬合過(guò)程中采樣點(diǎn)應(yīng)呈均勻隨機(jī)型分布。在達(dá)到一定的擬合精度后，采樣點(diǎn)數(shù)目對(duì)擬合精度的影響不明顯，而過(guò)多的采樣點(diǎn)數(shù)量會(huì)增加擬合運(yùn)算量，降低擬合效率。已有研究表明，傳統(tǒng)的Zernike多項(xiàng)式模型使用邊緣集中型采樣點(diǎn)分布更有利于提高面形擬合精度[3]。因此，采樣點(diǎn)分布方式對(duì)基于不同自由曲面表征模型的面形擬合精度的影響是不相同的。

(a)均勻隨機(jī)分布采樣點(diǎn)擬合后的面形矢高差

(b)邊緣集中分布采樣點(diǎn)擬合后的面形矢高差

(c)不同采樣點(diǎn)分布類型的面形擬合精度隨采樣點(diǎn)數(shù)目的變化

(c)Fitting precision versus number of sampling points for different distribution of sampling points

圖5使用不同采樣點(diǎn)分布的帶凸起的拋物面面形擬合結(jié)果

Fig.5Fitting results of bumpy paraboloid surface by using sampling points in different distributions

采樣點(diǎn)分布方式對(duì)Gaussian RBF-BS模型和Zernike多項(xiàng)式模型面形擬合精度的影響之所以出現(xiàn)完全不同的結(jié)果，其原因在于Zernike多項(xiàng)式模型是一種全局性表征模型，擬合過(guò)程中矢高差通常在口徑邊緣處顯著增大，因此邊緣集中型采樣點(diǎn)分布將增加口徑邊緣處的采樣點(diǎn)個(gè)數(shù)(如圖1(a)所示)，從而有利于面形擬合精度的提高。而Gaussian RBF-BS模型是一種具有局域性的表征模型，通過(guò)設(shè)置與面形斜率相關(guān)的形狀因子，使得擬合過(guò)程中的矢高差在口徑內(nèi)呈現(xiàn)相對(duì)均勻的分布特性，如圖4(a)、4(b)和圖5(a)、5(b) 所示。如果采用邊緣集中分布的采樣點(diǎn)，會(huì)造成擬合過(guò)程中的口徑內(nèi)欠采樣，反而降低了全口徑的擬合精度，而均勻隨機(jī)分布的采樣點(diǎn)則明顯有利于面形擬合精度的提高。隨著面形擬合精度的提高，擬合結(jié)果會(huì)越接近理想面形。因此，在光學(xué)系統(tǒng)設(shè)計(jì)過(guò)程中，采用均勻隨機(jī)分布方式獲得的初始結(jié)構(gòu)進(jìn)行系統(tǒng)像質(zhì)優(yōu)化可以得到更好的像質(zhì)。

5　設(shè)計(jì)實(shí)例

為了進(jìn)一步說(shuō)明不同采樣點(diǎn)分布方式對(duì)面形表征和系統(tǒng)設(shè)計(jì)中初始結(jié)構(gòu)求解的影響，本文以離軸三反系統(tǒng)為例，對(duì)系統(tǒng)中的三鏡分別采用由不同采樣點(diǎn)分布方式得到的兩種RBF-BS自由曲面進(jìn)行面形表征及系統(tǒng)優(yōu)化，比較系統(tǒng)設(shè)計(jì)結(jié)果的成像質(zhì)量。在面形表征過(guò)程中，基函數(shù)有64項(xiàng)，采樣點(diǎn)數(shù)目為4 000個(gè)。

離軸三反系統(tǒng)結(jié)構(gòu)如圖6所示。其中三鏡的半口徑為50 mm，系統(tǒng)的有效焦距為110 mm，設(shè)計(jì)工作波長(zhǎng)為632.8 nm，視場(chǎng)為±1.2°×±1°，這里在二維視場(chǎng)中選取14個(gè)視場(chǎng)點(diǎn)對(duì)系統(tǒng)像質(zhì)進(jìn)行綜合考察。

圖6　離軸三反系統(tǒng)結(jié)構(gòu)圖

Fig.6Schematic diagram of Three Mirror Anastigmat(TMA) system

本文采用光學(xué)設(shè)計(jì)軟件Code V對(duì)系統(tǒng)進(jìn)行像質(zhì)優(yōu)化及評(píng)估。系統(tǒng)三鏡的初始結(jié)構(gòu)是采用兩種不同采樣點(diǎn)分布方式通過(guò)面形擬合得到的兩個(gè)RBF-BS自由曲面表面，面形擬合矢高差如圖7所示，自由曲面的面形參數(shù)是由軟件面形自定義模塊導(dǎo)入的。這里選取全視場(chǎng)在100 lp/mm 處的平均MTF、最大畸變量和全視場(chǎng)平均波像差作為評(píng)估系統(tǒng)成像質(zhì)量的標(biāo)準(zhǔn)。在相同的優(yōu)化條件與環(huán)境下，經(jīng)過(guò)相同的優(yōu)化步驟得到的系統(tǒng)MTF曲線圖8所示，兩種采樣點(diǎn)分布方式得到的像質(zhì)如表1所示。

(a)均勻采樣

(b)邊緣集中采樣

DistributionMaxdistortion/%AverageMTFYRMSWFE(λ)Uniformlyrandomdistribution2.3672.6140.18877Edge-clustereddistribution2.7641.0390.41716

(a)采用均勻采樣方式得到的初始面形對(duì)應(yīng)系統(tǒng)MTF曲線

(a)MTF curve of system with initial surface which generated by uniformly random grid

(b)采用邊緣集中采樣方式得到的初始面形對(duì)應(yīng)系統(tǒng)MTF曲線

(b)MTF curve of system with initial surface which generated by edge-clustered grid

圖8兩種初始面形用于離軸三反系統(tǒng)的像質(zhì)對(duì)比

Fig.8Comparison of image quality of two initial surfaces used in TMA system

由設(shè)計(jì)結(jié)果可知，采用均勻隨機(jī)型采樣方式得到的初始面形在進(jìn)行系統(tǒng)優(yōu)化后，MTF達(dá)到0.72以上，最大畸變量為2.36%。由于采用均勻隨機(jī)型采樣點(diǎn)分布方式的面形擬合精度更高，矢高誤差相對(duì)較小，所以得到的初始面形更接近于理想面形。因此，對(duì)于口徑內(nèi)矢高變化較為劇烈的三鏡，采用均勻隨機(jī)型采樣點(diǎn)分布方式擬合得到的初始面形更容易獲得理想像質(zhì)。這與上述理論分析結(jié)果相一致。

6　結(jié)　論

本文針對(duì)基于Gaussian RBF-BS模型研究了兩種不同的采樣點(diǎn)分布類型，即均勻隨機(jī)型采樣點(diǎn)分布和邊緣集中型采樣點(diǎn)分布，對(duì)Gaussian BF-BS模型擬合自由曲面的面形擬合精度的影響。通過(guò)比較不同采樣點(diǎn)分布條件下離軸二次曲面和帶凸起的拋物面的擬合結(jié)果可以發(fā)現(xiàn)，均勻隨機(jī)型采樣點(diǎn)分布的面形擬合精度比邊緣集中型采樣點(diǎn)分布的面形擬合精度高，并且達(dá)到一定的擬合精度后，采樣點(diǎn)的數(shù)目不再影響擬合精度，而過(guò)多的采樣點(diǎn)數(shù)量則會(huì)增加擬合運(yùn)算量，降低擬合效率。這一結(jié)果不同于已有的關(guān)于Zernike多項(xiàng)式模型的研究報(bào)道，對(duì)Gaussian RBF-BS模型在自由曲面表征和設(shè)計(jì)中的應(yīng)用具有一定的理論指導(dǎo)意義。最后，以離軸三反系統(tǒng)設(shè)計(jì)為例，分析比較由不同采樣方式生成初始面形后系統(tǒng)的像質(zhì)優(yōu)化結(jié)果。結(jié)果顯示，采用均勻隨機(jī)型采樣方式得到的初始面形在進(jìn)行系統(tǒng)優(yōu)化后，最終全視場(chǎng)的平均調(diào)制傳遞函數(shù)(MTF)達(dá)到0.72以上，遠(yuǎn)高于由邊緣集中采樣方式生成初始面形后系統(tǒng)像質(zhì)的優(yōu)化結(jié)果。

[1]THOMPSON K P， ROLLAND J P. Freeform optical surfaces： a revolution in imaging optical design [J].Opt.Pho.News， 2012， 23(6)： 30-35.

[2]FUERSCHBACH K， JANNICK P R， KEVIN P T. A new family of optical systems employingφ-polynomial surfaces [J].Opt.Express， 2011，19(22)： 21919-21928.

[3]ILHAN K， KEVIN P T， JANNICK P R. Edge clustered fitting grids forφ-polynomial characterization of freeform optical surfaces[J].OpticalExpress， 2011，19(27)： 26962-26974.

[4]ILHAN K， KEVIN P T， JANNICK P R. Comparative assessment of freeform polynomials as optical surface descriptions[J].Opt.Express， 2012，20(20)： 22683-22691.

[5]王超.自由曲面表征函數(shù)及其應(yīng)用研究[D].長(zhǎng)春：中國(guó)科學(xué)院長(zhǎng)春光學(xué)精密機(jī)械與物理研究所， 2014.

WANG CH.ResearchonCharacterizationFunctionandApplicationofFree-formSurface[D]. Changchun： Changchun Institute of Optics， Fine Mechanics and Physics， Chinese Academy of Sciences， 2014.(in Chinese)

[6]OZAN C， ILHAN K， GREGORY E F，etal.. Application of radial basis functions to represent optical freeform surfaces [J].SPIE-OSA，2010， 7652：76520A-1.

[7]OZAN C， GREGORY E F， HASSAN F，etal.. Meshfree approximation methods for free-form surface representation in optical design with applications to head-worn displays [J].SPIE，7061： 70610D.

[8]OZAN C， BRENDAN M， HASSAN F，etal.. Optimal local shape description for rotationally non-symmetric optical surface design and analysis[J].OpticalExpress， 2008，16(3)： 1583-1589.

[9]OZAN C， SOPHIE V， HASSAN F，etal.. Application of radial basis functions to shape description in a dual-element off-axis magnifier[J].OpticalLetter， 2008，33(11)： 1237-1239.

[10]趙星，鄭義，張贊，等. 基于面形斜率的光學(xué)自由曲面表征[J].光學(xué) 精密工程，2015，23(7)：1954-1964.

ZHAO X， ZHENG Y， ZHANG Z，etal.. Characterization of freeform optical surfaces based on surface slope [J].Opt.PrecisionEng.， 2015，23(7)：1954-1964.(in Chinese)

[11]張博，王凌，常偉軍，等. 頭盔顯示器自由曲面光學(xué)組件的優(yōu)化設(shè)計(jì)[J].應(yīng)用光學(xué)，2014，35(2)：193-197.

ZHANG B， WANG L，CCHANG W J，etal.. Optimal design of free-form-surface optical component in helmet mounted display[J].JournalofAppliedOptics，2014，35(2)：193-197.(in Chinese)

[12]晏蕾，賈平，洪永豐.紅外雙視場(chǎng)光學(xué)系統(tǒng)的初始結(jié)構(gòu)求解[J].紅外，2011，32(1)：23-26.YANG L， JIA P， HONG Y F. Design of initial structure of infrared dual field-of-view optical system [J].Infrared，2011，32(1)：23-26.(in Chinese)

[13]崔恩坤，張葆，洪永豐，等.PW 法對(duì)連續(xù)變焦光學(xué)系統(tǒng)初始結(jié)構(gòu)的求解[J].應(yīng)用光學(xué)，2014(4)：586-591.

CUI E K， ZHANG B， HONG Y F，etal.. Design of initial structure of infrared zoom optical system with PW solution [J].JournalofAppliedOptics， 2014(4)：586-591.(in Chinese)

[14]馮婕，白瑜，邢廷文，等.Zernike多項(xiàng)式波面擬合精度研究[J].光電技術(shù)應(yīng)用，2011，26(2)：31-34.FENG J， BAI Y， XING T W，etal..Fitting accuracy of wavefront using Zernike polynomials [J].Electro-OpticTechnologyApplication， 2011，26(2)：31-34.(in Chinese)

[15]孫學(xué)真，蘇顯渝，荊海龍，等.抽樣點(diǎn)對(duì)基于Zernike多項(xiàng)式曲面擬合精度的影響[J].光學(xué)儀器，2008，30(4)：6-10.SUN X ZH， SU X Y， XING H L，etal.. The influence of sampling points on the precision of curved surface fitting based on Zernike polynomials [J].OpticalInstruments， 2008，30(4)：6-10.(in Chinese)

[16]謝蘇隆.Zernike多項(xiàng)式擬合曲面中擬合精度與采樣點(diǎn)數(shù)目研究[J].應(yīng)用光學(xué)，2010，31(6)：943-949.XIE S L. Sampling point number in curved surface fitting with Zernike polynomials [J].JournalofAppliedOptics， 2010，31(6)：943-949.(in Chinese)

[17]FORBES G W. Shape specification for axially symmetric optical surfaces [J].Opt.Express，2007，15：5218-5226.

[18]FORBES G W. Robust， efficient computational methods for axially symmetric optical aspheres [J].Opt.Express，2010，18： 19700-19712.

[19]OZAN C， BRENDAN M， HASSAN F，etal.. Optimal local shape description for rotationally non-symmetric optical surface design and analysis[J].OpticalExpress， 2008，16(3)： 1583-1589.

[20]鄢靜舟，雷凡，周必方，等.用Zernike多項(xiàng)式進(jìn)行波面擬合的幾種算法[J].光學(xué) 精密工程，1999，7(5)： 119-128.

YAN J ZH， LEI F， ZHOU B F，etal.. Several algorithms for wavefront fitting using Zernike polynomails[J].Opt.PrecisionEng.， 1999，7(5)： 119-128.(in Chinese)

武鵬飛(1968-)，男，遼寧興城人，博士，教授，1996年于南開(kāi)大學(xué)獲得博士學(xué)位，主要從事光學(xué)信息處理等方面的研究。E-mail： pwu@nankai.edu.cn

趙星(1980-)，男，北京人，博士，副教授，2007年于南開(kāi)大學(xué)獲得博士學(xué)位，主要從事光學(xué)信息處理以及光學(xué)系統(tǒng)設(shè)計(jì)等領(lǐng)域的研究工作。E-mail： zhaoxingtjnk@nankai.edu.cn

(版權(quán)所有未經(jīng)許可不得轉(zhuǎn)載)

Influence of sampling point distribution in freeform surfaces fitting with radial based function model

WU Peng-fei1， ZHANG Zan1，ZHENG Yi1， ZHANG Juan1， WANG Yao-yao1， ZHAO Xing1*， WANG Ling-jie2， WU Yan-xiong2

(1.InstituteofModernOptics，KeyLaboratoryofOpticalInformationScienceandTechnology，MinistryofEducation，NankaiUniversity，Tianjin300071，China；2.KeyLaboratoryofOpticalSystemAdvancedManufacturingTechnology，ChangchunInstituteofOptics，F(xiàn)ineMechanicsandPhysics，ChineseAcademyofSciences，Changchun130033，China)

*Correspondingauthor，E-mail：zhaoxingtjnk@nankai.edu.cn

In view of the importance of surface shape fitting precision of free-form surface model in researches such as free-form surface representation and selection of initiating structure of surface shape, the influence of different sampling point distribution patterns on surface shape fitting precision of the Gaussian Radial Basis Function Based on Slope( Gaussian RBF-BS) model was investigated. An off-axis quadric surface and an bumpy paraboloid were fitted with different sampling point distributions which were uniformly random points and edged-clustered points. The results show that the uniformly random grid is better for surface shape fitting due to the high precision. The influence of sampling point quantity on fitting precision remains constant when certain fitting precision is reached. The off-axis Three Mirror Anastigmat(TMA) system was optimized experimentally by using initial surface shapes generated with different sampling methods. The results show that, by using initial surface shape generated from uniformly random sampling method for the system optimization, the average Modulation Transfer Function (MTF) of whole field of view is more than 0.72, which is much higher than that by using surface shape generated by edge cluster sampling points. Thus the theoretical research results were verified.

optical design; freeform surface; fitting precision; distribution of sampling points

2016-03-04；

2016-04-29.

國(guó)家自然科學(xué)基金資助項(xiàng)目(No.11474169)；中科院科技創(chuàng)新基金重點(diǎn)項(xiàng)目；天津市應(yīng)用基礎(chǔ)與前沿技術(shù)研究計(jì)劃資助項(xiàng)目(No.15JCYBJC16900)

1004-924X(2016)07-1564-09

TH703

Adoi:10.3788/OPE.20162407.1564