CN105842848B - 哈特曼波前探测器电子倍增增益的自适应控制方法 - Google Patents

Abstract

本发明是自适应光学成像***中哈特曼波前探测器电子倍增增益的自适应控制方法，以使哈特曼波前探测器保持较高的信噪比同时避免亮度过饱和。针对被观测目标，探测到的光斑阵列中有一最高亮度值I_m，统计200帧～500帧中的I_m，计算其相对统计帧数的均值I_m相对的均方差值σ；设定理想最高亮度值比饱和亮度低3.0σ_m～3.2σ_m，其中σ_m为条件下的均方差值；再将光斑的亮度最大值相对一帧阵列中的光斑数计算平均值记为I_a，实时监测I_a；依据信号强度I与增益G的线性关系事先测得哈特曼波前探测器的技术特征参数k_a、k_m和理想最高亮度值即可以由当前帧探测数据I_a和增益G_a计算得到所需调整的增益值G_ms。增益自适应控制前的I_m和I_a曲线如图1(a)、控制后的I_m和I_a如图1(b)。

Description

哈特曼波前探测器电子倍增增益的自适应控制方法

技术领域

本发明属于自适应光学技术领域，是针对自适应光学成像***中哈特曼波前探测器电子倍增增益的控制方法。涉及自适应***对空间运动目标跟踪成像过程中目标亮度不断变化而对哈特曼波前探测器中的信号强度实时控制的技术，本质上是自适应控制哈特曼波前探测器背部的EMCCD的电子倍增增益。

背景技术

利用望远镜对天体目标进行观测时，由于大气湍流的随机干扰，成像光波前发生动态畸变，导致望远镜成像分辨率下降。连接于望远镜终端的自适应光学***可实时探测和校正大气湍流造成的畸变波前，以恢复望远镜理想的高分辨率成像。哈特曼波前探测器是当前自适应***中广泛应用的波前探测器。该探测器由前置的微透镜阵列与位于微透镜阵列焦平面上的相机组成，微透镜阵列将入射波前分割为子波前的阵列，使通过的光束在每个微透镜后的相机上会聚成点状光斑，计算光斑质心位置获得子波前的倾斜信息，进而重构出入射的整体畸变波前。图1展示了哈特曼波前探测器输出的一幅典型光斑阵列图像，微透镜会聚光斑位于离散的相机像素上，质心位置由像素上离散亮度分布的质心公式计算出来。

目标亮度越强信噪比越高，光斑的质心计算受噪声干扰越小，哈特曼波前探测器的探测精度越高；而目标亮度变弱时哈特曼波前探测器的探测精度下降，更严重的是大气湍流会造成接收口径内能量分布不均匀，某些光斑亮度会低于检出限，引起很大的波前探测误差。为了提高暗弱目标的波前探测精度，哈特曼波前探测器中的相机大多使用电子倍增型的EMCCD，通过设置适宜的电子倍增增益提高信噪比，减少弱暗光斑丢失的几率。

下面为叙述方便将“电子倍增增益”简称为“增益”。如何确定一个适宜的增益，是使用基于EMCCD的哈特曼波前探测器的关键问题，尤其运动目标在被跟踪成像过程中亮度发生较大变化时，设置固定的增益值不是造成局部信号亮度过饱和就是丢失弱暗光斑。例如，当观测围绕地球运动的人造卫星时，望远镜接收到的主要是人造卫星对太阳的反射光，反射光强随着人造卫星与太阳、地面望远镜的相对位置而变化【杨莉，在轨卫星星等建模与计算，应用光学，35(3)，365-370(2014)】，进而导致入射到哈特曼波前探测器中的信号光能量随之变化。如果设置的增益偏高有可能会在观测过程中出现EMCCD亮度过饱和而不能正常工作，如果设置的增益偏低又有可能会出现信噪比不足、局部光斑丢失甚至被观测目标丢失的情况。因此需随着被观测目标亮度的变化，自适应地控制哈特曼波前探测器的增益，以实现稳定的高信噪比，达到高精度波前探测的目的。

