引言

大口径光学设备主要用于外场远距离目标的监视和测量，其口径大、焦深短、调焦量程大，精准对焦相对较难。在使用中，常常会遇到目标与光学设备之间的距离未知、传感器靶面偏离理想像面位置较大、图像因离焦而模糊不清问题。针对大口径光学设备的离焦问题，传统的手动对焦方法依靠人眼确定对焦限度，精度低，很难达到图像质量最佳；目前通用的计算机视觉调焦方法^[1-2]需要在整个量程内来回调整调焦机构寻找最佳点，对于调焦量程大的大口径光学设备，调焦速度慢，调整过程中出现模糊图像，难以做到实时动态运行。随着一些图像复原方法的出现，建立光学设备的退化模型，利用传统图像复原方法^[3-4]对图像进行逆卷积处理，可以将模糊图像还原为清晰图像。目前，国内外还没有成像设备自我实时测量离焦量大小和方向的方法^[3]。外场大口径光学设备使用中，无法实时测量离焦量的大小和方向，因此，直接建立准确的离焦退化模型非常困难。

本文从光学镜头的成像原理出发，按照焦点附近光场变化规律和焦深变化范围，以一定的周期间隔选择某些假想的离焦量值，推算相应的离焦点扩散函数、相干光传递函数和非相干光传递函数，提出图像离焦程度判断方法，运用爬山法^[2]、逆滤波^[3-7]和离焦程度判断法，寻找大口径光学设备拍摄目标时刻的最佳离焦量值和方向，对模糊图像给予复原。依据计算得到的离焦量大小和方向，不需要在整个量程内来回调整调焦机构寻找最佳点，能很快将调焦机构调整到最佳位置点，调整过程中设备一直显示清晰图像，可以实现大口径光学设备的实时精准对焦。

1   光学设备的成像模型

无穷远距离目标的光源以发散光束经过光学设备镜头之后，入射到光学设备前，形成汇聚光波^[2]，在焦面上形成了点像。光波传播与成像过程见图1所示，从物到像经过了4个过程，物面$\Sigma 1 $产生的物光波，经自由空间传播到达光学系统前主面$\Sigma $2，光学系统的光瞳对光束限制，再经透镜传播之后到达后主面$\Sigma $3，再经自由传播到达实际像面$\Sigma $5，形成像面。

图  1  物平面光波经透镜到像的传播过程方框图

Figure  1.  Block diagram of light propagation from object to its image through the the lens

下载: 全尺寸图片 幻灯片

光学镜头的离焦成像的光路图见图2所示，图中：$\Sigma $1为物平面；$\Sigma $2为透镜前主面；$\Sigma $3为透镜后主面；$\Sigma $4为理想高斯像平面；$\Sigma $5为实际像面；$u$为物距；$v$为理想像距（即高斯像距）；${d_i}$为实际像距；$\delta $为平面$\Sigma $4与平面$\Sigma $5之间的距离，$\delta $定义为离焦量，一般$\delta $$\ll $$v$；$f$为镜头组等效焦距；$D$为光学系统的入射光瞳口径。

图  2  光传播与成像示意图

Figure  2.  Schematic diagram of light propagation and imaging

下载: 全尺寸图片 幻灯片

根据离焦量定义，有：

于是，由（1）式得：

物面$ \Sigma $1的物体点分布光波用函数${u_o}(\xi ,\eta )$表示，ξ为横坐标，η为纵坐标；光学系统的光瞳可通光瞳函数$p(\xi ',\eta ')$表示，ξ′为横坐标，η′为纵坐标；实际像平面$ \Sigma$5上形成的像点分布函数为Ψ_i(x, y)，x为横坐标，y为纵坐标。

2   离焦退化模型推导

光学设备的离焦退化模型就是光学成像镜头的点扩散函数和传递函数，成像波段为可见光，光学系统等效为透镜。

2.1   离焦退化模型的点扩散函数

设物体的平均光波波长为λ，光波的坡印廷矢量k=$\dfrac{{2 {\text{π}} }}{\lambda }$，像平面空间坐标为（x，y）。结合图2的成像过程，依据透镜的位相变换和夫琅禾费衍射理论公式^[8]可得：

