CN102607813B - 基于线光源的光学系统横向放大率测量方法 - Google Patents

基于线光源的光学系统横向放大率测量方法 Download PDF

Info

Publication number
CN102607813B
CN102607813B CN201210084537.1A CN201210084537A CN102607813B CN 102607813 B CN102607813 B CN 102607813B CN 201210084537 A CN201210084537 A CN 201210084537A CN 102607813 B CN102607813 B CN 102607813B
Authority
CN
China
Prior art keywords
spread function
function image
line
optical system
line source
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.)
Expired - Fee Related
Application number
CN201210084537.1A
Other languages
English (en)
Other versions
CN102607813A (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.)
Harbin Institute of Technology
Original Assignee
Harbin Institute of Technology
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 Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN201210084537.1A priority Critical patent/CN102607813B/zh
Publication of CN102607813A publication Critical patent/CN102607813A/zh
Application granted granted Critical
Publication of CN102607813B publication Critical patent/CN102607813B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Length Measuring Devices By Optical Means (AREA)

Abstract

基于线光源的光学系统横向放大率测量方法与装置属于以采用光学方法为特征的计量设备领域;本方法是以线光源为目标得到线状图像,在频域中利用实际测量得到的调制传递函数值第一次达到极小值时所对应的频率与理论截止频率相等,计算得到光学系统横向放大率;本装置在该装置光轴方向与图像传感器行或列方向所确定的平面内,线光源呈弯曲状,且所述的线光源上任意位置都准焦成像到图像传感器表面;采用本发明测量光学系统横向放大率,有利于减小单次测量结果之间的误差,进而提高测量结果重复性。

Description