发明内容

本发明针对运动目标观测中亮度发生变化的问题，对哈特曼波前探测器的增益进行自适应控制，目的是使哈特曼波前探测器保持较高的信噪比，同时避免亮度过饱和。

下面详述本发明的方法。

由哈特曼波前探测器的技术说明书查出饱和亮度值I_L；针对一个亮度恒定的被观测目标，首先将哈特曼波前探测器的增益G_EM设置为初始增益G₀＝10，统计哈特曼波前探测器采集的每一帧光斑阵列中像素上的最高亮度值I_m的起伏变化，逐渐增大哈特曼波前探测器的增益，使监测到的最高亮度值I_m逐渐逼近并小于哈特曼波前探测器的饱和亮度值I_L；在此条件下统计200帧～500帧中像素上的最高亮度值I_m，计算其相对统计帧数的均值I_m相对最高亮度均值的均方差值σ_m；为使哈特曼波前探测器在采样中不出现过饱和且接近饱和，需限定最高亮度均值等于理想最高亮度值令等于3.0σ_m～3.2σ_m，其中σ_m为哈特曼波前探测器的技术特征参数，只要测得σ_m就能使值被确定下来，因此理想最高亮度值也是哈特曼波前探测器的技术特征参数。

为了即时控制哈特曼波前探测器的信号亮度，还需要寻找一个判断每帧信号亮度的参量。因此再将一帧光斑阵列上每个光斑在像素上的亮度最大值相对阵列中的光斑数计算平均值记为I_a，简称为平均光斑最大亮度I_a。I_a相对I_m要稳定得多，作为即时判断每帧信号亮度的监测参量。进一步找出I_a与的关系即可以通过即时监测每帧的平均光斑最大亮度I_a来控制增益使哈特曼波前探测器的成立。依据CCD光强I的线性方程【StanislavVitek,Jiri Hozman,Image quality influenced by selected image sensorparameters,Proc.SPIE9273,92732J(2003).】

I＝k·N_aper·G_EM (1)

其中I对应本发明中的I_a、I_m；k为相对I的常量系数，即I_a、I_m相应的k为k_a、k_m，看出k也是哈特曼波前探测器的技术特征参数，需要说明的是各参量的平均值以及特征值，如I_m的统计平均值和特征值都与I_m是同一性质的参量，故对应同一常量系数k，由于参量平均值的离散误差小，(1)式左侧的I应尽可能使用参量平均值；(1)式中的N_aper为对应被观测目标亮度等级的参量，是在G_EM＝1的条件下光斑阵列上每一个微透镜子孔径内所有像素上的亮度信号之和相对子孔径数的平均值N、再相对所采集的所有光斑阵列的N求平均值得出N_aper，简称N为子孔径亮度、N_aper为子孔径亮度平均值。令k_a为当前帧平均光斑最大亮度I_a的常量系数，哈特曼波前探测器当前帧的增益为G_a，k_m为的常量系数，对应条件的哈特曼波前探测器调整后的增益G_ms须满足下式：

由(2)式可以看出：事先测得哈特曼波前探测器的技术特征参数k_a、k_m和理想最高亮度值即可以由当前帧探测数据I_a和G_a计算得到所需调整的增益值G_ms，施加这一新的增益，即可完成一次增益控制。

为测量哈特曼波前探测器的技术特征参数k_a、k_m常量系数和理想最高亮度值做哈特曼波前探测器输出亮度信号I与EMCCD增益关系的检测，搭建如图2所示的检测光学***。检测光学***由电源1、光纤光源2、小孔光阑3、第一透镜4、第二透镜5、第三透镜6、哈特曼波前探测器7、湍流屏8、计算机9组成。光路中各元件的作用为：电源1为可编程信号发生器，用以给光纤光源2供电，控制其亮度变化，模拟运动目标在被观测中信号亮度I的变化；小孔光阑3的直径为100μm，使通过小孔光阑3的光成为点光源，该点光源位于第一透镜4的前焦点上，使通过第一透镜4的光束为平行光束；第二透镜5与第三透镜6共焦点，使通过第一透镜4的平行光束缩束为与哈特曼波前探测器7口径匹配的平行光束而全部进入哈特曼波前探测器7，湍流屏8定量模拟产生一定强度的大气湍流，模拟信号采集过程中大气湍流对信号亮度I的扰动；计算机9中装有控制哈特曼波前探测器7的控制软件，且与哈特曼波前探测器7相连，实时读取、处理哈特曼波前探测器7的亮度信号I，进行增益的调整。

