CN112305737B - 自由曲面反射式红外成像系统的设计方法 - Google Patents
自由曲面反射式红外成像系统的设计方法 Download PDFInfo
- Publication number
- CN112305737B CN112305737B CN201910709139.6A CN201910709139A CN112305737B CN 112305737 B CN112305737 B CN 112305737B CN 201910709139 A CN201910709139 A CN 201910709139A CN 112305737 B CN112305737 B CN 112305737B
- Authority
- CN
- China
- Prior art keywords
- free
- field
- characteristic
- view
- point
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 58
- 238000003331 infrared imaging Methods 0.000 title claims abstract description 26
- 238000013461 design Methods 0.000 title claims abstract description 19
- 230000000007 visual effect Effects 0.000 claims abstract description 12
- 239000013598 vector Substances 0.000 claims description 40
- 238000010276 construction Methods 0.000 claims description 14
- 238000004364 calculation method Methods 0.000 claims description 8
- FGUUSXIOTUKUDN-IBGZPJMESA-N C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 Chemical compound C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 FGUUSXIOTUKUDN-IBGZPJMESA-N 0.000 claims description 2
- 238000005070 sampling Methods 0.000 claims description 2
- 238000003384 imaging method Methods 0.000 description 7
- 230000003287 optical effect Effects 0.000 description 5
- 238000005457 optimization Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000010187 selection method Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Images
Classifications
-
- 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
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N25/00—Circuitry of solid-state image sensors [SSIS]; Control thereof
- H04N25/40—Extracting pixel data from image sensors by controlling scanning circuits, e.g. by modifying the number of pixels sampled or to be sampled
- H04N25/41—Extracting pixel data from a plurality of image sensors simultaneously picking up an image, e.g. for increasing the field of view by combining the outputs of a plurality of sensors
-
- 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
-
- G—PHYSICS
- G02—OPTICS
- G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
- G02B27/00—Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00
- G02B27/0012—Optical design, e.g. procedures, algorithms, optimisation routines
-
- G—PHYSICS
- G02—OPTICS
- G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
- G02B27/00—Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00
- G02B27/09—Beam shaping, e.g. changing the cross-sectional area, not otherwise provided for
- G02B27/0938—Using specific optical elements
- G02B27/0977—Reflective elements
-
- G—PHYSICS
- G02—OPTICS
- G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
- G02B27/00—Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00
- G02B27/0025—Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00 for optical correction, e.g. distorsion, aberration
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Optics & Photonics (AREA)
- Engineering & Computer Science (AREA)
- Multimedia (AREA)
- Signal Processing (AREA)
- Image Analysis (AREA)
- Lenses (AREA)
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
- Transforming Light Signals Into Electric Signals (AREA)
Abstract
本发明涉及一种自由曲面反射式红外成像系统的设计方法,其具体包括:S1,选取已有且可扩展空间充足的小F数小视场系统作为初始系统,该初始系统的视场为X0×Y0;S2,根据待设计的系统的视场Xn×Yn和初始系统的视场X0×Y0,确定视场序列[X0,Y0],[X1,Y1],[X2,Y2],…,[Xn,Yn],其中X0<X1<X2<…<Xn,Y0<Y1<Y2<…<Yn;S3,根据视场X1×Y1,采用逐点求解的方法构建系统中的所有自由曲面;S4,以上一步得到的系统再次作为初始系统,根据视场X2×Y2,采用逐点求解的方法构建系统中所有的自由曲面;以及S5,以此类推,按照视场序列的顺序不断进行构建,直到达到最终视场Xn×Yn,从而得到待设计的自由曲面反射式红外成像系统。
Description
技术领域
本发明涉及光学设计领域,尤其涉及一种自由曲面反射式红外成像系统的设计方法。
背景技术
自由曲面是指无法用球面或非球面系数来表示的非传统曲面,通常是非回转对称的,结构灵活,变量较多,为光学设计提供了更多的自由度,可以大大降低光学系统的像差,减小系统的体积、重量与镜片数量,可以满足现代成像系统的需要,有着广阔的发展应用前景。由于自由曲面有非对称面并提供了更多的设计自由度,他们常被用在离轴非对称系统中。
现有的自由曲面离轴成像系统普遍采用直接设计方法得到,如微分方程法,多曲面同步设计方法和逐点构建与迭代方法,然而,现有的设计方法在设计过程中只能考虑有限数量的视场,在设计具有大相对孔径和大矩形视场角的自由曲面离轴成像系统时具有局限性。
发明内容
综上所述,确有必要提供一种自由曲面反射式红外成像系统的设计方法,采用该设计方法,可以得到同时具有小F数和大矩形视场的自由曲面反射式红外成像系统。
一种自由曲面反射式红外成像系统的设计方法,其具体包括以下步骤:
S1,选取已有且可扩展空间充足的小F数小视场系统作为初始系统,该初始系统的视场为X0×Y0;
S2,根据待设计的系统的视场Xn×Yn和初始系统的视场X0×Y0,确定视场序列[X0,Y0],[X1,Y1],[X2,Y2],…,[Xn,Yn],其中X0<X1<X2<…<Xn,Y0<Y1<Y2<…<Yn;
S3,根据视场X1×Y1,采用逐点求解的方法构建系统中的所有自由曲面;
S4,以上一步得到的系统再次作为初始系统,根据视场X2×Y2,采用逐点求解的方法构建系统中所有的自由曲面;以及
S5,以此类推,按照视场序列的顺序不断进行构建,直到达到最终视场Xn×Yn,从而得到待设计的自由曲面反射式红外成像系统。
相较于现有技术,本发明提供的自由曲面反射式红外成像系统的设计方法利用已有的小F数小矩形视场系统作为初始系统,等步长或不等步长地逐渐扩展视场,同时采用逐点求解的方法对系统进行构建,经过多次构建后得到同时具有小F数和大矩形视场的自由曲面反射式红外成像系统。
附图说明
图1为本发明实施例提供的自由曲面反射式红外成像系统的设计流程示意图。
图2为本发明实施例提供的初始系统的结构示意图。
图3为本发明实施例提供的在第一视场中对系统进行构建得到的结果。
图4为本发明实施例提供的在第二视场中对系统进行构建得到的结果。
图5为本发明实施例提供的在第四视场中对系统进行构建得到的结果。
图6为本发明实施例提供的在第六视场中对系统进行构建得到的结果。
图7为本发明实施例提供的在第九视场中对系统进行构建得到的结果。
主要元件符号说明
初始自由曲面离轴三反系统 100
主镜 102
次镜 104
三镜 106
红外光探测器 108
如下具体实施方式将结合上述附图进一步说明本发明。
具体实施方式
下面根据说明书附图并结合具体实施例对本发明的技术方案进一步详细表述。
本发明实施例提供一种自由曲面反射式红外成像系统的设计方法,参照图1所示,其具体包括以下步骤。
步骤一,选取已有且可扩展空间充足的小F数小视场系统作为初始系统,该初始系统的视场为X0×Y0。
首先选定一个具有小F数但视场不满足要求的自由曲面系统作为初始系统。该初始系统的视场为X0×Y0。该初始系统包括多个初始自由曲面,将初始系统中初始自由曲面定义为Lj(j=1,2,……L),即初始自由曲面的个数为L个。所述多个初始自由曲面的具体位置根据待设计的自由曲面反射式红外成像系统的实际需要进行选择。所述初始系统中初始自由曲面的数量根据实际需要进行设计。本实施例中,参照图2所示,所述初始系统为一初始自由曲面离轴三反系统100,所述初始自由曲面离轴三反系统100包括一初始主镜102、一初始次镜104、一初始三镜106、一红外光探测器108。所述初始主镜102、所述初始次镜104以及所述初始三镜106的反射面均为初始自由曲面,所述初始次镜104为一孔径光阑。该初始系统的F数为1.39,视场大小为带有偏置的8°×6°,中心视场为(0,-22°),经过缩放使得该系统的焦距为9.3毫米(mm)。
步骤二,根据待设计的系统的视场Xn×Yn和初始系统的视场X0×Y0,确定视场序列[X0,Y0],[X1,Y1],[X2,Y2],…,[Xn,Yn],其中X0<X1<X2<…<Xn,Y0<Y1<Y2<…<Yn。
通过对初始系统不断构建来满足成像要求,假设总共需要n步构建,最终待设计的系统的视场大小定义为Xn×Yn。则可以选择视场序列为[X0,Y0],[X1,Y1],[X2,Y2],…,[Xn,Yn],分别定义为初始视场、第一视场、第二视场、…、第n视场,其中X0<X1<X2<…<Xn,Y0<Y1<Y2<…<Yn。视场序列的选取方式可以按照等步长的方式选取,即X0、X1、X2、…、Xn中相邻数值之间的差值相等,Y0、Y1、Y2、…、Yn中相邻数值之间的差值相等;也可以按照不等步长的方式选取,即在最初的几次构建过程选用大步长,当接近最终视场时选用小步长。在视场较小时,系统构建过程稳定性高,选用大步长可以减少构建次数,节省算法时间;当接近最终视场时,系统性能较高,系统构建过程稳定性降低,选用小步长增加视场可以增加构建过程的稳定性。本实施例中,设定待设计的视场大小为40°×30°,需要9步构建,在最初的几次构建过程选用大步长,当接近最终视场时选用小步长,从而确定视场序列为初始视场[8°,6°],第一视场[16°,12°],第二视场[20°,15°],第三视场[24°,18°],第四视场[28°,21°],第五视场[32°,24°],第六视场[34°,25.5°],第七视场[36°,27°],第八视场[38°,28.5°]以及第九视场[40°,30°]。
步骤三,根据视场X1×Y1,采用逐点求解的方法构建系统中的所有自由曲面。
以步骤一选取的小F数小视场系统作为初始系统,依据根据视场X1×Y1,采用逐点求解的方法将初始自由曲面Lj(j=i=1,2…L)依次构建成自由曲面Nj(j=i=1,2…L)。所述自由曲面Nj的求解顺序不限,可以根据实际需要进行调换。本发明以自由曲面N1的构建过程为例详细说明自由曲面的构建。将初始自由曲面L1构建成自由曲面N1的过程包括:
S31,选取K条特征光线Ri(i=1,2…K);
S32,根据物像关系及斯涅尔定律逐点求解所述K条特征光线与自由曲面N1上的多个交点,进而计算得到所述自由曲面N1上的所有特征数据点Pi(i=1,2…K);以及
S33,将该特征数据点进行曲面拟合得到所述自由曲面N1的方程式。
在步骤S31中,所述K条特征光线的选取可以根据以下方法进行。
根据视场X1×Y1,选取M个特征视场:
特征视场的选取方法不限。离轴反射成像系统往往关于yOz平面对称,如果视场是偏置的,在(0,y0)为中心视场的情况下,可以按照如下方式选取特征视场:
本实施例中,中心视场为(0,-22°),因此,对于16°×12°的视场大小,选取特征视场为(0,-22°)、(0,-16°)、(0,-28°)、(8,-22°)、(8,-16°)、(8,-28°)。
在所述M个特征视场中选取K条特征光线:
将M个特征视场中的每个特征视场的孔径分成N等份,并从每一等份中选取不同孔径位置的P条特征光线,这样一共选取了K=M×N×P条对应不同视场不同孔径位置的特征光线。所述孔径可以为圆形、长方形、正方形、椭圆形或其他规则或不规则的形状。优选的,所述视场孔径为圆形孔径,将每个视场的圆形孔径等分成N个角度,间隔为因此有沿着每个角度的半径方向取P个不同的孔径位置,那么一共取K=M×N×P条对应不同视场不同孔径位置的特征光线。本实施例中,将每个特征视场的孔径分成14等份,并从每一等份中选取不同孔径位置的7条特征光线,这样每个特征视场采样98条特征光线,6个特征视场总共采样588条特征光线。
在步骤S32中,为了得到所述自由曲面N1上的所有特征数据点Pi(i=1,2…K),将借助特征光线Ri(i=1,2…K)与待求的自由曲面N1的前一个曲面N1'及后一个曲面N1″的交点。在求解每条特征光线Ri(i=1,2…K)对应的待求的自由曲面N1上的特征数据点Pi(i=1,2…K)时,将该特征光线Ri与前一个曲面N1'的交点定义为该特征光线的起点Si,特征光线Ri与后一个曲面N1″的交点定义为该特征光线的终点Ei。当待设计的系统与特征光线确定后,该特征光线Ri的起点Si是确定的,且易于通过光线追迹即物像关系得到,特征光线的终点Ei可通过物像关系求解。在理想状态下,特征光线Ri从N1'上的Si射出后,经过N1上的Pi,交于N1″上的Ei,并最终交目标面于其理想目标点Ti,ideal。如果N1″是目标面,特征光线的终点就是其目标点Ti;如果在待求自由曲面N1和目标面之间还有其他面,特征光线的终点是N1″上从特征数据点到其对应的目标点的光程的一阶变分为零的点,即
所述自由曲面N1上特征数据点Pi(i=1,2…K)可以通过以下两种计算方法获得。
第一种计算方法包括以下步骤:
步骤a,取定一条特征光线R1与所述自由曲面面N1对应的初始曲面的第一交点为特征数据点P1;
步骤c,过所述i(1≤i≤K-1)个特征数据点Pi分别做一第一切平面,得到i个第一切平面,该i个第一切平面与其余K-i条特征光线相交得到i×(K-i)个第二交点,从该i×(K-i)个第二交点中选取出与所述i(1≤i≤K-1)个特征数据点Pi距离最近的一个第二交点,作为所述自由曲面N1的下一个特征数据点Pi+1;
步骤d,重复步骤b和c,直到计算得到自由曲面N1上的所有特征数据点Pi(i=1,2…K),通过曲面拟合可以得到所述自由曲面N1的方程式。
第二种计算方法包括以下步骤:
步骤a',取定一第一条特征光线R1与所述自由曲面N1对应的初始曲面的第一交点为特征数据点P1;
步骤c',仅过所述第i(1≤i≤K-1)个特征数据点Pi做一第一切平面并与其余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对应的的方程式。
在步骤S33中,将所述自由曲面N1上的多个特征数据点Pi(i=1,2…K)进行曲面拟合包括以下步骤:
S331:将所述多个特征数据点Pi(i=1,2…K)在所述第一直角坐标系中拟合成一球面,得到所述球面的曲率c及其对应的曲率中心(xc,yc,zc);
S332:将中心采样视场主光线对应的特征数据点(xo,yo,zo)定义为球面的顶点,并以该球面的顶点为原点,过曲率中心与球面顶点的直线为z轴,建立一第二三维直角坐标系;
S333:将所述多个特征数据点Pi(i=1,2…K)在第一三维直角坐标系中的坐标(xi,yi,zi)及其对应的法向量(αi,βi,γi)变换为第二三维直角坐标系中的坐标(x'i,y'i,z'i)及其法向量(α'i,β'i,γ'i);
S334:根据所述多个特征数据点Pi(i=1,2…K)在第二三维直角坐标系中的坐标(x'i,y'i,z'i),以及步骤S331中求得的球面的曲率c,将特征数据点Pi(i=1,2…K)在第二三维直角坐标系中拟合成一个二次曲面,得到二次曲面系数k;以及
S335:将所述多个特征数据点Pi(i=1,2…K)在第二三维直角坐标系中的二次曲面上的坐标与法向量分别从坐标(x'i,y'i,z'i)和法向量(α'i,β'i,γ'i)中排除掉,得到残余坐标与残余法向量,将该残余坐标与残余法向量拟合成一个自由曲面,该自由曲面的方程式与所述二次曲面的方程式相加即可得到所述自由曲面N1的方程式。
通常的,光学系统关于三维直角坐标系的yoz面对称,因此,在步骤S333中,所述球面在第二三维直角坐标系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)分别为:
(xi″,yi″,zi″)=(xi′,yi′,zi′-zis′)和
在步骤S335中,所述将残余坐标与残余法向量拟合得到自由曲面的步骤包括:
S3351:在所述第二三维直角坐标系中,用排除掉二次曲面项的多项式曲面作为待构建自由曲面的面形描述,即
其中,gj(x,y)为多项式的某一项,P=(p1,p2,,pJ)T为系数集合;
S3352:获得残余坐标拟合误差d1(P),即所述残余坐标值(x″i,y″i,z″i)(i=1,2,…,n)与所述自由曲面在z'轴方向残余坐标差值的平方和,以及残余法向矢量拟合误差d2(P),即所述残余法向量Ni=(α″i,β″i,-1)(i=1,2,…,n)与所述自由曲面法向量的矢量差的模值的平方和,
其中,Z=(z1,z2,…,zI)T,U=(u1,u2,…,uI)T,V=(v1,v2,…,vI)T,
S3354:选择不同的权重w,并令所述评价函数f(P)的梯度▽f(P)=0,从而获得多组不同的P及其对应的多个自由曲面面形z=f(x,y;P);以及
S3355:获得具有最佳的成像质量的最终自由曲面面形。
可以理解,初始系统中的其他初始曲面构建成自由曲面的方法与构建自由曲面N1的方法相同。依据上述构建自由曲面的方法,依次构建每个曲面,从而得到所有的自由曲面。本实施例中,根据第一视场[16°,12°]对所述初始系统进行构建的结果如图3所示。
步骤四,以上一步得到的系统再次作为初始系统,根据视场X2×Y2,采用逐点求解的方法构建系统中所有的自由曲面。
以步骤三构建的系统作为初始系统,在第二视场X2×Y2中,再次构建系统中所有的自由曲面,其具体包括:选取第二视场X2×Y2中的K条特征光线;根据物像关系及斯涅尔定律逐点求解所述K条特征光线与待求解的自由曲面的多个交点,该多个交点作为待求解的自由曲面的特征数据点;将该特征数据点进行曲面拟合得到待求解的自由曲面的方程式。其中,选取特征光线、求解特征数据点以及进行曲面拟合的方法均与步骤三所提供的方法相同,在此不作详细叙述。本实施例中,以根据第一视场[16°,12°]构建的自由曲面系统再次作为初始系统,根据第二视场[20°,15°]对系统进行构建,结果如图4所示。
步骤五,以此类推,按照视场序列的顺序不断进行构建,直到达到最终视场Xn×Yn,从而得到待设计的自由曲面反射式红外成像系统。
以根据第二视场X2×Y2构建的自由曲面系统再次作为初始系统,在第三视场X3×Y3中,构建系统中所有的自由曲面。然后以根据第三视场X3×Y3构建的自由曲面系统再次作为初始系统,在第四视场X4×Y4中,构建系统中所有的自由曲面。以此类推,按照步骤二选择的视场序列逐步增加视场的大小,同时以上一步得到的系统作为初始系统进行系统构建,直到达到最终待设计的视场Xn×Yn。最后,在第n视场Xn×Yn中,以根第n-1视场Xn-1×Yn-1构建的自由曲面系统再次作为初始系统,构建系统中所有的自由曲面,从而得到一小F数大矩形视场的自由曲面反射式红外成像系统。本实施例中,在第四视场[28°,21°]和第六视场[34°,25.5°]中对系统进行构建,得到的结果分别如图5和图6所示。最终在第九视场[40°,30°]中进行构建得到的系统如图7所示,该系统的F数为1.39,视场达到40°×30°,该系统满足性能需求且已具备成像性能,可以作为初始结构用于后续优化。
所述自由曲面反射式红外成像系统的设计方法可进一步包括一对步骤五中得到的自由曲面反射式红外成像系统进行优化的步骤。具体地,将步骤五中得到的自由曲面反射式红外成像系统作为后续优化的初始结构。可以理解,该优化的步骤并不是必需的,可以根据实际需要设计。
本发明提供的自由曲面反射式红外成像系统的设计方法利用已有的小F数小矩形视场系统作为初始系统,等步长或不等步长地逐渐扩展视场,同时采用逐点求解的方法对系统进行构建,经过多次构建后得到同时具有小F数和大矩形视场的自由曲面反射式红外成像系统。
另外,本领域技术人员还可在本发明精神内作其它变化,当然这些依据本发明精神所作的变化,都应包含在本发明所要求保护的范围内。
Claims (9)
1.一种自由曲面反射式红外成像系统的设计方法,其具体包括以下步骤:
S1,选取已有且可扩展空间充足的小F数小视场系统作为初始系统,该初始系统的视场为X0×Y0;
S2,根据待设计的系统的视场Xn×Yn和初始系统的视场X0×Y0,确定视场序列[X0,Y0],[X1,Y1],[X2,Y2],…,[Xn,Yn],其中X0<X1<X2<…<Xn,Y0<Y1<Y2<…<Yn;
S3,根据视场X1×Y1,采用逐点求解的方法构建系统中的所有自由曲面;采用逐点求解的方法构建系统中的所有自由曲面为采用逐点求解的方法将初始系统中的初始曲面Lj依次构建成自由曲面Nj,其中,j=1,2…L,将初始曲面L1构建成自由曲面N1的方法包括:
S31,选取K条特征光线Ri,其中,i=1,2…K;
S32,根据物像关系及斯涅尔定律逐点求解所述K条特征光线与自由曲面N1上的多个交点,进而计算得到所述自由曲面N1上的所有特征数据点Pi,其中,i=1,2…K;以及
S33,将该特征数据点进行曲面拟合得到所述自由曲面N1的方程式;
依据构建自由曲面N1的方法,依次构建每个自由曲面,从而得到所有的自由曲面;
S4,以上一步得到的系统再次作为初始系统,根据视场X2×Y2,采用逐点求解的方法构建系统中所有的自由曲面;以及
S5,以此类推,按照视场序列的顺序不断进行构建,直到达到最终视场Xn×Yn,从而得到待设计的自由曲面反射式红外成像系统。
2.如权利要求1所述的自由曲面反射式红外成像系统的设计方法,其特征在于,在步骤S1中,所述初始系统为一初始自由曲面系统,该初始自由曲面系统包括多个初始自由曲面。
3.如权利要求1所述的自由曲面反射式红外成像系统的设计方法,其特征在于,在步骤S2中,视场序列的选取方式按照等步长的方式选取,X0、X1、X2、…、Xn中相邻数值之间的差值相等,Y0、Y1、Y2、…、Yn中相邻数值之间的差值相等。
4.如权利要求1所述的自由曲面反射式红外成像系统的设计方法,其特征在于,在步骤S2中,视场序列的选取方式按照不等步长的方式选取,在最初的几次构建过程选用大步长,当接近最终视场时选用小步长。
5.如权利要求1所述的自由曲面反射式红外成像系统的设计方法,其特征在于,在步骤S31中,选取K条特征光线Ri的方式为:根据视场X1×Y1,选取M个特征视场,并在所述M个特征视场中选取K条特征光线。
7.如权利要求5所述的自由曲面反射式红外成像系统的设计方法,其特征在于,所述K条特征光线的选取方式为:将所述M个特征视场中的每个特征视场的孔径分成N等份,并从每一等份中选取不同孔径位置的P条特征光线,一共选取K=M×N×P条对应不同视场不同孔径位置的特征光线。
8.如权利要求1所述的自由曲面反射式红外成像系统的设计方法,其特征在于,所述自由曲面N1上特征数据点Pi的计算方法包括以下步骤:
步骤a,取定一第一条特征光线R1与所述自由曲面N1对应的初始曲面的第一交点为特征数据点P1;
步骤c,仅过所述第i个特征数据点Pi做一第一切平面并与其余K-i条特征光线相交,得到K-i个第二交点,其中,1≤i≤K-1;从该K-i个第二交点中选取出与所述第i个特征数据点Pi距离最短的第二交点Qi+1,并将其对应的特征光线及与所述第i个特征数据点Pi的最短距离分别定义为Ri+1和Di;
步骤d,过特征数据点Pi之前已求得的i-1个第一特征数据点分别做一第二切平面,得到i-1个第二切平面,其中,1≤i≤K-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对应的方程式。
9.如权利要求1所述的自由曲面反射式红外成像系统的设计方法,其特征在于,在步骤S33中,将该特征数据点进行曲面拟合得到所述自由曲面N1的方程式包括以下步骤:
将所述多个特征数据点Pi在一第一三维直角坐标系中拟合成一球面,得到所述球面的曲率c及其对应的曲率中心(xc,yc,zc),其中,i=1,2…K;
将中心采样视场主光线对应的特征数据点(xo,yo,zo)定义为球面的顶点,并以该球面的顶点为原点,过曲率中心与球面顶点的直线为z轴,建立一第二三维直角坐标系;
将所述多个特征数据点Pi在第一三维直角坐标系中的坐标(xi,yi,zi)及其对应的法向量(αi,βi,γi)变换为第二三维直角坐标系中的坐标(x'i,y'i,z'i)及其法向量(α'i,β'i,γ'i),其中,i=1,2…K;
根据所述多个特征数据点Pi在第二三维直角坐标系中的坐标(x'i,y'i,z'i),以及上述求得的球面的曲率c,将特征数据点Pi(i=1,2…K)在第二三维直角坐标系中拟合成一个二次曲面,得到二次曲面系数k,其中,i=1,2…K;以及
将所述多个特征数据点Pi在第二三维直角坐标系中的二次曲面上的坐标与法向量分别从坐标(x'i,y'i,z'i)和法向量(α'i,β'i,γ'i)中排除掉,得到残余坐标与残余法向量,将该残余坐标与残余法向量拟合成一个自由曲面,该自由曲面的方程式与所述二次曲面的方程式相加即可得到所述自由曲面N1的方程式,其中,i=1,2…K。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910709139.6A CN112305737B (zh) | 2019-08-01 | 2019-08-01 | 自由曲面反射式红外成像系统的设计方法 |
TW108130927A TWI748237B (zh) | 2019-08-01 | 2019-08-28 | 自由曲面反射式紅外成像系統的設計方法 |
US16/910,370 US11025841B2 (en) | 2019-08-01 | 2020-06-24 | Method for designing a freeform surface reflective imaging system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910709139.6A CN112305737B (zh) | 2019-08-01 | 2019-08-01 | 自由曲面反射式红外成像系统的设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112305737A CN112305737A (zh) | 2021-02-02 |
CN112305737B true CN112305737B (zh) | 2022-03-18 |
Family
ID=74259510
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910709139.6A Active CN112305737B (zh) | 2019-08-01 | 2019-08-01 | 自由曲面反射式红外成像系统的设计方法 |
Country Status (3)
Country | Link |
---|---|
US (1) | US11025841B2 (zh) |
CN (1) | CN112305737B (zh) |
TW (1) | TWI748237B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113281903B (zh) * | 2021-06-15 | 2022-03-08 | 中国科学院长春光学精密机械与物理研究所 | 简单高效的光线与光学自由曲面交点位置迭代求解方法 |
CN113281902B (zh) * | 2021-06-15 | 2022-03-08 | 中国科学院长春光学精密机械与物理研究所 | 基于割线法的光线与光学自由曲面交点位置迭代求解方法 |
CN115079391B (zh) * | 2022-06-10 | 2024-01-30 | 清华大学 | 离轴反射系统和小f数多反离轴系统的设计方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105988213A (zh) * | 2015-02-05 | 2016-10-05 | 清华大学 | 自由曲面离轴光学系统的设计方法 |
CN107219626A (zh) * | 2017-07-01 | 2017-09-29 | 南京理工大学 | 结合面形和视场优化策略的自由曲面光学系统优化方法 |
CN108227184A (zh) * | 2016-12-22 | 2018-06-29 | 清华大学 | 自由曲面成像光学系统的设计方法 |
CN110031969A (zh) * | 2018-01-12 | 2019-07-19 | 清华大学 | 自由曲面离轴成像系统的设计方法 |
CN111487766A (zh) * | 2019-01-25 | 2020-08-04 | 清华大学 | 自由曲面离轴三反成像系统的设计方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104570340B (zh) * | 2013-10-24 | 2017-04-05 | 清华大学 | 自由曲面成像系统的设计方法 |
CN105739089B (zh) * | 2014-12-11 | 2018-04-03 | 清华大学 | 自由曲面离轴三反成像系统的设计方法 |
CN105988214B (zh) * | 2015-02-05 | 2018-07-17 | 清华大学 | 具有实出瞳的自由曲面离轴成像系统的设计方法 |
-
2019
- 2019-08-01 CN CN201910709139.6A patent/CN112305737B/zh active Active
- 2019-08-28 TW TW108130927A patent/TWI748237B/zh active
-
2020
- 2020-06-24 US US16/910,370 patent/US11025841B2/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105988213A (zh) * | 2015-02-05 | 2016-10-05 | 清华大学 | 自由曲面离轴光学系统的设计方法 |
CN108227184A (zh) * | 2016-12-22 | 2018-06-29 | 清华大学 | 自由曲面成像光学系统的设计方法 |
CN107219626A (zh) * | 2017-07-01 | 2017-09-29 | 南京理工大学 | 结合面形和视场优化策略的自由曲面光学系统优化方法 |
CN110031969A (zh) * | 2018-01-12 | 2019-07-19 | 清华大学 | 自由曲面离轴成像系统的设计方法 |
CN111487766A (zh) * | 2019-01-25 | 2020-08-04 | 清华大学 | 自由曲面离轴三反成像系统的设计方法 |
Also Published As
Publication number | Publication date |
---|---|
TW202113425A (zh) | 2021-04-01 |
TWI748237B (zh) | 2021-12-01 |
US20210037199A1 (en) | 2021-02-04 |
US11025841B2 (en) | 2021-06-01 |
CN112305737A (zh) | 2021-02-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110031957B (zh) | 自由曲面离轴三反成像系统 | |
CN108152948B (zh) | 离轴非球面光学系统的设计方法 | |
TWI557433B (zh) | 具有實出瞳的自由曲面離軸成像系統的設計方法 | |
TWI557436B (zh) | 自由曲面離軸光學系統的設計方法 | |
CN110031969B (zh) | 自由曲面离轴成像系统的设计方法 | |
CN112305737B (zh) | 自由曲面反射式红外成像系统的设计方法 | |
TWI544230B (zh) | 自由曲面離軸三反成像系統的設計方法 | |
CN111487766B (zh) | 自由曲面离轴三反成像系统的设计方法 | |
CN105988212B (zh) | 三维自由曲面的设计方法 | |
CN108152939B (zh) | 离轴非球面三反光学系统 | |
TWI628463B (zh) | 混合表面離軸三反光學系統 | |
TW201534974A (zh) | 應用於光學系統的自由曲面的設計方法 | |
CN108345106B (zh) | 混合表面光学系统的设计方法 | |
CN113126270B (zh) | 自由曲面离轴三反成像系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |