Research on focal length calibration method of oblique installation collimator

对飞行目标进行稳定跟踪并且准确解算其运行轨迹是光电经纬仪优良性能的重要体现^[1-3]。光电经纬仪焦距是影响轨迹解算准确性的主因之一 ^[4-6]，需要在出厂前对其焦距进行测试，测试方法为在平行光管焦面位置处放置双狭缝目标板，狭缝目标板经平行光管和待测经纬仪后用自带的CCD进行图像接收，利用双狭缝间距、平行光管焦距以及由CCD测得的狭缝像间距联合计算待测经纬仪焦距。由此可见，平行光管焦距是影响经纬仪焦距测试的关键因素，需要对其焦距进行精确标定。根据光电经纬仪布站方案，极端情况下系统工作最大俯仰角会达到60°，为了准确评价经纬仪焦距，在拱形检测架上需要安装不同倾角的焦面带有狭缝目标板的平行光管。通常情况下在标定焦距时都要将检测架上的平行光管拆卸下来水平放置，由于平行光管安装位置较高，频繁拆卸会给测试人员带来一定的安全风险，另外最重要的是平行光管水平放置和倾斜安装时由于受力状态不同其光学参数会有较大差异^[7-9]，因此需要对检测架上倾斜安装的平行光管焦距进行原位检测。段亚轩、赵磊、邵晶等采用放大率法测试倾斜状态下平行光管焦距时需要将基准平行光管和显微测量系统均倾斜架设并调整三者之间的相对空间位置关系从而保证光轴之间互相平行^[10-12]，这种情况下，不仅架设、调整困难而且由于重力的干扰，基准平行光管的光学参数也会发生变化，进而影响焦距的标定结果。吴玲玲、袁士东、姚震等采用精密测角法测试平行光管焦距时要求全站仪视轴和待标定平行光管光轴均处于水平状态^[13-17]，若平行光管倾斜架设，用全站仪测试两狭缝张角时会因投影误差的存在而导致两狭缝间的张角测试值和理论值存在较大差异，从而影响焦距标定结果。因此研究倾斜安装状态下平行光管焦距的标定方法具有极大的工程应用价值。

图1为采用精密测角法测试平行光管焦距的基本原理图，将平行光管放置于整平的光学平台上，将星点目标板安装于平行光管焦面位置处并用卤钨灯光源照亮目标板，在平行光管正前方放置全站仪并调整全站仪3个脚螺使其垂直轴垂直于水平面，旋转水平轴使得全站仪视轴处于水平状态。转动全站仪垂直轴轴系使得星点像在方位方向处于全站仪的视场中心，若此时星点像在俯仰方向偏离视场中心，则调整平行光管调节旋钮最终使得星点像在方位方向和俯仰方向均处于全站仪视场中心。取下星点目标板将间距为L的双狭缝目标板安装于平行光管焦面位置处，调整目标板相位使得通过全站仪观察的狭缝像和全站仪竖丝平行，转动全站仪垂直轴轴系依次瞄准对称的两狭缝像，全站仪测得的双狭缝对平行光管的张角记为ω，按公式（1）计算平行光管焦距f，则：

Figure 1.  Schematic diagram of focal length calibration based on precision goniometry

式中：$ f $为平行光管焦距，单位为mm；$ L $为狭缝间距，单位为mm；$ \omega $为双狭缝对平行光管的张角，单位为rad。

通常情况下在实验室采用图1的原理标定焦距即可，但为了避免重力因素对光管焦距产生影响以及准确评价光电经纬仪焦距，需要对平行光管进行原位检测。为了实现平行光管焦距的原位检测，根据光电经纬仪的实际工作角度范围，在检测架上安装不同倾斜角度的焦面带有狭缝目标板的平行光管，各光管光轴均相交于拱形检测架球心，在地基环上架设全站仪并调平，使用卤钨灯光源照亮待标定光管的双狭缝目标板，狭缝目标经平行光管和全站仪光学系统后用人眼接收并瞄准，标定示意图如图2所示，倾斜安装条件下全站仪测得的双狭缝对平行光管的张角记为$ \omega ' $。

Figure 2.  Schematic diagram of focal length calibration under inclined installation

当全站仪视轴处于非水平状态，此时转动垂直轴轴系就会在方位方向产生投影误差，即倾斜状态下测得的张角$ \omega ' $和水平状态下测得的张角$ \omega $不相等，倾斜角度越大两者相差越大，这是由全站仪自身的测角原理所决定的，属于系统性、原理性误差，因此需要对倾斜安装状态下测试的张角值进行修正，为此建立以下数学模型：