哈特曼波前探测器7的技术特征参数k_a、k_m常量系数和理想最高亮度值的测量步骤如下：

1)遮光条件下测量哈特曼波前探测器7的暗背底，测量500帧，计算平均每帧像素上的亮度值分布并记录为背底图像。之后获得的每帧探测信号都需扣除背底图像。

2)使得光纤光源2出射刚刚超过哈特曼波前探测器7检出限的暗弱光且保持恒定，增益G_EM＝1，采集200帧～500帧光斑阵列图像，将每帧探测信号都扣除背底图像后计算每一帧平均光斑最大亮度I_a、最高亮度值I_m和子孔径亮度N，按照探测帧数排序绘制I_m和I_a的变化曲线如图3所示，看出I_m的变化起伏程度高于I_a，而I_a相对稳定，通过I_a判断当前帧光斑阵列的亮度更为合适；再统计平均连续200帧～500帧的平均光斑最大亮度记为最高亮度均值和子孔径亮度平均值N_aper。

3)保持增益G_EM＝1，设置光纤光源2的亮度逐渐增加，每改变一次亮度后，重复一遍2)中步骤，并记录所得到的系列和子孔径亮度平均值N_aper；绘制所设置增益G_EM下的分别随子孔径亮度平均值N_aper的变化曲线。

4)设置光纤光源2的亮度为步骤3)中所使用的最低亮度，然后设置不同的增益G_EM，每改变一次G_EM后，重复一遍3)中步骤，得到和两组线性关系如图4；图4(a)为不同增益G_EM下的一组直线，用每一直线斜率除以相应增益G_EM得出多个k_a的特例值，此处为6个，然后对6个特例值求平均得出较准确的k_a；图4(b)为不同增益G_EM下的一组直线，同样用每一直线斜率除以相应增益G_EM得出多个k_m的特例值，此处也为6个，然后对6个特例值求平均得出较准确的k_m。

5)接下来确定哈特曼波前探测器7的理想最高亮度值

由哈特曼波前探测器7的技术手册查出其背部的EMCCD的饱和亮度值I_L；再由上一步骤测量直线关系时获得的最大增益G_EM下的和200帧～500帧最高亮度值I_m，这些I_m已接近并小于哈特曼波前探测器的饱和亮度值I_L，利用这200～500个最高亮度值I_m和最高亮度均值计算出I_m相对最高亮度均值的均方差值σ_m，再令等于3.0σ_m～3.2σ_m，即可以获得哈特曼波前探测器7的理想最高亮度值

将上述步骤求得的k_m、k_a以及监测量I_a和相应的G_a代入(2)式即可求出当前应给哈特曼波前探测器7施加的增益G_ms。

进一步设置增益控制的其它条件：哈特曼波前探测器7背部的EMCCD每一次增益改变需要延迟时间数十毫秒，记为Δt，限定平均光斑最大亮度I_a的检测时间间隔T∈[Δt,Δt+10ms]的范围；另外监测量即平均光斑最大亮度I_a自身存在起伏，按照图3所示的I_a曲线，统计对应条件的I_a曲线起伏范围，记为ΔI_a，故设置增益调整控制判断条件为：当前帧的平均光斑最大亮度I_a和上一帧的平均光斑最大亮度I_a ^*之差满足|I_a-I_a ^*|≥3ΔI_a，当该条件成立时，需对当前增益G_a进行调整，否则无需对增益G_a进行调整；设置初始增益G₀＝10。

综上所述，哈特曼波前探测器7的增益自适应控制程序的思路为：

