打拿极光电倍增管逐级增益自动测试系统

刘欣妍, 钱芸生, 魏静雯

刘欣妍, 钱芸生, 魏静雯. 打拿极光电倍增管逐级增益自动测试系统[J]. 应用光学, 2022, 43(6): 1117-1123. DOI: 10.5768/JAO202243.0604010
引用本文: 刘欣妍, 钱芸生, 魏静雯. 打拿极光电倍增管逐级增益自动测试系统[J]. 应用光学, 2022, 43(6): 1117-1123. DOI: 10.5768/JAO202243.0604010
shu
LIU Xinyan, QIAN Yunsheng, WEI Jingwen. Automatic test system for progressive gain of dynode photomultiplier tube[J]. Journal of Applied Optics, 2022, 43(6): 1117-1123. DOI: 10.5768/JAO202243.0604010
Citation: LIU Xinyan, QIAN Yunsheng, WEI Jingwen. Automatic test system for progressive gain of dynode photomultiplier tube[J]. Journal of Applied Optics, 2022, 43(6): 1117-1123. DOI: 10.5768/JAO202243.0604010
shu

打拿极光电倍增管逐级增益自动测试系统

  • 南京理工大学 电子工程与光电技术学院，江苏 南京 210094

基金项目: 国家重大科学仪器设备研制重点项目（2016YFF0100400）
详细信息
    作者简介:

    刘欣妍（1998—），女，硕士研究生，主要从事微光成像器件测试技术研究。E-mail：liuxinyann@yeah.net

    通讯作者:

    钱芸生（1968—），男，博士，教授，主要从事微光与红外成像器件与系统及相关测试技术、原位测试余分析、光电测试、图像处理和仿真研究。E-mail：yshqian@mail.njust.edu.cn

  • 中图分类号: TN206

Automatic test system for progressive gain of dynode photomultiplier tube

  • School of Electronic and Optical Engineering, Nanjing University of Science and Technology, Nanjing 210094, China