基于线光源的光学系统横向放大率测量方法
技术领域
基于线光源的光学系统横向放大率测量方法与装置属于以采用光学方法为特征的计量设备领域,尤其涉及一种以线光源为目标,在频域利用线光源像来测量光学系统横向放大率的方法与装置。 
背景技术
光学系统横向放大率是医学以及精密测量领域中非常重要的参数,它不仅标明光学系统的技术指标,同样可以利用这项技术指标开展其它参数的精密测量。然而,如何获得一个光学系统的横向放大率,是开展这项工作的首要问题。 
一、光学系统横向放大率测量方法问题 
1987年07月,《医学物理》发表文章《论显微镜中物镜的放大率》,发现了显微镜中物镜的横向放大率经验公式与实际测量过程中产生的矛盾,该文章虽然没有给出物镜横向放大率的测量方法,但是该矛盾却引出了光学系统横向放大率的测量问题。 
而后续的一些文章,均显现出光学系统横向放大率测量的必要性。 
1999年03月,《黄山高等专科学校学报》第1卷第2期发表文章《关于几何光学中横向放大率的讨论》,该文章讨论了光学系统横向放大率的数学表达式,该方法的适用条件是近轴条件下的理想光学系统成像,而当这些条件不满足时,本文所总结的公式与实际光学系统横向放大率之间的误差却没有说明,更缺少对于这种误差,如何测量光学系统横向放大率方法的说明。 
2000年05月,《华南师范大学学报(自然科学版)》第2期发表文章《关于理想光具组横向放大率曲线的分析与运用》,本文根据光学系统横向放大率的定义总结出理想光具组横向放大率的计算公式,并绘制出横向放大率-物距像距曲线,该方法适用的条件仍然是理想光学系统近轴光线,而对于非理想条件下,经验公式中指出的横向放大率与实际横向放大率之间的误差却没有说明,更说明了测量光学系统横向放大率方法的必要性。 
2002年06月,《江西教育学院学报(自然科学)》第23卷第3期发表文章《用位相变换函数导出傍轴条件下透镜的物像距公式和横向放大率公式》,该文章以傅里叶光学为基础,利用透镜的相位变换作用推导出了傍轴条件下的物像距离公式以及光学系统横向放大率公式,然而,这篇文章的适用条件仍然是傍轴近似条件下的理想光学系统成像,同样具有同之前两篇文章相同的问题。 
因为存在光学系统横向放大率测量的迫切需求,所以在医学领域和精密测量等领域,均有学者提出自己的测量方法。 
2010年09月,《医学影像技术》第26卷增刊1发表文章《数字X线机固有放大率的测定》提供了一种放大率的测定方法,这种测定方法首先将小钢球固定在X线探测器上,摄片后用机器自带的标尺测出小球投影的直径;打印出照片,在阅片灯下用分规量好照片上投影的钢球直径,并用游标卡尺精确测出其数据,对比两组数据有误差异。同样用游标卡尺测出相应钢球的实际直径,可得出两直径之比,即X线影线放大率。由于该文章并非由精密测量领域人员所写,所以文章所采用的测量方法比较古老,沿用的是标尺测量物高,这种标尺测量具有一定的主观性,对测量结果影响较大。 
2003年09月,《河北职业技术师范学院学报》第17卷第3期发表文章《比较板法测望远镜放大率》,该文章介绍了一种光学系统横向放大率的新方法,这种方法与现行的普通物理实验所用方法相比,不仅原理简单、数据准确,而且更具操作性。然而这种方法任然没有摆脱传统方法的束缚,对像高的判断仍然沿用刻度尺读取目标长度的方法,因此同样具有主观性的问题。 
然而,这个问题随着CCD的迅速发展并广泛应用到精密测量领域而解决,同时,光学系统横向放大率的测量精度也相应得到了提高。 
1998年06月,《光电工程》第25卷第3期发表文章《CCD测望远系统放大率》,该文章介绍的方法原理简单,直接利用像高物高比来测定望远系统的放大率,该文章所介绍的方法与传统方法相比,像高不再采用标尺进行测量,而是通过刻线所占CCD像素个数与像素间距的乘积来判断,这种方法减少了测量过程中的主观因素,使测量结果更加准确。 
2002年03月,《物理实验》第22卷第3期发表文章《横向放大率法确定复合光学系统的基点》,2006年08月,《大学物理》第25卷第8期发表文章《横向放大率法测定光具组的基点》,这两篇文章将横向放大率扩展到了一个新的应用领域,用它来确定复合光学系统的基点,并得出重要结论,基点是光学系统横向放大率的函数。这个结论说明基点确定的准确与否直接与光学系统横向放大率的准确程度相关,因此,有必要精确测量光学系统横向放大率。而本文仍然沿用横向放大率的定义,即像高与物高比值进行测量,其中,像高的测量仍然沿用上一篇文章的测量原理,根据双缝所横跨的像素个数与像素间距的乘积来确定。 
对现有技术方法的陈述可以总结出以下结论,对于光学系统横向放大率的测量问题,无非是采用两种方法: 
1)利用光学系统横向放大率的定义,即像高和物高的比值来直接测量; 
2)根据光学系统横向放大率与某图像高度在特定光学系统中的特定关系,通过图像高度的获取实现对光学系统横向放大率的间接测量。 
无论是哪一种方法,均需要对像高进行判断,而现阶段的判断方法具有相同的技术特征: 
利用图像所横跨像素的个数与像素间距的乘积得到图像的高度信息。 
虽然该技术特征可以避免传统方法中用刻度尺测量像高的过程中的主观因素,但是这种方法也有自身的问题,因为对于像素个数的判断,只能是整数判断,每一侧的判断最多存在±0.5个像素的误差,两个边缘就可能存在±1个像素的误差,图像的尺寸越小,误差就会越大。虽然在理论上可以增大线光源的长度,通过用更多的像素来均摊误差得到弥补,但是对于大畸变光学系统,即不同视场下放大倍率不同的光学系统,增大线光源的长度同样会带来新的问题: 
1)增大目标尺寸,可能会使图像在长度上发生严重形变,这种情况下,不仅不能均摊误差,反而会使像素个数的判断误差更大,因此对于大畸变光学系统,该方法不适合在大视场范围内进行测量; 
2)针对大畸变光学系统,理应在每一个小视场范围内,精确测量该视场范围下的横向放大率,最终得到不同视场下的横向放大率曲线,但由于背景技术所采用的测量方法在小视场范围内单次测量结果之间误差较大,因此大畸变光学系统横向放大率测量重复性低。 
二、光学系统横向放大率测量装置问题 
国际专利分类号G01M11/02光学性质的测试领域,有两项发明专利公开了动像调制传递函数测量装置的组成: 
专利号ZL200810137150.1,授权公告日2010年09月29日,发明专利《动态目标调制传递函数测量方法与装置》,公开了一种高精度多功能的动像调制传递函数测量装置,该装置中也具有光源、光学系统以及图像传感器的结构,并且同样是光源经过光学系统成像到图像传感器表面。 
专利号ZL201010252619.3,授权公告日2012年01月11日,发明专利《动像调制传递函数测量装置》,在上一个专利所公开装置的基础上,进一步限定了装置中光学镜头的耦合方式以及测量的同步方式。 
但是这两项发明的特点是光源的运动轨迹是垂直于光轴的直线,对于有场曲的光学系统,光源运动的过程中,必然会造成图像的离焦,如果将这两项发明所公开的测量装置直接应用到本发明中,无法克服离焦造成的图像模糊问题以及图像灰度值变化问题,该问题会造成截止频率位置上的偏移,使测量结果的准确性受到影响。 
发明内容
本发明就是针对上述现有测量方法针对大畸变光学系统,不适合大视场范围内测量,而在小视场范围内,又存在横向放大率测量重复性低的问题,以及现有测量装置存在离焦的问题,提出了一种光学系统横向放大率的测量方法与装置,该方法可以在小视场范围内提高测量结果重复性,更适合测量大畸变光学系统横向放大率;该装置可以消除离焦对测量结果的影响,进一步提高测量结果重复性。 
本发明的目的是这样实现的: 
基于线光源的光学系统横向放大率测量方法,步骤如下: 
a、在物方放置长度为d的线光源,方向与图像传感器的行或列方向平行; 
b、图像传感器对线光源成像,得到初始点扩展函数图像;保持图像传感器曝光时间不变,移除线光源,图像传感器对背景成像,得到干扰图像,并将干扰图像中灰度值的最大值作为阈值; 
c.将第b步得到的初始点扩展函数图像中,线光源像所在行或列的整行或整列信息提取出来,作为初始线扩展函数图像,并初始线扩展函数图像中灰度值小于第b步所得阈值的像素的灰度值修正为0,得到修正线扩展函数图像,该修正线扩展函数图像具有n个元素; 
或者: 
将第b步得到的初始点扩展函数图像中,灰度值小于第b步所得阈值的像素的灰度值修正为0,作为修正点扩展函数图像;并将修正点扩展函数图像中,线光源像所在行或列的整行或整列信息提取出来,得到修正线扩展函数图像,该修正线扩展函数图像具有n个元素; 
d.对第c步得到的修正线扩展函数图像进行离散傅里叶变换并取模,得到调制传递函数图像,该调制传递函数图像具有同第c步得到的修正线扩展函数图像相同的元素个数n,即n个离散频谱分量,按照空间频率从小到大的顺序分别为M0、M1、M2、...、Mn-1,在该顺序下,调制传递函数值第一次达到极小值所对应的调制传递函数值为Mi,其下脚标序号为i,结合图像传感器的像素间距l,得到Mi所对应的空间频率值为:f=i/(nl); 
e.根据调制传递函数模型MTF(f)=|sinc(πfd′)|,结合第d步得到的空间频率值f,得到线光源像长度为:d′=1/f=nl/i; 
f.根据第a步线光源长度d和第e步得到的线光源像长度d′,计算得到光学系统横向放大率为:β=d′/d=nl/(id)。 
基于线光源的光学系统横向放大率测量装置,包括线光源、光学系统、图像传感器,所 述的线光源经过光学系统成像到图像传感器表面,并且,在该装置光轴方向与图像传感器行或列方向所确定的平面内,线光源呈弯曲状,且所述的线光源上任意位置都准焦成像到图像传感器表面。 
本发明的有益效果是: 
1)本发明采用的测量方法区别于传统空域测量方法,该方法以线光源为目标得到线状图像,在频域中利用实际测量得到的调制传递函数值第一次达到极小值时所对应的频率与理论截止频率相等,计算得到光学系统横向放大率;该特征使得采用长度短小的线光源时,会得到更高的截止频率,从而均摊截止频率的误差,使得单次测量结果之间的误差更小,进而提高测量结果重复性; 
2)本发明采用的测量装置在该装置光轴方向与图像传感器行或列方向所确定的平面内,线光源呈弯曲状,且所述的线光源上任意位置都准焦成像到图像传感器表面;该特征使测量得到的调制传递函数曲线更接近真实曲线,实际测量得到的截止频率位置更加准确,可以进一步减小单次测量结果之间的误差,提高测量结果重复性。 
附图说明
图1是基于线光源的光学系统横向放大率测量装置结构示意图 
图2是基于线光源的光学系统横向放大率测量装置平面光路图 
图3是基于线光源的光学系统横向放大率测量方法流程图 
图4是初始线扩展函数图像 
图5是修正线扩展函数图像 
图中:1线光源  2光学系统  3图像传感器 
具体实施方式
下面结合附图对本发明具体实施例做进一步详细描述。 
图1为基于线光源的光学系统横向放大率测量装置结构示意图,其平面光路图如图2所示;该装置包括线光源1、光学系统2、图像传感器3,所述的线光源1经过光学系统2成像到图像传感器3表面,并且,在该装置光轴方向与图像传感器3行方向所确定的平面内,线光源1呈弯曲状,且所述的线光源1上任意位置都准焦成像到图像传感器3表面;其中,线光源1的横向长度为3mm,图像传感器3的像素间距为5.6μm。 
基于线光源的光学系统横向放大率测量方法,流程图如图3所示,该方法步骤如下: 
a、在物方放置长度为d=3mm的线光源1,方向与图像传感器3行方向平行; 
b、图像传感器3对线光源1成像,得到初始点扩展函数图像;保持图像传感器3曝光时间不变,移除线光源1,图像传感器3对背景成像,得到干扰图像,并将干扰图像中灰度值的最大值作为阈值,该阈值为10; 
c.将第b步得到的初始点扩展函数图像中,线光源像所在行的整行信息提取出来,作为初始线扩展函数图像,如图4所示,并初始线扩展函数图像中灰度值小于第b步所得阈值的像素的灰度值修正为0,得到修正线扩展函数图像,如图5所示,该修正线扩展函数图像具有n=1280个元素; 
或者: 
将第b步得到的初始点扩展函数图像中,灰度值小于第b步所得阈值的像素的灰度值修正为0,作为修正点扩展函数图像;并将修正点扩展函数图像中,线光源像所在行的整行信息提取出来,得到修正线扩展函数图像,如图5所示,该修正线扩展函数图像具有n=1280个元素; 
d.对第c步得到的修正线扩展函数图像进行离散傅里叶变换并取模,得到调制传递函数图像,该调制传递函数图像具有同第c步得到的修正线扩展函数图像相同的元素个数n=1280,即1280个离散频谱分量,按照空间频率从小到大的顺序分别为M0、M1、M2、...、M1279,在该顺序下,调制传递函数值第一次达到极小值所对应的调制传递函数值为M42,其下脚标序号为i=42,结合图像传感器3的像素间距l=5.6μm,得到M42所对应的空间频率值为:f=i/(nl)=42/(1280×5.6×10-3)=5.8594lp/mm; 
e.根据调制传递函数模型MTF(f)=|sinc(πfd′)|,结合第d步得到的空间频率值f=5.8594lp/mm,得到线光源像长度为:d′=1/f=nl/i=1280×5.6×10-3/42=0.1707mm; 
f.根据第a步线光源1长度d=3mm和第e步得到的线光源像长度d′=0.1707mm,计算得到光学系统3横向放大率为:β=d′/d=nl/(id)=1280×5.6×10-3/(42×3)=0.0569。 