常量参数：ΔI_a、σ_m、Δt、k_m、k_a、ΔG_a；

G₀＝10；

循环：每隔Δt检测I_a；

判断是否|I_a-I_a ^*|≥3ΔI_a；

是，则依据(2)式计算G_ms，将当前增益值G_a调整为G_ms；否，则保持G_a；

进入下一个循环。

附图说明

图1是哈特曼波前探测器输出的一幅典型光斑阵列图像，微透镜会聚光斑于离散的相机像素上。

图2是测量哈特曼波前探测器的亮度线性方程常量系数的检测光学***，由电源1、光纤光源2、小孔光阑3、第一透镜4、第二透镜5、第三透镜6、哈特曼波前探测器7、湍流屏8、计算机9组成。其中电源1为可编程信号发生器，用以给光纤光源2供电，控制其亮度变化；小孔光阑3的直径为100μm，使通过小孔光阑3的光成为点光源，该点光源位于第一透镜4的前焦点上，使通过第一透镜4的光束为平行光束；第二透镜5与第三透镜6共焦点，使通过第一透镜4的平行光束缩束为与哈特曼波前探测器7口径匹配的平行光束而全部进入哈特曼波前探测器7，湍流屏8定量模拟产生一定强度的大气湍流；计算机9中装有控制哈特曼波前探测器7的控制软件，且与哈特曼波前探测器7相连，实时读取、处理哈特曼波前探测器7的探测信号强度I，进行增益的调整。

图3是光纤光源2的亮度恒定、哈特曼波前探测器7在一固定增益G_EM下探测信号强度I的波动曲线。其中曲线11和12分别为每帧中平均光斑最大亮度I_a与最高亮度值I_m的波动曲线，13为最高亮度均值14为I_m相对的均方差值σ。

图4是哈特曼波前探测器7在不同增益设置下、光纤光源2的亮度逐渐增加探测信号强度I与子孔径亮度平均值N_aper的线性关系，图中每一条线的斜率为亮度线性方程系数k与设置增益G_EM的乘积。其中(a)是统计平均的平均光斑最大亮度与子孔径亮度平均值N_aper的线性关系，直线21、22、23、24、25、26分别对应增益G_EM＝100、200、300、400、500、600，由(a)组直线的平均斜率能够求出k_a；(b)是最高亮度均值与子孔径亮度平均值N_aper的线性关系，直线31、32、33、34、35、36分别对应增益G_EM＝100、200、300、400、500、600，由(b)组直线的平均斜率能够求出k_m。

图5是模拟探测信号强度以正弦形迅速变化时通过自适应增益控制后的稳定信号强度的效果。其中(a)是G_EM＝1时正弦形迅速变化的探测信号强度I，曲线41为波动幅值约40ADU、峰值只有70ADU的平均光斑最大亮度I_a的曲线，曲线42为波动幅值约70ADU、峰值只有120ADU的最高亮度值I_m的曲线；(b)是自适应增益控制后正弦波动几乎祛除的稳定信号强度I的效果，曲线51为增益控制后强度升高至约6000ADU的平均光斑最大亮度I_a的曲线，曲线52为增益控制后强度升高至约11500ADU的最高亮度值I_m的曲线。

图6是在正弦形光信号输出条件下设置不同固定增益G_EM＝1、100、200和采用自适应增益控制的哈特曼波前探测器7中对应一个正弦周期的信噪比变化曲线，其中曲线61为增益G_EM＝1时的信噪比SNR曲线，曲线64为采用自适应增益控制的信噪比SNR曲线，垂直线62为增益G_EM＝100条件下哈特曼波前探测器7由于过饱和采集信号的截止线，所采集到的部分信号信噪比SNR曲线与曲线64重合，垂直线63为增益G_EM＝200条件下哈特曼波前探测器7由于过饱和采集信号的截止线，所采集到的部分信号信噪比SNR曲线也与曲线64重合。

具体实施方式

模拟被观测目标的光强度发生大幅度的迅速变化，采用本发明方法对哈特曼波前探测器增益进行自适应控制，提高弱信号时的探测信噪比。实施步骤及效果如下：

