Infrared modulation detecting technology of electro-optic radar
-
摘要: 因红外调制检测技术用于检测扫描积分型光电雷达调制器的工作性能，为此设计一套与调制器固有属性匹配的光学传递和激光准直系统。运用标准伺服控制理论为调制脉码信号建立Ⅲ型条件稳定的轴系扫描模型，通过截获、跟踪视场切换及驱动电压调幅，使模拟红外汇聚的激光扫描光束与3°×3°、40′×40′十字靶标重合，计量扫描视场及零位精度偏差，采集标定反射镜振动位置和扫描频率的过零脉冲，用以实现红外信号调制性能的检测。测试结果表明，该红外调制检测技术对调制器扫描视场的精度测量优于0.5′(1σ)，零位误差的精度测量不大于5′(1σ)，系统检测的相对误差不大于±1%，满足光电雷达红外光轴瞄准线工艺装配的精度要求，可应用于线列扫积型调制器的性能检测。Abstract: Infrared modulation detection technique is used to detect performance of scanning integral type electro-optical radar modulator. An optical transmission and laser collimation system are designed in this paper, which matches inherent properties of the modulator. Based on standard servo control theory, axis scanning model of Ⅲ system with stable condition is established for modulation of pulse code signal. Through FOV's switching of intercepting and tracking, and modulation of driving voltage amplitude, simulated infrared-converged laser scanning beam is coincided with 3°×3°, 40′×40′ cross target. Deviation of scanning FOV and zero position precision are metered, calibration mirror vibration position and scanning frequency zero-crossing pulse are collected in order to achieve infrared signal modulation performance detection. Results show that accuracy of infrared modulation detection technique is better than 0.5′(1σ), accuracy of zero error is less than 5′(1σ), and relative error of the system is less than ± 1%.It meets accuracy requirements of electro-optical radar in infrared optical axis assembly process, and can be applied to linear array scanning integral type modulator performance testing.
-
引言
线阵CCD图像传感器具有视野极大、扫描频率高、分辨率高等特点，广泛应用于非接触式测量和表面检测等领域[1]。二维坐标测量系统主要是对基准面上的坐标进行测量，广泛应用于机床、路径规划、障碍定位等领域[2]。目前，坐标测量可以分为接触式测量和非接触式测量，接触式测量具有适用范围的局限性，而非接触式测量具有广泛的适用范围，非接触式测量方法主要有激光扫描、影像测量、照相测量和工业CT扫描等[3]。激光扫描方法由半导体激光器发出激光，聚光透镜在被测曲面上结成光电并反射，利用光敏元件接收其散射光，根据其在光敏元件上的位置，从而测出被测点的坐标位置[4]。如朱博源等[5]将多个角锥棱镜作为无源待测靶点，通过构建非球面反射接收模型实现了信号发射端与接收端的结合，建立了多靶点光信号匹配机制以实现无源多靶点交会定位；影像测量方法通过光学放大系统对被测物体放大，CCD图像采集系统采集图像特征并送入计算机，实现高效地检测[6]。如白玮本等[7]采用基于结构光的嵌入式图像采集系统，借助光栅发生器对被测物进行扫描，CCD相机抓取被测物光条纹信息，通过以太网传输至PC端；照相测量方法通过相机多角度采集被测对象得到图像信息，依据空间物体投影原理，利用被测对象表面的标志点特征计算被测对象的三维空间信息，多应用于三维重构领域[8]；工业CT扫描方法是将计算机断层成像技术应用于工业中，其基本原理是依据辐射在被检测物体中的减弱和吸收特性，广泛应用于医学和三维重构领域[9-10]；角度交会法通过两个角度的测量进而计算得到被测对象坐标，具有测量精度高，适用范围广等特点，广泛应用于工程项目[11]。
综上所述，基于影像测量法的高效优势和角度交会法的高精度特点，针对基准平面障碍物检测定位问题，本文设计了一种基于角度交会的嵌入式二维坐标图像测量系统，将图像检测技术应用于二维坐标测量，可以发挥线阵CCD图像传感器大视野、高扫描频率、高分辨率和嵌入式设备实时性高等优势，能够实时、准确、快速和高效地测量显示障碍物二维坐标，可以广泛应用于基准平面障碍物检测定位，具有一定的实用价值和意义。
1 系统测量原理
如图1所示，二维坐标测量系统主要由被测对象、光源、光路模块、传感模块、单片机（microcontroller unit, MCU）和显示模块组成。二维坐标测量系统将光电技术结合，在外部光源照射的条件下，光路模块发射光信号到被测物体表面，采集模块检测到物体表面反射的光线并将采集的数据传输给单片机，单片机对数据进行处理和算法运算，最终得到被测对象的二维坐标信息并在显示模块上面显示。
![]()
图 1 二维坐标测量系统总体框图Figure 1. Overall block diagram of two-dimensional coordinate measurement system二维坐标测量系统采用角度交会法对被测对象进行坐标计算。角度交会法是一种通过角度的测量，进而得到未知坐标的算法，具有测量算法公式简单、系统结构装置简易及测量精度较高的特点。如图2所示，A点放置位置为固定的步进电机带动激光光源，B点放置位置为固定的线阵CCD图像传感器，C点为被测对象的位置。假设点A、B和C的坐标依次为（XA,YA）、（XB,YB）和（XC,YC），A点和B点坐标已知，若要求解C点的坐标，则需要测量角$ \alpha $和$ \beta $的大小，然后应用式（1）得到C点坐标：
$$ \left\{\begin{array}{l} X_{C}=\dfrac{X_{B} \cdot \tan \beta}{\tan \alpha+\tan \beta} \\ Y_{C}=\dfrac{X_{B} \cdot \tan \beta \cdot \tan \alpha}{\tan \alpha+\tan \beta} \end{array} \right.$$ (1) ![]()
图 2 角度交会法求解坐标示意图Figure 2. Schematic diagram of solving coordinates by angle intersection method在图2中，AC为入射光线，CB为反射光线，$ \alpha $角度的大小可以通过步进电机的转动角来测得，$ \beta $角度的大小可以通过线阵CCD输出的信号处理得到，采用角度交会法便可以求解出三角模型中C点的坐标。
2 硬件设计
测量系统采用STM32F103ZET6单片机，核心芯片为STM32F103ZET6；传感模块采用TSL1401线阵CCD图像传感器，选用90°广角镜头；显示模块采用NT35310 TFT-LCD显示屏；同时采用28BYJ-48步进电机和ULN2003步进电机驱动，以满足测量系统设计精度的需求。系统电路原理如图3所示。
MCU驱动线阵CCD图像传感器和显示模块工作，ULN2003步进电机驱动步进电机工作。MCU利用STM32的定时器中断功能编写线阵CCD的驱动时序，无需外加AD转换模块就可以实现将线阵CCD图像传感器采集的信号从模拟信号转换为数字信号。STM32使用FSMC驱动TFT-LCD显示屏，实现显示模块的显示功能。因单片机引脚输出电压过低无法驱动步进电机，所以采用ULN2003步进电机驱动和外设电源的方案，通过MCU将电信号传输给ULN2003，间接实现驱动步进电机的功能。测量系统的各个模块相互配合，完成数据的采集、处理和运算，并对二维坐标进行实时显示。
3 软件设计
软件设计部分主要是MCU对步进电机转动、线阵CCD图像传感器数据采集和TFT-LCD显示进行控制，经运算处理得到相关参数，从而得到被测对象的坐标并显示。MCU控制步进电机转动设定角度，从而带动激光光源转动并扫描平面内障碍物；MCU驱动线阵CCD图像传感器实时采集视场范围内的状况，并将采集的信息转化为电信号传输到MCU进行处理得到相关参数；运算处理得到二维坐标，驱动TFT-LCD显示被测对象的坐标信息。在MCU和外部设备相互配合下，完成实验设定条件的控制、数据采集、数据处理、数据计算和结果显示。该系统的主程序流程图如图4所示。
由图4可知，系统上电后进行初始化，步进电机复位并开始运行程序；步进电机带动激光光源以单位脉冲周期$ 5.625 / 64^{\circ} $的步距角进行转动，并对脉冲数进行计数；当激光光源扫描到被测物时，实时返回脉冲数并计算此时的角度$ \alpha $，此时线阵CCD对被测对象进行数据采集并处理计算得到角度$ \beta $，若未扫描到被测对象，步进电机仍旧以单位脉冲步距角带动激光光源转动直至扫描到被测对象；当被测对象扫描完成时，单片机计算得到角度$ \alpha $和$ \beta $的期望，通过运算处理得到被测对象的坐标，最终在TFT-LCD显示屏上显示，若被测对象未完成扫描，步进电机继续带动激光光源进行扫描直至完成被测对象的扫描。
3.1 CCD驱动及采集
利用STM32定时器完成TSL1401线阵CCD图像传感器驱动时序。线阵CCD包含128个感光单元，tqt为采样周期间隔，驱动时序图如图5所示。采用时钟管脚PA7产生CLK信号，IO串口PA6作为SI信号输入，从而输出模拟信号AO[12]。
![]()
图 5 TSL1401线阵CCD驱动时序图Figure 5. Drive time-sequence diagram of TSL1401 linear array CCD利用STM32的定时器来触发ADC的启动转换[13]。系统使用单次单通道ADC1触发时序图如图6所示，tSTAB为上电延时；ADC_CLK为ADC1时钟；SET ADON为设置ADC_CR寄存器的ADON位；EOC为规则组转换结束事件标志。设置PA1为模拟输入，通过RCC_CFGR设置ADC1分频因子，最终完成ADC1的模式配置，实现单次转换模式、触发方式选择、数据对齐方式等功能。
3.2 步进电机驱动
如图7所示，系统选用四相八拍的28BYJ48步进电机，对步进电机施加一系列特定顺序并且连续不断地控制脉冲时，便可以连续不断地转动。每个脉冲信号使得步进电机的一相或两相绕组的通电状态改变一次，即转子转过一定角度，当通电状态的改变完成1个循环时，转子转过1个步距角。所以启动步进电机，只需要依次给各个相输入高电平信号。系统采用八拍方式驱动，STM32的4个GPIO串口依次与步进电机驱动的IN1-IN4端口连接，即STM32的PC6、PC7、PC8和PC9串口配置为I/O口，类似PC9、PC8、PC7和PC6这4个端口的输出信号，可以实现对步进电机启动、停止、正转和反转的控制。
3.3 CCD标定及$\;{\beta} $ 计算
采用双直线法标定线阵CCD。如图8所示，在线阵CCD图像传感器的视场范围内进行位置标定，在二维坐标测量扫描区的远近平行位置$Y_{{\rm{F}}}$和$Y_{{\rm{N}}}$处对线阵CCD标定采样，其中$Y_{{\rm{F}}}$和$Y_{{\rm{N}}}$为线阵CCD标定的2条采样直线，并且$ \angle M_{i} $和$ \angle N_{i} $ 与$ \beta $相等；假设被测物体在$ m_{i} $和$ n_{i} $点时，线阵CCD图像传感器检测$ m_{i} $和$ n_{i} $点的反射光线重合，则线阵CCD产生的输出信号也相同，可以据此对线阵CCD标定。在远位置$Y_{{\rm{F}}}$上取等间距的20个采样点，将被测对象依次放在采样点处，线阵CCD图像传感器进行数据采集，微控制单元对线阵CCD采集的数据进行处理得到采样点CCD的像元值$ {C_{{x_i}}} $。采用同样的方法对近位置$Y_{{\rm{N}}}$处进行标定采样，因为$ m_{i} $和$ n_{i} $点处进入线阵CCD的入射角$ \angle M_{i} $和$ \angle N_{i} $相等，即可以根据已得的像元值进一步找到$ n_{i} $点。相反，已知$ m_{i} $和$ n_{i} $点的坐标，可以得到$ M_{i} $和$ N_{i} $角度的大小，即可以得到像元值$ {C_{{x_i}}} $与$ \angle M_{i} $和$ \angle N_{i} $ 一一对应的关系。取像元值$ {C_{{x_i}}} $为横坐标, $ \angle M_{i} $和$ \angle N_{i} $为纵坐标，可以得到关于$ {C_{{x_i}}} $与$ \angle M_{i} $和$ \angle N_{i} $的散点图，采用最小二乘法对数据进行拟合，得到$ {C_{{x_i}}} $与$ \angle M_{i} $ 和$ \angle N_{i} $的关系式[14]。
如图9所示，标定图案为在白色平板上刻画黑色竖线，相邻竖线间隔可选。首先选用相邻竖线间隔10 mm的标定图案，将标定图案竖直放置于$Y_{{\rm{F}}}=10$ 的直线上，并将最左边边缘线与$ X=0 $对齐；然后，步进电机带动激光光源在设定的速度下转动，激光光源会在被照射的标定图像竖线上形成点$ m_{i} $；之后，线阵CCD图像传感器进行数据采集，并将信号输出给MCU；最后，MCU对输入的信号进行处理，得到点$ m_{i} $对应的平均像元值，并记录实验数据，如表1所示。
表 1 $ m_{i} $采样点与对应像元值测试数据Table 1. Sampling points $ m_{i} $ and corresponding pixel value test data${\boldsymbol{m} }_{ {\boldsymbol{i} } }\_{ {\boldsymbol{X} } }_{ {\boldsymbol{i} } }$ ${C_{{x_i}}}$ ${\boldsymbol{m} }_{ {\boldsymbol{i} } }\_{ {\boldsymbol{X} } }_{ {\boldsymbol{i} } }$ ${C_{{x_i}}}$ ${\boldsymbol{m} }_{ {\boldsymbol{i} } }\_{ {\boldsymbol{X} } }_{ {\boldsymbol{i} } }$ ${C_{{x_i}}}$ 1 114 8 72 15 26 2 109 9 66 16 20 3 103 10 62 17 16 4 98 11 54 18 9 5 91 12 45 19 5 6 84 13 40 20 0 7 78 14 32 选用相邻竖线间隔5 mm的标定图案，将标定图案竖直放置于$Y_{{\rm{N}}}=5$MT4 的直线上，重复上述操作，得到点$ n_{i} $对应的平均像元值，并进行记录，如表2所示。
表 2 $ n_{i} $采样点与对应像元值测试数据Table 2. Sampling points $ n_{i} $ and corresponding pixel value test data${\boldsymbol{n} }_{ {\boldsymbol{i} } }\_{ {\boldsymbol{X} } }_{ {\boldsymbol{i} } }$ ${C_{{x_i}}}$ ${\boldsymbol{n} }_{ {\boldsymbol{i} } }\_{ {\boldsymbol{X} } }_{ {\boldsymbol{i} } }$ ${C_{{x_i}}}$ ${\boldsymbol{n} }_{ {\boldsymbol{i} } }\_{ {\boldsymbol{X} } }_{ {\boldsymbol{i} } }$ ${C_{{x_i}}}$ 5.0 0 8.5 78 12.0 33 5.5 114 9.0 72 12.5 26 6.0 108 9.5 66 13.0 20 6.5 103 10.0 61 13.5 15 7.0 98 10.5 54 14.0 9 7.5 90 11.0 45 14.5 5 8.0 84 11.5 40 15.0 1 对比表1和表2可知，点$ n_{i} $对应的像元值几乎可以与点$ m_{i} $对应的像元值一一对应。利用Matlab将上述的2个表格中的点$ m_{i} $和$ n_{i} $描绘在直角坐标系中，如图10所示，并做相关的辅助线，线阵CCD采集$ m_{i} n_{i} $线段上点的像元值输出基本相同，即$ m_{i} n_{i} $线段上点对应的像元值相同，证明采用双直线标定线阵CCD是成立的。
![]()
图 10 $ m_{i} $和$ n_{i} $采样点展示图Figure 10. Display diagram of sampling points $ m_{i} $ and $ n_{i} $$Y_{{\rm{F}}}=10$和$Y_{{\rm{N}}}=5$ 的采样点$ m_{i} $和$ n_{i} $坐标已知，由式（2）计算得$ \angle M_{i} $和$\angle N_{i}$的大小，因$ \angle M_{i} $和$\angle N_{i}$与$ \beta_{i} $相等，即求得$ \beta_{i} $的大小。由表1和表2计算线段$ m_{i} n_{i} $与X轴的夹角$ \beta_{i} $，并将$ \beta_{i} $和相应的像元值记录如表3所示。
表 3 像元值与对应$ \beta_{i} $数据Table 3. Image value and the corresponding $ \beta_{i} $ data$ {C_{{x_i}}} $ $ {{\boldsymbol{\beta}} _{\boldsymbol{i}}} $ $ {C_{{x_i}}} $ $ {{\boldsymbol{\beta}} _{\boldsymbol{i}}} $ $ {C_{{x_i}}} $ $ {{\boldsymbol{\beta}} _{\boldsymbol{i}}} $ 0 45.0 78 73.3 32 111.8 114 48.0 72 78.8 26 116.6 108 51.3 66 84.3 20 121.0 103 55.0 61 90.0 15 125.0 98 59.0 54 95.7 9 128.7 90 63.4 45 101.3 5 132.0 84 68.2 40 106.7 1 135.0 $$ \begin{array}{l} \tan \left(\angle M_{i}\right)=\tan \left(\angle N_{i}\right)=-\left(\dfrac{Y_{m_i}-Y_{n_i}}{X_{m_i}-X_{n_i}}\right) \\ \tan \left(\beta_{i}\right)=-\left(\dfrac{Y_{m_i}}{X_{m_i}-10}\right) \end{array} $$ (2) 将像元值${C_{{x_i}}}$作为横坐标，$ \beta_{i} $作为纵坐标绘制相应的散点图，利用Matlab进行线性拟合如图11所示。设${\beta _i}=a \cdot {C_{{x_i}}}+b$，最小二乘法线性回归方程如式（3）所示，即$ \beta_{i} $与像元值的关系式为
$$ \beta_{i}=-0.79{C_{{x_i}}}+136.8 $$ (3) 式中：
$$ a=\frac{N \cdot \displaystyle\sum\limits_{i=1}^{N} x_{i} y_{i}-\left(\displaystyle\sum\limits_{i=1}^{N} x_{i}\right) \cdot\left(\displaystyle\sum\limits_{i=1}^{N} y_{i}\right)}{N \cdot \displaystyle\sum\limits_{i=1}^{N} x_{i}^{2}-\left(\displaystyle\sum\limits_{i=1}^{N} x_{i}\right)^{2}}=-0.79 $$ (4) $$ b=\frac{\left(\displaystyle\sum\limits_{i=1}^{N} x_{i}^{2}\right) \cdot\left(\displaystyle\sum\limits_{i=1}^{N} y_{i}\right)-\left(\displaystyle\sum\limits_{i=1}^{N} x_{i}\right) \cdot\left(\displaystyle\sum\limits_{i=1}^{N} x_{i} y_{i}\right)}{N \cdot \displaystyle\sum\limits_{i=1}^{N} x_{i}^{2}-\left(\displaystyle\sum\limits_{i=1}^{N} x_{i}\right)^{2}}=136.8 $$ (5) ![]()
图 11 $ \beta_{i} $与像元值的最小二乘法线性回归方程拟合图Figure 11. $ \beta_{i} $ and least square linear regression equation for pixel values3.4 坐标测量算法
被测对象在线阵CCD图像传感器的视场范围内，步进电机带动激光光源装置放置于原点(0,0)，线阵CCD图像传感器放置于点(10,0)处，MCU直接控制步进电机的转动，间接调整激光光源与X轴夹角的变化。步进电机的转角即就是激光光源的转角，选用型号为28BYJ-48、步距角为$ 5.625 / 64^{\circ} $、额定电压为5V的步进电机，步进电机转动的角度等于脉冲个数与步距角的乘积且与$ \alpha $互余，如式（6）所示：
$$ \alpha=90-N \cdot 5.625 / 64^{\circ} $$ (6) 由式（3）和式（6）可得入射光线和反射光线的方程式。
入射光线方程式：
$$ Y_{C}=X_{C} \cdot \tan (\alpha)=X_{C} \cdot \tan \left(90-N \cdot 5.625 / 64^{\circ}\right) $$ (7) 反射光线方程式：
$$ \begin{split} Y_{C}=&-\left(X_{C}-10\right) \cdot \tan \left(\beta\right)=-\left(X_{C}-10\right) \cdot\\ & \tan \left(-0.79 {C_{{x_i}}}+136.8\right) \end{split} $$ (8) 由式（6）和式（7）联合求解，得XC与像元值的关系式：
$$ \begin{split} X_{C}=&\frac{10 \cdot \tan \left(\beta_{i}\right)}{\tan (\alpha)+\tan \left(\beta_{i}\right)}=\\ &\frac{10 \cdot \tan \left(-0.79 {C_{{x_i}}}+136.8\right)}{\tan \left(90-N \cdot 5.625 / 64^{\circ}\right)+\tan \left(-0.79 {C_{{x_i}}}+136.8\right)} \end{split} $$ (9) 由式（6）和式（8）联合求解，得YC与像元值的关系式：
$$ \begin{split} &Y_{C}=\frac{10 \cdot \tan \left(\beta_{i}\right) \cdot \tan (\alpha)}{\tan (\alpha)+\tan \left(\beta_{i}\right)}=\\ &\frac{10 \cdot \tan \left(-0.79 {C_{{x_i}}}+136.8\right) \cdot \tan \left(90-N \cdot 5.625 / 64^{\circ}\right)}{\tan \left(90-N \cdot 5.625 / 64^{\circ}\right)+\tan \left(-0.79 {C_{{x_i}}}+136.8\right)} \end{split} $$ (10) 4 测试结果及分析
二维坐标图像测量系统如图12所示，采用控制变量法分别对X轴和Y轴坐标进行分析。将被测对象放置于Y=10处，在线阵CCD图像传感器的视野范围内以及X轴坐标准确的前提下，以1 cm为间隔移动被测对象，得到Y轴固定X轴变化的二维坐标，如表4所示。
表 4 Y=10时坐标测量数据Table 4. Coordinate test data when Y=10$ {\boldsymbol{Y}} $ $ {\boldsymbol{X}} $ ${\boldsymbol{Y} }_{ {{\rm{test}} } }$ $ {\boldsymbol{X}} $ ${\boldsymbol{Y} }_{ {{\rm{test}} } }$ $ {\boldsymbol{X}} $ ${\boldsymbol{Y} }_{ {{\rm{test}} } }$ 10 1 10.22 8 9.86 15 10.20 10 2 10.07 9 9.87 16 10.17 10 3 10.00 10 9.91 17 10.10 10 4 9.96 11 9.97 18 10.09 10 5 9.92 12 10.04 19 9.94 10 6 9.89 13 10.12 10 7 9.87 14 10.17 计算知Y=10处的测量值$\bar{Y}_{ {{\rm{test}} }}=10.019$，同上测量Y=2至Y=11间隔0.5 cm的直线上的坐标，并对数据进行处理，记录如表5所示。对Y轴坐标的测量误差进行多项式线性拟合，得到如图13所示的Y轴坐标测量值误差的拟合多项式误差曲线[15]。误差曲线的6次项和5次项系数比较小可以忽略，误差曲线可以表示为$y=0.01 x^{4}-0.09 x^{3}+0.38 x^{2}- 0.83 x+0.75$，即可以取Y轴坐标测量值误差的拟合曲线的相反数作为修正值，则被测对象的Y轴修正坐标为
$$ \begin{split} {Y_i} = & \frac{{10 \cdot \tan \left( { - 0.79{C_{{x_i}}} + 136.8} \right) \cdot \tan \left( {90 - N \cdot 5.625/{{64}^{\circ} }} \right)}}{{\tan \left( {90 - N \cdot 5.625/{{64}^{\circ} }} \right) + \tan \left( { - 0.79{C_{{x_i}}} + 136.8} \right)}} -\\ &0.01{x^4} + 0.09{x^3} - 0.38{x^2} + 0.83x - 0.75 \end{split}$$ 表 5 不同Y轴坐标测量数据Table 5. Different Y-axis coordinate measurement data$ {\boldsymbol{Y}} $ $\bar{\boldsymbol{Y} }_{ {{\rm{test}} } }$ $ \Delta {\boldsymbol{Y}} $ $ {\boldsymbol{Y}} $ $\bar{\boldsymbol{Y} }_{ {{\rm{test}} } }$ $ \Delta {\boldsymbol{Y}} $ 2.0 2.045 0.045 7.0 6.994 −0.006 2.5 2.515 0.015 7.5 7.508 0.008 3.0 3.000 0.000 8.0 8.024 0.024 3.5 3.491 −0.009 8.5 8.531 0.031 4.0 3.983 −0.017 9.0 9.038 0.038 4.5 4.481 −0.020 9.5 9.534 0.034 5.0 4.975 −0.025 10.0 10.020 0.020 5.5 5.474 −0.026 10.5 10.506 0.006 6.0 5.975 −0.025 11.0 10.988 −0.012 6.5 6.483 −0.017 ![]()
图 13 Y轴坐标的测量误差多项式拟合曲线图Figure 13. Polynomial fitting curve of measurement error for Y-axis coordinates同上操作，测量X=2至X=19间隔1 cm的直线上的坐标，对X轴坐标进行分析处理如表6所示。
表 6 不同X轴坐标测量数据Table 6. Different X-axis coordinate measurement data$ {{\boldsymbol{X}}} $ $\overline{ { {\boldsymbol{X} } } }_{{\rm{test}}}$ $ \Delta {{\boldsymbol{X}}} $ $ {{\boldsymbol{X}}} $ $\overline{ { {\boldsymbol{X} } } }_{{\rm{test}}}$ $ \Delta {{\boldsymbol{X}}} $ 2 1.982 −0.018 11 10.987 −0.013 3 2.984 −0.016 12 11.977 −0.023 4 3.999 −0.001 13 12.974 −0.026 5 5.013 0.013 14 13.976 −0.024 6 6.022 0.022 15 14.988 −0.012 7 7.026 0.026 16 16.001 0.001 8 8.021 0.021 17 17.013 0.013 9 9.012 0.012 18 18.026 0.026 10 10.00 0.000 19 19.034 0.034 同上操作，误差曲线的5次项系数比较小可以忽略，如图14所示，误差曲线可以表示为$y= 0.004\;6 x^{4}-0.006\;5 x^{3}+0.044 x^{2}-0.12 x+0.092$，即被测对象的X轴修正坐标为
$$ \begin{split} {X_i} = &\frac{{10 \cdot \tan \left( { - 0.79{C_{{x_i}}} + 136.8} \right)}}{{\tan \left( {90 - N \cdot 5.625/{{64}^{\circ} }} \right) + \tan \left( { - 0.79{C_{{x_i}}} + 136.8} \right)}} -\\ &0.000\;46{x^4} + 0.006\;5{x^3} - 0.044{x^2} + 0.12x - 0.092 \end{split}$$ ![]()
图 14 X轴坐标的测量误差多项式拟合曲线图Figure 14. Polynomial fitting curve of measurement error for X-axis coordinates将上述X轴和Y轴坐标修正公式加入系统程序，得到如图15所示的点(5,6)坐标的测量结果。与未修正的坐标测量结果对比，可知修正后的坐标测量结果误差减小了一些，X轴坐标误差由−0.006 0 cm减小到−0.001 3 cm，Y轴坐标误差由−0.039 9 cm减小到−0.010 4 cm，即这样的坐标修正方法是可行的。
5 结论
基于角度交会法设计了一种二维坐标图像测量系统，分析处理测量结果得到修正公式，进而修正坐标计算公式提高坐标测量精度。系统采用线阵CCD图像传感器为传感模块，发挥线阵CCD大视野、高扫描频率、高分辨率的特点，能够准确、快速、可靠的采集被测对象的信息；搭建三角测距光路模型，测量两个角度参数，采用角度交会法计算二维坐标；利用Matlab分别分析处理X轴坐标和Y轴坐标，对坐标测量误差进行多项式拟合，得到坐标修正公式，并将其应用于坐标计算，实验测量结果表明修正后的系统测量精度更高。
-
![]()
图 2 调制器在红外物镜组件工作原理示意图
Figure 2. Working schematic of modulator in infrared objective lens
表 1 测试结果
Table 1 Test results
测试项 3°×3°截获场 40′×40′跟踪场 激光理论偏转角度 1.5° 20′ 扫描反射镜理论振幅 0.75° 10′ 驱动工作电压 9~13 V 17~20 V 扫描视场精度 优于0.5′(1σ) 优于0.2′(1σ) 零位误差精度 不大于5′(1σ) 不大于5′(1σ) 
下载: 导出CSV
-
[1] 陈洁, 张若岚.应用于红外搜索跟踪和态势感知系统的全景成像技术[J].红外技术, 2016, 38(4):269-279. http://d.old.wanfangdata.com.cn/Periodical/hwjs201604001 Chen Jie, Zhang Ruolan. Panoramic imaging technology applied in IRST and status awareness system[J]. Infrared Technology, 2016, 38(4):269-279. http://d.old.wanfangdata.com.cn/Periodical/hwjs201604001
[2] 贾俊涛, 唐跃峰, 毛鑫, 等.红外全景扫描跟踪成像系统设计与实现[J].应用光学, 2013, 34(3):407-412. http://d.old.wanfangdata.com.cn/Periodical/yygx201303005 Jia Juntao, Tang Yuefeng, Mao Xin, et al. Design and implement of infrared panoramic scanning and tracking imaging system[J]. Journal of Applied Optics, 2013, 34(3):407-412. http://d.old.wanfangdata.com.cn/Periodical/yygx201303005
[3] 王卫华, 黄宗福, 何艳, 等.一种基于线阵扫描成像的红外搜索跟踪一体化系统设计[J].信号处理, 2010, 26(9):1312-1317. doi: 10.3969/j.issn.1003-0530.2010.09.006 Wang Weihua, Huang Zongfu, He Yan, et al. An incorporate system design of infrared search and track based on line-array-can imaging[J]. Signal Processing, 2010, 26(9):1312-1317. doi: 10.3969/j.issn.1003-0530.2010.09.006
[4] 范宏波.基于1152×6长波线列探测器的高性能红外搜索预警系统[J].红外技术, 2010, 32(1):20-24. doi: 10.3969/j.issn.1001-8891.2010.01.005 Fan Hongbo. A high performance IRST system based on 1152×6 LWIR detectors[J]. Infrared Technology, 2010, 32(1):20-24. doi: 10.3969/j.issn.1001-8891.2010.01.005
[5] 焦明印.采用480×6元长波红外探测器的搜索光学系统设计[J].应用光学, 2012, 33(6):1011-1013. http://d.old.wanfangdata.com.cn/Periodical/yygx201206001 Jiao Mingyin. Optical design of searching system using 480×6 LWIR detector[J]. Journal of Applied Optics, 2012, 33(6):1011-1013. http://d.old.wanfangdata.com.cn/Periodical/yygx201206001
[6] 罗军辉, 姬红兵, 王延光.全向红外搜索跟踪系统的设计与实现[J].红外与激光工程, 2007, 36(6):869-872. doi: 10.3969/j.issn.1007-2276.2007.06.026 Luo Junhui, Ji Hongbing, Wang Yanguang. System design and realization of the all-direction IRST[J]. Infrared and Laser Engineering, 2007, 36(6):869-872. doi: 10.3969/j.issn.1007-2276.2007.06.026
[7] 安永泉, 禹健. 576×6长波红外探测器成像系统设计[J].激光与红外, 2009, 39(2):173-177. doi: 10.3969/j.issn.1001-5078.2009.02.015 An Yongquan, YuJian. Infrared imaging system design for 576×6 linear detector[J].Laser & Infrared, 2009, 39(2):173-177. doi: 10.3969/j.issn.1001-5078.2009.02.015
[8] 杨乐, 孙强, 王健.长波红外连续变焦光学系统设计[J].红外与激光工程, 2012, 41(4):999-1004. http://www.cnki.com.cn/Article/CJFDTotal-HWYJ201204030.htm Yang Le, Sun Qiang, Wang Jian. Design of long-wave infrared continuous zoom optical system[J]. Infrared and Laser Engineering, 2012, 41(4):999-1004. http://www.cnki.com.cn/Article/CJFDTotal-HWYJ201204030.htm
[9] 骆守俊, 夏寅辉, 杨宁宁, 等.扫描型长波红外连续变焦光学系统[J].中国光学, 2015(1):107-113. doi: 10.3969/j.issn.1007-2276.2012.04.031 Luo Shoujun, Xia Yinhui, Yang Ningning, et al. Long-wavelength infrared continuous zoom scanning optical system[J]. Chinese Optics, 2015(1):107-113. doi: 10.3969/j.issn.1007-2276.2012.04.031
[10] 杨育周，张双垒，龚学艺，等. 基于面阵探测系统的扫描成像信息获取方式研究[J]. 红外技术，2014(6):479-484. http://d.old.wanfangdata.com.cn/Periodical/zggxyyygxwz201501013 Yang Yuzhou, Zhang Shuanglei, Gong Xueyi, et al. Research of scanning imaging method with focal plane array[J]. Infrared Technology, 2014(6):479-484. http://d.old.wanfangdata.com.cn/Periodical/zggxyyygxwz201501013
-
期刊类型引用(8)
1. 习啸天，韩军，张岩. 一种大视场高分辨率的复眼光学系统设计. 应用光学. 2024(02): 314-320 . 
本站查看
2. 陈炳旭，杨旭，张智强，栾晓宇，蒋锐. 宽波段折反射式航空相机设计. 红外技术. 2024(08): 864-871 . 百度学术
3. 王跃，刘海英，陈晓宇，朱海滨，杨旭，姜彦名，王英. 基于三线靶标的航空相机地面分辨率质量评价. 光电技术应用. 2024(05): 34-38 . 百度学术
4. 尹云飞. 航空相机像旋补偿双向控制中的内环补偿器设计. 机械管理开发. 2024(12): 99-101+104 . 百度学术
5. 刘金沧，黎治坤，王英谋，李翔，郭雨晴. 摆扫式倾斜影像数据处理关键问题研究. 科学技术创新. 2023(08): 17-21 . 百度学术
6. 吴雄雄，王惠林，宁飞，张文博，刘吉龙，王明超. 大气湍流对航空光电侦察图像质量影响因素分析. 应用光学. 2022(01): 1-8 . 本站查看
7. 张磊，许慧超，王云，付杨，赵越，汤晓峰，石钰，孙勇，张弦，贾金升. 大口径光纤倒像器制备工艺研究. 应用光学. 2021(04): 735-742 . 本站查看
8. 邹甜妹，韦晓孝，吕元，万新军，李雪宸，张昊. 面向光学装调微小距离测量的系统设计. 光学仪器. 2021(05): 48-53 . 百度学术
其他类型引用(5)
计量
- 文章访问数:
- HTML全文浏览量:
- PDF下载量:
- 被引次数: 13
陕公网安备 61011302001501号 