CN102620911B - 采用点目标像拼合的光学系统横向放大率测量方法与装置 - Google Patents

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)针对大畸变光学系统，理应在每一个小视场范围内，精确测量该视场范围下的横向放大率，最终得到不同视场下的横向放大率曲线，但由于背景技术所采用的测量方法在小视场范围内单次测量结果之间误差较大，因此大畸变光学系统横向放大率测量重复性低。

二、光学系统横向放大率测量装置问题

国际专利分类号G01M 11/02光学性质的测试领域，有两项发明专利公开了动像调制传递函数测量装置的组成：

专利号ZL200810137150.1，授权公告日2010年09月29日，发明专利《动态目标调制传递函数测量方法与装置》，公开了一种高精度多功能的动像调制传递函数测量装置，该装置中也具有光源、光学系统以及图像传感器的结构，并且同样是光源经过光学系统成像到图像传感器表面。

专利号ZL201010252619.3，授权公告日2012年01月11日，发明专利《动像调制传递函数测量装置》，在上一个专利所公开装置的基础上，进一步限定了装置中光学镜头的耦合方式以及测量的同步方式。

但是这两项发明的特点是光源的运动轨迹是垂直于光轴的直线，对于有场曲的光学系统，光源运动的过程中，必然会造成图像的离焦，如果将这两项发明所公开的测量装置直接应用到本发明中，无法克服离焦造成的图像模糊问题以及图像灰度值变化问题，该问题会造成截止频率位置上的偏移，使测量结果的准确性受到影响。

发明内容

本发明就是针对上述现有测量方法针对大畸变光学系统，不适合大视场范围内测量，而在小视场范围内，又存在横向放大率测量重复性低的问题，以及现有测量装置存在离焦的问题，提出了一种光学系统横向放大率的测量方法与装置，该方法可以在小视场范围内提高测量结果重复性，更适合测量大畸变光学系统横向放大率；该装置可以消除离焦对测量结果的影响，进一步提高测量结果重复性。

本发明的目的是这样实现的：

采用点目标像拼合的光学系统横向放大率测量方法，步骤如下：

a.图像传感器对静态点目标第一次成像，得到第一帧初始静态点目标图像，并提取点目标像所在像素坐标位置(x₁，y₁)；

b.使点目标沿图像传感器行或列方向运动，位移量为d，之后保持点目标为静止状态；

c.保持图像传感器曝光时间不变，图像传感器对静态点目标第二次成像，得到第二帧初始静态点目标图像，并提取点目标像所在像素坐标位置(x₂，y₂)；

d.移除点目标并保持图像传感器曝光时间不变，图像传感器对背景成像，得到干扰图像，并将干扰图像中灰度值的最大值作为阈值；

e.第a步得到的第一帧初始静态点目标图像，将灰度值小于第d步所得阈值的像素的灰度值修正为0，得到第一帧修正静态点目标图像；第c步得到的第二帧初始静态点目标图像，将灰度值小于第d步所得阈值的像素的灰度值修正为0，得到第二帧修正静态点目标图像；将第一帧修正静态点目标图像和第二帧修正静态点目标图像相叠加，叠加后新图像中两个点目标像所在行或列的所有像素灰度值相加并除以2，得到新灰度值；并将第a步得到的像素坐标位置(x₁，y₁)和第c步得到的像素坐标位置(x₂，y₂)连线所覆盖的像素的灰度值替换为新灰度值，得到构造点扩展函数图像；

或者：

将第a步得到的第一帧初始静态点目标图像和第c步得到的第二帧初始静态点目标图像相叠加，并将叠加后图像中灰度值小于第d步所得阈值2倍的像素的灰度值修正为0，得到修正叠加图像；将修正叠加图像中两个点目标像所在行或列的所有像素灰度值相加并除以2，得到新灰度值；并将第a步得到的像素坐标位置(x₁，y₁)和第c步得到的像素坐标位置(x₂，y₂)连线所覆盖的像素的灰度值替换为新灰度值，得到构造点扩展函数图像；

f.第e步得到的构造点扩展函数图像，将线状光斑所在行或列的整行或整列信息提取出来，作为构造线扩展函数图像，该构造线扩展函数图像具有n个元素；