1)搭建如图2所示的模拟光路***，由电源1、光纤光源2、小孔光阑3、第一透镜4、第二透镜5、第三透镜6、哈特曼波前探测器7、湍流屏8、计算机9组成。

电源1为可编程信号发生器，用以给光纤光源2供电，控制其亮度变化；小孔光阑3的直径为100μm，使通过小孔光阑3的光成为点光源，该点光源位于第一透镜4的前焦点上，使通过第一透镜4的光束为平行光束；第二透镜5与第三透镜6共焦点，使通过第一透镜4的平行光束缩束为与哈特曼波前探测器7口径匹配的平行光束而全部进入哈特曼波前探测器7，湍流屏8模拟大气湍流对哈特曼波前探测器7中光斑亮度分布影响；计算机9中安装控制哈特曼波前探测器7的控制软件及自适应增益控制程序，且与哈特曼波前探测器7相连，实时读取、处理哈特曼波前探测器7的信号强度I并进行自适应增益控制。

光路***中各元件的技术参数为：

电源1，美国海洋光学公司，型号HL-2000-HP-232R，最大输出功率20W；

光纤光源2，大恒新纪元科技公司，型号为GCI-020107；

小孔光阑3，大恒新纪元科技公司，型号为GCO-01100；

第一透镜4、第二透镜5、第三透镜6的焦距分别为200mm、200mm、100mm；

哈特曼波前探测器7中的微透镜阵列由20×20个呈正方形排列的微透镜组成，其背部的EMCCD是英国e2v Technologies公司的产品，型号为CCD220，其饱和亮度值I_L＝2¹⁴ADU≈16300ADU，在超过8％的像素上出现亮度过饱和时增益自动重置为1，增益改变需要延迟时间Δt＝50ms，具有120×120个像素、像素尺寸为48μm，每一个微透镜对应6×6像素的子窗口；微透镜孔径d＝288μm，焦距f＝19.35mm，光斑的爱里斑直径为2.44λf/d＝90μm，其中波长λ＝550nm；

湍流屏8是美国Lexitek.Inc.公司的产品，型号为Near-Index-Match^TM phaseplate，自带控制盒，能够控制产生符合要求的大气湍流，此处设置湍流强度对应大气相干长度10cm、望远镜接收口径2m、格林伍德频率为50Hz。

2)遮光条件下测量哈特曼波前探测器7的暗背底，测量500帧，计算平均每帧像素上的亮度值分布并记录为背底图像，存入计算机9中，之后获得的每帧探测信号都需扣除此背底图像。

3)设置电源1输出较低的恒定电流，使得光纤光源2出射刚刚超过哈特曼波前探测器7检出限的暗弱光，设置哈特曼波前探测器7的增益G_EM＝1，保持光纤光源2的亮度恒定，采集500帧哈特曼波前探测器7实时探测的光斑阵列图像，存入计算机9中，每帧探测信号都扣除背底图像，计算每一帧平均光斑最大亮度I_a、最高亮度值I_m和子孔径亮度N；再统计连续500帧的平均光斑最大亮度记为最高亮度均值和子孔径亮度平均值N_aper。

4)保持增益G_EM＝1，设置光纤光源2的亮度逐渐增加，每改变一次亮度后，重复一遍3)中步骤，并记录所得到的系列和子孔径亮度平均值N_aper；绘制增益G_EM＝1条件下的分别随子孔径亮度平均值N_aper的变化曲线。

5)设置不同的增益G_EM＝100、200、300、400、500、600，每改变一次G_EM后，重复一遍4)中步骤；绘制图4(a)所示的一组直线；绘制图4(b)所示的第二组直线；利用增益G_EM＝600条件下获得的500帧光斑阵列的平均光斑最大亮度I_a、最高亮度值I_m和最高亮度均值计算此条件下的I_a起伏范围ΔI_a、I_m相对值的均方差值σ_m，得到ΔI_a≈80ADU、σ_m≈1000ADU。

