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

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

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

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

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

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

发明内容

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

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

采用线光源的光学***横向放大率测量方法，步骤如下：

a、在物方放置长度为d的线光源，方向与图像传感器的行或列方向平行；

b、图像传感器对线光源成像，得到初始点扩展函数图像；保持图像传感器曝光时间不变，移除线光源，图像传感器对背景成像，得到干扰图像，并将干扰图像中灰度值的最大值作为阈值；

c.将第b步得到的初始点扩展函数图像中，线光源像所在行或列的整行或整列信息提取出来，作为初始线扩展函数图像，该初始线扩展函数图像具有n个元素；并将这n个元素中灰度值小于第b步所得阈值的像素的灰度值修正为0，得到修正线扩展函数图像，该修正线扩展函数图象具有n个元素；

或者：

将第b步得到的初始点扩展函数图像中，灰度值小于第b步所得阈值的像素的灰度值修正为0，作为修正点扩展函数图像；并将修正点扩展函数图像中，线光源像所在行或列的整行或整列信息提取出来，得到修正线扩展函数图像，该修正线扩展函数图象具有n个元素；

d.对第c步得到的修正线扩展函数图像进行离散傅里叶变换并取模，得到调制传递函数图像，该调制传递函数图像具有同第c步得到的修正线扩展函数图像相同的元素个数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)；

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

f.根据第a步线光源长度d和第e步得到的线光源像长度取值范围，计算得到光学***横向放大率取值范围为：β_min＝d_min′/d＝nl/((i+1)d)以及β_max＝d_max′/d＝nl/((i-1)d)；

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

采用线光源的光学***横向放大率测量装置，包括线光源、光学***、图像传感器，所述的线光源经过光学***成像到图像传感器表面，并且，在该装置光轴方向与图像传感器行或列方向所确定的平面内，线光源呈弯曲状，且所述的线光源上任意位置都准焦成像到图像传感器表面。

本发明的有益效果是：

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个离散频谱分量，按照空间频率从小到大的顺序分别为M₀、M₁、M₂、…、M₁₂₇₉，在该顺序下，调制传递函数值第一次达到极小值所对应的调制传递函数值为M₄₂，其下脚标序号为i＝42，结合图像传感器3的像素间距l＝5.6μm，得到M₄₁和M₄₃所对应的空间频率值为：f_min＝(i-1)/(nl)＝(42-1)/(1280×5.6×10^-3)＝5.71991p/mm以及f_max＝(i+1)/(nl)＝(42+1)/(1280×5.6×10^-3)＝5.99891p/mm；

e.根据调制传递函数模型MTF(f)＝|sinc(πfd′)|，结合第d步得到的空间频率范围f_min＝5.71991p/mm和f_max＝5.99891p/mm，得到线光源像长度取值范围：d_max′＝1/f_min＝nl/(i-1)＝1280×5.6×10^-3/(42-1)＝0.1748mm以及d_min′＝1/f_max＝nl/(i+1)＝1280×5.6×10^-3/(42+1)＝0.1667mm；

f.根据第a步线光源1长度d＝3mm和第e步得到的线光源像长度取值范围d_min′＝0.1667mm以及d_max′＝0.1748mm，计算得到光学***横向放大率取值范围为：β_min＝d_min′/d＝nl/((i+1)d)＝1280×5.6×10^-3/((42+1)×3)＝0.0556以及β_max＝d_max′/d＝nl/((i-1)d)＝1280×5.6×10^-3/((42-1)×3)＝0.0583；

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

Claims (1)

1.采用线光源的光学***横向放大率测量方法，其特征在于所述方法步骤如下：

a、在物方放置长度为d的线光源，方向与图像传感器的行或列方向平行；

b、图像传感器对线光源成像，得到初始点扩展函数图像；保持图像传感器曝光时间不变，移除线光源，图像传感器对背景成像，得到干扰图像，并将干扰图像中灰度值的最大值作为阈值；

c.将第b步得到的初始点扩展函数图像中，线光源像所在行或列的整行或整列信息提取出来，作为初始线扩展函数图像，该初始线扩展函数图像具有n个元素；并将这n个元素中灰度值小于第b步所得阈值的像素的灰度值修正为0，得到修正线扩展函数图像，该修正线扩展函数图象具有n个元素；

或者：

将第b步得到的初始点扩展函数图像中，灰度值小于第b步所得阈值的像素的灰度值修正为0，作为修正点扩展函数图像；并将修正点扩展函数图像中，线光源像所在行或列的整行或整列信息提取出来，得到修正线扩展函数图像，该修正线扩展函数图象具有n个元素；

d.对第c步得到的修正线扩展函数图像进行离散傅里叶变换并取模，得到调制传递函数图像，该调制传递函数图像具有同第c步得到的修正线扩展函数图像相同的元素个数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)；

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

f.根据第a步线光源长度d和第e步得到的线光源像长度取值范围，计算得到光学***横向放大率取值范围为：β_min＝d_min′/d＝nl/((i+1)d)以及β_max＝d_max′/d＝nl/((i-1)d)；

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