(高清版)GB∕T 44221-2024 光學(xué)系統(tǒng)波前像差的測定 夏克-哈特曼光電測量法

GB/T44221—2024光學(xué)系統(tǒng)波前像差的測定夏克-哈特曼光電測量法2024-07-24發(fā)布2025-02-01實施國家市場監(jiān)督管理總局國家標準化管理委員會前言 1范圍 13術(shù)語和定義 4測量原理及方法 24.1測量原理 24.2光學(xué)系統(tǒng)波前像差測量方法 2 3 45.1測量環(huán)境 45.2樣品 4 46.1測量儀 46.2輔助鏡頭 5 7.1測量前準備 5 67.3光路對準 6 68測量數(shù)據(jù)處理 6 7附錄A(資料性)波前復(fù)原方法 8 附錄C(資料性)測量報告 參考文獻 TⅢ1光學(xué)系統(tǒng)波前像差的測定夏克-哈特曼光電測量法本文件描述了采用夏克-哈特曼法測量光學(xué)系統(tǒng)波前像差的原理及方法、測量條件、設(shè)備及裝置、測量步驟以及測量數(shù)據(jù)處理。本文件適用于采用夏克-哈特曼法測量光學(xué)系統(tǒng)波前像差的測量，也適用于光學(xué)零件面形偏差的測量。2規(guī)范性引用文件本文件沒有規(guī)范性引用文件。3術(shù)語和定義下列術(shù)語和定義適用于本文件。光波傳播時的等相位面。Φ波前與理想波前的偏差。被測光學(xué)表面相對于參考光學(xué)表面的偏差。通過子孔徑的斜率計算得到入射波前的相位分布的過程。儀器能夠檢測的光學(xué)零件或系統(tǒng)的通光孔徑。使平行光管發(fā)出的平行光照射在試樣上，再由試樣反射回平行光管，根據(jù)焦點附近像的情況測定2ΦpvΦms32——被測望遠系統(tǒng)；3——輔助平面鏡。圖2自準直法測量光學(xué)系統(tǒng)波前像差示意圖采用直接法進行測量時，測量儀輸出的平行光通過被測光學(xué)系統(tǒng)后直接進入測量儀的夏克-哈特曼傳感器，測量結(jié)果即為光學(xué)系統(tǒng)的波前像差。圖3以望遠系統(tǒng)為例，給出了測量儀采用直接法測量光學(xué)系統(tǒng)波前像差的示意圖。1——測量儀的光源與系統(tǒng)光路部分；2——被測望遠系統(tǒng)；3——測量儀的夏克-哈特曼傳感器。圖3直接法測量光學(xué)系統(tǒng)波前像差示意圖4.3光學(xué)零件面形偏差的測量測量球面光學(xué)零件面形偏差采用自準直法。測量儀輸出的平行光經(jīng)輔助鏡頭匯聚到待測球面的球心后被待測球面反射，再經(jīng)輔助鏡頭平行入射到測量儀，得到待測球面的面形偏差，其結(jié)果是待測面形偏差的兩倍。測量凹球面光學(xué)零件面形偏差的示意圖見圖4。測量前，應(yīng)使用標準球面鏡標定帶有輔助鏡頭的測量儀。2——輔助鏡頭；圖4測量凹球面光學(xué)零件面形偏差的示意圖4b)相對濕度：不大于75%;5GB/T44221—202427898——微透鏡陣列；10——軟件與顯示單元；12——系統(tǒng)光路。56圖6測量儀示意圖包含輔助透鏡組和輔助反射鏡等附件。宜對帶有輔助鏡頭的測量儀進行自校；如果無法實現(xiàn)自輔助透鏡組的F數(shù)(透鏡組焦距與有效孔徑之比)應(yīng)不大于被測樣品的曲率半徑與通光孔徑的比值，且越接近越好。輔助透鏡組的后焦距應(yīng)大于被測凸球面的曲率半徑。輔助反射鏡口徑應(yīng)大于被測光學(xué)系統(tǒng)的出射口徑。7.1測量前準備測量前準備工作如下：6圖見圖7;8測量數(shù)據(jù)處理77GB/T44221—2024j——記錄數(shù)據(jù)的組數(shù)；9測量報告8GB/T44221—2024(資料性)波前復(fù)原方法A.1概述測量儀通過測量波前局部斜率來復(fù)原入射波前，波前復(fù)原包括質(zhì)心計算、斜率計算和波前重構(gòu)三個步驟。質(zhì)心計算是確定光斑的質(zhì)心位置坐標，斜率計算是由理論光斑質(zhì)心位置與實際光斑質(zhì)心位置的偏差，以及微透鏡陣列的子孔徑和焦距參數(shù)，計算出波前的局部斜率。波前重構(gòu)是根據(jù)子孔徑的斜率計算得到入射波前的相位分布。實現(xiàn)波前重構(gòu)的算法有多種，常用的有模式重構(gòu)法和區(qū)域重構(gòu)法。A.2模式重構(gòu)法模式重構(gòu)法將全孔徑內(nèi)的波前相位展開成不同的模式(例如：平移、傾斜、離焦、像散、彗差和球差等),用全孔徑內(nèi)的測量數(shù)據(jù)去求解各模式的系數(shù)，得到完整的波前展開式，重構(gòu)出入射波前。模式重構(gòu)法常用的展開多項式是Zernike正交多項式。當被測波前在圓域內(nèi)時按公式(A.1)用Zernike多項式描述。a?——平均相位波前；n——Zernike多項式系數(shù)的模式階數(shù)；ak——第k項Zernike多項式系數(shù)；Z(x,y)——第k項Zernike多項式；E——波前相位測量誤差。測量儀第i個子孔徑內(nèi)的平均斜率與Zernike多項式系數(shù)的關(guān)系可用公式(A.2)和公式(A.3)表示。式中：(A.2)(A.3)G(i)——第i個子孔徑在x方向的平均斜率；Z(i)——第k項Zernike多項式在x方向上的分量；e,(i)——波前相位第i個子孔徑內(nèi)在x方向上的測量誤差；G(i)——第i個子孔徑在y方向的平均斜率；Z(i)——第k項Zernike多項式在y方向上的分量；設(shè)測量儀有m個子孔徑，并取模式函數(shù)系列Z(x,y)的前n項進行波前重構(gòu)，則模式重構(gòu)的矩陣表達式為公式(A.4),可用公式(A.5)表示。GB/T44221—2024G=DA+ε……….(A.5)G——波前相位斜率向量，包括測量儀所有子孔徑測量的x和y方向平均斜率；A——Zernike系數(shù)矩陣。測量儀測量得到波前相位斜率向量G后，求出D的廣義逆D+,解出Zernike系數(shù)矩陣A的值，計將每個子孔徑上波前相位離散化，用若干個點取代連續(xù)面，一個完整波前被細分成一定數(shù)量的子區(qū)間(子孔徑)。區(qū)域重構(gòu)法利用各子孔徑邊界上測量的波前梯度或相位差數(shù)據(jù)，重構(gòu)整個波前相位。根據(jù)測量參數(shù)的性質(zhì)(梯度或相位差)、重構(gòu)波前相位的位置和重構(gòu)的算法的不同，有多種重構(gòu)波前的方法。按照相位測量點和重構(gòu)點相對位置不同，有休晉模型、弗雷德模型和紹契威爾模型三種重要重圖A.1休晉模型9圖A.2弗雷德模型圖A.3紹契威爾模型GB/T44221—2024(資料性)B.1概述Zernike多項式是由無窮數(shù)量的多項式完全集組成的，用極坐標表述，徑向坐標r和角度坐標θZernike多項式在單位圓內(nèi)部是連續(xù)正交的，但在單位圓內(nèi)部的離散的坐標上是不具備正交性質(zhì)的。Zernike多項式和光學(xué)檢測中觀測到的像差多項式的形式是一致的，常被用來描述波前特性。在B.2給出的48項Zernike多項式中，Z?項是平移項，是第0組系數(shù)；Z?到Z?項是第2組系數(shù)，代表波前的高Z4和Z?項代表了水平和豎直方向的像散，Z?和Z?項代表水平和豎直方向的彗差，Zg項代表了一階球差；Z?到Z??項是第6組系數(shù)；Z?到Z?項是第8組系數(shù)，Z?5到Z?5項是第10組系數(shù)，Z??到Z?8項是第12組系數(shù)。B.2前48項Zernike多項式數(shù)學(xué)表達式參見ISO/TR14999-2:2019,給出前48項Zernike多項式序列數(shù)學(xué)表達式如下：Z?(r,θ)=1Z?(r,θ)=rcosθZ?(r,θ)=rsinθZ?(r,θ)=2r2—1Z?(r,θ)=r2cos20Z?(r,θ)=(3r2—2)rcosθZ?(r,θ)=(3r2—2)rsinθZ?(r,θ)=6r?—6r2+1Z?(r,θ)=r3cos30Z(r,θ)=(4r2—3)r2cos20Z12(r,θ)=(4r2—3)r2sin20Z(r,θ)=(10r?—12r2+3)rcosθZ14(r,θ)=(10r?—12r2+3)rsinθZ?(r,θ)=20r?—30r?+12r2—1Z?(r,θ)=r?cos40Z??(r,θ)=r?sin40Z??(r,θ)=(5r2—4)r3cos30Z??(r,θ)=(5r2—4)r3sin30Z?(r,θ)=(15r?—20r2+6)r2cos20Z,(r,θ)=(15r?—20r2+6)r2sin20Z?2(r,θ)=(35r?—60r?+30r2—4)rcosθZ?3(r,θ)=(35r?—60r?+30r2—4)rsinθZ??(r,θ)=70r?—140r?+90r?—20r2+1Z??(r,θ)=r?cos50Z??(r,θ)=r?sin50Z??(r,θ)=(6r2—5)r?cos40Z??(r,θ)=(6r2—5)r?sin40Z??(r,θ)=(21r?—30r2+10)r3cos30Z??(r,θ)=(21r?—30r2+10)r3sin30Z?(r,θ)=(56r?—105r?+60r2—10)r2cos20Z?(r,θ)=(56r?—105r?+60r2—10)r2sin20Z??(r,θ)=(126r?—280r?+210r?—60r2+5)rcosθZ??(r,θ)=(126r?—280r?+210r?—60r2+5)rsinθZ?5(r,θ)=252r1?—630r?+560r?—210r?+30r2-1Z??(r,θ)=r?cos60Z?7(r,θ)=r?sin60Z??(r,θ)=(7r2—6)r?cos50Z?(r,0)=(7r2—6)r?sin50Z??(r,θ)=(28r?—42r2+15)r?cos40Z??(r,θ)=(28r?—42r2+15)r?sin40Z?(r,θ)=(84r?—168r?+105r2—20)r3cos30Z?3(r,θ)=(84r?—168r?+105r2—20)r3sin30Z(r,θ)=(210r?—504r?+420r?—140r2+15)r2cos20Z??(r,θ)=(210r?—504r?+420r?—140r2+15)r2sin20Z?6(r,θ)=(462r10—1260r?+1260r?—560r?+105r2—6)rcosθZ??(r,θ)=(462r10—1260r?+1260r?—560r?+105r2—6)rsinθZ??(r,θ)=924r12—2772r1?+3150r?—1680r?