TW201935071A - 自由曲面離軸三反成像系統 - Google Patents
自由曲面離軸三反成像系統 Download PDFInfo
- Publication number
- TW201935071A TW201935071A TW107103160A TW107103160A TW201935071A TW 201935071 A TW201935071 A TW 201935071A TW 107103160 A TW107103160 A TW 107103160A TW 107103160 A TW107103160 A TW 107103160A TW 201935071 A TW201935071 A TW 201935071A
- Authority
- TW
- Taiwan
- Prior art keywords
- axis
- free
- coordinate system
- dimensional rectangular
- rectangular coordinate
- Prior art date
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 87
- 229910052727 yttrium Inorganic materials 0.000 claims description 16
- 238000013519 translation Methods 0.000 claims description 3
- 210000001747 pupil Anatomy 0.000 claims description 2
- 239000013598 vector Substances 0.000 description 26
- 238000010276 construction Methods 0.000 description 22
- 238000000034 method Methods 0.000 description 21
- 238000010586 diagram Methods 0.000 description 8
- 230000003287 optical effect Effects 0.000 description 8
- 238000013461 design Methods 0.000 description 6
- 230000004075 alteration Effects 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 230000000717 retained effect Effects 0.000 description 3
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N Silicium dioxide Chemical compound O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 239000007769 metal material Substances 0.000 description 2
- 238000002310 reflectometry Methods 0.000 description 2
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 description 1
- 229910052782 aluminium Inorganic materials 0.000 description 1
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 1
- 229910052802 copper Inorganic materials 0.000 description 1
- 239000010949 copper Substances 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 1
- 239000010931 gold Substances 0.000 description 1
- 229910052737 gold Inorganic materials 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 229910052755 nonmetal Inorganic materials 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- HBMJWWWQQXIZIP-UHFFFAOYSA-N silicon carbide Chemical compound [Si+]#[C-] HBMJWWWQQXIZIP-UHFFFAOYSA-N 0.000 description 1
- 229910010271 silicon carbide Inorganic materials 0.000 description 1
- 235000012239 silicon dioxide Nutrition 0.000 description 1
- 239000000377 silicon dioxide Substances 0.000 description 1
- 238000001931 thermography Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H10—SEMICONDUCTOR DEVICES; ELECTRIC SOLID-STATE DEVICES NOT OTHERWISE PROVIDED FOR
- H10F—INORGANIC SEMICONDUCTOR DEVICES SENSITIVE TO INFRARED RADIATION, LIGHT, ELECTROMAGNETIC RADIATION OF SHORTER WAVELENGTH OR CORPUSCULAR RADIATION
- H10F39/00—Integrated devices, or assemblies of multiple devices, comprising at least one element covered by group H10F30/00, e.g. radiation detectors comprising photodiode arrays
- H10F39/80—Constructional details of image sensors
- H10F39/806—Optical elements or arrangements associated with the image sensors
- H10F39/8067—Reflectors
-
- G—PHYSICS
- G02—OPTICS
- G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
- G02B17/00—Systems with reflecting surfaces, with or without refracting elements
- G02B17/02—Catoptric systems, e.g. image erecting and reversing system
- G02B17/06—Catoptric systems, e.g. image erecting and reversing system using mirrors only, i.e. having only one curved mirror
- G02B17/0626—Catoptric systems, e.g. image erecting and reversing system using mirrors only, i.e. having only one curved mirror using three curved mirrors
-
- G—PHYSICS
- G02—OPTICS
- G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
- G02B17/00—Systems with reflecting surfaces, with or without refracting elements
- G02B17/02—Catoptric systems, e.g. image erecting and reversing system
- G02B17/06—Catoptric systems, e.g. image erecting and reversing system using mirrors only, i.e. having only one curved mirror
- G02B17/0626—Catoptric systems, e.g. image erecting and reversing system using mirrors only, i.e. having only one curved mirror using three curved mirrors
- G02B17/0642—Catoptric systems, e.g. image erecting and reversing system using mirrors only, i.e. having only one curved mirror using three curved mirrors off-axis or unobscured systems in which not all of the mirrors share a common axis of rotational symmetry, e.g. at least one of the mirrors is warped, tilted or decentered with respect to the other elements
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N25/00—Circuitry of solid-state image sensors [SSIS]; Control thereof
- H04N25/60—Noise processing, e.g. detecting, correcting, reducing or removing noise
- H04N25/61—Noise processing, e.g. detecting, correcting, reducing or removing noise the noise originating only from the lens unit, e.g. flare, shading, vignetting or "cos4"
- H04N25/615—Noise processing, e.g. detecting, correcting, reducing or removing noise the noise originating only from the lens unit, e.g. flare, shading, vignetting or "cos4" involving a transfer function modelling the optical system, e.g. optical transfer function [OTF], phase transfer function [PhTF] or modulation transfer function [MTF]
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Optics & Photonics (AREA)
- Engineering & Computer Science (AREA)
- Multimedia (AREA)
- Signal Processing (AREA)
- Lenses (AREA)
Abstract
本發明涉及一種自由曲面離軸三反成像系統包括:一主鏡,一次鏡,一三鏡,以及一探測器,光線依次經過所述主鏡、次鏡和三鏡的反射後,被所述探測器接收到並成像,該自由曲面離軸三反成像系統的視場角大於等於60°×1°,且F數小於等於2.5。
Description
本發明涉及光學設計領域,尤其涉及一種自由曲面離軸三反成像系統。
自由曲面是指無法用球面或非球面係數來表示的非傳統曲面,通常是非回轉對稱的,結構靈活,變數較多,為光學設計提供了更多的自由度,可以大大降低光學系統的像差,減小系統的體積、重量與鏡片數量,可以滿足現代成像系統的需要,有著廣闊的發展應用前景。由於自由曲面有非對稱面並提供了更多的設計自由度,他們常被用在離軸非對稱系統中。
先前的自由曲面離軸三反光學系統主要針對大視場或者小F數高解析度成像,無法同時實現大視場角和小F數,不能在大的成像範圍內同時得到高解析度圖像。
綜上所述,確有必要提供一種可以同時實現大視場角和小F數的自由曲面離軸三反成像系統。
一種自由曲面離軸三反成像系統,包括:一主鏡,以所述主鏡的頂點為第一原點定義一第一三維直角坐標系(X,Y,Z),在該第一三維直角坐標系(X,Y,Z)中,所述主鏡的反射面為一xy多項式自由曲面;一次鏡,以所述次鏡的頂點為第二原點定義一第二三維直角坐標系(X',Y',Z'),且該第二三維直角坐標系(X',Y',Z')由所述第一三維直角坐標系(X,Y,Z)沿Z軸負方向和Y軸負方向平移得到,在該第二三維直角坐標系(X',Y',Z')中,所述次鏡的反射面為一x'y'多項式自由曲面;一三鏡,以所述三鏡的頂點為第三原點定義一第三三維直角坐標系(X'',Y'',Z''),且該第三三維直角坐標系(X'',Y'',Z'')由所述第二三維直角坐標系(X',Y',Z')沿Z'軸正方向和Y'軸正方向平移得到,在該第三三維直角坐標系(X'',Y'',Z'')中,所述三鏡的反射面為一x''y''多項式自由曲面;以及一探測器,光線依次經過所述主鏡、次鏡和三鏡的反射後,被該探測器接收到並成像,其中,該自由曲面離軸三反成像系統的視場角大於等於60°×1°,且F數小於等於2.5。
相較於先前技術,本發明提供的自由曲面離軸三反光學系統的視場角大於等於60°×1°的同時F數小於等於2.5,因此,本發明提供的自由曲面離軸三反光學系統可以同時實現大視場角和小F數的成像,能在大的成像範圍內同時得到高解析度圖像。
下面將結合附圖及具體實施例對本發明作進一步的詳細說明。
本發明提供一種自由曲面離軸成像系統的設計方法,其包括以下步驟:
步驟一,建立初始系統並選取特徵視場。
所述初始系統包括複數個初始曲面,且該初始系統中的一個初始曲面對應待設計的自由曲面離軸成像系統中一個曲面。所述複數個初始曲面可以為平面、球面等。所述複數個初始曲面的具體位置根據待設計的自由曲面離軸成像系統的實際需要進行選擇。所述初始系統中初始曲面的數量根據實際需要進行設計。本實施例中,所述初始系統為一初始平面三反系統,該初始平面三反系統包括三個初始平面。
所述特徵視場的數量遠大於初始系統中初始曲面的數量。優選的,所述特徵視場數量為初始曲面數量的50倍以上。更優選的,所述特徵視場數量為初始曲面數量的50~60倍。所述特徵視場可以根據實際需要進行選擇。優選的,所述選取特徵視場可以在弧失方向和子午方向等距採樣M個特徵視場,分別記作φ1
,φ2
,……φM
,假設中心視場為(0,φy0
)。本實施例中,所述初始三反平面系統在子午方向工作在8°到9°視場,在弧失方向工作在−30°到30°視場。請參閱圖1,所述初始三反平面系統關於YOZ平面對稱,在弧失方向只考慮一半的視場,每隔0.5°採樣一個特徵視場,在子午方向採樣三個特徵視場,在弧失方向採樣61個特徵視場,一共選取183個特徵視場。
在所述M個特徵視場中的每個特徵視場中均選取K條特徵光線。所述K條特徵光線的選取可通過以下方法進行:
將M個視場中的每個視場的孔徑分成N等份,並從每一等份中選取不同孔徑位置的P條特徵光線,這樣一共選取了K=M×N×P條對應不同視場不同孔徑位置的特徵光線。所述孔徑可以為圓形、長方形、正方形、橢圓形或其他規則或不規則的形狀。優選的,所述視場孔徑為圓形孔徑,將每個視場的圓形孔徑等分成N個角度,間隔為φ,因此有N=2π/φ,沿著每個角度的半徑方向取P個不同的孔徑位置,那麼一共取K=M×N×P條對應不同視場不同孔徑位置的特徵光線。本實施例中,將每個視場的孔徑分成14等份,並從每一等份中選取不同孔徑位置的7條特徵光線,這樣每個特徵視場採樣98條特徵光線。
步驟二、逐漸拓展構建特徵視場,同時將所述初始系統構建為自由曲面系統。
通常自由曲面離軸成像系統在弧失方向對稱,因此構建時只需考慮弧失方向半視場的特徵光線。請參閱圖2,優選的,選取所述中心視場(0,φy0
)作為第一被構建的特徵視場φ1
,在後續的拓展構建特徵視場過程中,子午方向上構建特徵視場沿著正負兩個方向拓展,弧失方向上構建特徵視場沿弧失正方向拓展。本實施例中,由於初始平面三反系統在子午方向僅有1°視場,因此,在拓展構建特徵視場過程中,只對弧失方向進行拓展構建特徵視場,子午方向三個特徵視場在每步構建時同時考慮。
請參閱圖3,將初始系統中初始曲面定義為Lj
(j=1,2,……L),即初始曲面的個數為L個,所述逐漸拓展構建特徵視場,同時將所述初始系統構建為自由曲面系統具體包括以下步驟: S21:在所述M個特徵視場φi
(i=1,2,……M)中的至少一個特徵視場作為第一構建特徵視場,在該第一構建特徵視場中將所述初始曲面L1
構建成自由曲面N1
; S22,增加至少一個特徵視場對所述第一構建特徵視場進行拓展得到一第二構建特徵視場,在該第二構建特徵視場中將初始曲面L2
構建成自由曲面N2
;以及 S23,以此類推,每次增加至少一個特徵視場對構建特徵視場進行拓展,並使用拓展後的構建特徵視場將初始曲面構建成自由曲面,直到L個初始曲面構建成L個自由曲面。
步驟S21中,所述使用特徵視場φ1
將所述初始曲面L1
構建成自由曲面N1
包括:計算得到所述自由曲面N1
上的所有特徵數據點Pi
(i=1,2…K);以及將該特徵數據點進行曲面擬合得到所述自由曲面N1
的方程式。
為了得到自由曲面N1
上的所有特徵數據點Pi
(i=1,2…K),將借助特徵光線Ri
(i=1,2…K)與自由曲面N1
的前一個曲面及後一個曲面的交點。在求解每條特徵光線Ri
(i=1,2…K)對應的自由曲面N1
上的特徵數據點Pi
(i=1,2…K)時,將該特徵光線Ri
與前一個曲面的交點定義為該特徵光線的起點Si
,特徵光線Ri
與後一個曲面的交點定義為該特徵光線的終點Ei
。當待設計的系統與特徵光線確定後,該特徵光線Ri
的起點Si
是確定的,且易於通過光線追跡即物像關係得到,特徵光線的終點Ei
可通過物像關係求解。在理想狀態下,特徵光線Ri
從Si
射出後,經過Pi
,交於Ei
,並最終交目標面於其理想目標點Ti,ideal
。如果自由曲面N1
是目標面,特徵光線的終點就是其目標點Ti
;如果在自由曲面N1
和目標面之間還有其他面,特徵光線的終點是Ω''上從特徵數據點到其對應的目標點的光程的一階變分為零的點,即:。
所述自由曲面N1
上特徵數據點Pi
(i=1,2…K)可以通過以下計算方法獲得: 步驟a,取定第一構建特徵視場中的第一條特徵光線R1
與所述自由曲面N1
對應的初始曲面的第一交點為特徵數據點P1
; 步驟b,在得到第i(1≦i≦K−1)個特徵數據點Pi
後,根據斯涅爾定律的向量形式求解第i個特徵數據點Pi
處的單位法向量,進而求得Pi
處的單位切向量; 步驟c,僅過所述第i(1≦i≦K−1)個特徵數據點Pi做一第一切平面並與特徵視場φ1
中其餘K−i條特徵光線相交,得到K−i個第二交點,從該K−i個第二交點中選取出與所述第i個特徵數據點Pi
距離最短的第二交點Qi+1
,並將其對應的特徵光線及與所述第i個特徵數據點Pi
的最短距離分別定義為Ri+1
和D; 步驟d,過特徵數據點Pi
(1≦i≦K−1)之前已求得的i−1個第一特徵數據點分別做一第二切平面,得到i−1個第二切平面,該i−1個第二切平面與所述特徵光線Ri+1
相交得到i−1個第三交點,在每一第二切平面上每一第三交點與其所對應的特徵數據點Pi
形成一交點對,在所述交點對中,選出交點對中距離最短的一對,並將距離最短的交點對的第三交點和最短距離分別定義為Q(i+1)
'和Di
'; 步驟e,比較Di
與Di
',如果Di
≦Di
',則把Qi+1
取為下一個特徵數據點Pi+1
,反之,則把Q(i+1)
'取為下一個特徵數據點Pi+1
;以及 步驟f,重複步驟b到e,直到計算得到自由曲面N1
上的所有特徵數據點Pi
(i=1,2…K),通過曲面擬合可以得到所述自由曲面N1
對應的的方程式。
步驟b中,每個特徵數據點Pi處的單位法向量可以根據斯涅爾(Snell)定律的向量形式求解。當待求的自由曲面Ω為折射面時,則每個特徵數據點Pi (i=1,2…K)處的單位法向量滿足:(1), 其中,、分別是沿著特徵光線入射與出射方向的單位向量,n, n' 分別為自由曲面N1
前後兩種介質的折射率。
類似的,當自由曲面N1
為反射面時,則每個特徵數據點Pi
(i=1,2…K)處的單位法向量滿足:(2)。
由於,所述特徵數據點Pi
(i=1,2…K)處的單位法向量與所述特徵數據點Pi
(i=1,2…K) 處的切平面垂直。故,可以得到特徵數據點Pi
(i=1,2…K) 處的切平面。
將所述初始系統所在的空間定義一第一三維直角坐標系。優選的,將光束傳播方向定義為z軸,垂直於z軸方向的平面為xoy平面。
將所述自由曲面N1
上的複數個特徵數據點Pi
(i=1,2…K)進行曲面擬合包括以下步驟: S211:將所述複數個特徵數據點Pi
(i=1,2…K)在所述第一直角坐標系中擬合成一球面,得到所述球面的曲率c及其對應的曲率中心(xc
, yc
, zc
); S212:將中心採樣視場主光線對應的特徵數據點(xo
, yo
, zo
)定義為球面的頂點,並以該球面的頂點為原點,過曲率中心與球面頂點的直線為z軸,建立一第二三維直角坐標系; S213:將所述複數個特徵數據點Pi
(i=1,2…K)在第一三維直角坐標系中的坐標(xi
, yi
, zi
)及其對應的法向量(αi
, βi
, γi
)變換為第二三維直角坐標系中的坐標(x'i
, y'i
, z'i
)及其法向量(α'i
, β'i
, γ'i
); S214:根據所述複數個特徵數據點Pi
(i=1,2…K)在第二三維直角坐標系中的坐標(x'i
, y'i
, z'i
),以及步驟S42中求得的球面的曲率c,將特徵數據點Pi
(i=1,2…K)在第二三維直角坐標系中擬合成一個二次曲面,得到二次曲面係數k;以及 S215:將所述複數個特徵數據點Pi
(i=1,2…K)在第二三維直角坐標系中的二次曲面上的坐標與法向量分別從坐標(x'i
, y'i
, z'i
)和法向量(α'i
, β'i
, γ'i
)中排除掉,得到殘餘坐標與殘餘法向量,將該殘餘坐標與殘餘法向量擬合成一個自由曲面,該自由曲面的方程式與所述二次曲面的方程式相加即可得到所述自由曲面N1
的方程式。
通常的,成像系統關於三維直角坐標系的yoz面對稱,因此,步驟S211中,所述球面在第二三維直角坐標系yoz面內相對於在第一三維直角坐標系yoz面內的傾斜角θ為:。
所述複數個特徵數據點Pi
(i=1,2…K)在第二三維直角坐標系中的坐標(x'i
, y'i
, z'i
)與法向量(α'i
, β'i
, γ'i
)與在第一三維直角坐標系中的坐標(xi
, yi
, zi
)和法向量(αi
, βi
, γi
)的關係式分別為:和。
在第二三維直角坐標系中,將在二次曲面上的特徵數據點的坐標與法向量分別定義為(x'i
, y'i
, z'is
)和(α'is
, β'is
, γ'is
)。將法向的z'分量歸一化為-1,將在二次曲面上的特徵數據點坐標(x'i
, y'i
, z'is
)與法向量(α'is
, β'is
, γ'is
)分別從坐標(x'i
, y'i
, z'i
)和法向量(α'i
, β'i
, γ'i
)排除掉之後,得到殘餘坐標(x''i
, y''i
, z''i
)與殘餘法向量(α''s
, β''s
, -1)分別為:和。
步驟S215中,所述將殘餘坐標與殘餘法向量擬合得到自由曲面的步驟包括: 在所述第二三維直角坐標系中,用排除掉二次曲面項的多項式曲面作為待構建自由曲面的面形描述,即:, 其中,為多項式的某一項,為係數集合; 獲得殘餘坐標擬合誤差d1
(P),即所述殘餘坐標值(x''i
, y''i
, z''i
)(i = 1, 2, …, n)與所述自由曲面在z'軸方向殘餘坐標差值的平方和;以及殘餘法向向量擬合誤差d2
(P),即所述殘餘法向量Ni
= (α''i
, β''i
, -1)(i = 1, 2, …, n) 與所述自由曲面法向量的向量差的模值的平方和,;, 其中,,,,; 獲得評價函數,其中,w為權重且大於0; 選擇不同的權重w,並令所述評價函數的梯度,從而獲得多組不同的P及其對應的複數個自由曲面面形z = f (x, y; P);以及 獲得具有最佳的成像品質的最終自由曲面面形。
可以理解,初始系統中的其它初始曲面構建成自由曲面的方法與構建自由曲面N1
的方法基本相同,其不同盡在於選取的構建特徵視場不同。在構建方法中計算其他自由曲面上的特徵數據點,以及其他自由曲面上的特徵數據點進行曲面擬合的方法均與計算自由曲面N1
上特徵數據點的方法以及曲面擬合方法相同。
優選的,步驟S21中,在所述M個特徵視場φi
(i=1,2,……M)中選取一個特徵視場作為第一構建特徵視場;步驟S22和S23中每次增加一個特徵視場對構建特徵視場進行構建。即,步驟二的具體步驟優選包括:S21,使用特徵視場φ1
將所述初始曲面L1
構建成自由曲面N1
;S22:增加特徵視場φ2
,在特徵視場φ1
和特徵視場φ2
中將所述初始曲面L2
構建成自由曲面N2
;以及S23,以此類推,每次增加一個特徵視場對構建時的特徵視場進行拓展,使用拓展後的特徵視場將初始曲面構建成自由曲面,直到所有的初始曲面均構建成自由曲面。
步驟三、對所述自由曲面系統中的自由曲面的構建區域進行延展,並在延展後的構建區域中重新構建自由曲面,進而得到待設計的自由曲面離軸成像系統。
請參閱圖4,步驟三具體包括以下步驟: S31,增加至少一個特徵視場對步驟二中拓展後的構建視場進行延展得到一延展後的第I構建特徵視場,在該延展後的第I構建特徵視場中將自由曲面N1
進行構建得到自由曲面N'1
; S32,在延展後的第I構建特徵視場的基礎上增加至少一個特徵視場得到一延展後的第II構建特徵視場,在該延展後的第II構建特徵視場中將自由曲面N2
進行構建得到自由曲面N'2
; S33,以此類推,每次增加至少一個特徵視場對構建特徵視場進行延展,使用延展後的構建特徵視場將所述自由曲面系統中的L個自由曲面構建成L個新的自由曲面,然後多次迴圈該步驟,直到所述M個特徵視場均被使用完;以及 S34,所有M個特徵視場均被使用完之後,M個特徵視場中仍然有至少一個特徵視場在至少一個自由曲面上沒有被構建使用,使用所有M個特徵視場重新構建該至少一個自由曲面。
請參閱圖4,在某實施例中,當自由曲面系統中的自由曲面Nk
被重新構建成自由曲面N'k
時,所述M個特徵視場均被使用完,此時採用所有M個特徵視場重新構建該自由曲面N'k
之外的其它自由曲面。
步驟S31中,根據物像關係及斯涅爾定律逐點求解所述延展後的第I特徵視場中的特徵光線與自由曲面N1
的複數個交點,該複數個交點作為自由曲面N'1
的特徵數據點,依據物像關係分別求解每個特徵數據點處的法向量,將該複數個交點進行曲面擬合得到自由曲面N'1
的方程式。該將複數個交點進行曲面擬合得到自由曲面N'1
的方程式的擬合方法與步驟二中將特徵數據點Pi
(i=1,2…K)擬合得到自由曲面N1
的方法相同。
優選的,每次增加一個特徵視場對所述自由曲面系統中的自由曲面的構建區域進行延展。步驟三優選具體包括:S31,增加特徵視場φL+1
,在特徵視場φ1
、特徵視場φ1
、……以及特徵視場φL+1
中將所述自由曲面系統中的自由曲面N1
進行重新構建得到自由曲面N'1
;S32,增加特徵視場φL+2
,在特徵視場φ1
、特徵視場φ2
、……以及特徵視場φL+2
中將自由曲面系統中的自由曲面N2
進行重新構建得到自由曲面N'2
;以及,以此類推,每次增加一個特徵視場對構建區域進行延展,使用延展後的構建特徵視場將對自由曲面進行重新構建,直至選取的特徵視場均被使用完。
請參閱圖5,為所述自由曲面離軸成像系統的整個設計流程。
所述自由曲面離軸成像系統的設計方法可進一步包括一對步驟三中得到的自由曲面離軸成像系統進行優化的步驟。具體地,將步驟三中得到的自由曲面離軸成像系統作為後續優化的初始系統。可以理解,該優化的步驟並不是必需的,可以根據實際需要設計。
為了提高構建時數值計算穩定性,也可以首先採用上述方法構建一第一自由曲面離軸成像系統,在構建完成後再將該第一離軸三反成像系統放大一定的倍數得到所述自由曲面離軸成像系統。
所述自由曲面離軸成像系統中自由曲面的求解順序不限,可以根據實際需要進行調換。請參閱圖6,本實施例中,首先使用特徵視場φ1
=0.5°構建三鏡,然後逐漸增加構建特徵視場,按照三鏡-次鏡-主鏡的順序依次構建,構建完成後,進行自由曲面延展過程,延展過程的重構順序與構建順序一致。
本發明提供的自由曲面離軸成像系統的設計方法中在不斷拓展構建特徵視場同時輪流構建自由曲面,並通過延展構建區域重新構建自由曲面,該方法可以得到同時具有大相對孔徑和寬視場角的自由曲面離軸成像系統。而且,在拓展構建特徵視場的過程中可以對已構建特徵視場進行像差補償,進而提高最終得到的自由曲面離軸成像系統的成像品質。
請參閱圖7,本發明實施例提供一種採用上述方法設計得到的自由曲面離軸三反成像系統100。其中,該自由曲面離軸三反成像系統100包括相鄰且間隔設置的一主鏡102、一次鏡104、以及一三鏡106,該次鏡104為光闌面。所述主鏡102、次鏡104和三鏡106的面形均為自由曲面。從光源出射的光線依次經過所述主鏡102、次鏡104、以及三鏡106的反射後,被一探測器108接收到並成像。為了描述方便,將所述主鏡102所處的空間定義一第一三維直角坐標系(X,Y,Z)、次鏡104所處的空間定義一第二三維直角坐標系(X',Y',Z')、三鏡106所處的空間定義一第三三維直角坐標系(X'',Y'',Z'')。
所述主鏡102的頂點為所述第一三維直角坐標系(X,Y,Z)的原點,通過主鏡102頂點的一條水準方向的直線為Z軸,向左為負向右為正,Y軸在圖8所示的平面內,垂直於Z軸向上為正向下為負,X軸垂直於YZ平面,垂直YZ平面向裡為正向外為負。
在所述第一三維直角坐標系(X,Y,Z)中,所述主鏡102的反射面為xy的多項式自由曲面,該xy多項式自由曲面的方程式可表達為:, 其中,z為曲面矢高,c為曲面曲率,k為二次曲面係數,Ai
是多項式中第i項的係數。由於所述自由曲面離軸三反成像系統100關於YZ平面對稱,因此,可以僅保留X的偶次項。
優選的,所述主鏡102的反射面一xy多項式自由曲面,所述xy多項式為x的偶次多項式,x的最高次數為8次,該xy多項式自由曲面的方程式可表達為:。
本實施例中,所述主鏡102反射面的xy多項式中曲率c、二次曲面係數k以及各項係數Ai
的值請參見表1。可以理解,曲率c、二次曲面係數k以及各項係數Ai
的值也不限於表1中所述,本領域技術人員可以根據實際需要調整。 表1主鏡的反射面的xy多項式中的各係數的值
所述第二三維直角坐標系(X',Y',Z')的原點為所述次鏡104的頂點。所述第二三維直角坐標系(X',Y',Z')為所述第一三維直角坐標系(X,Y,Z)沿Z軸負方向和Y軸負方向平移得到。本實施例中,所述第二三維直角坐標系(X',Y',Z')由所述第一三維直角坐標系(X,Y,Z)沿Y軸負方向平移約272.306mm,再沿Z軸反方向平移約518.025mm,然後以X軸為旋轉軸沿逆時針方向的旋轉角度約為31.253°得到。所述第一三維直角坐標系(X,Y,Z)的原點到所述第二三維直角坐標系(X',Y',Z')的原點的距離約為585.235mm。
在所述第二三維直角坐標系(X',Y',Z')中,所述次鏡104的反射面為x'y'的多項式自由曲面,該x'y'多項式自由曲面的方程式可以表達為:, 其中,z'為曲面矢高,c'為曲面曲率,k'為二次曲面係數,Ai
'是多項式中第i項的係數。由於所述自由曲面離軸三反成像系統100關於Y'Z'平面對稱,因此,可以僅保留X'的偶次項。
優選的,所述次鏡104的反射面為一x'y'多項式自由曲面,所述x'y'多項式為x'的偶次多項式,x'的最高次數為8次,該x'y'多項式自由曲面的方程式可表達為:。
本實施例中,所述次鏡104反射面的x'y'多項式中曲率c'、二次曲面係數k'以及各項係數Ai
'的值請參見表2。可以理解,曲率c'、二次曲面係數k'以及各項係數Ai
'的值也不限於表2中所述,本領域技術人員可以根據實際需要調整。 表2次鏡的反射面的x'y'多項式中的各係數的值
所述第三三維直角坐標系(X'',Y'',Z'')的原點為所述第三反射鏡的頂點。該第三三維直角坐標系(X'',Y'',Z'')為所述第二三維直角坐標系(X',Y',Z')沿Z'軸正方向和Y'軸正方向平移得到。本實施例中,所述第三三維直角坐標系(X'',Y'',Z'')為所述第二三維直角坐標系(X',Y',Z')沿Z'軸正方向平移約346.467mm,再沿Y' 軸正方向平移約141.540,然後繞X軸逆時針旋轉角度約為20.079°得到。所述第二三維直角坐標系(X',Y',Z')的原點到所述第三三維直角坐標系(X'',Y'',Z'')的原點的距離約為374.263mm。
在所述第三三維直角坐標系(X'',Y'',Z'')中,所述三鏡106的反射面為x''y''的多項式自由曲面,該x''y''多項式曲面的方程式可以表達為:, 其中,z''為曲面矢高,c''為曲面曲率,k''為二次曲面係數,Ai
''是多項式中第i項的係數。由於所述自由曲面離軸三反成像系統100關於Y''Z''平面對稱,因此,可以僅保留x''的偶次項。
優選的,所述三鏡106的反射面一x''y''多項式自由曲面,所述x''y''多項式為x''的偶次多項式,x''的最高次數為8次,該x''y''多項式自由曲面的方程式可表達為:,
本實施例中,所述三鏡106反射面面形的x''y''多項式中,曲率c''、二次曲面係數k''以及各項係數Ai
''的值請參見表3。可以理解,曲率c''、二次曲面係數k''以及各項係數Ai
''的值也不限於表3中所述,本領域技術人員可以根據實際需要調整。 表3三鏡的反射面的x''y''多項式中的各係數的值
所述主鏡102、次鏡104和三鏡106的材料不限,只要保證其具有較高的反射率即可。所述主鏡102、次鏡104和三鏡106可選用鋁、銅等金屬材料,也可選用碳化矽、二氧化矽等無機非金屬材料。為了進一步增加所述主鏡102、次鏡104和三鏡106的反射率,可在其各自的反射面鍍一增反膜,該增反膜可為一金膜。所述主鏡102、次鏡104和三鏡106的尺寸不限。
在第一三維直角坐標系(X,Y,Z)中,所述探測器108的中心與所述三鏡106的頂點沿Z軸反方向的距離約為391.895mm,所述探測器108的中心沿Y軸正方向偏離Z軸,偏離量約為33.491mm。所述探測器108與XY平面沿順時針方向的角度約為7.435°。所述探測器108的尺寸根據實際需要進行設計。本實施例中,所述探測器108的尺寸為5μm×5μm。
所述自由曲面離軸三反成像系統100的等效入瞳直徑為60毫米。
所述自由曲面離軸三反成像系統100在子午方向上採用了離軸視場。所述自由曲面離軸三反成像系統100的視場角大於等於60°×1°。本實施例中,所述自由曲面離軸三反成像系統100的視場角為60°×1°,其中,在弧矢方向的角度為-30°至30°,在子午方向的角度為8°至9°。
所述自由曲面離軸三反成像系統100的工作波長範圍為380奈米到780奈米。當然,所述自由曲面離軸三反成像系統100的工作波長並不限於本實施例,本領域技術人員可以根據實際需要調整。
所述自由曲面離軸三反成像系統100的等效焦距為150mm。
所述自由曲面離軸三反成像系統100的F數小於等於2.5。本實施例中,所述自由曲面離軸三反成像系統100的F數為2.5,相對孔徑大小D/f為所述F數的倒數,即,所述自由曲面離軸三反成像系統100的相對孔徑大小D/f為0.4。
請參閱圖8,為自由曲面離軸三反成像系統100在可見光波段下部分視場角的調製傳遞函數MTF,從圖中可以看出,各視場MTF曲線都基本達到了衍射極限,表明該自由曲面離軸三反成像系統100具有很高的成像品質。
請參閱圖9,為所述自由曲面離軸三反成像系統100的波像差,平均值為0.089λ,其中λ=530.5 nm。說明該自由曲面離軸三反成像系統100的成像品質很好。
請參閱圖10,為所述自由曲面離軸三反成像系統100的相對畸變圖。從圖10中可以看出,所述自由曲面離軸三反成像系統100的最大相對畸變為8.1%。相比於其他大視場角的自由曲面離軸三反成像系統,該自由曲面離軸三反成像系統100的畸變較小,說明自由曲面離軸三反成像系統100的成像誤差較小,品質較高。
本發明提供的自由曲面離軸三反成像系統100採用離軸三反系統,沒有中心遮攔,視場角可以達到60°×1°,具有很大的視場角,成像範圍較大;該自由曲面離軸三反成像系統100的F數只有2.5,能夠獲得高解析度圖像;該自由曲面離軸三反成像系統100的結構緊湊。
本發明提供的自由曲面離軸三反成像系統100應用領域涉及到對地觀測、空間目標探測、天文觀測、多光譜熱成像、立體測繪等。本發明提供的自由曲面離軸三反成像系統100在可見光波段達到了衍射極限,可以在可見光下進行使用,也可以在紅外波段進行使用。
綜上所述,本發明確已符合發明專利之要件,遂依法提出專利申請。惟,以上所述者僅為本發明之較佳實施例,自不能以此限制本案之申請專利範圍。舉凡習知本案技藝之人士援依本發明之精神所作之等效修飾或變化,皆應涵蓋於以下申請專利範圍內。
100‧‧‧自由曲面離軸三反成像系統
102‧‧‧主鏡
104‧‧‧次鏡
106‧‧‧三鏡
108‧‧‧探測器
圖1為本發明實施例提供的特徵視場的選取方法。
圖2為本發明實施例提供的構建特徵視場的拓展方向示意圖。
圖3為本發明提供的拓展構建特徵視場並構建自由曲面系統的示意圖。
圖4為本發明實施例提供的延展構建區域並構建自由曲面離軸成像系統的示意圖。
圖5為本發明實施例提供的自由曲面離軸成像系統的設計流程示意圖。
圖6為本發明實施例提供的自由曲面離軸三反成像系統的構建過程示意圖。
圖7為本發明實施例提供的自由曲面離軸三反成像系統的光路圖。
圖8為本發明實施例提供的自由曲面離軸三反成像系統的MTF曲線。
圖9為本發明實施例提供的自由曲面離軸三反成像系統的平均RMS彌散斑直徑圖。
圖10為本發明實施例提供的自由曲面離軸三反成像系統的畸變網格圖。
無
Claims (10)
- 一種自由曲面離軸三反成像系統,包括: 一主鏡,以所述主鏡的頂點為第一原點定義一第一三維直角坐標系(X,Y,Z),在該第一三維直角坐標系(X,Y,Z)中,所述主鏡的反射面為一xy多項式自由曲面; 一次鏡,以所述次鏡的頂點為第二原點定義一第二三維直角坐標系(X',Y',Z'),且該第二三維直角坐標系(X',Y',Z')由所述第一三維直角坐標系(X,Y,Z)沿Z軸負方向和Y軸負方向平移得到,在該第二三維直角坐標系(X',Y',Z')中,所述次鏡的反射面為一x'y'多項式自由曲面; 一三鏡,以所述三鏡的頂點為第三原點定義一第三三維直角坐標系(X'',Y'',Z''),且該第三三維直角坐標系(X'',Y'',Z'')由所述第二三維直角坐標系(X',Y',Z')沿Z'軸正方向和Y'軸正方向平移得到,在該第三三維直角坐標系(X'',Y'',Z'')中,所述三鏡的反射面為一x''y''多項式自由曲面;以及 一探測器,光線依次經過所述主鏡、次鏡和三鏡的反射後,被該探測器接收到並成像; 其中,該自由曲面離軸三反成像系統的視場角大於等於60°×1°,且F數小於等於2.5。
- 如請求項1所述的自由曲面離軸三反成像系統,其中,所述第二三維直角坐標系(X',Y',Z')由所述第一三維直角坐標系(X,Y,Z)沿Y軸負方向平移272.306mm,再沿Z軸反方向平移518.025mm,然後以X軸為旋轉軸沿逆時針方向旋轉31.253°得到。
- 如請求項1所述的自由曲面離軸三反成像系統,其中,所述第三三維直角坐標系(X'',Y'',Z'')為所述第二三維直角坐標系(X',Y',Z')沿Z'軸正方向平移346.467mm,再沿Y' 軸正方向平移141.540,然後繞X'軸逆時針旋轉20.079°得到。
- 如請求項1所述的自由曲面離軸三反成像系統,其中,所述第一原點到所述第二原點的距離為585.235mm,所述第二原點到所述第三原點的距離為374.263mm。
- 如請求項1所述的自由曲面離軸三反成像系統,其中,在該第一三維直角坐標系(X,Y,Z)中,所述主鏡的反射面為一xy多項式自由曲面,所述xy多項式為x的偶次多項式,x的最高次數為8次,該xy多項式自由曲面的方程式為:, 其中,c=2.0905292558E-03,k=-1.21892397257173,A2=-7.4132168689E-01,A3=-1.0489399277E-03,A5=-1.0932187535E-03,A7=1.8886774078E-07,A9=3.9367755612E-08,A10=1.8215989925E-10,A12=2.0074956485E-10,A14=3.4948173329E-09,A16=-1.5663370553E-13,A18=-2.3270934790E-13,A20=-4.1540232365E-11,A21=-7.5480347775E-17,A23=4.3918989072E-16,A25=-8.2899003141E-15,A27=1.7435762237E-13,A29=1.3339713122E-19,A31=3.0368239305E-18,A33=8.0061751667E-18,A35=-4.7950132667E-17,A36=5.1808677261E-23,A40=-1.2808456638E-20,以及A42=6.5955002284E-19。
- 如請求項1所述的自由曲面離軸三反成像系統,其中,在第二三維直角坐標系(X',Y',Z')中,所述次鏡的反射面為一x'y'多項式自由曲面,所述x'y'多項式為x'的偶次多項式,x'的最高次數為8次,該x'y'多項式自由曲面的方程式為:, 其中,c'=1.2419254309E-03,k'=11.4470143760146,A2'=2.6039347329E-01,A3'=-1.5571495307E-03,A5'=-1.8499719195E-03,A7'=3.1123110562E-07,A9'=-7.2750566448E-08,A10'=-1.0565839855E-08,A12'=-2.4835983560E-08,A14'=-1.6584532798E-08,A16'=-2.1393403767E-12,A18'=-1.6810005234E-12,A20'=-1.3947835440E-11,A21'=-3.1895103407E-14,A23'=-8.4530987667E-13,A25'=-4.1377307296E-13,A27'=8.5456614986E-14,A29'=9.6174103349E-15,A31'=-7.1427426063E-15,A33'=-3.1310574623E-15,A35'=2.8784215520E-15,A36'=6.8582114820E-17,A40'=-1.5699572247E-15,以及A42'=6.5300451569E-17。
- 如請求項1所述的自由曲面離軸三反成像系統,其中,第三三維直角坐標系(X'',Y'',Z'')中,所述三鏡的反射面為一x''y''多項式自由曲面,所述x''y''多項式為x''的偶次多項式,x''的最高次數為8次,該x''y''多項式自由曲面的方程式為:, 其中,c''=-2.3117990132E-03,k''=-9.0353982823E-02,A2''=2.3151773248E-02,A3''=1.3411119803E-04,A5''=9.7954301258E-05,A7''=4.0553954355E-08,A9''=2.1746460869E-08,A10''=1.6438448222E-10,A12''=2.2410934162E-10,A14''=4.5429699280E-11,A16''=1.1872420633E-13,A18''=1.5798764234E-13,A20''=1.3425352595E-14,A21''=3.1961386864E-16,A23''=8.2640583723E-16,A25''=5.0212511084E-17,A27''=1.5492827151E-16,A29''=7.1429537548E-19,A31''=1.1507858347E-18,A33''=-6.5780484704E-19,A35''=4.2158929949E-18,A36''=7.6045895800E-22,A40''=-1.0675040865E-22,以及A42''=2.0087308226E-20。
- 如請求項1所述的自由曲面離軸三反成像系統,其中,在第一三維直角坐標系(X,Y,Z)中,所述探測器的中心與所述三鏡的頂點沿Z軸負方向的距離為391.895mm;所述探測器的中心沿Y軸正方向偏離Z軸,偏離量為33.491mm;且所述探測器與XY平面沿順時針方向的角度為7.435°。
- 如請求項1所述的自由曲面離軸三反成像系統,其中,所述自由曲面離軸三反成像系統的視場角為60°×1°,其中,在弧矢方向的角度為-30°至30°,在子午方向的角度為8°至9°。
- 如請求項1所述的自由曲面離軸三反成像系統,其中,所述自由曲面離軸三反成像系統的等效入瞳直徑為60毫米。
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| ??201810033376.0 | 2018-01-12 | ||
| CN201810033376.0 | 2018-01-12 | ||
| CN201810033376.0A CN110031957B (zh) | 2018-01-12 | 2018-01-12 | 自由曲面离轴三反成像系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| TW201935071A true TW201935071A (zh) | 2019-09-01 |
| TWI724275B TWI724275B (zh) | 2021-04-11 |
Family
ID=67212621
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| TW107103160A TWI724275B (zh) | 2018-01-12 | 2018-01-30 | 自由曲面離軸三反成像系統 |
Country Status (3)
| Country | Link |
|---|---|
| US (1) | US11043521B2 (zh) |
| CN (1) | CN110031957B (zh) |
| TW (1) | TWI724275B (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| TWI734389B (zh) * | 2020-01-15 | 2021-07-21 | 鴻海精密工業股份有限公司 | 自由曲面離軸三反成像系統 |
| TWI754295B (zh) * | 2020-05-29 | 2022-02-01 | 鴻海精密工業股份有限公司 | 自由曲面離軸三反光學系統 |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112305738B (zh) * | 2019-08-01 | 2022-02-08 | 清华大学 | 自由曲面反射式红外成像系统 |
| CN112882188A (zh) * | 2019-11-29 | 2021-06-01 | 浙江舜宇光学有限公司 | 摄像装置 |
| CN110794557B (zh) * | 2019-11-29 | 2025-05-02 | 浙江舜宇光学有限公司 | 摄像装置 |
| CN110927943A (zh) * | 2019-12-26 | 2020-03-27 | 中国科学院长春光学精密机械与物理研究所 | 一种离轴三反式漫反射板照明系统 |
| CN111190273B (zh) * | 2020-02-28 | 2021-10-15 | 莆田学院 | 用于空间遥感相机的大视场紧凑型光学系统 |
| CN112180576B (zh) * | 2020-09-02 | 2021-09-14 | 中国科学院西安光学精密机械研究所 | 一种制冷型自由曲面离轴三反光学系统 |
| CN114879351A (zh) * | 2021-02-05 | 2022-08-09 | 清华大学 | 无对称自由曲面光学系统 |
| CN113031296B (zh) * | 2021-03-24 | 2022-08-19 | 长春长光智欧科技有限公司 | 可快速装调的金属基自由曲面三反光学系统的装调方法 |
| CN115793189B (zh) * | 2022-10-14 | 2024-08-30 | 中国科学院西安光学精密机械研究所 | 一种全自由曲面成像系统 |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2000171717A (ja) * | 1998-12-07 | 2000-06-23 | Olympus Optical Co Ltd | 結像光学系 |
| JP4292609B2 (ja) * | 1999-01-06 | 2009-07-08 | 株式会社ニコン | 軸外し反射光学系 |
| JP5049464B2 (ja) * | 2004-12-28 | 2012-10-17 | 富士フイルム株式会社 | 反射型光学系およびこれを用いた投写型表示装置 |
| US20120081800A1 (en) * | 2009-04-20 | 2012-04-05 | Dewen Cheng | Optical see-through free-form head-mounted display |
| DE102012022967A1 (de) * | 2012-11-21 | 2014-05-22 | Carl Zeiss Meditec Ag | Reflexfreies, optisches System für eine Funduskamera |
| CN104898261B (zh) * | 2014-03-05 | 2017-10-24 | 清华大学 | 自由曲面离轴三反光学系统 |
| CN105739073B (zh) * | 2014-12-11 | 2019-02-12 | 清华大学 | 自由曲面离轴三反光学系统 |
| CN105334607B (zh) * | 2015-12-05 | 2017-11-14 | 中国航空工业集团公司洛阳电光设备研究所 | 一种无遮拦纯反射光学系统 |
-
2018
- 2018-01-12 CN CN201810033376.0A patent/CN110031957B/zh active Active
- 2018-01-30 TW TW107103160A patent/TWI724275B/zh active
- 2018-12-23 US US16/231,535 patent/US11043521B2/en active Active
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| TWI734389B (zh) * | 2020-01-15 | 2021-07-21 | 鴻海精密工業股份有限公司 | 自由曲面離軸三反成像系統 |
| TWI754295B (zh) * | 2020-05-29 | 2022-02-01 | 鴻海精密工業股份有限公司 | 自由曲面離軸三反光學系統 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110031957A (zh) | 2019-07-19 |
| TWI724275B (zh) | 2021-04-11 |
| US20190221599A1 (en) | 2019-07-18 |
| CN110031957B (zh) | 2021-01-05 |
| US11043521B2 (en) | 2021-06-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| TWI724275B (zh) | 自由曲面離軸三反成像系統 | |
| TWI677728B (zh) | 自由曲面離軸成像系統的設計方法 | |
| CN108241207B (zh) | 用作斜视相机镜头的光学系统的设计方法 | |
| TWI622813B (zh) | 離軸非球面光學系統的設計方法 | |
| CN111487766B (zh) | 自由曲面离轴三反成像系统的设计方法 | |
| TWI557433B (zh) | 具有實出瞳的自由曲面離軸成像系統的設計方法 | |
| CN108241204B (zh) | 用作斜视相机镜头的光学系统 | |
| TWI622800B (zh) | 離軸非球面三反光學系統 | |
| CN110727092A (zh) | 一种基于自由曲面的离轴反射式两镜扩束系统 | |
| CN108345094B (zh) | 混合表面离轴三反光学系统 | |
| JP3933439B2 (ja) | 光学系 | |
| TWI639930B (zh) | 混合表面光學系統的設計方法 | |
| TWI748237B (zh) | 自由曲面反射式紅外成像系統的設計方法 | |
| TW202107146A (zh) | 自由曲面反射式紅外成像系統 | |
| TWI734389B (zh) | 自由曲面離軸三反成像系統 | |
| CN111487754A (zh) | 自由曲面离轴三反成像系统 |