式中j为虚宗量。

根据高斯成像公式，理想像的垂直放大倍率为

像平面（x，y）处的场只与物平面上对应点近旁一个小区域有关，这时exp因子的宗量变化不大于几分之一弧度，因而可采用（5）式近似：

将（1）、（2）、（4）、（5）式代入（3）式，变量置换Mη、M ξ→η₀、ξ₀，整理后得：

其中：积分区域$\Sigma $₃为直径为D的圆面；$\Sigma $₁表示目标区域；$\Sigma $₄表示理想像区域(高斯像面)。令

作变量代换之后，（7）式可变形为

式中$h({\xi _0} - x,{\eta _0} - y)$为光学设备的离焦退化模型的点扩散函数$h\left( {x,y} \right)$的卷积表示，是实际像函数${\psi _i}(x,y)$与理想像函数${u_i}(x,y)$之间关系，表示从高斯像面理想像函数到实际物平面像函数的点扩散函数。令E₀为光波的幅度，E₀取为

则有：

引入新坐标，令

其中$r = \sqrt {{x^2} + {y^2}} $，表示像平面点到像面中心（光轴）的距离。代入（9）式并归一化得：

（10）式为推导得出的光学成像镜头离焦退化图像的点扩散函数形式，式中：D为光学设备入射孔径；v为理想像距；δ 为离焦量； λ为光的平均波长。对于物体上的点光源，在实际像面的光幅度分布函数为可以看出，其点像扩散函数具有圆对称性。（10）式中：

${J_0}\left( x \right)$为0阶贝塞尔函数^[9]，x为自变量。令光轴向变量 ω₁为

令${I_0} = {\left| {{E_0}} \right|^2}$，I₀表示光的强度。当像面点与中心的距离r=0时，光轴上的光强为

（15）式中，当$ - 3.2 \leqslant {\omega _1} \leqslant 3.2$时，光强I/I₀下降20%，于是焦深ε^[2]有：

离焦量不超过焦深（即$\delta \leqslant \varepsilon $）时，不影响像面清晰度；超出焦深（即$\delta \geqslant \varepsilon $）时就会影响像面的清晰度。由于点扩散函数具有圆对称性，因此我们取某一方向进行计算。设某一方向空间量化尺寸（即像元尺寸）为b，当前位置点的像元数为N。

令图像径向变量ω₀为

则(12)式、（13）式变为

当离焦量δ=0时，焦面上的点像衍射图像为

式中${J_1}\left( x \right)$为1阶贝塞尔函数^[2]。物体上的点光源，在像面的点像幅度扩散函数为

物体上的点光源，在像面的点像强度扩散函数为

2.2   离焦退化模型的传递函数

2.2.1   相干光照明传递函数

取像面空间坐标$ x$、$y$相对应的空间频谱为${F_x}$、${F_{\rm{y}}}$，根据传递函数定义^[10]，传递函数为点扩散函数的傅里叶变换：

将$r = \sqrt {{x^2} + {y^2}} $，$\rho = \sqrt {{F_x}^2 + {F_y}^2} $和（10）式代入（24）式得：

式中：

下标re表示实部，下标im 表示虚部。

2.2.2   非相干光照明下传递函数

存在波像差的非相干光照明下传递函数^[10]为

式中$\Delta \omega (x,y)$为光学成像系统的波像差，根据像差光学系统的频谱分析^[5]，离焦的波像差近似计算为^[5]

式中：${\rho _0} = \dfrac{D}{{\lambda \cdot f}}$，ρ₀为离焦量δ为0的光学系统的截止频率；x、y为积分变量；D为光学系统入射孔径；f 为光学镜头焦距。（27）式代入（26）式得：

式中$\alpha ({F_x},{F_y})$是两个错开$(\lambda {d_i}\rho )$的光学系统入射光瞳重合部分，图3中的阴影部分为积分区间。取${\rm tg}i = - {{{{F_y}} / F}_x}$后，将x、y坐标系顺时针旋转i角度到坐标系ζ、η。

图  3  积分区间$\Sigma $示意图

Figure  3.  Sketch map of integral interval $\Sigma $

下载: 全尺寸图片 幻灯片

由（29）式得：