设图3为平行光管狭缝目标板上的点与全站仪角度之间的映射关系，$ MOQ $面为双狭缝所在的平面，$ O $点为焦平面的中心，$ {O_1} $点为平行光管主点，$ N $点为目标板十字横线与狭缝的交点，$ ON $为中心点到狭缝的距离，$ OQ $平行于水平面，$ OM $与$ OQ $垂直，以$ OQ $和$ OM $分别为横轴和纵轴建立直角坐标系$ OXY $。$ {N_1}Q\parallel MO $，$ M{N_1}\parallel OQ $，$ M' {O_1}Q' $为水平面，$ {O_1}P $为$ MP $在水平面内的投影，$ M $、$ {N_1} $、$ O $、$ Q $在水平面内的投影点为$ M' $、$ N' $、$ O' $、$ Q' $，令$ OQ = x $，$ OM = y $，$ O{O_1} = f $，$ \angle O{O_1}P = \angle {E_o} $，$ \angle P{O_1}N' = \angle A $，$ \angle N{O_1}N' = \angle E $，由空间几何关系，可推导出如下关系式：

Figure 3.  Mapping between the point on the target board and the total station angle

点$\left( {x,y} \right)$距离平行光管焦面中心的距离为：

则由公式（2）~（4）联立解得平行光管焦距：

采用光刻工艺制做一块狭缝目标板，目标板形状为圆形，直径为30 mm，基底材料为K9玻璃，具体布局见图4，狭缝目标板中心制作有十字线，十字竖线一侧平行地刻划有三条狭缝，十字线和狭缝透光，其余部分均不透光，$ {N_1} $、$ {N_2} $、$ {N_3} $分别为狭缝1、狭缝2、狭缝3与十字横线的交点，$ O $为十字线交点。三条狭缝距离十字竖线的距离分别用$ O{N_1} $、$ O{N_2} $、$ O{N_3} $表示，使用德国马尔公司制造的10106 HA型三坐标测量仪对狭缝到十字竖线的距离进行测量，测量结果由小到大依次为2.9245、6.5432、9.4979 mm。

Figure 4.  Layout of slit target board

将制做好的目标板安装到检测架上平行光管的焦面位置处，在地基环上架设瑞士徕卡公司制造的TC2003型全站仪，如图5所示，使用全站仪首先瞄准目标板经平行光管出射的十字线交叉点$ O $，记方位角为$ \angle {A_o} $，俯仰角为$ \angle {E_o} $，其次瞄准十字横线与狭缝的交点，记方位角为$ \angle A $，俯仰角为$ \angle E $。

Figure 5.  Acquisition site of testing point angle of target board

转动目标板并通过全站仪目镜进行观察，使得线段$ {N_1}{N_3} $经平行光管出射后分别成像于全站仪分划板的第一象限、第二象限、第三象限、第四象限及与分划板的竖丝和横丝平行，按照上述测试流程依次瞄准交点$ O $、$ {N_1} $、$ {N_2} $、$ {N_3} $并采集对应的方位角和俯仰角，过程数据如表1所示。

将线段${N_1}、{N_3}$分别位于不同象限及与分划板横丝平行时所采集的$ \angle {E_o} $，$ \angle A $（另$ \angle {A_o} = 0 $，则$ \angle A - \angle {A_o} = \angle A $），$ \angle E $代入公式（5）可解得修正投影误差后的平行光管焦距，最终解算的焦距值如表2所示。由表2可知，经投影误差修正后由同名测试点解算的焦距最大差值为4.41 mm，测试点$ {N_1} $的焦距平均值为1982.19 mm，波动值为4.40 mm，测试点$ {N_2} $的焦距平均值为2009.94 mm，波动值为3.39 mm，测试点$ {N_3} $的焦距平均值为2039.15 mm，波动值为3.79 mm。各测试点焦距平均值出现较大差异是由畸变引入的原理性误差所致。平行光管半视场为15 mm，则$ {N_1} $点位于0.19视场位置处，$ {N_2} $点位于0.44视场位置处，$ {N_3} $点位于0.63视场位置处，由平行光管的光学设计结果可知0.19视场处系统畸变为0.02%，该值可忽略，0.44及0.63视场处系统畸变分别为−1.3%和−2.8%，根据畸变设计结果对$ O{N_2} $、$ O{N_3} $的长度进行修正重新代入公式（5）计算焦距，如表3所示，测试点$ {N_2} $的焦距平均值变为1983.81 mm、测试点$ {N_3} $的焦距平均值变为1982.06 mm，修正畸变后，测试点$ {N_2} $和$ {N_3} $的焦距平均值与测试点$ {N_1} $的焦距平均差值分别为1.62 mm和−0.13 mm。

