Research progress of narrow-linewidth mid-infrared laser
-
摘要: 3 μm~5 μm中红外激光处于大气窗口波段，对应着众多原子或分子的特征吸收峰，在医疗诊断、大气环境监测、空间通信以及光电对抗等诸多领域具有非常重要的应用价值。在这些应用领域，人们往往要求光源拥有窄谱宽和快速波长调谐功能，而窄谱宽激光具有较小的谱宽、能量集中，是满足这些应用的理想光源。总结了实现窄谱宽3 μm~5 μm中红外激光输出的Fe2+/Cr2+离子掺杂固体激光器和氟化物光纤激光器谱宽压缩技术，以量子级联激光器为例，展示了几种激光稳频的措施，重点阐述了结构紧凑、全固化的中红外光参量振荡器的调谐原理和压缩谱宽所采取的技术，对课题组在窄谱宽光参量振荡器方面的研究工作进行了介绍，并对窄谱宽中红外激光技术的研究前景进行了展望。Abstract: 3 μm~5 μm is one of the transmission windows of the atmosphere, including many atoms or molecules' characteristic absorption peaks. 3 μm~5 μm laser souses have wide applications, such as medical diagnosis, environmental monitoring, spectral analysis, military countermeasure, and so on. In these application fields, light sources are often required to have narrow linewidth and wavelength tuning capability. Narrow-linewidth laser is an ideal light source for these applications due to its small spectrum range and concentrative energy. The technology of narrow-linewidth of 3 μm~5 μm laser were summarized, including the linewidth compression method of the Fe2+/Cr2+ ion-doped solid-state laser and fluoride fiber laser, and several methods of quantum cascade laser frequency stabilization were discussed. The compact and all-solid-state mid-infrared optical parametric oscillator was mainly introduced. And the related research work of our group on the narrow linewidth of an optical parametric oscillator was introduced. Finally, the research prospect of the narrow-linewidth mid-infrared laser has been prospected.
-
Keywords:
- laser technique /
- narrow-linewidth /
- solid-state laser /
- fiber laser /
- optical parametric oscillator
-
0. 引言
利用图像全站仪技术进行自动工业摄影测量在技术上已成为可能。采用图像全站仪进行摄影测量具备高精度、无接触、面测量、自动化及无人为照准误差等特性，在天文测量、变形测量及精密工程测量等领域有广泛应用前景[1-7]。高精度的相机标定是图像全站仪进行精密摄影测量一个基础而又不可或缺的问题，也是目前图像全站仪及其应用亟待解决的关键技术之一，但目前相关的研究鲜见报道，相近的研究也主要集中在天文测量和变形测量领域，且只针对望远镜相机定焦测量展开研究，未对广角相机展开研究[8-14]。
EHRHART等人在使用徕卡MS50图像全站仪进行工程结构动态监测时，对该望远镜相机检校进行了研究，认为对于变形监测而言，可以不考虑相机的像主点坐标和畸变，只需要相机的角分辨率即可[11]。国内信息工程大学天文课题组在使用TS50i图像全站仪进行天文测量时，对相机的像平面倾斜、像主点坐标、主距及畸变系数等进行了检校[12]，提出的检校方法简单且方便操作，但相机成像模型并不严密，未考虑相机在全站仪设备中的偏心元素(包括3个位置偏移量x0, y0, z0和3个角度偏移量α, β, γ)，其本质是将全站仪坐标系统等价于相机坐标系统，这会引起相机标定误差。对于望远镜相机而言，相机与全站仪望远镜同轴安置，位置偏移量级为亚毫米级，角度偏移量级为角秒级，全站仪仪器中心近似为相机光心对相机标定结果影响较小，但对于广角相机而言，相机与全站仪望远镜显著偏心安置，仪器中心与相机光心的位置偏差为厘米级，将仪器中心当作相机光心会显著引起相机标定误差，使得相机标定参数在不同摄影距离下标定解算的结果不具备一致性，无法实现图像测量和全站仪角度测量的精确统一，从而限制了图像全站仪在精密工程测量领域中的应用。综上所述，对图像全站仪相机科学的标定方法展开研究具有应用意义。
偏心元素的求解是图像全站仪相机标定独有的问题，由于图像全站仪相机与全站仪刚性连接，偏心元素的各项数值均为常数，在一般标定方法中存在参数耦合问题，需用特殊的方法进行参数解耦。ZHANG提出使用检校条，通过多测站近距离旋转摄影的方式，构建成像共线方程，采用摄影测量前方交会的方式完成偏心元素检校[15]。ZHU提出相机参数逐级标定法，利用全站仪做小角旋转单独标定相机主距，将之作为已知量逐步求解其它参数，避免参数之间的耦合[16]。两者均从理论上指出图像全站仪可通过旋转摄影的方式构建标定场，但缺乏实验来验证旋转摄影标定方法的有效性。
本文作者基于旋转摄影理论，采用图像全站仪单目标旋转摄影标定的方法，推导并建立了图像全站仪相机标定模型，同时，开展实验对模型在不同摄影距离上进行了相机标定，均获得了高精度的标定结果。
1. 单目标旋转摄影标定方法
本文中提出的单目标旋转摄影标定方法主要是指首先通过计算机程序控制全站仪相机对单个靶标进行扫描式摄影，然后利用相机摄影时的姿态角及对应靶标影像的中心像素坐标来构建相机成像方程，最后利用构建的成像方程实现对相机标定参数的求解。该方法不仅充分发挥了图像全站仪平台特性，同时对目标对象的布设以及标定过程的自动化程度方面，相对于传统的方法，具备无可比拟的优势。
单目标旋转摄影标定方法中，人们借助图像全站仪的GEOCOM二次开发接口[17]，通过编写计算机程序实现靶标图像自动化数据采集，整个测量和拍摄作业自动化程度高。具体流程是：通过计算机控制指令驱动望远镜在水平和垂直方向上指向规划的角度并记录当前望远镜姿态，相机与望远镜固联，随全站仪旋转驱动；摄影标定目标不动，仪器每次旋转规划的角度后，对靶标进行拍摄，使目标近似均匀的成像在电荷耦合器件(charge-coupled device, CCD)不同区域。该过程可等价成相机不动而靶标在以相机为中心的球面上移动[18]，目标的每一个姿态角均可等价于一个方向控制点，且控制点精度仅取决于测角精度。在几十米范围内，图像全站仪测角精度可达到0.5″，远高于测距精度，这使得旋转摄影标定方法获得的精度比控制网标定方法更高。采用单靶标进行拍照，可将靶标做得相对较大，有利于提升靶标中心提取精度。用旋转摄影的方法进行点位扩展，相当于在空间中布设高精度的异面控制点，能近乎零成本和更高效地建立高精度且视域内均匀分布的空间标定场。
2. 单目标旋转摄影标定的理论模型
2.1 坐标系
相机标定过程中，涉及到不同坐标系间的坐标转换，主要有全站仪坐标系、望远镜坐标系、相机坐标系以及像素坐标系。为保持与全站仪观测量的一致性，均采用了与全站仪坐标系一致的左手坐标系。
(a) 全站仪坐标系(total station，TS)OTS-XTSYTSZTS，原点为全站仪中心，ZTS轴取全站仪竖轴，XTS轴平行水平度盘零方向，YTS轴由左手系确定。
(b) 望远镜坐标系(telescope，T)OT-XTYTZT，原点取全站仪中心，ZT轴取全站仪望远镜视准轴，YT轴取全站仪横轴，XT轴由左手系确定。当望远镜坐标系处于置零状态时与全站仪坐标系一致。
(c) 相机坐标系(camera，C)OC-XCYCZC，原点为相机光心，ZC轴取相机主光轴，XC轴与像素坐标系XP轴平行且同方向，YC轴由左手系确定。
(d) 像素坐标系(pixel，P)OP-XPYP，建立在相机坐标系ZC=-f的2维平面坐标系，f为相机等效主距。原点为图像左上角，取行方向为XP轴，列方向为YP轴。
图像全站仪相机与全站仪刚性连接，偏心元素为常量，相机坐标系与望远镜坐标系的空间关系如图 1所示。各坐标系的空间关系示意图如图 2所示，其中相机坐标系的ZC轴延长线与像素坐标系平面的交点为像主点(u0, v0)。
![图 1 T坐标系与C坐标系空间关系图]()
图 1 T坐标系与C坐标系空间关系图Figure 1. Spatial relationship between T-coordinate system and C-coordinate system2.2 相机成像模型
首先，建立T坐标系至S坐标系的位姿转换。设全站仪照准靶标P时，全站仪水平角、天顶角和距离读数分别为θ0, φ0, s，P点在全站仪坐标系下的坐标为：
\boldsymbol{P}_{\mathrm{TS}}=s\left[\begin{array}{c} \sin \varphi_{0} \cos \theta_{0} \\ \sin \varphi_{0} \sin \theta_{0} \\ \cos \varphi_{0} \end{array}\right] (1) 驱动全站仪旋转到某一规划指向，此时全站仪水平角、天顶角读数为分别为θ, φ。根据TS坐标系与T坐标系的空间旋转关系可知，TS坐标系绕ZTS轴旋转θ，再绕YTS轴旋转φ，可实现TS坐标系到T坐标系的转换，即：
\begin{equation*} \boldsymbol{P}_{\mathrm{T}}=\boldsymbol{R}_{Y_{\mathrm{TS}}}(\varphi) \cdot \boldsymbol{R}_{Z_{\mathrm{TS}}}(\theta) \cdot \boldsymbol{X}_{\mathrm{TS}} \end{equation*} (2) 式中: RZTS(θ)、RYTS(φ)为对应旋转矩阵，可表示为：
\begin{align*} & \boldsymbol{R}_{Z_{\mathrm{TS}}}(\theta)=\left[\begin{array}{ccc} \cos \theta & \sin \theta & 0 \\ -\sin \theta & \cos \theta & 0 \\ 0 & 0 & 1 \end{array}\right] \end{align*} (3) \begin{align*} & \boldsymbol{R}_{Y_{\mathrm{TS}}}(\varphi)=\left[\begin{array}{ccc} \cos \varphi & 0 & -\sin \varphi \\ 0 & 1 & 0 \\ \sin \varphi & 0 & \cos \varphi \end{array}\right] \end{align*} (4) 式(2)可表示为：
\begin{equation*} \boldsymbol{P}_{\mathrm{T}}=s \cdot \boldsymbol{n}_{\mathrm{T}} \end{equation*} (5) 式中: nT为靶标中心在望远镜坐标系下的单位方向向量，其计算表达式为：
\begin{gather*} \boldsymbol{n}_{\mathrm{T}}=\left[\begin{array}{l} s_{x} \\ s_{y} \\ s_{z} \end{array}\right]= \\ {\left[\begin{array}{c} \cos \varphi \sin \varphi_{0} \cos \left(\theta_{0}-\theta\right)-\sin \varphi \cos \varphi_{0} \\ \sin \varphi_{0} \sin \left(\theta_{0}-\theta\right) \\ \sin \varphi \sin \varphi_{0} \cos \left(\theta_{0}-\theta\right)+\cos \varphi \cos \varphi_{0} \end{array}\right]} \end{gather*} (6) 其次，建立T坐标系至C坐标系的位姿转换，根据T坐标系与C坐标系的空间关系有：
\boldsymbol{P}_{\mathrm{C}}=\boldsymbol{R}_{X_{\mathrm{T}}}(\alpha) \cdot \boldsymbol{R}_{Y_{\mathrm{T}}}(\beta) \cdot \boldsymbol{R}_{Z_{\mathrm{T}}}(\gamma) \cdot \boldsymbol{P}_{\mathrm{T}}+\left[\begin{array}{l} x_{0} \\ y_{0} \\ z_{0} \end{array}\right] $$ (7) 式中: α, β, γ为T坐标系与C坐标系的角度偏移量; x0, y0, z0为位置偏移量。
考虑到无论是广角相机或是望远镜相机，C坐标系至T坐标系的空间转换均为小角度旋转，可将旋转矩阵线性化并取一次项近似代替，有：
\boldsymbol{P}_{\mathrm{C}}=\left[\begin{array}{ccc} 1 & \gamma & -\beta \\ -\gamma & 1 & \alpha \\ \beta & -\alpha & 1 \end{array}\right] \boldsymbol{P}_{\mathrm{T}}+\left[\begin{array}{l} x_{0} \\ y_{0} \\ z_{0} \end{array}\right] (8) 式(8)也可写为：
\begin{equation*} \boldsymbol{P}_{\mathrm{C}}=s \cdot \boldsymbol{n}_{\mathrm{C}} \end{equation*} (9) 式中: nC为靶标中心在相机坐标系下的单位方向向量, 其计算表达式为：
\boldsymbol{n}_{\mathrm{C}}=\left[\begin{array}{c} c_{x} \\ c_{y} \\ c_{z} \end{array}\right]=\left[\begin{array}{c} s_{x}+\gamma s_{y}-\beta s_{z}+x_{0} / s \\ -\gamma s_{x}+s_{y}+\alpha s_{z}+y_{0} / s \\ \beta s_{x}-\alpha s_{y}+s_{z}+z_{0} / s \end{array}\right] (10) 最后，建立C坐标系至P坐标系的坐标转换，由相机成像共线方程有：
\begin{equation*} \frac{x-u_{0}}{c_{x}}=\frac{y-v_{0}}{c_{y}}=\frac{f}{c_{z}} \end{equation*} (11) 由式(11)可得图像全站仪成像模型为：
\left\{\begin{array}{l} x=f_{x} \frac{c_{x}}{c_{z}}+u_{0} \\ y=f_{y} \frac{c_{y}}{c_{z}}+v_{0} \end{array}\right. (12) 式中: x, y为靶标中心像素坐标; fx, fy为相机等效主距; u0, v0为像主点坐标。式(12)也可表示为：
\left\{\begin{array}{l} x=f_{x} \frac{s_{x}+\gamma s_{y}-\beta s_{z}+x_{0} / s}{\beta s_{x}-\alpha s_{y}+s_{z}+z_{0} / s}+u_{0} \\ y=f_{y} \frac{-\gamma s_{x}+s_{y}+\alpha s_{z}+y_{0} / s}{\beta s_{x}-\alpha s_{y}+s_{z}+z_{0} / s}+v_{0} \end{array}\right. (13) 2.3 相机标定方程
对式(13)进行线性化，可得线性化观测误差方程为：
\begin{equation*} \boldsymbol{V}=\boldsymbol{A} \boldsymbol{X}-\boldsymbol{L} \end{equation*} (14) 式中: V为观测残差矩阵; A为参数系数矩阵; X为待解的相机参数; L为观测矩阵。其中：
\begin{gather*} \boldsymbol{X}=\left[u_{0}, v_{0}, f_{x}, f_{y}, \alpha, \beta, \gamma, x_{0}, y_{0}, z_{0}\right]^{\mathrm{T}} \end{gather*} (15) \boldsymbol{L}=\left[\begin{array}{l} x \\ y \end{array}\right] (16) 相机参数系数矩阵的计算表达式为：
\boldsymbol{A}=\left[\begin{array}{l} \frac{\partial x}{\partial u_{0}} \frac{\partial x}{\partial v_{0}} \frac{\partial x}{\partial f_{x}} \frac{\partial x}{\partial f_{y}} \frac{\partial x}{\partial \alpha} \frac{\partial x}{\partial \beta} \frac{\partial x}{\partial \gamma} \frac{\partial x}{\partial x_{0}} \frac{\partial x}{\partial y_{0}} \frac{\partial x}{\partial z_{0}} \\ \frac{\partial y}{\partial u_{0}} \frac{\partial y}{\partial v_{0}} \frac{\partial y}{\partial f_{x}} \frac{\partial y}{\partial f_{y}} \frac{\partial y}{\partial \alpha} \frac{\partial y}{\partial \beta} \frac{\partial y}{\partial \gamma} \frac{\partial y}{\partial x_{0}} \frac{\partial y}{\partial y_{0}} \frac{\partial y}{\partial z_{0}} \end{array}\right] (17) 对式(17)求导取1阶项，简化整理可得：
\mathit{\boldsymbol{A = }}\left[ \begin{array}{l} 1\quad 0\quad\frac{{{c_x}}}{{{c_z}}}\quad 0\quad \quad \left( {x - {u_0}} \right)\frac{{{s_y}}}{{{c_z}}}\quad \quad \left( {{u_0} - x} \right)\left( {\frac{{{s_x}}}{{{c_z}}} + \frac{{{s_z}}}{{{c_x}}}} \right)\quad \left( {x - {u_0}} \right)\frac{{{s_y}}}{{{c_x}}}\quad \frac{{x - {u_0}}}{{s{c_x}}}\quad 0\quad \frac{{{u_0} - x}}{{s{c_z}}}\\ 0\quad 1\quad 0\quad \frac{{{c_y}}}{{{c_z}}}\quad \left( {{v_0} - y} \right)\left( {\frac{{{s_z}}}{{{c_y}}} + \frac{{{s_y}}}{{{c_z}}}} \right)\quad \quad \left( {y - {v_0}} \right)\frac{{{s_x}}}{{{c_z}}}\quad\quad \left( {{v_0} - y} \right)\frac{{{s_x}}}{{{c_y}}}\quad \quad 0\quad \frac{{y - {v_0}}}{{s{c_y}}}\quad \frac{{{v_0} - y}}{{s{c_z}}} \end{array} \right] (18) 采用最小二乘方法求解相机参数，有：
\begin{equation*} \boldsymbol{X}=\left(\boldsymbol{A}^{\mathrm{T}} \boldsymbol{P} \boldsymbol{A}\right)^{-1} \boldsymbol{A} \boldsymbol{P} \boldsymbol{L} \end{equation*} (19) 式中: P为观测量的权矩阵。由于标定参数间存在的强相关性，直接组建标定方程求解时标定方程可能奇异，无法正确求解，应当分析变量间的相关性，解决参数强相关性问题。
2.4 相关性分析
通过参数之间的偏导数关系，对参数的相关性展开分析，具体为：
\begin{gather*} {\left[\begin{array}{llll} \frac{\partial \beta}{\partial x_{0}} & \frac{\partial \beta}{\partial u_{0}} & \frac{\partial x_{0}}{\partial u_{0}} & \frac{\partial z_{0}}{\partial f_{x}} \\ \frac{\partial \alpha}{\partial y_{0}} & \frac{\partial \alpha}{\partial v_{0}} & \frac{\partial y_{0}}{\partial v_{0}} & \frac{\partial z_{0}}{\partial f_{y}} \end{array}\right]=} \\ {\left[\begin{array}{llll} \frac{\partial x / \partial x_{0}}{\partial x / \partial \beta} & \frac{\partial x / \partial u_{0}}{\partial x / \partial \beta} & \frac{\partial x / \partial u_{0}}{\partial x / \partial x_{0}} & \frac{\partial x / \partial f_{x}}{\partial x / \partial z_{0}} \\ \frac{\partial y / \partial y_{0}}{\partial y / \partial \alpha} & \frac{\partial y / \partial v_{0}}{\partial y / \partial \alpha} & \frac{\partial y / \partial v_{0}}{\partial y / \partial y_{0}} & \frac{\partial y / \partial f_{y}}{\partial y / \partial z_{0}} \end{array}\right]} \end{gather*} (20) 经计算整理有：
\begin{array}{c} {\left[\begin{array}{llll} \frac{\partial \beta}{\partial x_{0}} & \frac{\partial \beta}{\partial u_{0}} & \frac{\partial x_{0}}{\partial u_{0}} & \frac{\partial z_{0}}{\partial f_{x}} \\ \frac{\partial \alpha}{\partial y_{0}} & \frac{\partial \alpha}{\partial v_{0}} & \frac{\partial y_{0}}{\partial v_{0}} & \frac{\partial z_{0}}{\partial f_{y}} \end{array}\right]=} \\ {\left[\begin{array}{cccc} \frac{-c_{z} / s}{c_{x} s_{x}+c_{z} s_{z}} & \frac{-c_{z}^{2} / f_{x}}{c_{x} s_{x}+c_{z} s_{z}} & \frac{s c_{z}}{f_{x}} & \frac{s c_{z}}{-f_{x}} \\ \frac{c_{z} / s}{c_{y} s_{y}+c_{z} s_{z}} & \frac{c_{z}^{2} / f_{y}}{c_{y} s_{y}+c_{z} s_{z}} & \frac{s c_{z}}{f_{y}} & \frac{s c_{z}}{-f_{y}} \end{array}\right]} \end{array} (21) 式中: 相机主距fx和fy可视作常数，偏心元素的各项数值均为小量常数，可计算出sz, cz的值趋于1，远大于sx, sy, cx, cy。当考虑单距离旋转摄影标定时，即摄影距离s为常量时，式(21)中各项导数计算值均近似为常数，即β, x0, u0强相关，α、y0、v0强相关，fx, fy, z0强相关，故通过单距离观测进行相机标定时，仅可考虑变量u0, v0, fx, fy, γ进行标定解算，无法解算出参数α, β, x0, y0, z0。
当考虑多距离旋转摄影标定时，即摄影距离s为变量时，式(21)中∂β/∂u0, ∂α/∂v0的值近似不变，其余项的计算值线性变化，即x0, y0, z0与摄影距离s线性相关，故通过多距离观测整体解算可解耦x0, y0, z0，但无法解算α, β。经分析认为，需改变相机相对全站仪的位置才能解耦α, β，这可以通过切换全站仪度盘来改变相机相对全站仪的位置，采用双盘观测的方式解耦α, β。
2.5 双盘观测方程
将望远镜绕水平方向旋转180°，再绕竖直方向旋转180°可实现全站仪度盘切换。切换观测度盘会改变相机相对全站仪的位置，假设在盘左观测中偏心元素为x0, y0, z0, α, β, γ，在盘右观测中偏心元素则为-x0, -y0, z0, -α, -β, γ。令盘左观测时靶标像素坐标为(xl, yl)，盘右观测时靶标像素坐标为(xr, yr)。盘右图像成正像显示，观测数据需作镜像处理，将式(12)对应参数进行修正，可得盘右观测时对应的成像方程为：
\left\{\begin{array}{l} x_{\mathrm{r}}=f_{x} \frac{-s_{x}-\gamma s_{y}-\beta s_{z}+x_{0} / s}{-\beta s_{x}+\alpha s_{y}+s_{z}+z_{0} / s}+\left(R-u_{0}\right) \\ y_{\mathrm{r}}=f_{y} \frac{\gamma s_{x}-s_{y}+\alpha s_{z}+y_{0} / s}{-\beta s_{x}+\alpha s_{y}+s_{z}+z_{0} / s}+\left(C-v_{0}\right) \end{array}\right. (22) 式中: R为图像的行数; C为图像的列数。构建双盘观测方程为：
\left\{\begin{array}{l} x_{1}=f_{x} \frac{s_{x}+\gamma s_{y}-\beta s_{z}+x_{0} / s}{\beta s_{x}-\alpha s_{y}+s_{z}+z_{0} / s}+u_{0} \\ y_{1}=f_{y} \frac{-\gamma s_{x}+s_{y}+\alpha s_{z}+y_{0} / s}{\beta s_{x}-\alpha s_{y}+s_{z}+z_{0} / s}+v_{0} \\ x_{\mathrm{r}}=f_{x} \frac{-s_{x}-\gamma s_{y}-\beta s_{z}+x_{0} / s}{-\beta s_{x}+\alpha s_{y}+s_{z}+z_{0} / s}+\left(R-u_{0}\right) \\ y_{\mathrm{r}}=f_{y} \frac{\gamma s_{x}-s_{y}+\alpha s_{z}+y_{0} / s}{-\beta s_{x}+\alpha s_{y}+s_{z}+z_{0} / s}+\left(C-v_{0}\right) \end{array}\right. (23) 利用式(23)，在多个距离上扫描目标的盘左和盘右成像坐标，通过整体平差解算可求解出全部相机参数。借鉴分步求解思路来降低参数耦合的影响，先假设参数α、β值为0°，通过整体平差解出参数fx, fy, u0, v0, γ, x0, y0, z0，此时α、β耦合在u0、v0中，后通过某种方式单独解算α、β，再对u0、v0进行修正。
由式(23)可构建双盘观测量：
\left\{\begin{align*} \frac{x_{1}+x_{\mathrm{r}}}{2}= & \frac{f_{x}}{2}\left(\frac{s_{x}+\gamma s_{y}-\beta s_{z}+x_{0} / s}{\beta s_{x}-\alpha s_{y}+s_{z}+z_{0} / s}+\right. \\ & \left.\frac{-s_{x}-\gamma s_{y}-\beta s_{z}+x_{0} / s}{-\beta s_{x}+\alpha s_{y}+s_{z}+z_{0} / s}\right)+\frac{R}{2} \\ \frac{y_{1}+y_{\mathrm{r}}}{2}= & \frac{f_{y}}{2}\left(\frac{-\gamma s_{x}+s_{y}+\alpha s_{z}+y_{0} / s}{\beta s_{x}-\alpha s_{y}+s_{z}+z_{0} / s}+\right. \\ & \left.\frac{\gamma s_{x}-s_{y}+\alpha s_{z}+y_{0} / s}{-\beta s_{x}+\alpha s_{y}+s_{z}+z_{0} / s}\right)+\frac{C}{2} \end{align*}\right. (24) 忽略小量βsx, αsy，整理化简可得：
\left\{\begin{array}{l} \frac{x_{1}+x_{\mathrm{r}}}{2}=-\beta f_{x}+f_{x} \frac{x_{0}}{s s_{z}}+\frac{R}{2} \\ \frac{y_{1}+y_{\mathrm{r}}}{2}=\alpha f_{y}+f_{y} \frac{y_{0}}{s s_{z}}+\frac{C}{2} \end{array}\right. (25) 可得到α, β的计算式为:
\left\{\begin{array}{l} \beta=\frac{x_{0}}{s s_{z}+z_{0}}+\frac{R}{2 f_{x}}-\frac{x_{1}+x_{\mathrm{r}}}{2 f_{x}} \\ \alpha=-\frac{y_{0}}{s s_{z}+z_{0}}-\frac{C}{2 f_{y}}+\frac{y_{1}+y_{\mathrm{r}}}{2 f_{y}} \end{array}\right. (26) 在解得α, β的基础上，对像主点坐标做相应修正，可获得实际像主点坐标值u0, v0：
\left\{\begin{array}{l} \bar{u}_{0}=u_{0}-f_{x} \beta \\ \bar{v}_{0}=v_{0}+f_{y} \alpha \end{array}\right. (27) 2.6 相机畸变模型
由于镜头设计的复杂性和工艺水平等因素的影响，镜头折射率不一致，光学器件存在偏移，导致实际成像点(x′, y′)与理论像点(x, y)之间存在像素偏差(δx, δy)，即：
\left[\begin{array}{l} x \\ y \end{array}\right]=\left[\begin{array}{l} x^{\prime} \\ y^{\prime} \end{array}\right]+\left[\begin{array}{l} \delta_{x} \\ \delta_{y} \end{array}\right] (28) Brown-Conrady畸变理论认为存在3种畸变类型，分别为径向畸变、偏心畸变和薄棱镜畸变[19]。同时考虑这3种畸变，相机畸变模型可表示为：
\left\{\begin{align*} \delta_{x}= & \left(x-u_{0}\right)\left(k_{1} r^{2}+k_{2} r^{4}+k_{3} r^{6}\right)+ \\ & 2 \rho_{1}\left(x-u_{0}\right)\left(y-v_{0}\right)+\rho_{2}\left(r^{2}+2 x^{2}\right)+s_{1} r^{2} \\ \delta_{y}= & \left(y-v_{0}\right)\left(k_{1} r^{2}+k_{2} r^{4}+k_{3} r^{6}\right)+ \\ & 2 \rho_{2}\left(x-u_{0}\right)\left(y-v_{0}\right)+\rho_{1}\left(r^{2}+2 y^{2}\right)+s_{2} r^{2} \end{align*}\right. (29) 式中: r=\sqrt{\left(x-u_0\right)^2+\left(y-v_0\right)^2}，表示像点到透视中心的像素距离；k1, k2, k3为径向畸变系数；ρ1和ρ2为偏心畸变系数；s1和s2为薄棱镜畸变系数。
3. 标定实验及分析
3.1 实验描述
基于上述理论，本文作者对TS60图像全站仪进行相机标定实验。TS60同时具有望远镜相机和广角相机，望远镜相机为变焦相机，视场角为1.5°，广角相机为定焦相机，视场角为19.4°，成像规格为2560 pixel×1920 pixel(R×C)。标定靶标采用同心圆环靶标，设计有十字丝便于全站仪初始精确照准。设计在8 m、12 m、16 m、20 m、24 m、30 m的距离进行望远镜相机和广角相机的图像采集。设定合适的摄影间隔使得标定靶标在视域内均匀分布成像，本文中设计旋转摄影照片数量为26×24，对应广角相机旋转摄影方位角和高度角摄影间隔分别为35.77′和29.25′，对应望远镜相机旋转摄影方位角和高度角摄影间隔分别为2.77′和2.25′，在单个距离采集624张靶标图像。基于C#实现计算机与全站仪的蓝牙通讯，参照手册开发计算机程序控制全站仪相机驱动到规划的角度姿态对靶标进行拍照，实现靶标图像自动化数据采集。实验室现场环境如图 3所示。望远镜相机获得的部分靶标图像集如图 4所示。
相机标定流程如下：(a)通过中心提取算法[20-21]获取靶标中心像素坐标(x, y)，并记录对应的摄影方向(θ, φ)；(b)仅考虑参数u0, v0, fx, fy, γ进行解算，作为后续求解的初值；(c)采用多距离双盘观测整体解算模型解算变量u0, v0, fx, fy, γ, x0, y0, z0；(d)构建双盘观测量解算α, β，在此基础上对像主点坐标进行改正；(e)广角相机考虑相机畸变参数k1、k2，望远镜相机考虑畸变参数k1，求解相机畸变参数；(f)通过畸变参数对靶标像素坐标做相应改正，重复步骤(b)~(f)，直至求解的参数变化量小于容许值，获得最终相机标定参数。
3.2 实验结果及分析
表 1中为采用单个距离的望远镜相机盘左和盘右观测数据进行相机标定的结果。表中, 全站仪观测度盘面值(face)为0表示盘左观测，度盘面值为1表示盘右观测，mr表示平均重投影点位误差，为中心提取算法获取的靶标中心像素坐标与通过模型计算得到的靶标中心像素坐标的平均差值。表 2中为联合所有观测距离的望远镜相机盘左和盘右观测数据分别进行标定解算的结果。
表 1 望远镜相机单距离观测标定解算Table 1. Calibration solution of telescope camera for single distance observationface s/m u0/pixel v0/pixel fx /pixel fy /pixel γ /(″) mr/pixel 0 8.11 962.47 1280.89 105588.46 105588.45 487.45 0.12 1 8.11 958.52 1279.37 105565.96 105592.29 491.36 0.10 0 12.27 962.05 1281.44 105484.74 105526.00 490.01 0.10 1 12.27 961.19 1280.47 105492.32 105529.98 478.04 0.12 0 15.83 962.02 1281.02 105512.10 105550.65 492.63 0.08 1 15.83 958.84 1278.62 105512.43 105542.27 477.48 0.07 0 19.85 962.37 1283.05 105522.13 105499.94 461.41 0.37 1 19.85 959.45 1279.93 105501.19 105508.14 489.65 0.14 0 23.93 960.35 1283.33 105465.85 105487.35 481.27 0.12 1 23.93 961.68 1280.78 105471.87 105584.78 522.80 0.18 0 29.70 961.57 1283.24 105462.67 105444.87 475.50 0.26 1 29.70 958.84 1277.76 105447.33 105467.91 450.16 0.19 表 2 望远镜相机多距离观测标定解算Table 2. Calibration solution of telescope camera for multiple distance observationface u0/pixel v0/pixel fx/pixel fy/pixel α/(″) β/(″) γ/(″) x0/mm y0/mm z0/mm k1/10-10 mr/pixel 0 960.63 1282.91 105392.33 105382.56 -6.30 2.30 496.28 0.16 -0.15 -15.77 -3.13 0.88 1 959.46 1279.69 105412.70 105452.90 -6.30 2.30 493.64 0.01 -0.01 -12.13 -4.76 1.15 由表 1与表 2可知：望远镜相机不同距离同盘观测的像主点坐标较为稳定，同一距离不同度盘观测的像主点坐标平均差异为2 pixel；单距离观测相机主距有随摄影距离增加而减小的趋势，与多距离观测主距最大差异为205 pixel，对应引起图像边缘的角度最大变化约7.0″；单盘解算平均重投影误差为0.15 pixel，对应图像分辨率0.3″，多距离解算平均重投影误差为1.02 pixel，对应的图像角分辨率约2.0″。
表 3中为采用单个距离的广角相机盘左和盘右观测数据进行相机标定的结果。表 4中为联合所有观测距离的广角相机盘左和盘右观测数据分别进行标定解算的结果。由表 3与表 4可知：广角相机单距离观测的像主点坐标存在明显变化，其中u0分量有随摄影距离增加而减小的趋势，v0分量有随摄影距离增加而增大的趋势，体现了偏心元素与相机参数耦合对相机标定的影响；单距离观测相机主距与多距离观测主距最大差异为59 pixel，对应引起图像边缘的角度最大变化约23″；单盘解算平均重投影误差为0.72 pixel，对应图像分辨率15.5″，多距离解算平均重投影误差为0.6 pixel，对应的图像角分辨率约13.3″。
表 3 广角相机单距离观测标定解算Table 3. Calibration solution of wide-angle camera for single distance observationface s/m u0/pixel v0/pixel fx/pixel fy/pixel γ/(″) mr/pixel 0 8.11 1033.24 1267.07 9438.46 9433.63 -64.21 0.91 1 8.11 1030.93 1265.02 9439.43 9434.14 -71.99 0.86 0 12.27 1009.08 1273.66 9414.37 9409.81 -37.40 0.92 1 12.27 1007.02 1271.27 9414.72 9410.60 -46.64 0.88 0 15.83 998.50 1276.43 9407.02 9402.95 -46.28 0.76 1 15.83 996.38 1273.87 9405.88 9401.83 -52.03 0.80 0 19.85 991.33 1278.15 9405.75 9402.92 -55.70 0.57 1 19.85 989.26 1276.18 9407.07 9401.63 -51.05 0.61 0 23.93 986.43 1279.06 9412.55 9411.97 -15.19 0.61 1 23.93 984.45 1277.25 9409.08 9406.26 -49.37 0.47 0 29.70 981.64 1280.41 9396.30 9401.08 -19.98 0.80 1 29.70 979.66 1278.31 9400.80 9401.36 -41.03 0.48 表 4 广角相机多距离观测标定解算Table 4. Calibration solution of wide-angle camera for single distance observationface u0/pixel v0/pixel fx/pixel fy/pixel α/(″) β/(″) γ/(″) x0/mm y0/mm z0/mm k1/10-9 k2/10-15 mr/pixel 0 962.48 1285.71 9380.27 9379.64 -50.16 40.97 -45.15 60.58 -15.88 -48.09 3.33 -3.26 0.61 1 960.61 1283.43 9384.54 9383.22 -50.16 40.97 -50.53 60.22 -15.88 -45.34 2.32 -2.32 0.58 望远镜相机多距离观测数据进行整体观测的重投影误差图如图 5所示。广角相机多距离观测数据进行整体观测的重投影误差图如图 6所示。图 5、图 6为利用相机多距离观测标定结果作为相机参数计算的靶标理论像素坐标与中心提取算法获得的靶标中心像素坐标的误差矢量图。从图中可发现，图像中心视域重投影误差较小，视域边缘重投影误差相对较大，说明相机中心视域内的图像测量精度更高。
![图 5 望远镜相机整体解算重投影误差图]()
图 5 望远镜相机整体解算重投影误差图Figure 5. Overall solution of reprojection error diagram of telescope camera![图 6 广角相机整体解算重投影误差图]()
图 6 广角相机整体解算重投影误差图Figure 6. Overall solution of reprojection error diagram of wide-angle camera从实验结果来看，有：(a)望远镜单距离简化模型中，望远镜相机可视作同轴安置，不同摄影距离相机像主点位置较为稳定，变化量主要体现在分量引起的主距变化；广角相机应视为偏心安置，偏心元素对相机标定解算的影响极大，主要表现为分量引起相机像主点位置在不同摄影距离下明显变化；(b)通过多距离双盘观测联合解算模型可以有效解决标定参数的耦合问题，由此得到的相机参数适用于不同摄影距离的图像测量，其中望远镜相机多距离变焦摄影图像的角度分辨率可达2.0″，单个距离定焦摄影时的角度分辨率可达0.3″；广角相机多距离摄影图像角度分辨率达13.3″；只利用相机的中心视域进行图像测量有更高的角度分辨率，相机标定结果可满足高精度的图像测量要求；(c)望远镜相机多距离多参数模型的标定解算精度低于单距离简化模型标定解算精度，这主要是因为望远镜相机采用了对焦系统，不同距离拍摄过程中需要重新调焦，由此会带来相机焦距和主光轴方向的变化。
在传统相机标定法中，标定场是由全站仪人工测量的系列非共面点的绝对坐标组成，采用这种方式整个控制网的布设过程费时耗力且成本高昂，且控制点的布设还存在通视上的限制；有时囿于加工精度的制约，标定板和标定块不能无限制增大，在广角相机标定时和在较远距离进行相机标定时，标定板和标定块在成像视域内所占图幅较小，控制点难以在成像视域内均匀分布，这会降低传统相机标定解算精度。同时，传统相机标定法要求对控制场进行多距离多视角的摄影，拍摄过程中更多的是使用人工拍摄，拍摄作业的自动化程度低。这些存在的问题，在本文中都得到了克服，很大程度上保证了方法的优势。
4. 结论
采用单目标旋转摄影的图像全站仪相机标定方法，建立了多距离双盘观测的联合解算模型，可以实现广角相机和望远镜相机的精确标定。
(a) 与传统的摄影测量标定方法相比，采用基于旋转摄影的图像全站仪标定方法，无需建立精密的标定控制网，可实现免控制点的精密标定，标定过程可做到自动化。
(b) 多距离双盘观测的联合解算模型能很好地克服参数相关问题，标定解算得到的相机参数可用于不同摄影距离的图像测量，最终望远镜相机图像角度分辨率达2.0″，定焦摄影时图像角度分辨率达0.3″；广角相机图像角度分辨率达13.3″，这个结果优于传统的方法。
(c) 可进一步推广，将相机搭载在全站仪上，采用本文中提出的旋转摄影相机标定方法和解算模型，实现其它相机免控制网的全自动高精度标定。
本文中的研究结果有望提高全站仪的测量精度。
-
![]()
![]()
![]()
![]()
![]()
图 13 窄谱宽MgO ∶PPLN OPO实验原理示意图
Figure 13. Schematic setup of narrow-linewidth MgO ∶PPLN OPO
![]()
表 1 窄谱宽中红外固体激光器研究进展
Table 1 Progress of narrow-linewidth solid laser
年份 晶体 抽运源 输出特性 谱宽 波长/μm 功率/W 能量/mJ 2019[10] Cr∶ZnSe 掺Tm光纤激光 2.455~2.65 0.475 — ≤900 MHz 2019[11] Fe∶ZnSe Ho, Pr∶LLF激光 3.957 — 0.016 23.2 nm 2019[12] Cr∶ZnSe Tm∶YLF激光 2.338~2.572 0.92 — 50 MHz 2018[13] Fe∶ZnSe Er∶YAG激光 4.24 — 8.89 10 nm 2017[14] Fe∶ZnSe 掺Er光纤激光 4.05 >3.5 — < 1 nm 2017[15] Fe∶CdMnTe Fe∶ZnSe激光 5.223 0.81 — 1 nm 2016[16] Cr∶ZnSe 掺Tm光纤激光 2.3~2.7 1.33 — 0.14 nm 2016[17] Cr∶ZnSe Ho∶YAG激光 2.450~2.570 5.5 — < 1 nm 2015[18] Fe∶ZnSe Er∶YAG激光 4.122 0.076 — 6 nm 2015[19] Cr∶ZnSe 掺Tm光纤激光 2.077~2.777 0.120 — < 0.1 nm 
下载: 导出CSV
表 2 3 μm波段窄谱宽光纤振荡器研究进展
Table 2 Progress of narrow-linewidth fiber lasers at 3 μm
年份 光纤基质 输出特性 谱宽 波长/μm 功率/W 2021[22] Er∶ZBLAN 2.7 约0.2 0.4 nm 2017[23] Er∶ZBLAN 2.71~2.82 1.5 1.5 nm 2015[24] Ho, Pr∶ZBLAN 2.825~2.975 7.2 < 0.14 nm 2015[25] Er∶ZBLAN 2.794 0.012 < 20 kHz 2014[26] Er∶ZBLAN 2.8 0.98 0.9 nm 2013[27] Ho, Pr∶ZBLAN 2.914 0.011 < 0.4 nm 2007[28] Er∶ZBLAN 2.7~2.83 2 1.27 GHz
下载: 导出CSV
表 3 窄谱宽中红外OPO研究进展
Table 3 Progress of narrow-linewidth mid-infrared OPO
年份 抽运源 工作物质 输出特性 谱宽 波长/μm 功率/W 2020[37] 光纤激光 MgO∶PPLN 3.7 0.0735 0.271 nm 2018[38] Nd∶YVO4 MgO∶PPLN 2.9~4.1 1.1 — 2017[39] 光纤激光 MgO∶PPLN 3.017 2.5 0.42 nm 2017[40] 掺Yb光纤激光 MgO∶PPLN 2.85~3.05 2.67 0.35 nm 2014[41] Nd∶YAG MgO∶PPLN 2.1 7.1 0.6 nm 2014[42] 掺Yb光纤激光 MgO∶PPLN 3.4 11 < 1.65 MHz 2013[43] Nd∶YAG MgO∶PPLN 3.0~3.5 — < 5.4 GHz 2012[44] 光纤激光 MgO∶PPLN 2.7~4.2 1 80 kHz
下载: 导出CSV
-
[1] 范晋祥. 美国弹道导弹防御系统的红外系统与技术的发展[J]. 红外与激光程, 2006, 35(5): 536-540. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ200605007.htm FAN J X. Status quo and trend of infrared system and technologies for America's ballistic missile defense system[J]. Infrared and Laser Engineering, 2006, 35(5): 536-540(in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ200605007.htm
[2] MOULTON P, DERGACHEV A, ISYANOVA Y, et al. Recent advances in solid state lasers and nonlinear optics for remote sensing[C]// Conference on Lidar Remote Sensing for Industry and Environment Monitoring Ⅲ. Bellingham, USA: International Society for Optical Engineering, 2003: 193-202.
[3] GUO B J, WANG Y, PENG C, et al. Laser-based mid-infrared reflectance imaging of biological tissues[J]. Optics Express, 2004, 12(1): 208-219. DOI: 10.1364/OPEX.12.000208
[4] VAN HERPEN M, TE LINTEL HEKKERT S, BISSON S E, et al. Development of a powerful continuously tunable mid-infrared CW PPLN OPO for trace gas detection[C]// ALT'01 International Conference on Advanced Laser Technologies. Bellingham, USA: International Society for Optical Engineering, 2002: 16-21.
[5] VAINIO M, SILTANEN M, PELTOLA J, et al. Grating-cavity continuous-wave optical parametric oscillators for high-resolution mid-infrared spectroscopy[J]. Applied Optics, 2011, 50(4): A1-A10. DOI: 10.1364/AO.50.0000A1
[6] RICHTER D, FRIED A, WERT B P, et al. Development of a tunable mid-IR difference frequency laser source for highly sensitive airborne trace gas detection[J]. Applied Physics, 2002, B75(2/3): 281-288.
[7] KRZEMPEK K, JAHJAH M, LEWICKI R, et al. CW DFB RT diode laser based sensor for trace-gas detection of ethane using novel compact multipass gas absorption cell[J]. Applied Physics, 2013, B112(4): 461-465.
[8] ELVIN R, HOTH G W, WRIGHT M, et al. Cold-atom clock based on a diffractive optic[J]. Optics Express, 2019, 27(26): 38359-38366. DOI: 10.1364/OE.378632
[9] REN T, WU C, YU Y, et al. Development progress of 3-5 μm mid-infrared lasers: OPO, solid-state and fiber laser[J]. Applied Sciences, 2021, 11(23): 11451. DOI: 10.3390/app112311451
[10] TURNER E J, McDANIEL S A, TABIRYAN N, et al. Rapidly tunable HIP treated Cr ∶ZnSe narrow-linewidth laser[J]. Optics Express, 2019, 27(9): 12282-12288. DOI: 10.1364/OE.27.012282
[11] LI Y Y, JU Y L, DAI T Y, et al. A gain-switched Fe ∶ZnSe laser pumped by a pulsed Ho, Pr: LLF laser[J]. Chinese Physics Letters, 2019, 36(4): 24-26.
[12] WANG Q, LIU C, QI L, et al. Wavelength tunable single-frequency Cr ∶ZnSe laser[C]// 2019 International Conference on Optical Instruments and Technology: Advanced Laser Technology and Applications. Bellingham, USA: International Society for Optical Engineering, 2019: 114370H.
[13] DAI S, FENG G, HONG Z, et al. 4.24 μm mid-infrared laser based on a single Fe2+-doped ZnSe microcrystal[J]. Optics Letters, 2018, 43(3): 411-414. DOI: 10.1364/OL.43.000411
[14] EVANS J W, STITES R W, HARRIS T R. Increasing the performance of an Fe ∶ZnSe laser using a hot isostatic press[J]. Optical Materials Express, 2017, 7(12): 4296-4303. DOI: 10.1364/OME.7.004296
[15] EVANS J W, DOLASINSKO B D, HARRIS T R, et al. Demonstration and power scaling of an Fe ∶CdMnTe laser at 5.2 microns[J]. Optical Materials Express, 2017, 7(3): 860-867. DOI: 10.1364/OME.7.000860
[16] STITES R W, McDANIEL S A, BARNES J O, et al. Hot isostatic pressing of transition metal ions into chalcogenide laser host crystals[J]. Optical Materials Express, 2016, 6(10): 3339-3353. DOI: 10.1364/OME.6.003339
[17] YUAN J H, CHEN Y, YANG H Y, et al. Investigation of a gain-switched Cr2+ ∶ZnSe laser pumped by an acousto-optic Q-switched Ho ∶YAG laser[J]. Quantum Electronics, 2016, 46(9): 772-776. DOI: 10.1070/QEL16157
[18] LANCASTER A, COOK G, McDANIEL S A, et al. Mid-infrared laser emission from Fe ∶ZnSe cladding waveguides[J]. Applied Physics Letters, 2015, 107(3): 885-895.
[19] MAcDONALD J R, BEECHER S J, LANCASTER A, et al. Ultrabroad mid-infrared tunable Cr ∶ZnSe channel waveguide laser[J]. IEEE Journal of Selected Topics in Quantum Electronics, 2014, 21(1): 375-379.
[20] CHEN M, CUI H, LI W, et al. Reparative effect of diffusion process on host defects in Cr2+ doped ZnS/ZnSe[J]. Journal of Alloys and Compounds, 2014, 597(17): 124-128.
[21] McDANIEL S A, BERRY P A, SCHEPLER K L, et al. Gain-switched operation of ultrafast laser inscribed waveguides in Cr ∶ZnSe[C]// Solid State Lasers Ⅹ Ⅹ Ⅳ: Technology and Devices. Bellingham, USA: International Society for Optical Engineering, 2015: 93420E.
[22] SHEN Y, WANG Y, ZHU F, et al. 200 μJ, 13 ns Er ∶ZBLAN mid-infrared fiber laser actively Q-switched by an electro-optic modulator[J]. Optics Letters, 2021, 46(5): 1141-1144. DOI: 10.1364/OL.418950
[23] SHEN Y, WANG Y, LUAN K, et al. High peak power actively Q-switched mid-infrared fiber lasers at 3 μm[J]. Applied Physics, 2017, B123(4): 105-111.
[24] CRAWFORD S, HUDSON D D, JACKSON S D. High-power broadly tunable 3 μm fiber laser for the measurement of optical fiber loss[J]. IEEE Photonics Journal, 2015, 7(3): 1502309.
[25] BERNIER M, MICHAUD-BELLEAU V, LEVASSEUR S, et al. All-fiber DFB laser operating at 2.8 μm[J]. Optics Letters, 2015, 40(1): 81-84. DOI: 10.1364/OL.40.000081
[26] 沈炎龙, 黄珂, 朱峰, 等. LD泵浦瓦级单模高掺铒中红外光纤激光器[J]. 光子学报, 2014, 43(3): 0314002. https://www.cnki.com.cn/Article/CJFDTOTAL-GZXB201403016.htm SHEN Y L, HUANG K, ZHU F, et al. Laser diode-pumped watt-level single mode heavily erbium-doped mid-infrared fiber laser[J]. Acta Photonica Sinica, 2014, 43(3): 0314002(in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-GZXB201403016.htm
[27] HUDSON D D, WILLIAMS R J, WITHFORD M J, et al. Single-frequency fiber laser operating at 2.9 μm[J]. Optics Letters, 2013, 38(14): 2388-2390. DOI: 10.1364/OL.38.002388
[28] ZHU X Sh, JAIN R. Compact 2 W wavelength-tunable Er ∶ZBLAN mid-infrared fiber laser[J]. Optics Letters, 2007, 32(16): 2381-2383. DOI: 10.1364/OL.32.002381
[29] BAYRAKLI I. Frequency-stabilized narrow-linewidth double-mode quantum cascade laser[J]. Optical and Quantum Electronics, 2022, 54(1): 22 (2022).
[30] CAPPELLI F, GALLI I, BORRI S, et al. Subkilohertz linewidth room-temperature mid-infrared quantum cascade laser using a molecular sub-doppler reference[J]. Optics Letters, 2012, 37(23): 4811-4813. DOI: 10.1364/OL.37.004811
[31] SHEHZAD A, BROCHARD P, MATTHEY R, et al. 10 kHz linewidth mid-infrared quantum cascade laser by stabilization to an optical delay line[J]. Optics Letters, 2019, 44(14): 3470-3473. DOI: 10.1364/OL.44.003470
[32] BORRI S, GALLI I, CAPPELLI F, et al. Direct link of a mid-infrared QCL to a frequency comb by optical injection[J]. Optics Letters, 2012, 37(6): 1011-1013. DOI: 10.1364/OL.37.001011
[33] ZHAO B, WANG X, WANG C. Strong optical feedback stabilized quantum cascade laser[J]. ACS Photonics, 2020, 7(5): 1255-1261. DOI: 10.1021/acsphotonics.0c00189
[34] 聂鸿坤, 宁建, 张百涛, 等. 光学超晶格中红外光参量振荡器研究进展[J]. 中国激光, 2021, 48(5): 0501008. https://www.cnki.com.cn/Article/CJFDTOTAL-JJZZ202105009.htm NIE H K, NING J, ZHANG B T, et al. Recent progress of optical-superlattice-based mid-infrared optical parametric oscillators[J]. Chinese Journal of Lasers, 2021, 48(5): 0501008 (in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-JJZZ202105009.htm
[35] WANG X C, WANG Y H, ZHENG H, et al. Wide-tunable mid infrared intra-cavity optical parametric oscillator based on multi-period MgO ∶PPLN[J]. Current Optics and Photonics, 2021, 5(1): 59-65.
[36] RICCIARDI I, MOSCA S, PARISI M, et al. Sub-kHz-linewidth mid-infrared optical parametric oscillator[C]// Conference on Lasers and Electro-Optics. New York, USA: IEEE, 2014: STh1N. 3.
[37] FENG J, CHENG X, LI X, et al. Highly efficient mid-infrared generation from low-power single-frequency fiber laser using phase-matched intracavity difference frequency mixing[J]. Applied Sciences-Basel, 2020, 10(21): 7454-7461. DOI: 10.3390/app10217454
[38] ZHAO J, CHENG P, XU F, et al. Watt-level continuous-wave single-frequency mid-infrared optical parametric oscillator based on MgO ∶PPLN at 3.68 μm[J]. Applied Sciences-Basel, 2018, 8(8): 1345-1352. DOI: 10.3390/app8081345
[39] 邢廷伦, 王礼, 胡舒武, 等. 3μm低阈值MgO ∶PPLN-OPO布拉格体光栅腔谱宽压窄研究[J]. 中国激光, 2017, 44(1): 10101006. https://www.cnki.com.cn/Article/CJFDTOTAL-JJZZ201701024.htm XING Y L, WANG L, HU Sh W, et al. Cavity-linewidth narrowing of 3 μm low threshold MgO ∶PPLN-OPO by volume Bragg grating[J]. Chinese Journal of Lasers, 2017, 44(1): 0101006 (in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-JJZZ201701024.htm
[40] XING T, WANG L, HU S, et al. Widely tunable and narrow-bandwidth pulsed mid-IR PPMgLN-OPO by self-seeding dual etalon-coupled cavities[J]. Optics Express, 2017, 25(25): 31810-31815. DOI: 10.1364/OE.25.031810
[41] JIAO Z, GUO J, HE G, et al. Narrowband intracavity MgO ∶PPLN optical parametric oscillator near degeneracy with a volume Bragg grating[J]. Optics and Laser Technology, 2014, 56: 230-233. DOI: 10.1016/j.optlastec.2013.08.023
[42] ZEIL P, THILMANN N, PASISKEVICIUS V, et al. High-power, single-frequency, continuous-wave optical parametric oscillator employing a variable reflectivity volume Bragg grating[J]. Optics Express, 2014, 22(24): 29907-29913. DOI: 10.1364/OE.22.029907
[43] DOLASINSKI B, POWERS P. Narrow bandwidth tunable optical parametric generator[C]// Nonlinear Frequency Generation and Conversion. Bellingham, USA: International Society for Optical Engineering, 2013: 8604H.
[44] RICCIARDI I, DE TOMMASI E, MADDALONI P, et al. A narrow-linewidth optical parametric oscillator for mid-infrared high-resolution spectroscopy[J]. Molecular Physics, 2012, 110(17): 2103-2109. DOI: 10.1080/00268976.2012.699640
[45] JACOBSSON B, CANALIAS C, PASISKEVICIUS V, et al. Narrowband and tunable ring optical parametric oscillator with a volume Bragg grating[J]. Optics Letters, 2007, 32(22): 3278-3280. DOI: 10.1364/OL.32.003278
[46] PENG Y, WEI X, NIE Z, et al. High-power, narrow-bandwidth mid-infrared PPMgLN optical parametric oscillator with a volume Bragg grating[J]. Optics Express, 2015, 23(24): 30827-30832. DOI: 10.1364/OE.23.030827
[47] HE G, GUO J, JIAO Z, et al. High-efficiency near-degenerate PPMgLN optical parametric oscillator with a volume Bragg grating[J]. Optics Letters, 2012, 37(8): 1364-1366. DOI: 10.1364/OL.37.001364
[48] LI K, YANG S, WANG X, et al. Frequency chirped intensity modulated mid-infrared light source based on optical parametric oscillation[J]. IEEE Photonics Journal, 2020, 12(1): 1500409.
[49] ERUSHIN E, NYUSHKOV B, IVANENKO A, et al. Tunable injection-seeded fan-out-PPLN optical parametric oscillator for high-sensitivity gas detection[J]. Laser Physics Letters, 2021, 18(11): 116201-116207. DOI: 10.1088/1612-202X/ac2585
[50] 卞进田, 叶庆, 孙晓泉. ZnGeP2 OPO产生4.3μm波段窄线宽激光实验研究[J]. 国防科技大学学报, 2018, 40(4): 9-14. https://www.cnki.com.cn/Article/CJFDTOTAL-GFKJ201804002.htm BIAN J T, YE Q, SUN X Q. ZnGeP2 optical parametric oscillator 4.3 μm laser with narrow line-width[J]. Journal of National University of Defense Technology, 2018, 40(4): 9-14(in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-GFKJ201804002.htm
[51] BIAN J T, KONG H, YE Q, et al. Narrow-linewidth BaGa4Se7 optical parametric oscillator[J]. Chinese Optics Letters, 2022, 20(4): 041901.
计量
- 文章访问数: 9
- HTML全文浏览量: 0
- PDF下载量: 16