令$\alpha ' = {\text{π}} {\left( {\dfrac{D}{f}} \right)^2}\left( {\dfrac{\delta }{\lambda }} \right)\dfrac{{{F_x}}}{{{\rho _0}}}$，$\; \beta ' = {\text{π}} {\left( {\dfrac{D}{f}} \right)^2}\left( {\dfrac{\delta }{\lambda }} \right)\dfrac{{{F_y}}}{{{\rho _0}}}$，$\; \rho = $$ \sqrt {F_x^2 + F_y^2} $，$\gamma = \dfrac{\rho }{{{\rho _0}}}$，$\alpha = \dfrac{1}{2}{\omega _1}\gamma $，${\theta _0} = \arccos \left( {\dfrac{\rho }{{{\rho _0}}}} \right)$。ω₁为光轴向变量，${\omega _1} = \dfrac{{{\text{π}} \cdot {D^2} \cdot \delta }}{{2{f^2} \cdot \lambda }}$（对于无穷远目标像距v≈f），ρ为像面空间的频率，ρ<ρ₀。于是有：$\alpha '\cos i + $$\alpha '\cos i + \beta '\sin i = 0,\; \alpha '\cos i - \beta '\sin i = {\text{π}} {\left( {\dfrac{D}{f}} \right)^2}\left( {\dfrac{\delta }{\lambda }} \right)\dfrac{\rho }{{{\rho _0}}}$。

将（30）式代入（29）式，推导得：

推导（31）式，经过积分计算得光学系统的离焦退化模型的传递函数：

当γ为0时，H（0，ω₁）=1；当离焦量δ为0时，ω₁=0，（32）式为

光学系统成像要求H(γ)≥0，H(γ)=0空间频率成份的光将被截止。（32）式方程H(γ)=0存在多个γ解，第一个解（即最小值解）用γ_max表示，截止时γ_max相对应的空间频率ρ_max为截止频率，用ρ_max表示。H(γ)<0为传函反相，图像上表现为“黑白颠倒”。

可以证明：ω₁=0时，γ_max=1；ω₁越大，γ_max越小；在γ<γ_max情况下，（32）式中正弦函数中ω₁·γ·(γ−x)<4。故（32）式中的正弦函数泰勒近似展开，积分得：

其中：

3   离焦模糊图像的复原方法

计算图像清晰度，以此作为评判图像的依据，找到图像清晰度的极大值点，把模糊图像复原。对于给定焦距和入射孔径的摄影物镜，采用迭代逼近法寻找离焦衍射模糊图像的离焦量。

3.1   离焦图像的维纳逆滤波法复原算法

光学设备成像过程中，因为存在离焦导致图像退化，图像退化过程^[11-16]可用下式表示：

（35）式中，g(x,y,δ)为模糊图象，f(x,y)为理想图像，h(x,y,δ)为点扩散函数，n(x,y)为噪声图像。对（35）式两边取傅里叶变化得：

(36)式中，$H({F_x},{F_y},\delta )$为传递函数，$G\left( {{F_x},{F_y}} \right)$为模糊图像的傅里叶频谱，$N({F_x},{F_y})$为噪声的傅里叶频谱，$U({F_x},{F_y})$为理想图像的傅里叶频谱，F_x和F_y为空间频率。将（36）式转换为矩阵形式为

（37）式中[N]、[G]、[F]、[H]为$N({F_x},{F_y})$、$G\left( {{F_x},{F_y}} \right)$、$U({F_x},{F_y})$、$H({F_x},{F_y})$的矩阵形式。利用最小二乘法^[5],求得[$U$]的最小二乘估计^[14-17][$\hat U$]为

取傅里叶逆变换得复原后的理想图像矩阵为

3.2   图像的离焦程度分析模型

已知光学设备成像时刻的离焦量，建立离焦传递函数，由（39）式，可以得到复原后的原始理想图像。然而，光学设备拍摄时刻的离焦量无法实时测得，我们假想一系列离焦量值，用这些值估计真实离焦量，进行图像复原，会取得一系列复原图像，计算这些复原图像的清晰度值，清晰度值反映了假想离焦量与真实离焦量的逼近程度，也就是离焦程度。将一系列复原图像中，清晰度值最大的图相对应的估计离焦量值，我们称之为“最佳离焦量值”。