为了解耦投影误差并获得焦距的真值，根据全站仪工作原理，转动水平轴轴系进行俯仰方向的角度测量时不存在投影误差，因此可以利用表1中线段${N_1},{N_3}$与分划板竖丝平行时的过程数据解算焦距并作为真值来验证模型的正确性，需要注意的是，公式中的$ \omega $需以$ 2\left| {\angle {E_o} - \angle E} \right| $代替，最终解算的焦距值如表4所示。

同样由于畸变的存在，各测试点解算的焦距值差异较大，测试点$ {N_1} $、$ {N_2} $、$ {N_3} $的焦距值分别为1980.03、2009.58、2039.91 mm，进行畸变修正后对应的焦距值分别为1980.03、1983.45、1982.79 mm，波动值为3.43 mm，各测试点焦距平均值即真值为1982.09 mm。线段${N_1}、{N_3}$位于不同象限及与分划板横丝平行时的各测试点在修正投影误差和畸变后焦距平均值为1982.69 mm，此值与真值差值为0.6 mm，验证了模型的正确性及鲁棒性。

修正畸变后将线段${N_1}、{N_3}$与分划板横丝平行时所采集的$ \angle {A_o} $、$ \angle {E_o} $、$ \angle A $及$ \angle E $代入公式（2）中可解得传统意义下的焦距，需要注意的是公式中的$ \omega $需以$ 2\left| {\angle A - \angle {A_o}} \right| $代替，最终解算的焦距值如表5所示。

由表5可知，修正畸变但未修正投影误差的情况下由测试点$ {N_1} $、$ {N_2} $、$ {N_3} $解算的焦距结果分别为996.42、995.23、995.22 mm，焦距平均值与真值之差为986.47 mm，相对误差为50.2%。

令公式（5）中的

则焦距f的表达式可变换为：

$\angle A$、$\angle E$、$\angle {E_\textit{o}}$及$O{N_1}$对焦距的灵敏度系数分别如下式所示：

设$\angle A$、$\angle E$、$\angle {E_o}$及$O{N_1}$的测量标准不确定度分别为$ {\mu _A} $、${\;\mu _E}$、$ {\;\mu _{{E_o}}} $、$ {\;\mu _{O{N_1}}} $且相互独立。则线段${N_1}、{N_3}$位于不同象限及与分划板横丝平行时焦距测量结果的不确定度应是所有不确定度分量的合成，用合成标准不确定度${\;\mu _{{c}}}$来表示，如下式所示：

线段${N_1}、{N_3}$与分划板竖丝平行时焦距的解算公式可表示为：

合成标准不确定度用${\;\mu _{{c}}}^\prime$来表示，如下式所示：

通过查验全站仪和三坐标测量仪计量证书可知$ \;{\mu }_{A}\text{=0}\text{.3}″ $，$ \;{\mu }_{E}\text{=}{\mu }_{{E}_{o}}\text{=}0.4″ $，$\; {\mu _{O{N_1}}}{\text{ = }}1.0\;{\text{μ}} {\rm{m}}$。修正畸变后将采集的过程数据及相应的测量标准不确定度代入公式（11）和公式（13）中，可得各测试点对应焦距的合成标准不确定度，取包含因子k=2，则展伸不确定度U=2${\;\mu _{{c}}}$（2$\;{\mu _{{c}}}'$），如表6所示。

由表6可知：线段${N_1}、{N_3}$位于不同象限及与分划板横丝平行时，各测试点的焦距测量结果展伸不确定度均大于与竖丝平行时的焦距测量结果展伸不确定度，测量结果的展伸不确定度与焦距真值的最大相对误差为0.36%，该值远小于GB/T 9917.1—2002 照相镜头中实测焦距对名义焦距的相对误差不超过±5%的规定。

为了精准评价系统焦距，建立了倾斜安装状态下光管焦距标定的数学模型，结合各测试工况下所采集的过程数据证实了目标狭缝在分划板中的相位可以是随机值，不需要和全站仪竖丝严格平行。线段${N_1}、{N_3}$与分划板横丝平行时的各测试点在修正畸变但未修正投影误差的情况下焦距平均值与真值相对误差为50.2%。线段${N_1}、{N_3}$位于不同象限及与分划板横丝平行时的各测试点在修正投影误差和畸变后焦距平均值与真值差值为0.6 mm。测量结果的展伸不确定度与焦距真值的最大相对误差为0.36%，该值远小于GB/T 9917.1—2002 照相镜头中实测焦距对名义焦距的相对误差不超过±5%的规定。实验结果表明：该模型具有普适性和较高的解算精度，对于倾斜安装条件下平行光管焦距的原位检测具有极大的工程应用价值。