Claims (1)

1.基于线光源的光学系统横向放大率测量方法,其特征在于所述方法步骤如下:
a、在物方放置长度为d的线光源,方向与图像传感器的行或列方向平行;
b、图像传感器对线光源成像,得到初始点扩展函数图像;保持图像传感器曝光时间不变,移除线光源,图像传感器对背景成像,得到干扰图像,并将干扰图像中灰度值的最大值作为阈值;
c.将第b步得到的初始点扩展函数图像中,线光源像所在行或列的整行或整列信息提取出来,作为初始线扩展函数图像,并将初始线扩展函数图像中灰度值小于第b步所得阈值的像素的灰度值修正为0,得到修正线扩展函数图像,该修正线扩展函数图像具有n个元素;
或者:
将第b步得到的初始点扩展函数图像中,灰度值小于第b步所得阈值的像素的灰度值修正为0,作为修正点扩展函数图像;并将修正点扩展函数图像中,线光源像所在行或列的整行或整列信息提取出来,得到修正线扩展函数图像,该修正线扩展函数图像具有n个元素;
d.对第c步得到的修正线扩展函数图像进行离散傅里叶变换并取模,得到调制传递函数图像,该调制传递函数图像具有同第c步得到的修正线扩展函数图像相同的元素个数n,即n个离散频谱分量,按照空间频率从小到大的顺序分别为M0、M1、M2、...、Mn-1,在该顺序下,调制传递函数值第一次达到极小值所对应的调制传递函数值为Mi,其下脚标序号为i,结合图像传感器的像素间距l,得到Mi所对应的空间频率值为:f=i/(nl);
e.根据调制传递函数模型MTF(f)=|sinc(πfd′)|,结合第d步得到的空间频率值f,得到线光源像长度为:d′=1/f=nl/i;
f.根据第a步线光源长度d和第e步得到的线光源像长度d′,计算得到光学系统横向放大率为:β=d′/d=nl/(id)。
CN201210084537.1A 2012-03-17 2012-03-17 基于线光源的光学系统横向放大率测量方法 Expired - Fee Related CN102607813B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210084537.1A CN102607813B (zh) 2012-03-17 2012-03-17 基于线光源的光学系统横向放大率测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210084537.1A CN102607813B (zh) 2012-03-17 2012-03-17 基于线光源的光学系统横向放大率测量方法