1） 清晰度定义

点像的强度分布曲线越陡，光学系统对于一般的目标物体所成的像边缘越清晰，依据光学系统锐度的定义^[14-16]，具有离焦量δ的实际像面点像的清晰度表示为

式中：a为传感器靶面的最大尺寸；r为空间水平或垂直方向的坐标变量；ρ为水平或垂直方向的空间频率；ρ_max为截止频率。

2） 离焦量的估计值与真实值之差越小，复原图像清晰度越大

根据（36）式和（40）式估计复原图像的清晰度：

ρ_max对应x和y方向的最大截止空间频率为${F_{x0}}$和${F_{y0}}$,代入（41）式，离散化为

空间频率高的部分反映图像的细节，选择高频F_{x_high}、F_{y_high}成份分析清晰度。考虑x方向，取实部

同理，考虑y方向，取实部

式中：${u_x}$、${u_y}$为常数；${K^2}({F_{y\_\rm high}})$为${F_{y\_\rm high}}$的多项式；${K^2}({F_{x\_\rm high}})$为${F_{x\_\rm high}}$的多项式；${U^2}({F_{x\_\rm high}})$和${U^2}({F_{y\_\rm high}})$为理想图像的高频强度值，计算可得：当$\delta - \delta ' = 0$时，（44）式中${A_c}({F_{x\_\rm high}},\delta ')$取最大值，表示复原图像的清晰度为最大值。

4   基于离焦复原的视频自动调焦方法

4.1   光学设备的调焦范围

外场光学设备用于观测遥远的物体，在满足焦深条件下，当物体物距u从无穷远u_∞到一定距离u₊内变化时，焦面的焦深范围内就为不离焦像面，不需要调整。${u_\infty } \leftrightarrow {u_ + }$视为无穷远成像景深距离，此时，像距v=焦距f。

如图4所示，当目标离开无穷远，在近距离u₊与u变化时，即u₊↔u之间，像面离开传感器靶面位置，超出焦深,需要调整传感器像面位置，使此接近理想相面位置v，像面位置变化为：$f + \varepsilon \leftrightarrow $$ (1 + \dfrac{f}{{u - f}})f - \varepsilon $。ε为焦深，f 为焦距，$\Delta z$为调焦量程。

图  4  外场光学设备调焦范围示意图

Figure  4.  Diagram of optical equipment focusing field range

下载: 全尺寸图片 幻灯片

于是，像面调整量程$\Delta z$为

大口径光学设备外场使用时，物距远大于焦距，即u>>f，取u=N·f，N为一个大整数，对于大口径光学设备N>>600。于是有：

从式（46）可以看出：大口径设备的焦距长而焦深短，调焦步进电机的调整量程大。

4.2   图像离焦程度的爬山法判断方法

一幅图像是否存在离焦，需要通过图像处理判断，应用爬山法^[2-5]判断图像离焦程度。

4.2.1   概略离焦量值寻找

以2倍的焦深为步长，寻找最大可能离焦量值，也就是清晰度值最大时对应的假设离焦量，即概略离焦量值，过程如下：

1） 取一些假想离焦量值${\omega _1}[{i}] = \left\{ {0,2{\text{π}} ,{\rm{4}}{\text{π}} ,{\rm{6}}{\text{π}} , \cdots , }\right.$$ \left. {2i{\text{π}} , \cdots ,{\rm{20}}{\text{π}} } \right\}$，找到这幅图像的像元尺寸a。结合给定的设备的光学D、焦距f参数，对照这幅图像的像元尺寸a，代入相干传递函数计算式（25），分别得出11个不同离焦量的相干传递函数H_c(ρ)，ρ为空间频率，并保存于计算机系统内存中。

2） 取实时图像u(i,j)，作傅立叶变换变换，得到傅立叶谱U(F_x,F_y)。

3） 将H_c(ρ)和U(ρ)代入式（42）、（43）、（44），计算得出11个A(ω₁(i))。找出11个A(ω₁(i))值中的最大值A_max(ω₁(i)),对应的2iπ定义为2i₀π，由2i₀π计算出离焦量δ₀。如图5所示。

