CN114019677A - 无限制像方场分布的表征方法 - Google Patents

无限制像方场分布的表征方法 Download PDF

Info

Publication number
CN114019677A
CN114019677A CN202111371574.6A CN202111371574A CN114019677A CN 114019677 A CN114019677 A CN 114019677A CN 202111371574 A CN202111371574 A CN 202111371574A CN 114019677 A CN114019677 A CN 114019677A
Authority
CN
China
Prior art keywords
plane
field distribution
lens
calculating
aperture diaphragm
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.)
Granted
Application number
CN202111371574.6A
Other languages
English (en)
Other versions
CN114019677B (zh
Inventor
陈涛
金思宇
余毅
刘�英
陈小林
王成龙
张晓宁
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Changchun Institute of Optics Fine Mechanics and Physics of CAS
Original Assignee
Changchun Institute of Optics Fine Mechanics and Physics of CAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Changchun Institute of Optics Fine Mechanics and Physics of CAS filed Critical Changchun Institute of Optics Fine Mechanics and Physics of CAS
Priority to CN202111371574.6A priority Critical patent/CN114019677B/zh
Publication of CN114019677A publication Critical patent/CN114019677A/zh
Application granted granted Critical
Publication of CN114019677B publication Critical patent/CN114019677B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Optics & Photonics (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明提供无限制像方场分布的表征方法，用以计算包括孔径光阑和傅里叶透镜的光学***的光强分布，光学***的光线经过孔径光阑，通过傅里叶透镜达到焦平面和离焦面，计算方法包括如下步骤：计算孔径光阑面平面场分布的复振幅透过率；计算待测离焦面处的离焦因子；频域内计算复振幅透过率与离焦因子的傅里叶变换；计算待测离焦面的光强分布。在给定***焦距、工作波长条件下，通过本发明中采用的像方弥散分布表征方法，即可求出任意形状、任意离轴量下的孔径光阑在任意离焦面的光强分布，从而分析任意孔径光阑，在任意离焦面下的强度变化情况。

Description

无限制像方场分布的表征方法
技术领域
本发明涉及光学原理领域，尤其涉及无限制像方场分布的表征方法。
背景技术
准确计算像方空间光强分布对于光学设计的意义十分重大，通常孔径光阑成像遵守的是菲涅尔衍射理论。而传统的菲涅尔衍射理论仅仅适用于轴对称孔径光阑成像光学***，例如圆形孔径、矩形孔径等，而且仅能计算焦平面处的能量分布，难以满足当前科技发展对普适性像方空间能量分布的计算，如自动对焦***、超分辨***、景深扩展***等。
发明内容
本发明为了解决上述问题，提出了一种新的基于标量衍射理论的任意孔径形状和任意离焦量下离焦像面光场强度分布的数值计算方法。为实现上述目的，本发明采用以下具体技术方案：
无限制像方场分布的表征方法，用以计算光学***的光强分布，光学***包括孔径光阑和傅里叶透镜，光学***中发射的光线经过孔径光阑，通过傅里叶透镜到达焦平面和离焦面，计算方法包括步骤：
S1、计算孔径光阑面平面场分布的复振幅透过率t(x,y)；
S2、计算待测离焦面处的离焦因子s(x,y)，离焦因子s(x,y)用下式表达为：
Figure BDA0003362525630000011
其中，f为焦距，Δz为离焦量，k为波矢，x为，y为；
S3、频域内计算复振幅透过率t(x,y)与离焦因子s(x,y)的傅里叶变换f(fx,fy)，傅里叶变换f(fx,fy)用下式表达为：
f(fx,fy)＝F{t(x,y)×s(x,y)} (2)；
S4、计算待测离焦面的光强分布Iz(xz,yz)，光强分布Iz(xz,yz)用下式表达为：
Figure BDA0003362525630000021
其中，A为入射光波振幅，λ为入射波波长。
进一步的，在步骤S1之前还包括步骤S0，步骤S0包括：
S0、获取光学***的光学参数值：
光学***包括垂直于光轴的第一平面和第二平面，第一平面经过傅里叶透镜的入射面与光轴的交点，第二平面经过傅里叶透镜的出射面与光轴的交点；
光学参数值至少包括：入射光波振幅A，入射波波长λ，焦距f，离焦量Δz，光阑面平面场分布U0(x0,y0)，复振幅透过率t(x0,y0)，孔径光阑与傅里叶透镜的距离d0，第一平面的平面场分布Ul(x,y)，第二表面的平面场分布Ul′(x,y)，待测离焦面的场分布Uz(xz,yz)。
进一步的，根据角谱理论得到第一平面的平面场分布Ul(x,y)满足下式的条件：
Figure BDA0003362525630000022
(fx,fy)为频域坐标，并分别可以用下式表达为:
Figure BDA0003362525630000023
Figure BDA0003362525630000024
其中，x，y分别为透镜的第一平面的坐标系下的坐标值，T(fx,fy)为t(x0,y0)的傅里叶变换；
令孔径光阑面的坐标轴平行于傅里叶透镜的第一平面的坐标轴，两个坐标轴的坐标原点位于光轴上，且各自坐标系下的单位长度相等，则可用t(x0,y0)表示t(x,y)。
本发明能够取得以下技术效果：
本发明的方法为分析任意孔径光阑，在任意离焦面下的强度变化情况，为应用离焦面的光强变化趋势进行科学研究奠定理论基础，并解决了数值计算中出现的频域压缩的问题，该方法有着极为广泛的应用范围。
附图说明
图1为本发明公开的计算方法的流程图；
图2为本发明公开的光学***的结构示意图。
附图标记：
1、孔径光阑；2、傅里叶透镜；3、第一平面；4、第二平面。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及具体实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，而不构成对本发明的限制。
本发明提供了一种无限制像方弥散分布表征方法，即对任意孔径光阑任意离焦面处的光强分布进行准确计算的方法。图2示出了光学***的结构，光学***包括孔径光阑1和透镜。透镜选用在误差允许范围内，像差较小的透镜，如显微物镜、光刻物镜等。优选计算结果准确的傅里叶透镜2。光线经过任意形状、任意位置的孔径光阑，通过傅里叶透镜2达到焦平面和离焦面，计算方法包括如下步骤：
S0、获取光学***的光学参数值；
光学***包括垂直于光轴的第一平3面和第二平面4，第一平面3经过傅里叶类型的透镜2的入射面与光轴的交点，第二平面4经过透镜2的出射面与光轴的交点；
光学参数包括：入射光波振幅A，入射波波长λ，焦距f，离焦量Δz，光阑面平面场分布U0(x0,y0)，第一平面上的复振幅透过率t(x,y)，孔径光阑1与透镜2的第一平面3的距离d0，第一平面3、即紧靠透镜前表面的平面场分布Ul(x，y)，第二表面4、即紧靠透镜后表面的平面场分布Ul′(x,y)，任意的待测离焦面的场分布Uz(xz,yz)。特别的，当Δz＝0时，离焦面场分布等于焦平面场分布。
S1、计算孔径光阑面平面场分布的复振幅透过率t(x0,y0)；
根据角谱理论，其原理是计算得到光波传播到紧靠透镜前的平面上的场分布满足条件为：
Figure BDA0003362525630000041
因待测离焦面的场分布Uz(xz,yz)可以表示为第二平面4的平面场分布U′l(x，y)经过菲涅尔衍射形成的，(fx，fy)为频域坐标，并分别用下式表达为:
Figure BDA0003362525630000042
Figure BDA0003362525630000043
公式(2)-(3)表示频率取值与待测的离焦面的坐标的关系。其中，xz，yz为离焦面坐标系坐标值。T(fx，fy)为t(x0，y0)的傅里叶变换；光阑面的坐标系是(x0，y0)，第一表面的坐标系是(x,y)，t(x0,y0)和t(x,y)是同一个物理量复振幅透过率在不同坐标下的表示。通过坐标系的转换函数可以实现互相转换，实现t(x0,y0)到t(x,y)的坐标转换。后文中在选择相同单位长度的情况下，可用t(x,y)直接表示t(x0,y0)，不需要添加额外系数。
具体的，令光阑面的坐标轴平行于傅里叶透镜2的出射表面的坐标轴，两个坐标轴的坐标原点位于***光轴上，且各自坐标系下的单位长度相等，则t(x,y)可表示为t(x0,y0)。因此可以根据坐标位置可以计算出t(x，y)。计算后的结果可以代入后文的步骤S3中。
S2、计算待测离焦面处的离焦因子s(x，y)，离焦因子s(x，y)用下式表达为：
Figure BDA0003362525630000044
其中，k为波矢，x、y为第一表面的坐标系的坐标值。将计算后的结果可以代入后文的步骤S3中。
S3、将计算后的离焦因子s(x，y)和t(x，y)结果，代入公式(5)进行傅里叶变换计算，得到f(fx,fy)：
Figure BDA0003362525630000051
S4、计算待测离焦面的光强分布Iz(xz，yz)，光强分布Iz(xz，yz)用下式表达为：
Figure BDA0003362525630000052
具体地，离焦面的场分布Uz(xz，yz)可以表示为透镜后表面场分布U′l(x，y)经过菲涅尔衍射形成，运用菲涅尔衍射公式，化简得：
Figure BDA0003362525630000053
令
Figure BDA0003362525630000054
整理得：
Figure BDA0003362525630000055
根据频域卷积定理，则公式(8)变为：
Figure BDA0003362525630000056
如将公式(1)带入公式(9)，会出现右边的公式将同时包含(x,y)和(x0,y0)两种坐标系表示的量的现象，且这些量之间是非独立的，不能直接相乘，所以必须转换成统一的形式。所以需要通过两个坐标系设定相同的单位长度，使t(x0，y0)能用t(x,y)来表达，将变量形式统一到与(x,y)相关。
运用步骤S1中提到的t(x0，y0)到t(x，y)的坐标转换，代入S1中的公式(1)，用
Figure BDA0003362525630000057
表示位相因子，得到：
Figure BDA0003362525630000058
再反向运用频域卷积定理，整理得离焦面场分布函数为：
Figure BDA0003362525630000059
根据光强分布是场分布模的平方，又有位相因子|E|＝1，由此可得到离焦面光强表达式
Figure BDA0003362525630000061
根据光强分布是场分布模的平方，Iz(xz,yz)即为任意光阑任意离焦面处的光强分布，从公式(6)中可以看出，离焦面的光强分布正比于光阑面的复振幅透过率函数t(x,y)与离焦因子s(x,y)乘积的傅里叶变换的平方。复振幅透过率t(x，y)随光阑形状和位置的变化而改变，离焦因子s(x，y)代表离焦引起的位相分布变化，在***工作波长和焦距固定的情况下，其只与离焦量Δz有关。
在给定***焦距、工作波长条件下，通过本发明中采用的像方弥散分布表征方法，即可求出任意形状、任意离轴量下的孔径光阑，在任意离焦面的光强分布，从而为分析任意孔径光阑，在任意离焦面下的强度变化情况，为应用离焦面的光强变化趋势进行科学研究奠定理论基础，例如光学自动对焦***、球差补偿、景深扩展和超分辨技术等。
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
尽管上面已经示出和描述了本发明的实施例，可以理解的是，上述实施例是示例性的，不能理解为对本发明的限制，本领域的普通技术人员在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。
以上本发明的具体实施方式，并不构成对本发明保护范围的限定。任何根据本发明的技术构思所作出的各种其他相应的改变与变形，均应包含在本发明权利要求的保护范围内。

Claims (4)

1.无限制像方场分布的表征方法，用以计算光学***的光强分布；所述光学***包括孔径光阑和透镜，所述光学***中发射的光线经过所述孔径光阑，通过所述透镜到达焦平面和离焦面；所述光学***包括垂直于光轴的第一平面和第二平面，所述第一平面经过所述透镜的入射面与所述光轴的交点，所述第二平面经过所述透镜的出射面与所述光轴的交点；其特征在于，计算方法包括步骤：
S1、计算孔径光阑面平面场分布的复振幅透过率t(x0，y0)；
S2、计算待测离焦面处的离焦因子s(x，y)，所述离焦因子s(x，y)用下式表达为：
Figure FDA0003362525620000011
其中，f为焦距，Δz为离焦量，k为波矢，x、y为所述第一表面的坐标系上的坐标值；
S3、频域内计算所述复振幅透过率t(x0，y0)与所述离焦因子s(x，y)的傅里叶变换f(fx，fy)，所述傅里叶变换f(fx，fy)用下式表达为：
f(fx，fy)＝F{t(x0，y0)×s(x，y)} (2)；
S4、计算所述待测离焦面的光强分布Iz(xz，yz)，所述光强分布Iz(xz，yz)用下式表达为：
Figure FDA0003362525620000012
其中，A为入射光波振幅，λ为入射波波长。
2.根据权利要求1所述的无限制像方场分布的表征方法，其特征在于，在所述步骤S1之前还包括步骤S0，所述步骤S0包括：
S0、获取所述光学***的光学参数值：
所述光学参数值至少包括：所述入射光波振幅A，所述入射波波长λ，所述焦距f，所述离焦量Δz，光阑面平面场分布U0(x0，y0)，第一平面上的复振幅透过率t(x，y)，所述孔径光阑与所述透镜的距离d0，所述第一平面的平面场分布Ul(x，y)，所述第二表面的平面场分布U′l(x，y)，所述待测离焦面的场分布Uz(xz，yz)。
3.根据权利要求2所述的无限制像方场分布的表征方法，其特征在于，根据角谱理论得到所述第一平面的平面场分布Ul(x，y)满足下式的条件：
Figure FDA0003362525620000021
所述(fx，fy)为频域坐标，并分别可以用下式表达为：
Figure FDA0003362525620000022
Figure FDA0003362525620000023
其中，所述T(fx，fy)为所述t(x0，y0)的傅里叶变换；F{Ul(x，y)}为所述Ul(x，y)的傅里叶变换；
令孔径光阑面的坐标轴平行于所述透镜的第一平面的坐标轴，两个所述坐标轴的坐标原点位于所述光轴上，且各自坐标系下的单位长度相等，则可用所述t(x，y)表示所述t(x0，y0)。
4.根据权利要求1所述的无限制像方场分布的表征方法，其特征在于，所述透镜为傅里叶透镜。
CN202111371574.6A 2021-11-18 2021-11-18 无限制像方场分布的表征方法 Active CN114019677B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111371574.6A CN114019677B (zh) 2021-11-18 2021-11-18 无限制像方场分布的表征方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111371574.6A CN114019677B (zh) 2021-11-18 2021-11-18 无限制像方场分布的表征方法

Publications (2)

Publication Number Publication Date
CN114019677A true CN114019677A (zh) 2022-02-08
CN114019677B CN114019677B (zh) 2022-07-26

Family

ID=80065291

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111371574.6A Active CN114019677B (zh) 2021-11-18 2021-11-18 无限制像方场分布的表征方法

Country Status (1)

Country Link
CN (1) CN114019677B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102252763A (zh) * 2011-04-15 2011-11-23 中国科学院长春光学精密机械与物理研究所 光学成像***光学波前的单焦面高精度测试方法
CN102721477A (zh) * 2012-06-04 2012-10-10 中国科学院光电技术研究所 一种用于相位差波前传感器的离焦量测量方法
CN108426699A (zh) * 2018-01-29 2018-08-21 中国科学院长春光学精密机械与物理研究所 一种光学***像面离焦量的检测方法
CN111580283A (zh) * 2020-05-09 2020-08-25 南京理工大学 一种基于相位恢复的单透镜计算成像方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102252763A (zh) * 2011-04-15 2011-11-23 中国科学院长春光学精密机械与物理研究所 光学成像***光学波前的单焦面高精度测试方法
CN102721477A (zh) * 2012-06-04 2012-10-10 中国科学院光电技术研究所 一种用于相位差波前传感器的离焦量测量方法
CN108426699A (zh) * 2018-01-29 2018-08-21 中国科学院长春光学精密机械与物理研究所 一种光学***像面离焦量的检测方法
CN111580283A (zh) * 2020-05-09 2020-08-25 南京理工大学 一种基于相位恢复的单透镜计算成像方法

Also Published As

Publication number Publication date
CN114019677B (zh) 2022-07-26

Similar Documents

Publication Publication Date Title
Jin et al. Dielectric metasurfaces for distance measurements and three-dimensional imaging
US6330109B1 (en) Optical system comprising a diffractive optical element, and method of designing the same
Sheppard et al. Optical microscopy with extended depth of field
Claus High resolution digital holographic synthetic aperture applied to deformation measurement and extended depth of field method
JP4917404B2 (ja) 位相物体の可視化方法とその顕微鏡システム
TW201821856A (zh) 離軸非球面三反光學系統
Yumoto et al. Stitching interferometry for ellipsoidal x-ray mirrors
CN112926192A (zh) 一种离轴自由曲面光学***偏振像差分析方法
CN116953923B (zh) 一种超透镜设计方法及超透镜
US11671706B2 (en) Optical device comprising a multi-order diffractive Fresnel lens (MOD-DFL) and an achromatizing compensation mechanism, and a method for enhancing images captured using the MOD-DFL
Tao et al. The three‐dimensional point spread functions of a microscope objective in image and object space
CN114019677B (zh) 无限制像方场分布的表征方法
CN110196491A (zh) 一种双视场红外成像***散焦的建模仿真方法
TWI778526B (zh) 自由曲面凹面光柵成像光譜儀的設計方法
WO2015059862A1 (ja) 回折格子レンズ、それを有する光学系の設計方法、画像算出プログラムおよび回折格子レンズの製造方法
CN115421297A (zh) 光学***最佳焦面位置的确定方法
Xie et al. An extended depth-of-field imaging system with a non-rotationally symmetric phase mask
Zhuang et al. Chromatic focal shift of optical system expressed by related wavelength formulas
Burge An easy way to relate optical element motion to system pointing stability
Zhao et al. Angles-centroids fitting calibration and the centroid algorithm applied to reverse Hartmann test
Harvey et al. Performance predictions of a Schwarzschild imaging microscope for soft x-ray applications
CN211454009U (zh) 一种具有稳定轴上光强分布的新型长焦深光学装置
Lü et al. Fraunhofer diffraction of a Bessel beam focused by an aperture lens
Katkovnik et al. A novel binary and multilevel phase masks for enhanced depth-of-focus infrared imaging
CN115586642B (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