Publications (2)

Publication Number Publication Date
CN102607813A CN102607813A (zh) 2012-07-25
CN102607813B true CN102607813B (zh) 2014-10-22

Family

ID=46525381

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210084537.1A Expired - Fee Related CN102607813B (zh) 2012-03-17 2012-03-17 基于线光源的光学系统横向放大率测量方法

Country Status (1)

Country Link
CN (1) CN102607813B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102607817B (zh) * 2012-03-17 2014-12-03 哈尔滨工业大学 利用线光源的光学系统横向放大率测量方法
CN102620914B (zh) * 2012-03-17 2014-12-03 哈尔滨工业大学 采用线光源的光学系统横向放大率测量方法
CN108106818B (zh) * 2017-12-11 2020-04-03 中国科学院上海光学精密机械研究所 光学成像系统倍率与畸变高精度测量装置及测量方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3938892A (en) * 1973-04-24 1976-02-17 The United States Of America As Represented By The United States National Aeronautics And Space Administration Electronic optical transfer function analyzer
JP2004037410A (ja) * 2002-07-08 2004-02-05 Yucaly Optical Laboratory Inc 変調伝達関数測定装置および変調伝達関数測定方法
CN101268353A (zh) * 2005-09-21 2008-09-17 诺瓦提斯公司 光学观察系统
CN101354307A (zh) * 2008-09-22 2009-01-28 哈尔滨工业大学 动态目标调制传递函数测量方法与装置
CN102607817A (zh) * 2012-03-17 2012-07-25 哈尔滨工业大学 利用线光源的光学系统横向放大率测量方法与装置
CN102620914A (zh) * 2012-03-17 2012-08-01 哈尔滨工业大学 采用线光源的光学系统横向放大率测量方法与装置

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3938892A (en) * 1973-04-24 1976-02-17 The United States Of America As Represented By The United States National Aeronautics And Space Administration Electronic optical transfer function analyzer
JP2004037410A (ja) * 2002-07-08 2004-02-05 Yucaly Optical Laboratory Inc 変調伝達関数測定装置および変調伝達関数測定方法
CN101268353A (zh) * 2005-09-21 2008-09-17 诺瓦提斯公司 光学观察系统
CN101354307A (zh) * 2008-09-22 2009-01-28 哈尔滨工业大学 动态目标调制传递函数测量方法与装置
CN102607817A (zh) * 2012-03-17 2012-07-25 哈尔滨工业大学 利用线光源的光学系统横向放大率测量方法与装置
CN102620914A (zh) * 2012-03-17 2012-08-01 哈尔滨工业大学 采用线光源的光学系统横向放大率测量方法与装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
于谦.光学检查镜头的MTF测试方法研究.《中国优秀硕士学位论文全文数据库·信息科技辑》.2010,(第08期),全文.
光学检查镜头的MTF测试方法研究;于谦;《中国优秀硕士学位论文全文数据库·信息科技辑》;20100815(第08期);全文 *