图  5  图像离焦程度的自动判断方框图

Figure  5.  Block diagram of automatic judgement for image defocus degree

下载: 全尺寸图片 幻灯片

4.2.2   精细离焦量值寻找

以小于1倍的焦深为步长，寻找最大可能离焦量值。此时，找到的清晰度值最大时对应的假想离焦量，即精细离焦量值。

1） 计算邻域离焦量清晰度值

以2i₀π±π分别为假想离焦量值，代入（43）和（44）式，分别计算$A(2{i_0}{\text{π}} + {\text{π}} )$和$A(2{i_0}{\text{π}} - {\text{π}} )$两个清晰度值。

2） 比较邻域离焦量清晰度值

比较$A(2{i_0}{\text{π}} + {\text{π}} )$和$A(2{i_0}{\text{π}} - {\text{π}} )$的大小。若$A(2{i_0}{\text{π}} + $$ {\text{π}} )$>$A(2{i_0}{\text{π}} - {\text{π}} )$，则在2i₀π～2i₀π+π之间选择7个点，计算$A(2{i_0}{\text{π}} + {\text{π}} /8)$、$A(2{i_0}{\text{π}} + {\text{π}} /4)$、$A(2{i_0}{\text{π}} + 3{\text{π}} /8)$、$A(2{i_0}{\text{π}} + $$ {\text{π}} /2)$、$A(2{i_0}{\text{π}} + 5{\text{π}} /8)$、$A(2{i_0}{\text{π}} + 3{\text{π}} /4)$、$A(2{i_0}{\text{π}} + $$7{\text{π}} /8) $中的最大值，找到对应的离焦量值；若$A(2{i_0}{\text{π}} + {\text{π}} )$<$A(2{i_0}{\text{π}} - {\text{π}} )$，则在2i₀π～2i₀π−π之间选择7个点，计算$A(2{i_0}{\text{π}} - {\text{π}} /8)$、$A(2{i_0}{\text{π}} - {\text{π}} /4)$、$A(2{i_0}{\text{π}} - 3{\text{π}} /8)$、$A(2{i_0}{\text{π}} - $${\text{π}} /2) $、$A(2{i_0}{\text{π}} - {\text{π}} /8)$、$A(2{i_0}{\text{π}} - 3{\text{π}} /4)$、$A(2{i_0}{\text{π}} - $$7{\text{π}} /8) $中的最大值，找出对应的离焦量值。

3） 若清晰度最大值对应的离焦量δ₀≤Δz（焦深），则图像不存在离焦问题。

4.3   离焦方向的判断方法

从光轴上光场分布（12）式可以得出，随着离焦量ω₁的变化，存在多个极大值点，ω₁=0处为0级极大点，ω₁=±6π处为1级极大点，级数越大，峰值越小。因此，在当前的像面位置的邻域寻找到的极大值点，不一定是最清晰位置。

为了兼顾实时性、调焦精度，我们利用图像复原信息判断图像离焦的方向，如图6所示，针对具有一定尺度的目标离焦模糊图像，依据图像复原处理前后的散斑面积大小变化来分析，复原后，目标区域面积变大，说明为负离焦，向远离镜头的方向运动；复原后，目标区域面积变小，说明为正离焦，向靠近镜头的方向运动。

5   实验结果

以某光学成像设备的外场实验图片复原为例进行实验，说明复原过程成像效果。

5.1   输入参数

选择光学摄像机的镜头口径D=27.5 mm，焦距为f=143 mm，感光波段为可见光（平均波长为λ=0.55 μm），单像元尺度b=10 μm，像距v≈f=143 mm。将已知参数代入（18）式，得ω₀=3.5π；代入（16）式，得焦深ε=±55λ=±30.29 μm；代入（27）式，得衍射极限频率ρ₀=349 lp/mm。

5.2   计算原图的傅里叶变化图

光学摄像机在离焦状态下拍摄的图片，如图7（a）所示，取像元数为400×400，作傅里叶变换，得到傅里叶频谱图见图7（b）所示。

图  7  原图与傅里叶频谱图

Figure  7.  Original image and its Fourier spectrum image

下载: 全尺寸图片 幻灯片