g.对第f步得到的构造线扩展函数图像进行离散傅里叶变换并取模，得到调制传递函数图像，该调制传递函数图像具有同第f步得到的构造线扩展函数图像相同的元素个数n，即n个离散频谱分量，按照空间频率从小到大的顺序分别为M₀、M₁、M₂、...、M_n-1，在该顺序下，调制传递函数值第一次达到极小值所对应的调制传递函数值为M_i，其下脚标序号为i，结合图像传感器的像素间距l，得到M_i-1和M_i+1所对应的空间频率值分别为：f_min＝(i-1)/(nl)以及f_max＝(i+1)/(nl)；

h.根据调制传递函数模型MTF(f)＝|sinc(πfd′)|，结合第g步得到的空间频率范围f_min和f_max，得到线状光斑长度取值范围：d_max′＝1/f_min＝nl/(i-1)以及d_min′＝1/f_max＝nl/(i+1)；

i.根据第b步点目标位移量d和第h步得到的线状光斑长度取值范围，计算得到光学系统横向放大率取值范围为：β_min＝d_min′/d＝nl/((i+1)d)以及β_max＝d_max′/d＝nl/((i-1)d)；

j.将β为变量，并以第i步得到的β_min和β_max为取值范围，在第g步得到的n个调制传递函数值中选取K个作为对比数据，这K个调制传递函数值分别是M_K1、M_K2、...、M_KK，采用遗传算法寻找到以下公式的最小值：该最小值所对应的光学系统横向放大率β即为所求。

采用点目标像拼合的光学系统横向放大率测量装置，包括点目标、光学系统、图像传感器、滑块以及垂直光轴方向的第一导轨，所述的点目标经过光学系统成像到图像传感器表面；并且，该装置还包括沿光轴方向的第二导轨，承载点目标的滑块安装在第一导轨和第二导轨上，滑块在第一导轨上的运动与滑块在第二导轨上的运动相配合，使点目标在任意视场位置都准焦成像到图像传感器表面。

本发明的有益效果是：

1)本发明采用的测量方法区别于传统空域测量方法，该方法以点光源为目标，使点目标处于不同视场下并对其两次成像，根据两幅点目标像构造出线状图像，在频域中寻找像素间距的取值范围，并根据与像素间距相关的实际调制传递函数曲线与理论调制传递函数曲线在最小二乘条件下重合度最好，利用遗传算法计算得到光学系统横向放大率；该特征使得采用长度短小的线光源时，会得到更高的截止频率，从而均摊截止频率的误差，使得单次测量结果之间的误差更小，进而提高测量结果重复性；

2)本发明采用的测量装置包括沿光轴方向的第二导轨，承载点目标的滑块安装在第一导轨和第二导轨上，滑块在第一导轨上的运动与滑块在第二导轨上的运动相配合，使点目标在任意视场位置都准焦成像到图像传感器表面；该特征使测量得到的调制传递函数曲线更接近真实曲线，实际测量得到的截止频率位置更加准确，可以进一步减小单次测量结果之间的误差，提高测量结果重复性。

附图说明

图1是采用点目标像拼合的光学系统横向放大率测量装置结构示意图

图2是采用点目标像拼合的光学系统横向放大率测量方法流程图

图3是第一帧初始静态点目标图像局部图

图4是第二帧初始静态点目标图像局部图

图5是构造点扩展函数图像局部图

图6是构造线扩展函数图像

图中：1点目标 2光学系统 3图像传感器 4滑块 5第一导轨 6第二导轨

具体实施方式

下面结合附图对本发明具体实施例做进一步详细描述。

图1为采用点目标像拼合的光学系统横向放大率测量装置结构示意图；该装置包括点目标1、光学系统2、图像传感器3、滑块4以及垂直光轴方向的第一导轨5，所述的点目标1经过光学系统2成像到图像传感器3表面；并且，该装置还包括沿光轴方向的第二导轨6，承载点目标1的滑块4安装在第一导轨5和第二导轨6上，滑块4在第一导轨5上的运动与滑块4在第二导轨6上的运动相配合，使点目标1在任意视场位置都准焦成像到图像传感器3表面；其中，点目标1为直径15μm的针孔，图像传感器3的像素间距为5.6μm。

采用点目标像拼合的光学系统横向放大率测量方法，流程图如图2所示，该方法步骤如下：