Also Published As

Publication number Publication date
CN102607813A (zh) 2012-07-25

Similar Documents

Publication Publication Date Title
Kumar et al. COSMOGRAIL: the COSmological MOnitoring of GRAvItational Lenses-XIV. Time delay of the doubly lensed quasar SDSS J1001+ 5027
CN102607815B (zh) 基于点目标像拼合的光学系统横向放大率测量方法与装置
CN102620910B (zh) 基于匀速动点目标的光学系统横向放大率测量方法
CN102620912B (zh) 利用点目标像拼合的光学系统横向放大率测量方法与装置
CN103575227A (zh) 一种基于数字散斑的视觉引伸计实现方法
CN103606147B (zh) 多台不共视场量测相机的坐标系转换标定方法
CN102607441B (zh) 利用匀速动点目标的图像传感器像素间距测量方法与装置
Hui et al. Line-scan camera calibration in close-range photogrammetry
CN102620911B (zh) 采用点目标像拼合的光学系统横向放大率测量方法与装置
CN102620668B (zh) 采用点目标像拼合的图像传感器像素间距测量方法与装置
CN102620913B (zh) 采用匀速动点目标的光学系统横向放大率测量方法与装置
CN102607813B (zh) 基于线光源的光学系统横向放大率测量方法
Gao et al. Stereo camera calibration for large field of view digital image correlation using zoom lens
CN108362226A (zh) 提高图像过曝区域相位测量精度的双四步相移法
CN102620914B (zh) 采用线光源的光学系统横向放大率测量方法
CN102607816B (zh) 利用匀速动点目标的光学系统横向放大率测量方法
Huang et al. Distance calibration between reference plane and screen in direct phase measuring deflectometry
CN102607442B (zh) 基于匀速动点目标的图像传感器像素间距测量方法与装置
Liu et al. Full-field calibration of color camera chromatic aberration using absolute phase maps
CN102607817B (zh) 利用线光源的光学系统横向放大率测量方法
CN102620667B (zh) 利用点目标像拼合的图像传感器像素间距测量方法与装置
CN102607444A (zh) 采用线光源的图像传感器像素间距测量方法与装置
Deng et al. A calibration method for system parameters in direct phase measuring deflectometry
CN102620671B (zh) 利用线光源的图像传感器像素间距测量方法与装置
CN102620670A (zh) 基于线光源的图像传感器像素间距测量方法与装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20141022