5.3   计算光学系统不同离焦量的点扩散函数曲线

选择不同的假想离焦量，代入（10）式，计算光学系统不同离焦量的点扩散函数曲线。如图8所示：（a）为离焦量δ=0时的焦面上的点像衍射图像的径向分布曲线；（b）为离焦量δ=ε（1倍焦深），ω₁=π时的焦面上的点像衍射图像的径向分布曲线；（c）为离焦量δ=2ε（2倍焦深），当ω₁=2π时的焦面上的点像衍射图像的径向分布曲线；（d）为离焦量δ=4ε（4倍焦深），ω₁=4π时的焦面上的点像衍射图像的径向分布曲线。可以看出，离焦量越大，点像弥散越大。离焦量接近4倍焦深时，图像中心点光强接近0，图像反相分布。

图  8  不同离焦量的点扩散函数曲线

Figure  8.  Point-spread function curves of varied defocus distance

下载: 全尺寸图片 幻灯片

5.4   计算传递函数函数曲线

以相对孔径为D/f=0.2选择不同离焦量δ，代入（33）式，计算传递函数函数曲线，图9为不同离焦量的传递函数图：（a）为离焦量为0时的传递函数图；(b)为离焦量为400个波长时的传递函数；(c)为离焦量为1 600个波长时的传递函数；(d)为离焦量为4 000个波长时的传递函数。可以看出，离焦量越大，曲线越抖，截止频率越低。

图  6  离焦方向的判断示意图

Figure  6.  Diagram of judgement for defocus direction

下载: 全尺寸图片 幻灯片

5.5   判断图像的离焦程度

依据4.2节的爬山法，选择在1.8π～2.2π之间每隔0.02π变化假想离焦量值，得出相干传递函数的实部与虚部，代入（43）式，对图像相干复原处理得到的x方向清晰度A_c值见图10（a）曲线，发现在1.86π～1.94π之间A_c值比较大。

为了进一步得到精确离焦量值，在1.86π～1.94π之间每隔0.02π变化假想离焦量值，得出相干传递函数的实部与虚部，代入（43）式，对图像复原处理得到的x方向清晰度A_c值见图10（b）曲线，可以看出，离焦量在1.87π~1.88π之间，清晰度取得最大值。

5.6   图像复原

取假想离焦量ω₁=1.875π，代入（33）式建立离焦传递函数，依据（38）式、（39）式，对离焦模糊图像图11(a)进行复原，得到复原图见图11（b）所示。

分别计算复原前后的图11(a)和图11（b）的清晰度值、最大灰度等级、最小灰度等级，并进行比较。将图11(a)和图11（b）的傅里叶变化图代入（40）式中计算，图像的清晰度值结果见表1所示。

从表1中可以看出，以离焦量ω₁=1.875π建立点扩散函数和传递函数，对模糊图像进行复原，复原前后图像统计结果为：复原前图像的清晰度计算结果为221，复原后图像的清晰度计算结果为598；复原前图像的最大灰度等级结果为111，复原后图像的最大灰度等级结果为181；复原前图像的最小灰度等级结果为13，复原后图像的最小灰度等级结果为5；由此可得出以下结论：

1） 离焦图像复原之后，清晰度得到提升；

2） 离焦图像复原之后，调制对比度得到提升。

图  9  不同离焦量的传递函数曲线

Figure  9.  Optical transfer function curves of varied defocus distance

下载: 全尺寸图片 幻灯片

图  10  不同离焦量的清晰度值曲线

Figure  10.  Definition value curves of varied defocus distance

下载: 全尺寸图片 幻灯片

5.7   调焦方向及调焦量

将ω₁=1.875π代入(14)式得到离焦量δ=55.7 μm。从图11可以得出，图像复原后，目标区域面积减小了，说明原图像为正离焦。故系统对焦时，传感器靶面应该向靠近镜头的方向运动55.7 μm。

5.8   对焦实时性分析比较

大口径光学设备的调焦全量程（≥12 mm）步进码不低于40 000个。若采用手动全量程调整对焦，每手动改变几个码，需要留出时间以便人眼观察比较几幅图像，手动单向调整至少需要240 s，来回调整至少需要480 s。若采用计算机视觉调焦方法全量程调整对焦，步进电机每调整几个码，需要留出时间给计算机来比较几幅图像清晰度值，全量程单向调整至少需要60 s，来回调整至少需要120 s。