6)利用图4(a)中的一组直线，用每一直线斜率除以相应增益G_EM得出增益G_EM＝100、200、300、400、500、600下的6个k_a的特例值，然后对6个特例值求平均得出较准确的k_a＝0.384；利用图4(b)的一组直线，同样用每一直线斜率除以相应增益G_EM得出G_EM＝100、200、300、400、500、600下的6个k_m的特例值，然后对6个特例值求平均得出较准确的k_m＝0.781。

7)接下来确定哈特曼波前探测器7的理想最高亮度值

依据哈特曼波前探测器7的技术手册，其背部的EMCCD的饱和亮度值I_L＝16300ADU；再由5)步骤求得的σ_m≈1000ADU，令即可以获得哈特曼波前探测器7的理想最高亮度值

将上述步骤求得的k_m、k_a以及监测量I_a和相应的G_a代入(2)式即可求出当前应给哈特曼波前探测器7施加的增益G_ms。

8)进一步设置增益控制的其它条件：哈特曼波前探测器7背部的EMCCD每一次增益改变需要延迟时间Δt＝50ms，限定平均光斑最大亮度I_a的检测时间间隔T＝55ms；由第5)步骤知ΔI_a≈80ADU，故设置增益调整控制判断条件为：当前帧的平均光斑最大亮度I_a和上一帧的平均光斑最大亮度I_a ^*之差满足|I_a-I_a ^*|≥240ADU；设置初始增益G₀＝10。

9)将所求得的哈特曼波前探测器7的技术特征参数

k_m＝0.781、k_a＝0.384、ΔI_a≈80ADU、T＝50ms，

代入哈特曼波前探测器7的增益自适应控制程序中，

设置初始增益G₀＝10；

自适应循环：每隔T检测一次I_a；

判断是否|I_a-I_a ^*|≥3ΔI_a；

是，则依据(2)式计算G_ms，将当前增益值G_a调整为G_ms；否，则保持G_a；

进入下一个循环。

10)设置电源1输出正弦形状变化电流，使光纤光源2输出的信号光强度发生大幅度的迅速变化，远超过实际被观测目标的光强变化范围与速度；设置固定增益G_EM＝1，记录哈特曼波前探测器7实时探测的光斑阵列信号，计算每帧中平均光斑最大亮度I_a和最高亮度值I_m，按帧序列排布的I_a和I_m均呈现为正弦形状的曲线，周期与电源1输出的正弦电流周期一致，G_EM＝1条件下的I_a和I_m曲线如图5(a)所示，其中曲线41为波动幅值约40ADU、峰值只有70ADU的平均光斑最大亮度I_a的曲线，曲线42为波动幅值约70ADU、峰值只有120ADU的最高亮度值I_m的曲线；设置G_EM＝100、200，分别在每个增益条件下重复上述的I_a和I_m的检测，但这两个条件下的最高亮度值I_m均在正弦信号的强度上升过程中出现过饱和，增益自动重置为1，未完成完整曲线的测量。

11)保持电源1输出正弦形状变化电流，记录哈特曼波前探测器7实时探测的光斑阵列信号，开启增益自适应控制程序，每隔Δt＝50ms计算一帧的平均光斑最大亮度I_a和最高亮度值I_m，再由I_a计算出理想最高亮度值对应的增益G_ms，自适应控制增益，使最高亮度值I_m和平均光斑最大亮度I_a曲线中的正弦波动几乎祛除，如图5(b)所示，其中曲线51为增益控制后平均强度升高至6000ADU的平均光斑最大亮度I_a的曲线，曲线52为增益控制后平均强度升高至11500ADU的最高亮度值I_m的曲线。