HTML全文

    摘要
  • 摘要:

    逐级增益是评价打拿极光电倍增管（photomultiplier tube, PMT）性能的重要参数之一，目前国内仅有一台逐级增益测试系统，依靠人工操作，测试效率不高。为提高逐级增益测试效率，设计了一套打拿极PMT逐级增益自动测试系统。该系统使用由可控光阑、电动挡板等组成的自动光源实现光源输出强度连续可调，并通过控制高压模块、分压器模块和继电器模块实现逐级电压的通断。利用精度为0.01 nA的电流计模块完成信号的采集、处理和传输，实现对逐级增益的自动测量。实验证明，该系统可有效测量打拿极PMT的逐级增益特性，测试重复性在2%以内，满足测试需求。

    Abstract:

    Progressive gain is one of the important parameters for evaluating the performance of dynode photomultiplier tube (PMT), but there is only one related test system at home, which relies on manual operation and the test efficiency is not high. In order to improve the efficiency of testing progressive gain, an automatic test system for the progressive gain of dynode PMT was designed. The continuous adjustment of the light source output intensity was achieved by an automatic light system composed of controllable disphragm, electric baffle and etc, and the on-off of the progressive voltage was controlled by the high-voltage module, the voltage divider module and the relay module. Finally, based on the galvanometer module with an accuracy of 0.01 nA, the acquisition, processing and transmission of the signal was completed, so as to realize the automatic measurement of the progressive gain. The experimental results show that the proposed system can effectively measure the progressive gain characteristics of the dynode PMT, and the test repeatability is within 2%, which meets the test requirements.

    • HTML全文

  • 引言

    利用相机对位置已知的控制点进行拍摄，通过拍摄得到的二维图像解算相机在世界坐标系下的位置和姿态，称为位姿估计问题，是机器人导航、计算机视觉等领域[1-3]的核心问题之一。位姿估计根据计算方法分为2种方式[4]：第一种一般通过解算一组与旋转参数和平移参数相关的多项式，进而获得位姿估计结果[5]，这类方法对于噪声敏感[3]，估计性能有限；另一种通常称为迭代类算法，这种算法通过建立代价函数[6]，将位姿估计问题转换为非线性最小二乘优化问题，然后利用Gauss-Newton[7]、Levenberg-Marquardt等方法进行位姿参数的估计，这类算法能够综合利用多点的冗余信息[8]，增加了估计的鲁棒性，迭代算法中还有一种称为正交迭代(orthogonal iteration, OI)算法，该算法将物空间共线性误差的最小化视为求解绝对定向问题，然后利用迭代方式进行参数估计[9]，这种方法计算量较小，速度较快[10]，并保证了全局收敛性，但是性能略逊于传统迭代算法。

    传统迭代算法是从信号的测量误差出发[11]，利用误差最小化的原则建立代价函数，从而得到估计值。这一类算法是一种最小二乘估计，试图使采样得到的信号和无噪声情况下的数据之差的平方达到最小，这一类方法对采样数据不做任何的统计假设，性能取决于噪声的特性，并且往往不是最佳估计，而由于没有对信号做任何统计假设，估计的统计性能是无法评价的[12]。由于相机拍摄的不确定性，图片采样所得的像素信号应视为随机信号，本文从随机信号的统计特性出发，根据采样信号的概率密度推导了最大似然位姿估计的一般形式。从理论上证明了利用单幅图像，在各向同性高斯噪声情况下传统迭代算法等效于最大似然估计，因此在理论上传统迭代算法只能应用于各向同性高斯噪声，对于复杂噪声，将会产生模型误差，导致算法失效。本文的最大似然算法将是推导复杂噪声下的位姿估计的基础；另外，根据信号的Fisher信息阵推导了位姿估计的克拉美-罗界(Cramér-Rao bound, CRB)，作为任何无偏估计的方差下界，克拉美-罗界可以作为位姿估计方法有效性的评价标准，并且根据所得结果可知，通过适当增加采样个数可以有效提高估计的性能。这是由于通过采样信号个数的增加，可以更好地估计随机噪声的统计特性，从而减小随机噪声对于估计的影响。

    1   信号模型

    假设有n个控制点，其在世界坐标系中的坐标为piw，i=1, 2, …, n，它们在相机坐标系下的坐标为pic，i=1, 2, …, n，相机对控制点进行拍摄，可以得到像平面上的n个像点，像点在图像坐标系下的物理坐标为pip，i=1, 2, …, n，在图像坐标系下的像素坐标为(ui, vi)T，如图 1所示。在k时刻对控制点进行拍摄，考虑采样过程的随机噪声影响，则观测量为如下随机向量：

    图  1  位姿估计问题的几何示意图
    Figure  1.  Geometric sketch of pose estimation
    下载: 全尺寸图片 幻灯片
    $$ \begin{array}{l} s\left( k \right) = \left[ {{{\hat u}_1}\left( k \right),{{\hat v}_1}\left( k \right),{{\hat u}_2}\left( k \right),{{\hat v}_2}\left( k \right), \cdots ,} \right.\\ \;\;\;\;\;\;\;\;\;\;\;\;{\left. {{{\hat u}_n}\left( k \right),{{\hat v}_n}\left( k \right)} \right]^{\rm{T}}} = \mathit{\boldsymbol{m}}\left( k \right) + \mathit{\boldsymbol{n}}\left( k \right) \end{array} $$ (1)

    式中：k=1, 2, ..., K, 其中:

    $$ \begin{array}{l} \mathit{\boldsymbol{m}}\left( k \right) = \left[ {{u_1}\left( k \right),{v_1}\left( k \right),{u_2}\left( k \right),{v_2}\left( k \right), \cdots } \right.\\ \;\;\;\;\;\;\;\;\;\;\;\;\;{\left. {{u_n}\left( k \right),{v_n}\left( k \right)} \right]^{\rm{T}}} \end{array} $$ (2)

    世界坐标与相机坐标存在以下关系：

    $$ \mathit{\boldsymbol{p}}_i^c\left( k \right) = \mathit{\boldsymbol{R}}\left( k \right)\mathit{\boldsymbol{p}}_i^w\left( k \right) + \mathit{\boldsymbol{t}}\left( k \right) $$ (3)

    由共线方程可以得到：

    $$ \left\{ \begin{array}{l} \frac{{x_i^p\left( k \right) - x_o^p}}{f} = \frac{{x_i^c\left( k \right)}}{{z_i^c\left( k \right)}}\\ \frac{{y_i^p\left( k \right) - y_o^p}}{f} = \frac{{y_i^c\left( k \right)}}{{z_i^c\left( k \right)}} \end{array} \right. $$ (4)

    式中(xop, yop)T为图像主点I在图像坐标系下的物理坐标，则像素坐标系与世界坐标系的转换关系可以表示为

    $$ \left\{ \begin{array}{l} {u_i}\left( k \right) = \frac{f}{{{d_x}}}\;\frac{{{\mathit{\boldsymbol{r}}_1}\left( k \right)\mathit{\boldsymbol{p}}_i^w\left( k \right) + {t_x}\left( k \right)}}{{{\mathit{\boldsymbol{r}}_3}\left( k \right)\mathit{\boldsymbol{p}}_i^w\left( k \right) + {t_z}\left( k \right)}} + \frac{{x_o^p}}{{{d_x}}}\\ {v_i}\left( k \right) = \frac{f}{{{d_y}}}\;\frac{{{\mathit{\boldsymbol{r}}_2}\left( k \right)\mathit{\boldsymbol{p}}_i^w\left( k \right) + {t_y}\left( k \right)}}{{{\mathit{\boldsymbol{r}}_3}\left( k \right)\mathit{\boldsymbol{p}}_i^w\left( k \right) + {t_z}\left( k \right)}} + \frac{{y_o^p}}{{{d_y}}} \end{array} \right. $$ (5)

    式中：tx(k)、ty(k)、tz(k)分别为t (k)从上到下的3个元素；r1(k)、r2(k)、r3(k)分别为R (k)的从上到下的3个横向量；R可以用3个欧拉角表示。假设在K次拍摄过程中，相机的姿态、控制点的世界坐标保持不变，即m不随时间变化，则采样信号可以表示为

    $$ \mathit{\boldsymbol{s}}\left( k \right) = \mathit{\boldsymbol{m}} + \mathit{\boldsymbol{n}}\left( k \right) $$ (6)

    2   最大似然估计

    最大似然估计是一种使用了信号概率模型的方法，其基本思想是参数的估计值应是使观测信号概率最大的值。假设信号为高斯随机变量，噪声的均值为零，则m即为观测量s(k)的均值，根据概率论可知，k时刻观测量的概率密度为[13]

    $$ \begin{array}{l} {p_{s\left( k \right)}} = \frac{1}{{\det \left( {{\rm{ \mathit{ π} }}{\mathit{\boldsymbol{Q}}_s}} \right)}}\exp \left\{ { - \left[ {{\mathit{\boldsymbol{s}}^H}\left( k \right) - } \right.} \right.\\ \;\;\;\;\;\;\;\;\left. {\left. {{\mathit{\boldsymbol{m}}^H}} \right]\mathit{\boldsymbol{Q}}_s^{ - 1}\left[ {\mathit{\boldsymbol{s}}\left( k \right) - \mathit{\boldsymbol{m}}} \right]} \right\} \end{array} $$ (7)

    式中：det(·)表示矩阵的行列式；Qs为随机信号s的方差。与待估计参数无关，待估计参数表示为

    $$ \mathit{\boldsymbol{\theta }} = {\left[ {\alpha ,\beta ,\gamma ,{t_x},{t_y},{t_z}} \right]^{\rm{T}}} $$ (8)

    则m(θ)表示为θ的函数，假设K次拍摄是统计独立的，则K次拍摄观测值的联合密度函数为

    $$ \begin{array}{l} {p_{s\left( 1 \right),s\left( 2 \right), \cdots ,s\left( k \right)}} = \mathop \Pi \limits_{k = 1}^k \frac{1}{{\det \left( {{\rm{ \mathit{ π} }}{\mathit{\boldsymbol{Q}}_s}} \right)}}\exp \left\{ { - \left[ {{s^H}\left( k \right) - } \right.} \right.\\ \;\;\;\;\;\;\;\;\left. {\left. {{\mathit{\boldsymbol{m}}^{\rm{H}}}\left( \mathit{\boldsymbol{\theta }} \right)} \right]\mathit{\boldsymbol{Q}}_s^{ - 1}\left[ {\mathit{\boldsymbol{s}}\left( k \right) - \mathit{\boldsymbol{m}}\left( \mathit{\boldsymbol{\theta }} \right)} \right]} \right\} \end{array} $$ (9)

    对等号两边求自然对数，去掉常数项并除以K，得到似然函数为

    $$ \begin{array}{l} L\left( \mathit{\boldsymbol{\theta }} \right) = - \ln \left[ {\det \left( {{\mathit{\boldsymbol{Q}}_s}} \right)} \right] - \frac{1}{K}\sum\limits_{k = 1}^K {\left[ {{\mathit{\boldsymbol{s}}^{\rm{H}}}\left( k \right) - } \right.} \\ \;\;\;\;\;\;\;\;\;\;\;\left. {{\mathit{\boldsymbol{m}}^{\rm{H}}}\left( \mathit{\boldsymbol{\theta }} \right)} \right]\left. {\mathit{\boldsymbol{Q}}_s^{ - 1}\left[ {\mathit{\boldsymbol{s}}\left( k \right) - \mathit{\boldsymbol{m}}\left( \mathit{\boldsymbol{\theta }} \right)} \right]} \right\} \end{array} $$ (10)

    (10) 式的第二项可以表示为

    $$ \begin{array}{l} \frac{1}{K}\sum\limits_{k = 1}^K {\left[ {{\mathit{\boldsymbol{s}}^{\rm{H}}}\left( k \right) - {\mathit{\boldsymbol{m}}^{\rm{H}}}\left( \theta \right)} \right]} \mathit{\boldsymbol{Q}}_s^{ - 1}\left[ {\mathit{\boldsymbol{s}}\left( k \right) - \mathit{\boldsymbol{m}}\left( \mathit{\boldsymbol{\theta }} \right)} \right] = \\ \rm{tr}\left\{ {\frac{1}{K}\sum\limits_{k = 1}^K {\left[ {{\mathit{\boldsymbol{s}}^{\rm{H}}}\left( k \right) - {\mathit{\boldsymbol{m}}^{\rm{H}}}\left( \theta \right)} \right]} \mathit{\boldsymbol{Q}}_s^{ - 1}\left[ {s\left( k \right) - \mathit{\boldsymbol{m}}\left( \mathit{\boldsymbol{\theta }} \right)} \right]} \right\} = \\ \rm{tr}\left\{ {\mathit{\boldsymbol{Q}}_s^{ - 1}\frac{1}{K}\sum\limits_{k = 1}^K {\left[ {\mathit{\boldsymbol{s}}\left( k \right) - \mathit{\boldsymbol{m}}\left( \mathit{\boldsymbol{\theta }} \right)} \right]\left[ {{\mathit{\boldsymbol{s}}^{\rm{H}}}\left( k \right) - {\mathit{\boldsymbol{m}}^{\rm{H}}}\left( \mathit{\boldsymbol{\theta }} \right)} \right]} } \right\} = \\ \rm{tr}\left\{ {\mathit{\boldsymbol{Q}}_s^{ - 1}{\mathit{\boldsymbol{C}}_s}\left( \mathit{\boldsymbol{\theta }} \right)} \right\} \end{array} $$ (11)

    式中：tr{·}表示取矩阵的迹；Cs(θ)定义为信号相关矩阵：

    $$ {\mathit{\boldsymbol{C}}_s}\left( \theta \right) \buildrel \Delta \over = \frac{1}{K}\sum\limits_{k = 1}^K {\left[ {\mathit{\boldsymbol{s}}\left( k \right) - \mathit{\boldsymbol{m}}\left( \mathit{\boldsymbol{\theta }} \right)} \right]\left[ {{\mathit{\boldsymbol{s}}^{\rm{H}}}\left( k \right) - {\mathit{\boldsymbol{m}}^{\rm{H}}}\left( \mathit{\boldsymbol{\theta }} \right)} \right]} $$ (12)

    于是得到似然函数：

    $$ L\left( \theta \right) = - \ln \left[ {\det \left( {{\mathit{\boldsymbol{Q}}_s}} \right)} \right] - \rm{tr}\left\{ {\mathit{\boldsymbol{Q}}_s^{ - 1}{\mathit{\boldsymbol{C}}_s}\left( \mathit{\boldsymbol{\theta }} \right)} \right\} $$ (13)

    (13) 式为最大似然估计的一般形式，使似然函数最大的参数θ即为位姿估计的最大似然估计，可以通过搜索或者更高效的迭代算法计算。

    当噪声为各向同性的高斯噪声时，Qs=δw2I，δw2为噪声功率，似然函数可以表示为

    $$ L\left( \theta \right) = - 2n\ln \delta _w^2 - \frac{1}{{\delta _w^2}}{\rm{tr}}\left\{ {{C_s}\left( \mathit{\boldsymbol{\theta }} \right)} \right\} $$ (14)

    由于似然函数取最大值的必要非充分条件为

    $$ \frac{\partial }{{\partial \delta _w^2}}L\left( \mathit{\boldsymbol{\theta }} \right) = 0 $$ (15)
    $$ \delta _w^2 - \frac{{\rm{tr}\left\{ {{\mathit{\boldsymbol{C}}_s}\left( \mathit{\boldsymbol{\theta }} \right)} \right\}}}{{2n}} $$ (16)

    将上式代入似然函数，去除常数项并除以2n，似然函数化简为

    $$ L\left( \mathit{\boldsymbol{\theta }} \right) = - {\rm{lntr}}\left\{ {{\mathit{\boldsymbol{C}}_s}\left( \mathit{\boldsymbol{\theta }} \right)} \right\} $$ (17)

    根据Cs(θ)的定义，求似然函数的最大值等效于求下式的最小值：

    $$ \begin{array}{l} \rm{tr}\left\{ {{\mathit{\boldsymbol{C}}_s}\left( \mathit{\boldsymbol{\theta }} \right)} \right\} = \frac{1}{K}\sum\limits_{k = 1}^K {\left[ {{s^{\rm{H}}}\left( k \right) - {\mathit{\boldsymbol{m}}^{\rm{H}}}\left( \mathit{\boldsymbol{\theta }} \right)} \right]\left[ {s\left( k \right) - \mathit{\boldsymbol{m}}\left( \mathit{\boldsymbol{\theta }} \right)} \right]} = \\ \frac{1}{K}{\sum\limits_{k = 1}^K {\sum\limits_{i = 1}^n {\left[ {\left( {{{\hat u}_i}\left( k \right) - \frac{f}{{{d_x}}}\;\frac{{{\mathit{\boldsymbol{r}}_1}\left( k \right)\mathit{\boldsymbol{p}}_i^w\left( k \right) + {t_x}\left( k \right)}}{{{\mathit{\boldsymbol{r}}_3}\left( k \right)\mathit{\boldsymbol{p}}_i^w\left( k \right) + {t_z}\left( k \right)}} + \frac{{x_o^p}}{{{d_x}}}} \right)} \right.} } ^2} + \\ \left. {{{\left( {{{\hat v}_i}\left( k \right) - \frac{f}{{{d_y}}}\;\frac{{{\mathit{\boldsymbol{r}}_2}\left( k \right)\mathit{\boldsymbol{p}}_i^w\left( k \right) + {t_y}\left( k \right)}}{{{\mathit{\boldsymbol{r}}_3}\left( k \right)\mathit{\boldsymbol{p}}_i^w\left( k \right) + {t_z}\left( k \right)}} + \frac{{y_o^p}}{{{d_y}}}} \right)}^2}} \right] \end{array} $$ (18)

    利用单幅图像时，上式与传统的迭代算法等效，因此，利用单幅图像时，在各向同性高斯噪声情况下传统迭代算法与最大似然算法等效。值得注意的是，上述结果是在(16)式的情况下得到的，即位姿估计是在噪声功率估计的基础上实现的，所以增加信号的采样，提高对噪声统计特性的估计，有利于提高参数估计的性能。

    3   克拉美-罗界分析

    为了了解参数估计的潜在性能，可以利用克拉美-罗界作为评价标准。克拉美-罗界给出了参数所有无偏估计的方差极限[14]，将估计误差的方差表示为C(θ)，对于所有无偏估计，有关系式：

    $$ \mathit{\boldsymbol{C}}\left( \mathit{\boldsymbol{\theta }} \right) \ge {\mathit{\boldsymbol{C}}_{{\rm{CR}}}}\left( \mathit{\boldsymbol{\theta }} \right) $$ (19)

    式中CCR(θ)表示克拉美-罗界，克拉美-罗界可以通过求解Fisher信息阵J的逆矩阵求得：

    $$ {\mathit{\boldsymbol{C}}_{CR}}\left( \mathit{\boldsymbol{\theta }} \right) = {\mathit{\boldsymbol{J}}^{ - 1}} $$ (20)

    Fisher信息阵J中的元素为

    $$ {\left[ \mathit{\boldsymbol{J}} \right]_{i,j}} = E\left[ {\frac{{\partial L\left( \mathit{\boldsymbol{\theta }} \right)}}{{\partial {\theta _i}}}\;\frac{{\partial L\left( \mathit{\boldsymbol{\theta }} \right)}}{{\partial {\theta _j}}}} \right] $$ (21)

    通过对(13)式求偏导并取期望可以得到[J]i, j通用的形式(推导可参见文献[15])：

    $$ \begin{array}{l} {\left[ \mathit{\boldsymbol{J}} \right]_{i,j}} = K\rm{tr}\left[ {\mathit{\boldsymbol{Q}}_s^{ - 1}\frac{{\partial {\mathit{\boldsymbol{Q}}_s}}}{{\partial {\theta _i}}}\mathit{\boldsymbol{Q}}_s^{ - 1}\;\frac{{\partial {\mathit{\boldsymbol{Q}}_s}}}{{\partial {\theta _j}}}} \right] + \\ \;\;\;\;\;\;\;\;\;\;\;\;2K{\mathop{\rm Re}\nolimits} \left[ {\frac{{\partial {\mathit{\boldsymbol{m}}^H}\left( \theta \right)}}{{\partial {\theta _i}}}\mathit{\boldsymbol{Q}}_s^{ - 1}\;\frac{{\partial \mathit{\boldsymbol{m}}\left( \theta \right)}}{{\partial {\theta _j}}}} \right] \end{array} $$ (22)

    (22) 式为Fisher信息阵的一般化形式，是许多克拉美-罗界推导的起点，需要注意的是，通常方差Qs是一个未知矩阵，因此上式应该包含方差的未知参数。

    在位姿估计问题中，由于(22)式的第1项只与方差所含未知参数有关，第2项只与待估计参数有关，可以将Fisher信息阵分块：

    $$ \mathit{\boldsymbol{J = }}\left[ \begin{array}{l} {\mathit{\boldsymbol{J}}_{{Q_s}{Q_s}}}\;\;\;\;0\\ 0\;\;\;\;\;\;\;\;{\mathit{\boldsymbol{J}}_{00}} \end{array} \right] $$ (23)
    $$ {\left[ {{\mathit{\boldsymbol{J}}_{{Q_s}{Q_s}}}} \right]_{i,j = }}K{\rm{tr}}\left[ {\mathit{\boldsymbol{Q}}_s^{ - 1}\frac{{\partial {\mathit{\boldsymbol{Q}}_s}}}{{\partial {q_i}}}\mathit{\boldsymbol{Q}}_s^{ - 1}\;\frac{{\partial {Q_s}}}{{\partial {q_j}}}} \right] $$ (24)
    $$ {\left[ {{\mathit{\boldsymbol{J}}_{00}}} \right]_{i,j}} = 2K{\mathop{\rm Re}\nolimits} \left[ {\frac{{\partial {\mathit{\boldsymbol{m}}^{\rm{H}}}\left( \mathit{\boldsymbol{\theta }} \right)}}{{\partial {\theta _i}}}Q_s^{ - 1}\;\frac{{\partial \mathit{\boldsymbol{m}}\left( \mathit{\boldsymbol{\theta }} \right)}}{{\partial {\theta _j}}}} \right] $$ (25)

    式中qi为方差所含的未知参数，当噪声为各点无关的各向同性高斯噪声时，Qs=δw2I，方差所含未知参数为δw2，则有：

    $$ {\mathit{\boldsymbol{J}}_{{Q_s}{Q_s}}} = {\mathit{\boldsymbol{J}}_{\delta _w^2\delta _w^2}} = \frac{{2nK}}{{{{\left( {\delta _w^2} \right)}^2}}} $$ (26)

    根据矩阵分块求逆原理，噪声功率估计的克拉美-罗界为

    $$ {\mathit{\boldsymbol{C}}_{CR}}\left( {\delta _w^2} \right) = \frac{{{{\left( {\delta _w^2} \right)}^2}}}{{2nK}} $$ (27)

    可以看出噪声功率估计的克拉美-罗界随着采样数的增加而减小，说明增加信号的采样，可以提高对噪声统计特性的估计。将Qs=δw2I代入(25)式，得到：

    $$ {\left[ {{\mathit{\boldsymbol{J}}_{\theta \theta }}} \right]_{i,j}} = 2K{\mathop{\rm Re}\nolimits} \left[ {\frac{{\partial {\mathit{\boldsymbol{m}}^{\rm{H}}}\left( \mathit{\boldsymbol{\theta }} \right)}}{{\partial {\theta _i}}}Q_s^{ - 1}\;\frac{{\partial \mathit{\boldsymbol{m}}\left( \mathit{\boldsymbol{\theta }} \right)}}{{\partial {\theta _j}}}} \right] $$ (28)

    定义

    $$ D \buildrel \Delta \over = \frac{{\partial {\mathit{\boldsymbol{m}}^{\rm{H}}}\left( \mathit{\boldsymbol{\theta }} \right)}}{{\partial \theta }} $$ (29)

    由于θ=[α, β, γ, tx, ty, tz]T，待估计参数的克拉美-罗界为

    $$ {\mathit{\boldsymbol{C}}_{CR}}\left( \theta \right) = \frac{{\delta _w^2}}{{2K}}{\left\{ {{\mathop{\rm Re}\nolimits} \left[ {\mathit{\boldsymbol{D}}{\mathit{\boldsymbol{D}}^{\rm{H}}}} \right]} \right\}^{ - 1}} $$ (30)

    从(30)式中可以看出，参数估计的克拉美-罗界随拍摄数量的增加线性下降，随噪声功率的增加线性增加。

    4   仿真

    利用计算机仿真对最大似然估计和克拉美-罗界进行分析，并与传统迭代算法以及OI算法进行比较。

    仿真参数：相机焦距f=35mm，单个像素尺寸dx=dy=62.5 μm，主点在图像中心，控制点从相机坐标系的[-2 2]×[-2 2]×[4 9] m3的长方形区域内随机选取，相机坐标系相对于世界坐标系的欧拉角为[20 10 30]°，平移矩阵为[2 3 10]m。图 2(a)是随机选取的控制点在相机坐标系下的坐标，图 2(b)中圆圈表示相机拍摄控制点时实际成像位置(像素坐标系)，圆圈旁边的10个点为噪声功率为1dB时，由拍摄的10张图片提取的坐标点。

    图  2  控制点在相机坐标系和成像平面上的位置
    Figure  2.  Positions of control points in camera coordinate system and image plane
    下载: 全尺寸图片 幻灯片

    图 3是传统迭代算法、OI算法、最大似然方法在不同的噪声功率下的参数估计均方根误差(root-mean-square error, RMSE)，在各噪声级下完成500次独立实验，其中最大似然方法利用了10张图片。从仿真结果中可以看出，3种方法在噪声功率较小的情况下都趋向于克拉美-罗界，传统迭代算法性能略优于OI算法，而最大似然方法的性能明显优于其他两种方法。

    图  3  不同噪声功率下各算法的参数估计均方根误差
    Figure  3.  RMSEs of different methods with different noise powers
    下载: 全尺寸图片 幻灯片

    图 4是最大似然估计在不同的图片数下的参数估计均方根误差，噪声功率为2 dB，每个图片数完成500次独立实验，可以看出参数估计的RMSE随着图片数的增加而减小，说明适当地增加拍摄数量可有效提高估计的性能。

    图  4  不同图片数下参数估计的均方根误差
    Figure  4.  RMSEs with different snapshot numbers
    下载: 全尺寸图片 幻灯片

    5   结论

    给出机器视觉位姿估计的信号模型，并从随机信号的角度出发，推导了基于机器视觉的最大似然位姿估计以及相应的克拉美-罗界，从理论上证明了利用单幅图像时，在各向同性高斯噪声情况下传统迭代算法与最大似然算法等效，通过仿真分析，可知在图像数量适当增加的情况下，估计性能得到了明显改进，适用于更高精度的位姿估计。并且推导的位姿估计克拉美-罗界作为任何无偏估计的方差下界，是位姿估计的性能极限，可以作为位姿估计方法有效性的评价标准；下一步，针对复杂的噪声情况，应在此基础上讨论各向异性噪声以及相关噪声情况下的位姿估计方法及其克拉美-罗界。

  • 图  1   盒栅状打拿极PMT的基本结构示意图

    Figure  1.   Schematic diagram of basic structure for box-grid type of dynode PMT

    下载: 全尺寸图片 幻灯片

    图  2   测试系统原理框图

    Figure  2.   Block diagram of test system principle

    下载: 全尺寸图片 幻灯片

    图  3   测试系统示意图

    Figure  3.   Schematic diagram of test system

    下载: 全尺寸图片 幻灯片

    图  4   继电器阵列原理图

    Figure  4.   Schematic diagram of relay array

    下载: 全尺寸图片 幻灯片

    图  5   测试系统软件界面

    Figure  5.   Interface of test system software

    下载: 全尺寸图片 幻灯片

    图  6   增益特性曲线

    Figure  6.   Fitting curves of gain characteristics

    下载: 全尺寸图片 幻灯片

    表  1   不同工作电压下3个PMT的逐级增益

    Table  1   Progressive gain of three PMT at different operating voltages

    工作电压/V不同管号PMT的增益
    91669171R09289
    3003.67×1024.84×1021.48×102
    4005.73×1038.01×1031.89×103
    5003.60×1046.87×1041.17×104
    6001.18×1064.07E×1056.43×104
    7001.21×1061.90×1062.27×105
    8005.15×1068.00×1069.54×105
    9001.44×1072.70×1072.69×106
    10003.50×1077.92×1077.20×106
    11001.10×1082.22×1081.58×107
    12002.82×1086.76×1084.77×107
    13008.82×1081.30×1098.76×107
    下载: 导出CSV

    表  2   PPH号PMT逐级增益测量结果

    Table  2   Progressive gain measurement results for PMT of PPH

    级数测试次数重复性/%
    123456
    16.2726.1866.3226.1446.0676.2671.53
    26.4836.5276.4846.4826.5416.5350.44
    35.2535.2895.2615.2675.2795.2690.24
    45.7355.6795.6245.5955.6495.6290.87
    52.9893.0123.0043.0063.0113.0050.28
    62.7192.7392.7412.7482.7442.7450.38
    74.9844.9764.9584.954.9814.9680.27
    85.3345.3635.3475.3655.3795.3670.30
    93.8793.9043.9093.8823.9083.9010.34
    103.1383.1383.1253.1503.1303.1400.28
    整管增益3.23×1063.27×1063.23×1063.15×1063.20×1063.7×1061.45
    下载: 导出CSV

    表  3   8160号PMT逐级增益测量结果

    Table  3   Progressive gain measurement results for PMT of 8160

    级数测试次数重复性/%
    123456
    111.08011.34211.30211.31611.30011.1490.95
    22.6262.6312.6302.6282.6352.6330.12
    32.6682.6372.6472.6422.6532.6500.40
    42.1402.0692.0882.0732.0842.0871.23
    51.9471.9281.9391.9351.9401.9430.34
    62.0822.0882.0852.0902.0832.0900.17
    72.4252.4102.4142.4172.4152.4180.21
    81.9531.9501.9481.9511.9451.9530.16
    92.0051.9972.0001.9961.9941.9990.19
    101.9051.9051.8991.9041.9001.9000.15
    111.8691.8561.8621.8651.8591.8670.27
    整管增益2.27×1042.18×1042.21×1042.20×1042.20×1042.21×1041.38%
    下载: 导出CSV

    表  4   PJ号PMT逐级增益测量结果

    Table  4   Progressive gain measurement results for PMT of PJ

    级数测试次数重复性/%
    123456
    13.8213.8433.8173.9983.8443.8331.78
    24.8844.8565.0804.8895.0414.8821.94
    37.0527.1407.0907.0997.0867.0480.48
    47.1467.1327.0257.0797.1827.1280.78
    58.5118.5828.4998.4688.4778.5020.47
    66.5776.5816.5416.6276.5446.5710.47
    78.858.8458.8508.8688.7558.8420.46
    86.7776.7356.7446.7576.6866.7380.45
    99.4009.4509.4369.3439.3059.3980.59
    整管增益2.97×1073.02×1073.0×1073.09×1072.98×1072.94×1071.71
    下载: 导出CSV

    表  5   高压电源输出和电流计输出6次测量结果

    Table  5   Six measurement results output by high-voltage power supply and galvanometer

    序号电压/V电流/μA序号电压/V电流/μA
    1996.01.58154999.81.5807
    21 001.01.579651 000.01.5823
    3998.51.580061 000.41.5803
    下载: 导出CSV
    • 参考文献(17)
  • [1] 张珮. 基于光电倍增管的数据采集系统设计[D]. 长春: 长春理工大学, 2017.

    ZHANG Pei. Design of data acquisition system based on photomultiplier tube[D]. Changchun: Changchun University of Science and Technology, 2017.
    [2] 赵艳, 潘超, 赵一鸣, 等. 光电倍增管寿命试验研究与系统设计[J]. 遥测遥控,2021,42(6):107-112.

    ZHAO Yan, PAN Chao, ZHAO Yiming, et al. Research on life test and system design for PMT photomultiplier tube[J]. Journal of Telemetry, Tracking and Command,2021,42(6):107-112.
    [3] 杨云开. 用于生物毒性的光电倍增管检测技术[J]. 数码世界,2019(7):47.

    YANG Yunkai. Photomultiplier tube detection technology for biotoxicity[J]. Digital World,2019(7):47.
    [4] T. Sekine, G. Delso, K. G. Zeimpekis, 等. 应用硅光电倍增检测器可以降低临床PET/MRI上使用18F-FDG的剂量[J]. 国际医学放射学杂志,2018,41(2):220.

    T. SEKINE, et al. Reduction of 18F-FDG dose in clinical PET/MR imaging by using silicon photomultiplier detectors[J]. International Journal of Medical Radiology,2018,41(2):220.
    [5] 赵富宽, 张波. PET/MRI研制及应用进展[J]. 中国医疗设备,2014,29(8):66-69. doi: 10.3969/j.issn.1674-1633.2014.08.020

    ZHAO Fukuan, ZHANG Bo. Development and application progress of PET/MRI[J]. China Medical Devices,2014,29(8):66-69. doi: 10.3969/j.issn.1674-1633.2014.08.020
    [6] 黄臻成, 沈韩, 唐健. 中微子探测中的光学应用[J]. 物理与工程,2022,32(1):33-36.

    HUANG Zhencheng, SHEN Han, TANG Jian. Application of optics in neutrino detections[J]. Physics and Engineering,2022,32(1):33-36.
    [7] 刘茂元, 陈鑫, 念聪. 光电倍增管在地基粒子天体物理实验中的应用[J]. 知识文库,2019(7):37.

    LIU Maoyuan, CHEN Xin, NIAN Cong. Application of photomultiplier tubes in ground-based particle astrophysics experiments[J]. Knowledge Base,2019(7):37.
    [8] 王浩东. 真空紫外光电倍增管现状及发展趋势分析[J]. 真空电子技术,2022(2):23-28.

    WANG Haodong. Status and development trend of UV-photomultiplier tubes[J]. Vacuum Electronics,2022(2):23-28.
    [9] 孙建宁, 任玲, 丛晓庆, 等. 一种大尺寸微通道板型光电倍增管[J]. 红外与激光工程,2017,46(4):18-22.

    SUN Jianning, REN Ling, CONG Xiaoqing, et al. Large-area micro-channel plate photomultiplier tube[J]. Infrared and Laser Engineering,2017,46(4):18-22.
    [10] 郭乐慧, 陈萍, 李立立, 等. 光电倍增管关键技术研究进展[J]. 真空电子技术,2020(4):1-13. doi: 10.16540/j.cnki.cn11-2485/tn.2020.04.01

    GUO Lehui, CHEN Ping, LI Lili, et al. Research progress on key technologies of photomultiplier tubes[J]. Vacuum Electronics,2020(4):1-13. doi: 10.16540/j.cnki.cn11-2485/tn.2020.04.01
    [11] 汪贵华. 光电子器件[M]. 2版. 北京: 国防工业出版社, 2014.

    WANG Guihua. Optoelectronic devices[M]. 2nd. Beijing: National Defense Industry Press, 2014.
    [12] 赵文锦. 光电倍增管的技术发展状态[J]. 光电子技术,2011,31(3):145-148. doi: 10.3969/j.issn.1005-488X.2011.03.001

    ZHAO Wenjin. Developments in technology of photomultipliers[J]. Optoelectronic Technology,2011,31(3):145-148. doi: 10.3969/j.issn.1005-488X.2011.03.001
    [13] 宋登元, 孙同文. 光电倍增管的特性、结构及选取考虑[J]. 中国仪器仪表,1999(5):33-35.

    SONG Dengyuan, SUN Tongwen. Properties, structures and choosing consideration for photomultiplier tubes[J]. China Instrmentation,1999(5):33-35.
    [14] 中国国家标准化管理委员会. 光电倍增管总规范: GB/T12564-2008[EB]. 北京: 中国标准出版社, 2008: 16-17.

    Standardization Administration of China. General Specification for Photomultiplier Tubes: GB/T12564-2008[EB]. Beijing: China Standard Press, 2008: 16-17.
    [15] 隋成华, 杜春年, 徐丹阳. 基于卤钨灯和LED的复合光源设计与实现[J]. 浙江工业大学学报,2017,45(3):351-354. doi: 10.3969/j.issn.1006-4303.2017.03.023

    SUI Chenghua, DU Chunnian, XU Danyang. Design and implementation of compound light source based on halogen lamp and LEDs[J]. Journal of Zhejiang University of Technology,2017,45(3):351-354. doi: 10.3969/j.issn.1006-4303.2017.03.023
    [16] 张文超, 王宇松, 卢可义. 光电倍增管(PMT)的有源偏置电路[J]. 生命科学仪器,2005,3(1):35-37. doi: 10.3969/j.issn.1671-7929.2005.01.010

    ZHANG Wenchao, WANG Yusong, LU Keyi. Active voltage biasing circuit for PhotoMultiplier tube(PMT)[J]. Life Science Instruments,2005,3(1):35-37. doi: 10.3969/j.issn.1671-7929.2005.01.010
    [17] 胡孟春, 叶文英, 周殿忠, 等. 两种光电倍增管增益与总工作电压的关系研究[J]. 核电子学与探测技术,2004,24(3):239-241. doi: 10.3969/j.issn.0258-0934.2004.03.006

    HU Mengchun, YE Wenying, ZHOU Dianzhong, et al. Study of the relationship photomultiplier tube gain with sum work voltage[J]. Nuclear Electronics & Detection Technology,2004,24(3):239-241. doi: 10.3969/j.issn.0258-0934.2004.03.006
    • 相关文章
    • 施引文献(3)

  • 期刊类型引用(0)

    其他类型引用(3)

    • 资源附件(0)
图(6)  /  表(5)
计量
  • 文章访问数:  155
  • HTML全文浏览量:  53
  • PDF下载量:  45
  • 被引次数: 3
出版历程
  • 收稿日期:  2022-08-11
  • 修回日期:  2022-10-09
  • 网络出版日期:  2022-11-13
  • 刊出日期:  2022-11-14

目录

摘要
HTML全文

  • 引言

  • 1   信号模型

  • 2   最大似然估计

  • 3   克拉美-罗界分析

  • 4   仿真

  • 5   结论

参考文献(17)
相关文章
施引文献(3)
资源附件(0)

/

返回文章
返回

相关数据库

广告合作

友情链接

帮助中心

扫描关注本刊微信平台

您是第  2913802位访问者，欢迎！

地址: 西安市电子三路西段9号 邮政编码: 710065 联系电话: 029-88288172

版权所有 © 《应用光学》编辑部

陕公网安备 61011302001501号       陕ICP备19002862号-1

北京仁和汇智信息技术有限公司 开发  技术支持: info@rhhz.net   百度统计