图  11  原图与复原图

Figure  11.  Original image and its restore image

下载: 全尺寸图片 幻灯片

表  1  图像复原前后的清晰度值

Table  1.  Definition values of original image and its restore image

项目	复原前	复原后
	图11(a)	图11（b）
清晰度值AC	221	598
最大灰度等级Hmax	111	181
最小灰度等级Hmin	13	5
调制对比度M=(Hmax−Hmin)/(Hmax+Hmin)	0.79	0.96

下载: 导出CSV 

| 显示表格

本文的基于离焦图像复原的光学设备实时视频对焦方法，在主频2.8 GHz的调焦系统工控机上试验，40多个不同假设离焦量的点扩散函数线阵、相干光传递函数线阵、相干光传递函数线阵事前预先计算并存处于计算机内存中，计算实时调用用时不到100 μs；取图像目标区域的线阵（1×100）进行离散傅立叶变换，不到300 μs时间；采用爬山法计算全量程的40多个不同假设离焦量对目标区域的线阵的复原清晰度值，不到900 μs时间；找出清晰度最大值对应的离焦量值需要200 μs；采用维纳逆滤波对原图像复原处理，提取原图像和复原图像的目标，分别计算目标面积并比较，得到离焦方向，需要4 ms；利用计算得到的离焦量大小和符号，驱动调焦机构直接到位，最多需要9 ms，如表2所示。从表2可以看出，实验总计耗时不超过15 ms。

表  2  基于离焦图像复原的光学设备实时视频对焦方法的上机试验时间

Table  2.  Time-consuming of real-time video focusing method based on defocus image restoration

过程	调用不同假设的离焦模型	线阵子图傅里叶变换	爬山法判断离焦程度	找出离焦量值	逆滤波找方向	机械对焦	总计
耗时	100 μs	300 μs	900 μs	200 μs	4 ms	9 ms	14.5 ms

下载: 导出CSV 

| 显示表格

成像系统实时采集图像和输出图像频率为25 fps，留给计算机图像处理和实时调焦时间为20 ms。比较手动对焦法、计算机视觉调焦法、基于离焦图像复原的光学设备实时视频对焦法3种方法的实时性结果，如表3所示。可以看出:大口径光学设备若采用手动对焦法、计算机视觉调焦法，调焦过程中，图像会出现模糊现象，不能满足实时性要求。采用基于离焦图像复原的光学设备实时视频对焦法，调焦机构一步到位，调焦过程中图像一直保持清晰，能够满足实时性精准对焦要求。

表  3  不同调焦方法的实时性结果比较

Table  3.  Comparison of real-time results of different focusing methods

不同方法	成像系统允许的图像 实时调焦时间	全量程调整 所用时间	成像效果	结果
手动对焦方法	20 ms	480 s	调整过程中，显示模糊图像	不实时，不精准
计算机视觉调焦方法	20 ms	120 s	调整过程中，显示模糊图像	不实时，不精准
基于离焦图像复原的光学设备实时视频对焦方法	20 ms	14.5 ms	调整过程中，显示清晰图像	实时且精准

下载: 导出CSV 

| 显示表格

6   结论

从光学镜头的成像原理出发，按照焦点附近光场变化规律和焦深变化范围，以一定的周期间隔选择某些离焦量值，推算离焦点扩散函数、相干光传递函数和非相干光传递函数，利用爬山法、逆滤波和离焦程度判断法，给出了寻找大口径光学设备拍摄目标时刻的最佳离焦量值和方向的方法，并对模糊图像给予了复原，图像的清晰度和调制对比度得到了提升。依据离焦量大小和方向，通过机械结构调整光学镜头焦面位置，实现大口径光学设备的实时精准对焦。最后，运用离焦模糊图像复原实验结果验证了方法的可行性，为解决大口径光学设备的外场使用中遇到的离焦模糊图像问题提供技术方法。在光学成像、计算机视觉及航天光学图像精密测量技术领域中具有重要的推广应用价值。