一种应用于数字全息重构过程的对焦方法
【技术领域】
[0001] 本发明设及一种对焦方法，尤其设及一种应用于数字全息重构过程的对焦方法。
【背景技术】
[0002] 数字全息术是基于光学全息术的基本原理，采用数字图像传感器(如:CCD和CMOS) 代替传统需要化学显影的光学全息板来记录全息图，将记录的强度信息转变为数字信号存 入计算机，再利用计算机模拟参考光对全息图进行衍射重构成像，从而实现数字重构的一 种新的全息术(详见参考专利:江竹青，王酷，陈依菲，等.数字全息成像在线重构显示系统 及方法：，CN103186090A[P].2013.)，现在已有基于数字全息的数字全息显微镜和数字全息 测量仪器设备。在数字重构的过程中需要根据实际物体与数字图像传感器之间的距离，通 过模拟光学衍射来实现对物体的再现重构，此时即实现了对焦。如果偏离了运个距离，重构 出来的像就会离焦变模糊。
[0003]典型的对焦方法有四大类:基于光场强度赌的方法，基于光场强度方差的方法，基 于光场强度梯度的方法和基于光场强度空间频率的方法(详见参考专利:蒋宏.空间频率对 比法作为成像光学系统的自动对焦判据，CN1303025A[P].2001参考文献:WangH，QinA， ZhaoB.Autofocusmethodforamplitudeandpurephaseobjectsreconstruction indigitalhologr曰phy[C]//Intern曰tion曰1ConferenceonOpticalInstrument曰tion andTechnology.InternationalSocietyforOpticsandPhotonics,2009:751312-75131Z-5.)。它们最初为相机成像系统对焦设计产生的，也可W应用于数字全息重构过程 的对焦，运些典型的方法存在对不同类型物体的适应能力相对较弱和计算量相对较大等方 面的不足，不太适合移动设备或者实时对焦处理等对计算要求严格的情况。另外还有多波 长的数字全息对焦方法（详见参考文献:GaoP,化0B,RuppR,etal.Autofocusingbased onwavelengthdependenceofdiffractionintwo-wavelengthdigitalholographic microscopy[J].Opticsletters，2012,37(7): 1172-1174.)，利用不同波长在衍射时会存 在差异，但在对焦面上多种波长差异最小的原理来判断是否对焦，运种方法需要多种不同 波长的激光，相比于普通的单波长数字全息术，会增加额外的硬件设备(如几台不同波长的 激光器），增加设备成本和实际操作难度。
【发明内容】
[0004]本发明是针对现有技术的不足，提出一种应用于数字全息重构过程的对焦方法。
[0005]本发明的目的是通过W下技术方案实现的，一种应用于数字全息重构过程的对焦 方法，该方法包括W下步骤：
[0006]步骤(1)采用数字全息方法得到全息图后，处理全息图得到物光场矩阵Uo;
[0007]步骤（2)用肉眼或者用常见长度测量工具粗略测量后估算出一个对焦距离区间 [Zl，Zk]，确保对焦距离在该区间内，此处对焦距离为成像对象与数字图像传感器的距离。设 置区间采样数为k，计算出采样间隔距湾
，每隔距离delta选取一个距离，得 I^z2=zi+delta，z3 =zi+2delta，z4 =zi+3delta，-..，zk=zi+(k-l)*delta;
[000引步骤(3)对反向传播距离Zl，Z2，Z3，...，Zk依次采用光学衍射算法得到光场矩阵化， 化，化，···，化；
[0009]步骤(4)在每个光场矩阵Ux中选相同位置需要对焦的目标区域，得到新矩阵ROIx， 即对焦窗口，矩阵ROIx大小为m*n，其中，x=l，2,…，k;
[0010] 步骤(5)通过矩阵ROIx计算实数R(zx)，公式如下；
[0011]
[001^ 其中，R0Ix(i，j)为矩阵ROIx第i行、觀列的元素，
[0013] 或者通过公式如下：
[0014]
[001引其中，R0Ix(i，j)为矩阵ROIx第i行、觀列的元素；
[0016] 步骤(6)计算数列y(zi)，y(Z2)，…，y(zk-i)，公式如下：
[0017] y(zx)= |R(Zx)-R(Zx+i) |2,(x= 1,2, . . . ,k-l),
[0018]步骤(7)对数列y(zi)，y(Z2)，…，y(zk-i)，直接遍历数列y的所有项或者采用对焦评 价函数捜索算法找出数列y的最小值y(zf)，最小值y(zf)所对应的Z即为对焦面与数字图像 传感器的距离Zf，距离数字图像传感器Zf处且垂直光传播方向的平面即为对焦面。
[0019] 本发明的明显效益在于:通过光学衍射算法在数字全息重构过程中重构出不同距 离的光场，对相邻的光场的振幅或者强度作差平方求和，找出最小值的方法，实现了在数字 全息仪器装备的重构过程中的准确、快速对焦，从而让数字全息仪器获取精确的物体位置 和清晰的图像，并且同时没有增加任何的硬件设备。
【附图说明】
[0020] 图1是具体实施例中的离轴数字全息的装置结构示意图。1为平行激光光束，2为平 面反射镜，3为分束镜，4为平面反射镜，5为用来成像的对象物体，6为分束镜，7为数字图像 传感器(在【具体实施方式】中为CCD)，，8为计算机，9为数字图像传感器到实际物体的距离，也 是对焦距离Zf。
[0021] 图2是本发明一种应用于数字全息重构过程的对焦方法的流程图
[0022] 图3是具体实施例中采用说明书
【发明内容】
里步骤(5)公式（1)后本发明得到的关于 y与Z的函数曲线图
[0023] 图4是具体实施例中采用说明书
【发明内容】
里步骤(5)公式(2)后本发明得到的关于 y与Z的函数曲线图
[0024] 图5是具体实施例对采集到的全息图衍射重构出来的，在对焦面出的物体强度图
【具体实施方式】
[0025] 下面结合附图详细说明本发明的数字全息重构过程的对焦方法的典型实施例。应 当说明的是，此处所描述的具体实施例仅用于解释本发明，并不用于限定本发明。
[0026] 图1为离轴数字全息的装置结构图。平面波激光束1由一个单波长的激光光源提 供，经过分束镜3，分为两束激光，一束经过反射镜2反射后射透过物体5，再通过分束镜6出 射到CCD上；另一束激光通过分束镜4W小倾斜角通过分束镜6入射到CCD上，两束激光在CCD 上发生干设，CCD记录下干设图（即全息图），CCD分辨率为2750巧200，并把数字化后的干设 图输入计算机。计算机经过傅里叶变换滤波处理滤除0级直流项和-1级的李生共辆项后得 到物光光场矩阵化，大小为2200巧750个元素。
[0027] 用直尺粗略测量后估算出聚焦区间[181mm，187mm]。设置区间采样数为201，故采 样间隔delta= 0.03mm，可W得到zi= 181mm，Z2 =zi+delta= 181.03mm，Z3 =zi+2delta= 181.06mm，Z4=zi+3delta= 181.09mm，...，Z2〇i=zi+(201-1 )*del1:a= 187.00mm。
[002引对反向传播距离Ζ1，Ζ2，Ζ3，···，Ζ201依次采用光学衍射算法的典型算法角谱传播算 法反向衍射得到光场矩阵化，化，化，…，化01;
[0029] 在每个光场矩阵Ux中选相同位置需要对焦的目标区域，得到新矩阵ROIx，即对焦窗 口，矩阵ROIx大小为64*64,其中，x=l，2,…，k;
[0030] 通过矩阵ROIx计算实数R(zx)，公式如下；
[0031]
[00创其中，R01x(i，j)为矩阵ROIx第i行、第j列的元素，
[0033] 或者通过公式如下：
[0034]
[00对其中，R01x(i，j)为矩阵ROIx第i行、觀列的元素；
[0036] 通过公式计算得到数列y(zi)，y(Z2)，…，y(Z2〇())，公式如下：
[0037] y(zx)=IR(zx)-R(zx+i)P,(x=l,2,... ,200),
[0038] 对实数数列y(zi)，y(Z2)，…，y(Z2〇o)，采用对焦评价函数捜索算法里的典型算法 Fibonacci捜索法找出数列y的最小值y(zo)，最小值y(zo)所对应的Z即为对焦面与CCD的距 离zf= 184.1mm，从而确定距离数字图像传感器zf= 184.1mm处且垂直光传播方向的平面即 为对焦面。
[0039] 关于y与Z的函数曲线图如图3和图4所示，其中图3为采用说明书
【发明内容】
步骤(4) 公式（1)得到的结果，图4为采用说明书
【发明内容】
步骤(4)公式(2)得到的结果，在所测得的 对焦面处重构的强度图如图5所示。
[0040] 表1为本方法和其他四种主流方法(基于光场强度赌的方法，基于光场强度方差的 方法，基于光场强度梯度的方法和基于光场强度傅里叶频谱权重的方法)在分别选取100 个、1000个等间隔点和64*64像素、512巧12像素、2048*2048像素对焦窗口条大小件下对焦 耗时对比，可W看出本发明在同等条件下，耗时要明显小于其他四种主流方法。
[0041] 表 1
[0042]
[0043]尽管特定参考实施例详细的描述了本发明，但在此描述的本发明实施例的意图不 是详尽或者局限于所公开的具体形式。相反，所选用的用于说明问题的实施例是为了使本 技术领域内的技术人员实施本发明而选择的。在不脱离本发明的权利要求所描述和限定的 本发明的实质范围的情况下，存在变型例或修改例。
【主权项】
1. 一种应用于数字全息重构过程的对焦方法，其特征在于，该方法包括以下步骤： 步骤(1)采用数字全息方法得到全息图后，处理全息图得到物光场矩阵u〇; 步骤（2)用肉眼或者用常见长度测量工具粗略测量后估算出一个对焦距离区间[Z1， zk]，确保对焦距离在该区间内，此处对焦距离为成像对象与数字图像传感器的距离。设置 区间采样数为k，计算出采样间隔距离，每隔采样间隔距离delta选取一个反 向传播距离，依次得到22 = 21+(161七&，23 = 21 + 2(161七&，24=21 + 3(161七3，."，21{ = 21+(1<-1)* delta； 步骤(3)对反向传播距离21，22,23，一，&依次采用光学衍射算法反向衍射得到光场矩阵 Ul，U2，U3，···，Uk; 步骤(4)在每个光场矩阵Ux中选相同位置需要对焦的目标区域，得到新矩阵R0Ix，即对 焦窗口，矩阵R〇Ix大小为m*n，其中，x=l，2,…，k; 步骤(5)通过矩阵R0Ix计算实数R(zx)，公式如下；其中，R〇Ix(i，j)为矩阵R〇Ix第i行、第j列的元素， 或者通过公式如下：其中，R〇Ix(i，j)为矩阵R〇Ix第i行、第j列的元素； 步骤(6)计算数列y(zi)，y(Z2)，…，y(zk-i)，公式如下： y(zx)= |R(zx)-R(zx+i) |2,(x=1,2, . . . ,k-l), 步骤(7)对数列y(zi) ,y(Z2),…，y(zk-i)进行处理找到对焦面。2. 根据权利要求1所述的方法，其特征在于，所述步骤(7)找到对焦面，具体如下： 直接遍历数列y的所有项或者采用对焦评价函数搜索算法找出数列y的最小值y(Zf)，最 小值y(Zf)所对应的Zf即为数字图像传感器到成像目标物体的距离，距离数字图像传感器z f 处且垂直光传播方向的平面，即为对焦面。
【专利摘要】本发明公开了一种应用于数字全息重构过程的对焦方法。该方法是在数字全息计算处理得到物光场后，物光场反向轴向传播，对两相邻截面光场强度(或振幅)的同一对焦区域矩阵元素作差，然后对所得到的矩阵的元素平方，再对平方后的矩阵所有元素求和得到一个函数值，依次对相邻等间距的轴向截面图像作上述处理即可以得到一组函数值，找到这组函数值的最小值，其对应的距离处即为对焦面的位置。该方法可应用于数字全息仪器装置的重构过程中的对焦，具有计算简单，适应范围广，对焦准确，不需要额外增加硬件装置等优点。
【IPC分类】G03H1/22
【公开号】CN105446111
【申请号】CN201610005898
【发明人】吕萌, 司徒国海
【申请人】中国科学院上海光学精密机械研究所
【公开日】2016年3月30日
【申请日】2016年1月6日