计算了在上述正弦形光信号输出条件下设置不同固定增益G_EM＝1、100、200和采用自适应增益控制的哈特曼波前探测器7中的信噪比SNR，如图6所示，其中曲线61为增益G_EM＝1时的信噪比SNR曲线，曲线64为采用自适应增益控制的信噪比SNR曲线，垂直线62为增益G_EM＝100条件下哈特曼波前探测器7由于过饱和采集信号的截止线，所采集到的部分信号信噪比SNR曲线与曲线64重合，垂直线63为增益G_EM＝200条件下哈特曼波前探测器7由于过饱和采集信号的截止线，所采集到的部分信号信噪比SNR曲线也与曲线64重合。图6的结果说明，对应光纤光源2的最低输出信号I_a＝35ADU时，G_EM＝1增益下的信噪比SNR只有12，而G_EM达到100时信噪比SNR达到上限22，G_EM再增至200信噪比SNR仍为22；随着光纤光源2的输出信号增强，G_EM＝1时的信噪比SNR最高增至22，而G_EM为100和200固定增益的信号强度先后出现过饱和而使信噪比SNR曲线垂直落下与G_EM＝1的SNR曲线重合，只有自适应增益控制的信噪比SNR升至30，可见自适应增益控制不仅使弱亮度目标的探测信噪比有效提高，使整个探测过程保持于信噪比的上限，同时避免了基于EMCCD的哈特曼波前探测器的信号过饱和。

Claims (2)

1.哈特曼波前探测器电子倍增增益的自适应控制方法，其特征是：

由哈特曼波前探测器的技术说明书查出饱和亮度值I_L；针对一个亮度恒定的被观测目标，首先将哈特曼波前探测器的增益G_EM设置为初始增益G₀＝10，统计哈特曼波前探测器采集的每一帧光斑阵列中像素上的最高亮度值I_m的起伏变化，逐渐增大哈特曼波前探测器的增益，使监测到的最高亮度值I_m逐渐逼近并小于哈特曼波前探测器的饱和亮度值I_L；在此条件下统计200帧～500帧中的最高亮度值I_m，计算其相对统计帧数的均值I_m相对最高亮度均值的均方差值σ_m；

设定最高亮度均值的理想值为令等于3.0σ_m～3.2σ_m；

再将一帧光斑阵列上每个光斑的亮度最大值相对阵列中的光斑数计算平均值记为I_a，简称为平均光斑最大亮度I_a，作为即时判断每帧信号亮度的监测参量；

进一步找出I_a与的关系即可以通过监测每帧的平均光斑最大亮度I_a来控制增益使哈特曼波前探测器的信号强度成立；依据CCD光强I的线性方程：

I＝k·N_aper·G_EM (1)

其中k为相对I的常量系数，I对应本发明中的I_a、I_m、σ，相应的k为k_a、k_m、k_σ，G_EM为设置的增益，N_aper为对应被观测目标亮度等级的参量，是在G_EM＝1的条件下一个微透镜子孔径内所有像素上的亮度信号之和相对一帧所有子孔径的平均值，简称N_aper为子孔径亮度平均值；在G_EM＝1的条件下，可以测量直线关系，此直线的斜率即为相应的常量系数k_a，也可以测量直线关系，直线的斜率即为相应的常量系数k_m；

令当前帧平均光斑最大亮度为I_a、哈特曼波前探测器当前帧的增益为G_a，对应条件的哈特曼波前探测器调整后的增益G_ms须满足下式：

<mrow> <mfrac> <msub> <mi>I</mi> <mi>a</mi> </msub> <msub> <mover> <mi>I</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mi>m</mi> <mi>s</mi> </mrow> </msub> </mfrac> <mo>=</mo> <mfrac> <msub> <mi>k</mi> <mi>a</mi> </msub> <msub> <mi>k</mi> <mi>m</mi> </msub> </mfrac> <mo>&amp;CenterDot;</mo> <mfrac> <msub> <mi>G</mi> <mi>a</mi> </msub> <msub> <mi>G</mi> <mrow> <mi>m</mi> <mi>s</mi> </mrow> </msub> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>

事先测得哈特曼波前探测器的技术特征参数k_a、k_m和理想最高亮度值即可以由当前帧探测数据I_a和G_a计算得到所需调整的增益值G_ms；