a.图像传感器3对静态点目标1第一次成像，得到第一帧初始静态点目标图像，其局部图如图3所示，并提取点目标像所在像素坐标位置(x₁，y₁)＝(371，539)；

b.使点目标1沿图像传感器3行方向运动，位移量为d＝6.104mm，之后保持点目标1为静止状态；

c.保持图像传感器3曝光时间不变，图像传感器3对静态点目标1第二次成像，得到第二帧初始静态点目标图像，其局部图如图4所示，并提取点目标像所在像素坐标位置(x₂，y₂)＝(371，601)；

d.移除点目标1并保持图像传感器3曝光时间不变，图像传感器3对背景成像，得到干扰图像，并将干扰图像中灰度值的最大值作为阈值，该阈值为10；

e.第a步得到的第一帧初始静态点目标图像，将灰度值小于第d步所得阈值的像素的灰度值修正为0，得到第一帧修正静态点目标图像；第c步得到的第二帧初始静态点目标图像，将灰度值小于第d步所得阈值的像素的灰度值修正为0，得到第二帧修正静态点目标图像；将第一帧修正静态点目标图像和第二帧修正静态点目标图像相叠加，叠加后新图像中两个点目标像所在行或列的所有像素灰度值相加并除以2，得到新灰度值；并将第a步得到的像素坐标位置(x₁，y₁)＝(371，539)和第c步得到的像素坐标位置(x₂，y₂)＝(371，601)连线所覆盖的像素的灰度值替换为新灰度值，得到构造点扩展函数图像，其局部图如图5所示；

或者：

将第a步得到的第一帧初始静态点目标图像和第c步得到的第二帧初始静态点目标图像相叠加，并将叠加后图像中灰度值小于第d步所得阈值2倍的像素的灰度值修正为0，得到修正叠加图像；将修正叠加图像中两个点目标像所在行或列的所有像素灰度值相加并除以2，得到新灰度值；并将第a步得到的像素坐标位置(x₁，y₁)＝(371，539)和第c步得到的像素坐标位置(x₂，y₂)＝(371，601)连线所覆盖的像素的灰度值替换为新灰度值，得到构造点扩展函数图像，其局部图如图5所示；

f.第e步得到的构造点扩展函数图像，将线状光斑所在行的整行信息提取出来，作为构造线扩展函数图像，如图6所示，该构造线扩展函数图像具有n＝1280个元素；

g.对第f步得到的构造线扩展函数图像进行离散傅里叶变换并取模，得到调制传递函数图像，该调制传递函数图像具有同第f步得到的构造线扩展函数图像相同的元素个数n＝1280，即1280个离散频谱分量，按照空间频率从小到大的顺序分别为M₀、M₁、M₂、...、M₁₂₇₉，在该顺序下，调制传递函数值第一次达到极小值所对应的调制传递函数值为M₂₁，其下脚标序号为i＝21，结合图像传感器3的像素间距l＝5.6μm，得到M₂₀和M₂₂所对应的空间频率值分别为：f_min＝(i-1)/(nl)＝(21-1)/(1280×5.6×10^-3)＝2.7902lp/mm以及f_max＝(i+1)/(nl)＝(21+1)/(1280×5.6×10^-3)＝3.0692lp/mm；

h.根据调制传递函数模型MTF(f)＝|sinc(πfd′)|，结合第g步得到的空间频率范围f_min＝2.7902lp/mm和f_max＝3.0692lp/mm，得到线状光斑长度取值范围：d_max′＝1/f_min＝nl/(i-1)＝1280×5.6×10^-3/(21-1)＝0.3584mm以及d_min′＝1/f_max＝nl/(i+1)＝1280×5.6×10^-3/(21+1)＝0.3258mm；

i.根据第b步点目标位移量d＝6.104mm和第h步得到的线状光斑长度取值范围d_min′＝0.3258mm以及d_max′＝0.3584mm，计算得到光学系统3横向放大率取值范围为：β_min＝d_min′/d＝nl/((i+1)d)＝1280×5.6×10^-3/((21+1)×6.104)＝0.0534以及β_max＝d_max′/d＝nl/((i-1)d)＝1280×5.6×10^-3/((21-1)×6.104)＝0.0587；

