CN101644606A - 光学成像系统波相差的测试方法 - Google Patents
光学成像系统波相差的测试方法 Download PDFInfo
- Publication number
- CN101644606A CN101644606A CN200910067352A CN200910067352A CN101644606A CN 101644606 A CN101644606 A CN 101644606A CN 200910067352 A CN200910067352 A CN 200910067352A CN 200910067352 A CN200910067352 A CN 200910067352A CN 101644606 A CN101644606 A CN 101644606A
- Authority
- CN
- China
- Prior art keywords
- optical system
- image
- wave
- phase
- tested
- 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.)
- Pending
Links
Images
Landscapes
- Testing Of Optical Devices Or Fibers (AREA)
Abstract
本发明涉及对光学成像系统的成像质量进行精确测评的检测方法,特别是一种光学成像系统波相差的测试方法。被照明光源照亮的星点目标板发出的光经平行光管后成为平行光,被测试光学系统接收此平行光后在其光学系统的像面上形成星点目标板的图像,由所述的设有显微物镜摄像头的CCD像机接收此图像后,经所述带有图像采集卡的计算机对图像进行采集,通过傅里叶变换计算和迭代算法实现由光学系统的星点图像获得波相差的数学分析与计算,使得光学系统的波相差复原,从而获得被测试系统的波相差图。其测试结果准确性高、所需测试技术条件易于实现、测试成本低。
Description
技术领域
本发明涉及对光学成像系统的成像质量进行精确测评的一种测试光学系统波相差的方法。
背景技术
所谓光学系统的波相差是指光学成像系统的光学波面与理想波面的位相误差。光学系统的波相差可以真实与直观地表述光学成像系统的质量状态。通常情况下,通过对光学成像系统波相差的测试可对光学成像系统的成像质量进行精确测评,掌握光学成像系统内部部件所处状态,用于指导光学系统的加工与安装、调试。
对无穷远距离点目标成像的光学系统,在其焦平面(像面上)上所成的像表述了光学系统成像质量,这个点目标像在光学上称为点扩散函数(PSF)。在光学成像质量测评领域中,可以通过对PSF的采集与测试而获得光学系统的传递函数(OTF),进而实现对光学成像系统分辨能力的测评。
光学系统的波相差与点扩散函数、光学传递函数之间存在着数学变换关系,可通过傅立叶变换或自相关计算实现从光学系统波相差到点扩散函数和光学传递函数的变换。本发明就是通过测试获得的点扩散函数图像,经过数学计算分析而得到被测试光学系统的波相差,从而掌握光学系统调试状态,判断光学系统是否装调到位。
目前测试光学系统的波相差是通过采用干涉方法,由干涉仪、基准平面镜和被测试光学成像系统构成测试光路,通过自准方式实现对光学系统波相差的测试。
参见图6,干涉仪测试光学成像系统的波相差测试过程,是由干涉仪12发出球面光波,球面光波的球心调试到被测试光学成像系统5像面上,球面光波通过被测试光学成像系统后转为平面波(平行光),此光束经平面反射镜11反射后,按原路返回,进入被测试光学成像系统5,在被测试光学成像系统5后形成球面波,在其像面上会聚后进入干涉仪12,与干涉仪12的基准光学波面干涉,产生干涉条纹,经过其内部图像采集、测试与计算系统处理后,给出被测试光学成像系统光学波像差。
这种方法的使用条件有一定的局限性。除了需要干涉仪外还要有相应口径面形良好的平面反射镜作为测试基准,因此,平面镜的面形质量将直接影响测试结果的准确性。对于口径1米~2米、面形误差均方值小于λ/40(λ=632.8nm)的平面基准镜在研制方面需要高成本,而且加工风险与检测难度极大,不易实现。
发明内容
本发明的目的是为克服目前惯用的干涉方法测试光学系统波相差存在的上述缺点,提出一种新型的光学成像系统波相差的测试方法,其测试结果准确性高、所需测试技术条件易于实现、测试成本低。
本发明光学成像系统波相差的测试方法,是通过在被测试光学系统的两端、沿该被测试光学系统光轴设置的平行光管和设有显微物镜摄像头的CCD像机,设置在平行光管的焦面上的星点目标板及与CCD像机相连接的带有图像采集卡的计算机按以下方式测得被测试光学系统的波相差:
被照明光源照亮的星点目标板发出的光经平行光管后成为平行光(模拟无穷远距离目标);被测试光学系统接收此平行光后在其光学系统的像面上形成星点目标板的图像;由所述的设有显微物镜摄像头的CCD像机接收此图像后,经所述带有图像采集卡的计算机对图像进行采集;通过傅里叶变换计算和迭代算法实现由光学系统的星点图像获得波相差的数学分析与计算,使得光学系统的波相差复原,从而获得被测试系统的波相差图。
本发明是通过采集被测试光学成像系统对点目标所成图像,基于信息光学理论和位相复原方法经计算分析而获得光学系统的波相差。本测试方法的具有以下优点:
1.测试成本低。在测试大口径光学成像系统波相差时,不需要相应口径的平面基准镜,所用设备为相应口径的平行光管。关键部件为口径与被测试光学系统口径匹配的主镜,其加工与检测与大口径的平面镜相比,在造价和难度上都大大降低,易于实现;同时一台相移干涉仪的产品价格(约100万元)远远高于星点图象采集分析系统的价格(约30万元)。
2.测试效率高。采用本发明方法测试系统,在测试光学系统波相差的同时,还可以实现光学系统分辨率、焦距和传递函数等光学系统主要技术指标的测试。当把平行光管焦面上的点目标板更换为分辨率目标板时可测试光学系统的分辨率;当把平行光管焦面上的点目标板更换为分划板时可测试光学系统的焦距;而用星点目标或狭缝目标时可以测试光学系统的传递函数,因此,大幅度提高检测效率。
3.测试波段广。由于干涉仪所采用的测量原理限制所使用的测试光源波段只能为某一单波段,以便产生干涉图;而本发明测试波段可通过在平行光管目标板附近加用不同波段的滤光片而实现多波段的选择,这种功能更能满足测评光学成像系统实际使用情况的需要。
附图说明
图1是本发明测试方法的原理示意图;
图2是CCD像机采集到的星点(PSF)原始图片;
图3是经平滑、去噪声处理后的星点图片;
图4是由星点图经过位相复原方法获得的光学系统波相差图片;
图5是由光学系统波相差图获得的干涉图片;
图6是目前干涉法测试光学系统波相差的原理示意图。
具体实施方式
以下结合附图给出的实施例对本发明方法作进一步详细说明。
对一种光学成像系统波相差的测试方法。
参照图1,将被测试光学成像系统5放置在一光学稳定平台10上,通过在被测试光学系统5的两端、沿该被测试光学系统光轴设置的平行光管4和设有显微物镜摄像头6的CCD像机7,在平行光管4的焦面上设置一由照明光源1照亮的星点目标板2,在星点目标板2和平行光管4之间设置一与欲测试波段相应的滤光片3,通过数据信号传输线8与CCD像机7相连接的带有图像采集卡的计算机9构成的测试系统,按以下方式测得被测试光学系统5的波相差:
被照明光源1照亮的星点目标板2发出的光经滤光片3、平行光管4后成为平行光(模拟无穷远距离目标)射出;被测试光学系统5接收此平行光后在其光学系统的像面上形成星点目标板的图像;所述的设有显微物镜摄像头6的CCD像机7接收此图像后把光信号转变为电信号,经由数据信号传输线8把电信号传输给带有图像采集卡的计算机9,计算机9对原始图像(如图2所示)通过相应的图像处理得到星点图像(如图3所示);然后经过位相复原计算软件对该星点图像的信息进行分析计算,最终给出被测试光学系统的波相差图象和数据信息(如图4、图5所示)。
所述位相复原计算软件采用公式:
光学系统的光瞳函数与振幅点扩散函数是傅立叶变换关系,一般用如下公式表示:
ASF(x,y)=∫∫P(x′,y′)exp[-2πi(xx′+yy′)]dx′dy′......................(1)
其中ASF(x,y)为振幅点扩散函数,P(x′,y′)为光瞳函数。
而P(x′,y′)=A(x′,y′)exp[ikW(x′,y′)]...............................(2)
A(x,y)为光瞳函数的振幅部分,是光瞳形状函数,在通光孔径内取1,通光孔径外为零;k=2π/λ,W(x′,y′)为光学系统的波相差。
而PSF(x,y)=ASF(x,y)×ASF*(x,y)=|ASF(x,y)|2...............(3)
因此,可以根据已知的PSF(x,y)即星点图和(1)、(2)、(3)式采用傅里叶变换方法和迭代方法实现由PSF(x,y)→ASF(x,y)→P(x′,y′)→W(x′,y′)的数学分析与计算,使得光学系统的波相差W(x′,y′)复原。
Claims (1)
1.一种光学成像系统波相差的测试方法,其特征在于,是通过在被测试光学系统5的两端、沿该被测试光学系统光轴设置的平行光管4和设有显微物镜摄像头6的CCD像机7,设置在平行光管4的焦面上的星点目标板2及与CCD像机7相连接的带有图像采集卡的计算机9按以下方式测得被测试光学系统5的波相差:
被照明光源1照亮的星点目标板2发出的光经平行光管4后成为平行光;被测试光学系统5接收此平行光后在其光学系统的像面上形成星点目标板的图像;由所述的设有显微物镜摄像头6的CCD像机7接收此图像后,经所述带有图像采集卡的计算机9对图像进行采集;通过傅里叶变换计算和迭代算法实现由光学系统的星点图像获得波相差的数学分析与计算,使得光学系统的波相差复原,从而获得被测试系统的波相差图。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200910067352A CN101644606A (zh) | 2009-07-28 | 2009-07-28 | 光学成像系统波相差的测试方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200910067352A CN101644606A (zh) | 2009-07-28 | 2009-07-28 | 光学成像系统波相差的测试方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN101644606A true CN101644606A (zh) | 2010-02-10 |
Family
ID=41656586
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200910067352A Pending CN101644606A (zh) | 2009-07-28 | 2009-07-28 | 光学成像系统波相差的测试方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101644606A (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101893485A (zh) * | 2010-07-16 | 2010-11-24 | 中国科学院长春光学精密机械与物理研究所 | 星点衍射波面相差的检测方法 |
CN102252763A (zh) * | 2011-04-15 | 2011-11-23 | 中国科学院长春光学精密机械与物理研究所 | 光学成像系统光学波前的单焦面高精度测试方法 |
CN102840964A (zh) * | 2012-08-27 | 2012-12-26 | 中国科学院长春光学精密机械与物理研究所 | 大口径、长焦距平行光管焦点实时监测系统 |
CN103207022A (zh) * | 2013-03-18 | 2013-07-17 | 中国科学院长春光学精密机械与物理研究所 | 相位复原测试过程中消除旋转非对称误差的标定方法 |
CN103697808A (zh) * | 2014-01-14 | 2014-04-02 | 中国科学院长春光学精密机械与物理研究所 | 点衍射移相式干涉仪高度传递函数测试装置及其测试方法 |
CN103698108A (zh) * | 2014-01-14 | 2014-04-02 | 中国科学院长春光学精密机械与物理研究所 | 球面波前测试干涉仪高度传递函数测试装置及其测试方法 |
CN109297683A (zh) * | 2018-10-11 | 2019-02-01 | 广州博冠光电科技股份有限公司 | 一种基于数字激光全息快速检测透镜中心偏的装置及方法 |
WO2020113855A1 (zh) * | 2018-12-03 | 2020-06-11 | 上海慧希电子科技有限公司 | 光学检测方法、系统及光学器件制造系统 |
CN115641376A (zh) * | 2022-10-17 | 2023-01-24 | 中国科学院长春光学精密机械与物理研究所 | 望远镜在轨位姿失调检测方法、装置、设备及介质 |
CN117148572A (zh) * | 2023-11-01 | 2023-12-01 | 中国科学院长春光学精密机械与物理研究所 | 一种用于模拟跟踪物体的装置、方法以及介质 |
-
2009
- 2009-07-28 CN CN200910067352A patent/CN101644606A/zh active Pending
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101893485A (zh) * | 2010-07-16 | 2010-11-24 | 中国科学院长春光学精密机械与物理研究所 | 星点衍射波面相差的检测方法 |
CN102252763A (zh) * | 2011-04-15 | 2011-11-23 | 中国科学院长春光学精密机械与物理研究所 | 光学成像系统光学波前的单焦面高精度测试方法 |
CN102252763B (zh) * | 2011-04-15 | 2012-09-19 | 中国科学院长春光学精密机械与物理研究所 | 光学成像系统光学波前的单焦面高精度测试方法 |
CN102840964A (zh) * | 2012-08-27 | 2012-12-26 | 中国科学院长春光学精密机械与物理研究所 | 大口径、长焦距平行光管焦点实时监测系统 |
CN103207022A (zh) * | 2013-03-18 | 2013-07-17 | 中国科学院长春光学精密机械与物理研究所 | 相位复原测试过程中消除旋转非对称误差的标定方法 |
CN103698108A (zh) * | 2014-01-14 | 2014-04-02 | 中国科学院长春光学精密机械与物理研究所 | 球面波前测试干涉仪高度传递函数测试装置及其测试方法 |
CN103697808A (zh) * | 2014-01-14 | 2014-04-02 | 中国科学院长春光学精密机械与物理研究所 | 点衍射移相式干涉仪高度传递函数测试装置及其测试方法 |
CN103697808B (zh) * | 2014-01-14 | 2016-07-06 | 中国科学院长春光学精密机械与物理研究所 | 点衍射移相式干涉仪高度传递函数测试装置及其测试方法 |
CN109297683A (zh) * | 2018-10-11 | 2019-02-01 | 广州博冠光电科技股份有限公司 | 一种基于数字激光全息快速检测透镜中心偏的装置及方法 |
CN109297683B (zh) * | 2018-10-11 | 2024-02-20 | 广州博冠光电科技股份有限公司 | 一种基于数字激光全息快速检测透镜中心偏的装置及方法 |
WO2020113855A1 (zh) * | 2018-12-03 | 2020-06-11 | 上海慧希电子科技有限公司 | 光学检测方法、系统及光学器件制造系统 |
US11841288B2 (en) | 2018-12-03 | 2023-12-12 | Shanghai Intelight Electronic Technology Co., Ltd. | Optical measurement method and system and optical device manufacturing system |
CN115641376A (zh) * | 2022-10-17 | 2023-01-24 | 中国科学院长春光学精密机械与物理研究所 | 望远镜在轨位姿失调检测方法、装置、设备及介质 |
CN117148572A (zh) * | 2023-11-01 | 2023-12-01 | 中国科学院长春光学精密机械与物理研究所 | 一种用于模拟跟踪物体的装置、方法以及介质 |
CN117148572B (zh) * | 2023-11-01 | 2024-01-02 | 中国科学院长春光学精密机械与物理研究所 | 一种用于模拟跟踪物体的装置、方法以及介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101644606A (zh) | 光学成像系统波相差的测试方法 | |
CN103926197B (zh) | 高空间分辨双轴差动共焦图谱显微成像方法与装置 | |
CN102499648B (zh) | 一种基于菲涅尔光谱仪的频域光学相干层析成像系统 | |
CN111751012B (zh) | 动态高分辨光学波前相位测量方法 | |
CN105588519A (zh) | 利用相位差异相位恢复技术检测大口径望远镜面形的方法 | |
CN103592108A (zh) | Ccd芯片调制传递函数测试装置及方法 | |
CN102506759A (zh) | 一种大口径非球面朗奇检测方法 | |
JP2014016253A (ja) | 屈折率分布計測方法、光学素子の製造方法、および、屈折率分布計測装置 | |
CN105318847A (zh) | 基于系统建模的非球面非零位环形子孔径拼接方法 | |
US20110112785A1 (en) | Measuring method and measuring apparatus | |
Wu et al. | An accurate method for shape retrieval and displacement measurement using bi-prism-based single lens 3D digital image correlation | |
CN103024427B (zh) | 相机调制传递函数的测试方法及其测试装置 | |
CN102889980A (zh) | 一种基于光栅剪切干涉检测系统的微透镜定焦检测方法 | |
CN110017791B (zh) | 光纤连接器端面参数测量装置及测量方法 | |
CN205482829U (zh) | 一种浅度非球面动态干涉检测装置 | |
CN203216701U (zh) | 传像光纤束像差检测装置 | |
CN111397634B (zh) | 星敏感器固定端面热变形的高分辨干涉检测装置及方法 | |
CN110631510B (zh) | 一种基于迈克尔逊结构的高精度测角装置及测角方法 | |
CN105571516A (zh) | 一种全视场低频外差干涉仪 | |
CN104634275A (zh) | 一种基于牛顿环的非球面实时干涉测量装置及方法 | |
CN101893485A (zh) | 星点衍射波面相差的检测方法 | |
CN104359654A (zh) | 一种光纤传像束两端面像元偏移量的测量装置及测量方法 | |
CN105716521B (zh) | 增大频域低相干光干涉测量范围的装置及其方法 | |
CN113063367A (zh) | 一种基于斜光轴数字图像相关方法的全场挠度实时测量系统及测量方法 | |
KR20140132725A (ko) | 광학 컴포넌트의 기하 구조 측정 방법 및 장치 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Open date: 20100210 |