增益控制的其它条件：哈特曼波前探测器(7)背部的EMCCD每一次增益改变需要延迟时间数十毫秒，记为Δt，限定平均光斑最大亮度I_a的检测时间间隔T∈[Δt,Δt+10ms]的范围；另外监测量即平均光斑最大亮度I_a自身存在起伏，记为ΔI_a，设置增益调整控制判断条件为：当前帧的平均光斑最大亮度I_a和上一帧的平均光斑最大亮度I_a ^*之差满足|I_a-I_a ^*|≥3ΔI_a，当该条件成立时，需对当前增益G_a进行调整，否则无需对增益G_a进行调整。

2.根据权利要求1所述的哈特曼波前探测器电子倍增增益的自适应控制方法，其特征是所述的哈特曼波前探测器的技术特征参数k_a、k_m常量系数和理想最高亮度值的测量方法如下：

搭建检测光学***；检测光学***由电源(1)、光纤光源(2)、小孔光阑(3)、第一透镜(4)、第二透镜(5)、第三透镜(6)、哈特曼波前探测器(7)、湍流屏(8)、计算机(9)组成；光路中各元件的作用为：电源(1)为可编程信号发生器，用以给光纤光源(2)供电，控制其亮度变化，模拟运动目标在被观测中信号亮度I的变化；小孔光阑(3)的直径为100μm，使通过小孔光阑(3)的光成为点光源，该点光源位于第一透镜(4)的前焦点上，使通过第一透镜(4)的光束为平行光束；第二透镜(5)与第三透镜(6)共焦点，使通过第一透镜(4)的平行光束缩束为与哈特曼波前探测器(7)口径匹配的平行光束而全部进入哈特曼波前探测器(7)，湍流屏(8)定量模拟产生一定强度的大气湍流，模拟信号采集过程中大气湍流对信号亮度I的扰动；计算机(9)中装有控制哈特曼波前探测器(7)的控制软件，且与哈特曼波前探测器(7)相连，实时读取、处理哈特曼波前探测器(7)的亮度信号I，进行增益的调整；

基于所搭建的检测***，哈特曼波前探测器(7)的技术特征参数k_a、k_m常量系数和理想最高亮度值的测量步骤如下：

1)遮光条件下测量哈特曼波前探测器(7)的暗背底，测量500帧，计算平均每帧像素上的亮度值分布并记录为背底图像，之后获得的每帧探测信号都需扣除背底图像；

2)使得光纤光源(2)出射刚刚超过哈特曼波前探测器(7)检出限的暗弱光且保持恒定，增益G_EM＝1，采集200帧～500帧光斑阵列图像，将每帧探测信号都扣除背底图像后计算每一帧平均光斑最大亮度I_a、最高亮度值I_m和子孔径亮度N，再统计平均连续200帧～500帧的平均光斑最大亮度记为最高亮度均值和子孔径亮度平均值N_aper；

3)保持增益G_EM＝1，设置光纤光源(2)的亮度逐渐增加，每改变一次亮度后，重复一遍2)中步骤，并记录所得到的系列和子孔径亮度平均值N_aper；绘制所设置增益G_EM下的分别随子孔径亮度平均值N_aper的变化曲线；

4)设置光纤光源(2)的亮度为步骤3)中所使用的最低亮度，然后设置不同的增益G_EM，每改变一次G_EM后，重复一遍3)中步骤，得到两组线性关系；先针对组线性关系，用每一直线斜率除以相应增益G_EM得出一个k_a的特例值，然后对全组的特例值求平均得出较准确的k_a；再用同样方法从组线性关系中求得k_m；

5)由哈特曼波前探测器(7)的技术手册查出其背部的EMCCD的饱和亮度值I_L；再由上一步骤测量直线关系时获得的最大增益G_EM下的和200帧～500帧最高亮度值I_m，这些I_m已接近并小于哈特曼波前探测器的饱和亮度值I_L，利用这200～500个最高亮度值I_m和最高亮度均值计算出I_m相对最高亮度均值的均方差值σ_m，再令等于3.0σ_m～3.2σ_m，即可以获得哈特曼波前探测器(7)的理想最高亮度值