j.将β为变量，并以第i步得到的β_min＝0.0534和β_max＝0.0587为取值范围，按照空间频率从小到大的顺序，将第g步得到的n＝1280个调制传递函数值绘制成一条曲线，选取这条曲线上从M₀开始到首个极大值，且不包括第g步得到的M₂₁，共K个作为比对数据，这K个调制传递函数值分别是M_K1、M_K2、...、M_KK，采用遗传算法寻找到以下公式的最小值：该最小值所对应的光学系统横向放大率β即为所求，经计算，β＝0.0558。

Claims (1)

1.采用点目标像拼合的光学系统横向放大率测量方法，其特征在于所述方法步骤如下：

a.图像传感器对静态点目标第一次成像，得到第一帧初始静态点目标图像，并提取点目标像所在像素坐标位置(x₁，y₁)；

b.使点目标沿图像传感器行或列方向运动，位移量为d，之后保持点目标为静止状态；

c.保持图像传感器曝光时间不变，图像传感器对静态点目标第二次成像，得到第二帧初始静态点目标图像，并提取点目标像所在像素坐标位置(x₂，y₂)；

d.移除点目标并保持图像传感器曝光时间不变，图像传感器对背景成像，得到干扰图像，并将干扰图像中灰度值的最大值作为阈值；

e.第a步得到的第一帧初始静态点目标图像，将灰度值小于第d步所得阈值的像素的灰度值修正为0，得到第一帧修正静态点目标图像；第c步得到的第二帧初始静态点目标图像，将灰度值小于第d步所得阈值的像素的灰度值修正为0，得到第二帧修正静态点目标图像；将第一帧修正静态点目标图像和第二帧修正静态点目标图像相叠加，叠加后新图像中两个点目标像所在行或列的所有像素灰度值相加并除以2，得到新灰度值；并将第a步得到的像素坐标位置(x₁，y₁)和第c步得到的像素坐标位置(x₂，y₂)连线所覆盖的像素的灰度值替换为新灰度值，得到构造点扩展函数图像；

或者：

将第a步得到的第一帧初始静态点目标图像和第c步得到的第二帧初始静态点目标图像相叠加，并将叠加后图像中灰度值小于第d步所得阈值2倍的像素的灰度值修正为0，得到修正叠加图像；将修正叠加图像中两个点目标像所在行或列的所有像素灰度值相加并除以2，得到新灰度值；并将第a步得到的像素坐标位置(x₁，y₁)和第c步得到的像素坐标位置(x₂，y₂)连线所覆盖的像素的灰度值替换为新灰度值，得到构造点扩展函数图像；

f.第e步得到的构造点扩展函数图像，将线状光斑所在行或列的整行或整列信息提取出来，作为构造线扩展函数图像，该构造线扩展函数图像具有n个元素；

g.对第f步得到的构造线扩展函数图像进行离散傅里叶变换并取模，得到调制传递函数图像，该调制传递函数图像具有同第f步得到的构造线扩展函数图像相同的元素个数n，即n个离散频谱分量，按照空间频率从小到大的顺序分别为M₀、M₁、M₂、…、M_n-1，在该顺序下，调制传递函数值第一次达到极小值所对应的调制传递函数值为M_i，其下脚标序号为i，结合图像传感器的像素间距l，得到M_i-1和M_i+1所对应的空间频率值分别为：f_min=(i-1)/(nl)以及f_max=(i+1)/(nl)；

h.根据调制传递函数模型MTF(f)=|sinc(πfd')|，结合第g步得到的空间频率范围f_min和f_max，得到线状光斑长度取值范围：d_max’=1/f_min=nl/(i-1)以及d_min’=1/f_max＝nl/(i+1)；

i.根据第b步点目标位移量d和第h步得到的线状光斑长度取值范围，计算得到光学系统横向放大率取值范围为：β_min=d_min'/d=nl/((i+1)d)以及β_max=d_max'/d＝nl/((i-1)d)；

j.将β为变量，并以第i步得到的β_min和β_max为取值范围，在第g步得到的n个调制传递函数值中选取K个作为对比数据，这K个调制传递函数值分别是M_K1、M_K2、…、M_KK，采用遗传算法寻找到以下公式的最小值：该最小值所对应的光学系统横向放大率β即为所求。