本发明涉及时间相关单光子计数深度成像
技术领域：
，具体为一种光子计数目标高效率高分辨力深度重建方法。
背景技术：
：近年来，光子计数型伪随机码幅度调制二维成像方法的研究已经广泛展开，成像过程中，存在以下两个问题：1)通过两路光子到达时间点的最佳检测运算，得到扩频光子计数值分布，这就是传统的基于最佳检测的相关接收技术。为了得到高精度的扩频光子计数值分布，往往通过延长积分时间来提高信噪比，以至于增大了采集的数据量、加重数据处理负担、降低成像速度和信号恢复效率。在特殊环境下，存在后向散射噪声光的干扰，而信号光衰减加剧，光子计数率很低，通常为每秒几个或更少，由于信息量的不足，难以恢复具有尖锐峰值的扩频光子计数值分布。2)构建系统的单光子探测器、时间到达记录仪等设备引入系统的瞬时干扰，该瞬时干扰可采用系统的冲激响应来描述。以卷积积分为理论依据，即输入信号通过线性系统冲激响应的移位、倍乘、叠加运算。信号通过系统后计算获得扩频光子计数值分布，该波形包含系统的瞬时干扰项(tamerf.refaat，syedismail，m.nurulabedin，scottm.spuler，shaned.mayor，andupendran.singh.lidarbackscattersignalrecoveryfromphototransistorsystematiceffectbydeconvolution[j].appl.opt，2008，47:5281-5288.joelf.campbell，binglin，aminr.nehrir，f.wallaceharrison，andmichaeld.obland.super-resolutiontechniqueforcwlidarusingfouriertransformreorderingandrichardson–lucydeconvolution[j].opt.lett，2014，39:6981-6984.)，导致失真的扩频光子计数值分布，使得距离分辨率往往无法达到伪随机码发送速率决定的距离分辨率，例如，码发送速率为1ghz，则最小距离单元为一个码元宽度1ns，对应距离为15厘米。实际探测中，可分辨的最小距离达不到15厘米。针对复杂的目标场景，单个光斑内往往包含多个距离的反射光子，当目标在横向和纵向距离上较为接近时，无法分辨多个距离，造成目标探测的遗漏，目标成像的细节性不高等问题，如图1所示，图1-(a)表示在激光光斑内，纵向方向上存在多个距离分布，探测器接收到两种距离的返回回波光子，纵向间隔距离分布比较密集时，目标的光子统计直方图会产生重叠如图1-(c)上方，采用峰值搜索法提取扩频光子计数值分布的峰值时，往往会遗漏部分被模糊的距离；光斑大小不变，在纵向分辨能力的限制下，图1-(b)中，当对横向方向上两个距离分布比较密集的目标成像时，无法分辨两个目标，如图1-(c)右下方所示。技术实现要素：本发明的目的在于提供一种光子计数目标高效率高分辨力深度重建方法，基于一体化光纤式伪随机码幅度调制偏移校正装置(申请号为201710325775.x)的成像系统，后期数据处理中，采用压缩感知的方法高校恢复扩频光子计数值分布，省去单独的噪声去除步骤，以最快的速度方向，加速估计下一次迭代的数值，运算速度快于线性搜索的迭代算法。为达到上述目的，根据本发明的一个方面，本发明提供如下技术方案：1.一种光子计数目标高效率高分辨力深度重建方法，基于一体化光纤式伪随机码幅度调制深度矫正装置的成像系统，其特征在于，包括以下步骤：s1、将成像系统的准直器固定在二维导轨平台上，去掉a单光子探测器、固定衰减器和分路器，伪随机序列发生器产生一路同步高电平脉冲，该脉冲产生时间和开始发送伪随机序列的时间一致，将同步高电平脉冲连接光子时间到达记录仪，作为开始计时的触发信号，伪随机序列发生器驱动激光器发送激光脉冲，激光脉冲信号通过成像系统后得到单点扩频光子计数值分布，表达为系统冲激响应和携带目标特征的目标光子计数值分布的卷积积分；s2、改变s1中的成像系统结构，将两个准直器面对面发射并接收光子，获取平均的系统冲激响应并拟合；s3、恢复成像系统结构到初始状态，基于最佳检测的运算机制，通过稀疏基表示方法，高效反演目标的扩频光子计数值分布；s4、采用基于贝叶斯理论和最大似然估计的迭代重建算法，引入矢量外推加速算法，估计最终携带目标特征的目标光子计数值分布；s5、根据最终携带目标特征的目标光子计数值分布，通过残差去除﹑峰值搜索和质心拟合方法提取距离数值；s6、逐点重构距离数值，得到重构的二维深度图像。本发明进一步设置为，所述步骤s1将准直器固定在二维导轨平台上，去掉a单光子探测器、固定衰减器和分路器，伪随机序列发生器产生一路同步高电平脉冲，该脉冲产生时间和开始发送伪随机序列的时间一致，将同步高电平脉冲连接光子时间到达记录仪，作为开始计时的触发信号，伪随机序列发生器驱动激光器发送激光脉冲，激光脉冲信号通过成像系统后得到的单点扩频光子计数值分布，表达为系统冲激响应和携带目标特征的目标光子计数值分布的卷积积分，具体为，将准直器固定在二维导轨平台上，去掉a单光子探测器、固定衰减器和分路器，伪随机序列发生器产生一路同步高电平脉冲，该脉冲产生时间和开始发送伪随机序列的时间一致，将同步高电平脉冲连接光子时间到达记录仪，作为开始计时的触发信号，伪随机序列发生器驱动激光器发送激光脉冲，激光脉冲发射后，先后通过大气介质、目标和接收系统，在短距离下在到达接收系统前，脉冲回波信号s(n)由目标冲激响应决定；通过接收系统后，目标冲激响应ho(n)和探测器的冲激响应hd(n)是造成接收的目标信号产生瞬时失真的主要因素，则伪随机光脉冲碰到目标后，往返一次，回波到达接收系统，探测到的单点扩频光子计数值分布r(n)为，式中，r(n)为携带目标信息并且被系统冲激响应模糊的单点扩频光子计数值分布，c(n)为互相关函数，为系统冲激响应，表征系统自身产生的瞬时干扰，f(n)＝ho(n)为携带目标特征的目标光子计数值分布，b为噪声光子。本发明进一步设置为，所述步骤s2改变s1中的成像系统结构，将两个准直器面对面发射并接收光子，获取平均的系统冲激响应并拟合，具体为，s2.1、在成像系统上，断开可调衰减器和光环行器的连接，将准直器b和可调衰减器连接，作为收准直器端，将连接光环行器的准直器a作为发准直器端，将准直器a和准直器b面对面放置，通过微调和粗调准直器的位置，使得从准直器a中发出的光信号导入接收的准直器b中，再由b单光子探测器接收后，得到接收路光子时间到达点；s2.2、将matlab的序列重构方法嵌入labview，基于labview，获取接收路的光子时间到达点，并重构0﹑1序列，表示为y(n)，伪随机序列发生器中存储的模板序列p(n)与y(n)进行快速傅里叶变换，得到扩频光子计数值分布，即系统冲激响应h，其代表系统的瞬时干扰，式中，表示傅里叶变换，为傅里叶逆变换；s2.3、测量若干次的系统冲激响应并作平均，得到平均的系统冲激响应；s2.4、采用高斯拟合法拟合平均的系统冲激响应，设平均的系统冲激响应数据(oi,hi)i＝1,2,3,…,n，用高斯函数描述为，式中，hmax、omax和s为待估参数，hmax代表高斯曲线的峰高，omax代表高斯曲线峰位置，s代表半高全宽；将(3)式两边取自然对数，化为，令则(4)式化为二次多项式拟合函数，考虑全部数据和测量误差εi，并以矩阵形式表示如下，简记为zn×1＝on×3b3×1 en×1，根据最小二乘原理，求得拟合常数b0﹑b1﹑b2构成的矩阵b的广义最小二乘解为，b＝(oto)-1otz(8)进而根据(5)式，求出待估参数hmax﹑omax和s，代入(3)式，从而得到拟合的高斯函数。本发明进一步设置为，所述步骤s3恢复成像系统结构并成像，基于最佳检测的运算机制，通过稀疏基表示方法，高效反演目标的扩频光子计数值分布，具体为，恢复成像系统结构，对目标成像，获得接收路光子到达点时间，s3.1、根据步骤s2.2的序列重构方法，重构接收路光子到达时间点，一个周期的序列长度为na，记为s3.2、生成数据压缩比α，令m＝αna，以此为参量，随机生成测量矩阵s3.3、基于最佳检测的运算机制，生成稀疏基ψ，基于最佳检测的运算机制，采用一路伪随机序列p(n)作为发送的参考码型，另一路为接收的经过延时后的光子到达时间点，重构后得到序列x(n)，两个序列最佳检测的运算表示为，其中，na为一个周期的序列长度，m为时间单元，扩频光子计数值分布r(m)的峰值对应横坐标即为时间延时，其代表目标所处的距离信息，基于最佳检测的运算方法，以参考序列为模板序列，拟采用包含时间位移的hankel矩阵结构，构造稀疏基ψ，该稀疏基的列元素包含所有可能的参考序列的时间位移，分以下两种情况建立，1)如果积分的时间ta和序列的周期tp相等，即模板序列和接收序列的长度相等，则构建的稀疏基表示为，其中，ψ(:,i)代表稀疏基的第i列向量，zeros代表零向量，[]t代表矩阵的转置，p(i:j)为从第i个元素到第j个元素的参考向量；3)当积分时间ta大于序列的周期tp，向量p右侧补零，直至p和x长度相等；s3.4、计算稀疏基的相干性δ(ψ)采用δ(ψ)表示稀疏基的相干性，最劣的相干性即为两个不同稀疏基元素之间内积的绝对值的最大值，δ(ψ)的取值范围为[0，1]，计算式(11)以防止出现相干的稀疏基元素对，式中，ψi为第i列元素构成的向量，ψj为第j列元素构成的向量；s3.5、重构扩频光子计数值分布设测量信号为为接收的经过延时后的光子到达时间信号，为变换系数，即待恢复的扩频光子计数值分布，其描述光子在时间轴上的分布特征，其具有尖锐的峰值，在时间分布上具有稀疏特性，即||r||0≤k，是k稀疏的，则有在稀疏基的表示下是k稀疏的，基于压缩感知的测量过程理解为对经过延时后的光子到达时间信号的m次独立测量，每次以线性函数的形式测量接收数据的线性投影m个测量向量构成的测量矩阵为其将高维信号x映射到低维，则测量过程表述为y＝φx＝φ(ψr)(12)设v表示噪声，带噪观测情况下，式(12)表示为，y＝ar v(13)其中a＝φψ；s3.5.1式(13)转化为以下lasso标准问题的求解，λ为代价参数，s3.5.2计算rk＝(ata ρi)-1[aty ρ(zk-1-uk-1)]，其中rk为第k次迭代计算得到的扩频光子计数值分布，即式(1)中的r(n)，ρ为惩罚系数，i为单位矩阵，()-1为矩阵的逆，()t为矩阵的转置，zk-1和uk-1均为第k-1次迭代的长度为2na-1的向量；s3.5.3计算其中prox为关于变量z的l1范数近端算子，具有显式解；s3.5.4计算uk＝uk-1 rk-zk；s3.6、判断δ(ψ)＜ε是否小于ε，如果小于ε，则k＝k 1，继续采用s3.5的步骤迭代更新r，如果大于设定数值，则停止迭代，输出r。本发明进一步设置为，所述步骤s4采用基于贝叶斯理论和最大似然估计的迭代重建算法，引入矢量外推加速算法，估计最终携带目标特征的目标光子计数值分布，具体为，s4.1、初始化拟合的冲激响应﹑最大迭代次数kmax和加速参数α；s4.2、设携带目标特征的目标光子计数值分布为f，简称为目标光子计数值，由贝叶斯模型可得，式(14)中，p(f/r)为已知扩频光子计数值分布情况下，目标光子计数值的概率分布情况，p(r/f)为已知目标光子计数值的情况下，扩频光子计数值概率分布情况，p(f)为目标光子计数概率分布；当p(f/r)取得最大数值，则恢复的目标信号接近理想的波形，式(14)的分母部分为常量，通过最大似然估计法求解其最大值，ml(f)＝max[p(r/f)](15)式(15)中，ml(f)表示采用最大似然估计的方法估计f，max为最大值，设待恢复的波形，即目标光子计数值f，其与系统的冲激响应h的卷积表示为t表示每个像素的积分时间，则联合概率密度分布l(r/f)为，式(16)等号两边求对数得到，对等式(17)两边求导，并令其为0，根据退化函数归一化的性质，可得，引入迭代运算，则有，式(19)是迭代计算公式的一般形式，其中，ψ(fk)表示线性搜索的迭代算法，对于一维扩频光子计数值分布，fk表示第k次迭代计算时产生的目标光子计数值的估计值，fk 1表示第k 1次迭代计算时产生的目标光子计数值的估计值；s4.3、基于差分表示的估计值变化引入根据式(19)从现有的估计值fk-1产生估计值fk的迭代出发，给出一阶泰勒矢量，fk＝ψ(fk-1)＝fk-1 gk(20)设fk逐渐收敛到第s次迭代计算时产生的估计值fs，在第k次迭代中引入的变化gk朝着解的梯度方向▽m(f)生成一条通向多维空间的解，目标光子计数值的求解过程中，需要最大化的函数可局部近似为多维二次函数m(f)，则有，m(f)＝-(f-fs)ta(f-fs)(21)▽m(f)＝-2a(f-fs)(22)式(22)中，()t表示转置，a表示平方、对称、正定矩阵，根据式(21)和(22)，采用差分形式给出第k次迭代引入的变化gk，如式(23)所示，gk＝ε▽m(fk)＝-2εa(fk-fs)(23)式(19)中，根据矩阵对角化原理，通过矩阵a的特征值构成的对角矩阵λa和特征向量s给出矩阵a，即，a＝sλas-1(24)设i为单位矩阵，λb是矩阵b的特征值构成的对角矩阵，b＝i-2εa,λb＝i-2ελa，f0表示估计值的初值，λbk为λb的k次方，ε表示一个用来防止λb出现负值的向量，根据式(20)、式(23)和式(24)，fk为，根据式(25)得到第k-1次迭代引入的变化gk-1，k-δ-1次迭代引入的变化gk-1-δ，k-δ次迭代的估计值fk-δ，gk-1＝sλbk(i-λb-1)s-1(f0-fs) fs(26)fk-δ＝sλbk-δs-1(f0-fs) fs(27)gk-δ-1＝sλbk-δ(i-λb-1)s-1(f0-fs)(28)式(28)中，λbk-δ为λb的k-δ次方，且有i-λb-δ≈δ(i-λb-1)；s4.4、计算估计值之差估计值fk和估计值fk-δ之差hk,-δ表示为，hk,-δ＝fk-fk-δ＝sλbk(i-λb-δ)s-1(f0-fs)(29)对于第k δ次迭代的估计值fk δ的预测为，fk δ＝sλbk δs-1(f0-fs) fs(30)第k δ-1次迭代引入的变化gk δ-1为，gk δ-1＝sλbk δ(i-λb-1)s-1(f0-fs)(31)估计值fk δ和估计值fk之差hk,δ为，hk,δ＝fk δ-fk＝sλbk δ(i-λb-δ)s-1(f0-fs)(32)式(32)中，λbk δ为λb的k δ次方，λb-δ为λb的-δ次方，λb-1为λb的-1次方，s-1为s的-1次方；s4.5、估计加速参数基于泰勒级数的一阶矢量外推法，通过最小二乘投影αkhk,-δ预测hk,δ，得，hk,δ＝αkhk,-δ(33)(hk,δ-αkhk,-δ)thk,-δ＝0(34)式(34)中，αk为加速参数，求解式(34)，得到，又因为i-λb-δ≈δ(i-λb-1),将i-λb-δ≈δ(i-λb-1)代入式(32)可得，hk,δ≈δsλbk(i-λb-1)s-1(f0-fs)＝δgk-1(36)以此类推得到，hk,-δ≈δgk-δ-1(37)将式(35)和式(37)代入式(35)，可得加速参数αk为，s4.6、估计携带目标特征的目标光子计数值设第k次待估计的目标光子计数值为lk，根据第k次估计中已经得到的fk、fk-1、αk，代入式lk＝fk αk(fk-fk-1)，得到的lk，k 1后，将lk代入得到fk 1，同理通过式(20)得到gk 1，通过式(38)得到αk 1，再代入lk＝fk αk(fk-fk-1)的变换式lk 1＝fk 1 αk 1(fk 1-fk)，得到lk 1，以此类推，对lk进行循环更新，直到达到最大迭代次数kmax，获得最终携带目标特征的目标光子计数值lkmax。本发明进一步设置为，所述步骤s5根据最终携带目标特征的目标光子计数值分布，通过残差去除﹑峰值搜索和质心拟合方法提取距离数值，具体为，s5.1、通过最终携带目标特征的目标光子计数值lkmax，最高幅度值选取阈值β，提取出满足lkmax＜β对应的所有k值，对应k的lkmax幅度值置为0，得到的数据记为lkmax_opt，即第k个深度单元对应的目标光子计数值；s5.2、搜索lkmax_opt中的峰值，并且提取峰值对应的k值；s5.3、在搜索到的峰值附近提取l个非零数值点，采用质心拟合算法计算距离数值d，式(39)中，τ为最小的距离分辨单元。本发明进一步设置为，所述步骤s6逐点重构距离数值，得到重构的二维深度图像，具体为，将准直器固定在二维导轨平台上，扫描获得每个像素点的光子到达时间点，采用步骤s3～步骤s5所述的方法，逐点重构距离数值，得到重构的二维深度图像。与现有技术相比，本发明具有的优势是：提出将压缩感知信号恢复方法引入基于最佳检测的单光子计数成像方法中，确立稀疏表示的信息高效恢复技术。通过稀疏基稀疏表达隐藏的延时信息，利用最佳检测的运算核心结构来构建稀疏基，并提出双相干性约束，进一步保证信息恢复的效率。该方法突破成像系统中前端接收系统的物理性能极限，省去冗余数据的存取，降低计算负担。丰富了压缩感知理论和最佳检测理论的科学内涵，为信息的高效恢复提供新的应用价值。并且本发明省去单独的噪声去除步骤，以最快的速度方向，加速估计下一次迭代的数值，运算速度快于线性搜索的迭代算法，具备节约硬件资源、降低能源损耗、提高成像效率的优势；提高了深度成像的横向和纵向分辨能力。附图说明图1为本发明成像目标结构示意图；图2为实时数据处理流程图，图2-(a)labview接收序列数据存取，图2-(b)matlab实时接收序列重构算法；图3为高效率高分辨力成像方法流程；图4为系统冲激响应图，其中实线为测量的平均的系统冲激响应，虚线为拟合的冲激响应；图5为最小均方误差对应最大迭代次数的确定图。具体实施方式下面结合说明书附图，对本发明作进一步的说明。与现有技术相比，本发明具有的优势是：提出将压缩感知信号恢复方法引入基于最佳检测的单光子计数成像方法中，确立结构化稀疏表示的信息高效恢复技术。利用最佳检测的算法核心结构来构建稀疏基，通过稀疏基稀疏表达隐藏的延时信息，并提出双相干性约束，进一步保证信息恢复的效率。该方法突破成像系统中前端接收系统的物理性能极限，省去冗余数据的存取，降低计算负担。丰富了压缩感知理论和最佳检测的理论的科学内涵，为信息的高效恢复提供新的应用价值。并且本发明省去单独的噪声去除步骤，以最快的速度方向，加速估计下一次迭代的数值，运算速度快于线性搜索的迭代算法，具备节约硬件资源、降低能源损耗、提高成像效率的优势；提高了深度成像的横向和纵向分辨能力。一种光子计数目标高效率高分辨力深度重建方法，基于一体化光纤式伪随机码幅度调制偏移校正装置(申请号为201710325775.x)的成像系统，包括以下步骤：s1、将成像系统的准直器固定在二维导轨平台上，去掉a单光子探测器、固定衰减器和分路器，伪随机序列发生器产生一路同步高电平脉冲，该脉冲产生时间和开始发送伪随机序列的时间一致，将同步高电平脉冲连接光子时间到达记录仪，作为开始计时的触发信号，伪随机序列发生器驱动激光器发送激光脉冲，激光脉冲信号通过成像系统后得到单点扩频光子计数值分布，表达为系统的冲激响应和携带目标特征的目标光子计数值分布的卷积积分，具体为，将准直器固定在二维导轨平台上，去掉a单光子探测器、固定衰减器和分路器，伪随机序列发生器产生一路同步高电平脉冲，该脉冲产生时间和开始发送伪随机序列的时间一致，将同步高电平脉冲连接光子时间到达记录仪，作为开始计时的触发信号，伪随机序列发生器驱动激光器发送激光脉冲，激光脉冲发射后，先后通过大气介质、目标和接收系统，在短距离下在到达接收系统前，脉冲回波信号s(n)由目标冲激响应决定；通过接收系统后，目标冲激响应ho(n)和探测器的冲激响应hd(n)是造成接收的目标信号产生瞬时失真的主要因素，则伪随机光脉冲碰到目标后，往返一次，回波到达接收系统，探测到的单点扩频光子计数值分布r(n)为，式中，r(n)为携带目标信息并且被系统的冲激响应模糊的单点扩频光子计数值分布，c(n)为互相关函数，为系统的冲激响应，表征系统自身产生的瞬时干扰，f(n)＝ho(n)为携带目标特征的目标光子计数值分布，b为噪声光子。s2、改变s1中的成像系统结构，将两个准直器面对面发射并接收光子，获取平均的系统冲激响应并拟合，具体为，s2.1、在一体化光纤式伪随机码幅度调制深度矫正装置的成像系统上，断开可调衰减器和光环行器的连接，将准直器b和可调衰减器连接，作为收准直器端，将连接光环行器的准直器a作为发准直器端，将准直器a和准直器b面对面放置，通过微调和粗调准直器的位置，使得从准直器a中发出的光信号导入接收的准直器b中，再由b单光子探测器接收后，得到接收路光子时间到达点；s2.2、基于labview，获取接收路的光子时间到达点如图2(a)，并采用matlab重构0﹑1序列如图2(b)，从获取的时间到达点中重构发序列码型，根据码发送速率和系统最小深度分辨率，比如码发送速率为2.5ghz，则码元宽度为400ps，时间相关记录仪分辨率为4ps，则时间分辨因子为400ps/4ps＝100，根据该通道的第一个时间到达点，通过时间分辨因子确定第一个‘1’码元的位置，再依次确定第二个、第三个的‘1’码元。接着，判断两个‘1’码元，即两个时间到达点之间序列零的个数，最后通过序列的对比和补零，得到长度一致的收发序列，基于matlab实现重构算法，获得一定积分时间下的序列码型。通过循环发送序列，一路是已经生成好的模板序列，一路是首尾相接的接收序列码型。伪随机码是循环发送的，通过接收路的光子时间到达点重构0﹑1序列，表示为y(n)，模板序列p(n)与y(n)进行快速傅里叶变换，得到扩频光子计数值分布，即系统冲激响应h，其代表系统的瞬时干扰，式中，表示傅里叶变换，为傅里叶逆变换；s2.3、测量50次的系统冲激响应并作平均，得到平均的系统冲激响应；s2.4、采用高斯拟合法拟合平均的系统冲激响应，设平均的系统冲激响应数据(oi,hi)i＝1,2,3,…,n，用高斯函数描述为，式中，hmax、omax和s为待估参数，hmax代表高斯曲线的峰高，omax代表高斯曲线峰位置，s代表半高全宽；将(3)式两边取自然对数，化为，令则(4)式化为二次多项式拟合函数，考虑全部数据和测量误差εi，并以矩阵形式表示如下，简记为zn×1＝on×3b3×1 en×1，根据最小二乘原理，求得拟合常数b0﹑b1﹑b2构成的矩阵b的广义最小二乘解为，b＝(oto)-1otz(8)进而根据(5)式，求出待估参数hmax﹑omax和s，代入(3)式，从而得到拟合的高斯函数。图4实线为测量得到的平均的系统冲激响应，采用高斯拟合法拟合平均的系统冲激响应如图4虚线，拟合曲线的半高全宽(fwhm)约为1.4ns，意味着仅仅采用峰值搜索法可分辨的目标的最小纵向间隔约为21厘米。再加上背景噪声、复杂的目标属性等外界条件，使得两个距离峰值不断融合，距离波形被模糊，有必要引入信号处理算法，高效去除系统的瞬时干扰，提高距离分辨率。获得了拟合的冲激响应后，二维图像的重构流程如图3所示。s3、恢复成像系统结构并成像，基于最佳检测的运算机制，通过稀疏基表示方法，高效反演目标的单点扩频光子计数值分布，具体为，恢复成像系统结构，对目标成像，获得接收路光子到达点时间，s3.1、根据步骤s2.2的序列重构方法，重构接收路光子到达时间点，一个周期的序列长度为na，记为s3.2、生成数据压缩比α，令m＝αna，以此为参量，随机生成测量矩阵s3.3、基于最佳检测的运算机制核心机制，生成稀疏基ψ，基于最佳检测的运算机制，采用一路伪随机序列p(n)作为发送的参考码型，另一路为接收的经过延时后的光子到达时间点，重构后得到序列x(n)，两个序列最佳检测的运算如下，其中，na为一个周期的序列长度，m为时间单元，扩频光子计数值分布r(m)的峰值对应横坐标即为时间延时，其代表目标所处的距离信息，其核心机制是接收序列和参考序列的各次位移相乘再求和的过程。基于最佳检测的运算方法，以参考序列为模板序列，拟采用包含时间位移的hankel矩阵结构，构造参考序列结构化的稀疏基，参考序列结构化的稀疏基的列元素包含所有可能的参考序列的时间位移，其分以下两种情况建立：1)如果积分时间ta和序列的周期tp相等，即模板序列和接收序列的长度相等，则构建的稀疏基表示为，其中，ψ(:,i)代表第i列向量，zeros代表零向量，[]t代表矩阵的转置，p(i:j)为从第i个元素到第j个元素的参考向量；4)当积分时间ta大于序列的周期tp，向量p右侧补零，直至p和x长度相等；s3.4、计算稀疏基的相干性δ(ψ)采用δ(ψ)表示稀疏基的相干性，最劣的相干性即为两个不同稀疏基元素之间内积的绝对值的最大值，δ(ψ)的取值范围为[0，1]，计算式(11)以防止出现相干的稀疏基元素对，式中，ψi为第i列元素构成的向量，ψj为第j列元素构成的向量；s3.5、重构扩频光子计数值分布设测量信号为为接收的经过延时后的光子到达时间信号，为变换系数，即待恢复的扩频光子计数值分布，其描述光子在时间轴上的分布特征，其具有尖锐的峰值，在时间分布上具有稀疏特性，即||r||0≤k，是k稀疏的，则有在稀疏基的表示下是k稀疏的，基于压缩感知的测量过程理解为对经过延时后的光子到达时间信号的m次独立测量，每次以线性函数的形式测量接收数据的线性投影m个测量向量构成的测量矩阵为其将高维信号x映射到低维，则测量过程表述为y＝φx＝φ(ψr)(12)设v表示噪声，带噪观测情况下，式(12)表示为，y＝ar v(13)其中a＝φψ；s3.5.1式(13)转化为以下lasso标准问题的求解，λ为代价参数，s3.5.2计算rk＝(ata ρi)-1[aty ρ(zk-1-uk-1)]其中rk为第k次迭代计算得到的扩频光子计数值分布，即式(1)中的r(n)，ρ为惩罚系数，i为单位矩阵，()-1为矩阵的逆，()t为矩阵的转置，zk-1和uk-1均为第k-1次迭代的长度为2na-1的向量；s3.5.3计算其中prox为关于变量z的l1范数近端算子，具有显式解；s3.5.4计算uk＝uk-1 rk-zk；s3.6、判断δ(ψ)＜ε是否小于ε，如果小于ε，则k＝k 1，继续采用s3.5的步骤迭代更新r，如果大于设定数值，则停止迭代，输出r。s4、采用基于贝叶斯理论和最大似然估计的迭代重建算法，引入矢量外推加速算法，估计最终携带目标特征的目标光子计数值分布，具体为，s4.1、初始化拟合的冲激响应﹑最大迭代次数kmax和加速参数α；s4.2、设携带目标特征的目标光子计数值分布为f，简称为目标光子计数值，由贝叶斯模型可得，式(14)中，p(f/r)为已知扩频光子计数值分布情况下，目标光子计数值的概率分布情况，p(r/f)为已知目标光子计数值的情况下，扩频光子计数值概率分布情况，p(f)为目标光子计数概率分布；当p(f/r)取得最大数值，则恢复的目标信号接近理想的波形，式(14)的分母部分为常量，通过最大似然估计法求解其最大值，ml(f)＝max[p(r/f)](15)式(15)中，ml(f)表示采用最大似然估计的方法估计f，max为最大值，设待恢复的波形，即目标光子计数值f，其与系统的冲激响应h的卷积表示为t表示每个像素的积分时间，则联合概率密度分布l(r/f)为，式(16)等号两边求对数得到，对等式(17)两边求导，并令其为0，根据退化函数归一化的性质，可得，引入迭代运算，则有，式(19)是迭代计算公式的一般形式，其中，ψ(fk)表示线性搜索的迭代算法，对于一维扩频光子计数值分布，fk表示第k次迭代计算时产生的目标光子计数值的估计值，fk 1表示第k 1次迭代计算时产生的目标光子计数值的估计值；s4.3、基于差分表示的估计值变化引入根据式(19)从现有的估计值fk-1产生估计值fk的迭代出发，给出一阶泰勒矢量，fk＝ψ(fk-1)＝fk-1 gk(20)设fk逐渐收敛到第s次迭代计算时产生的估计值fs，在第k次迭代中引入的变化gk朝着解的梯度方向▽m(f)生成一条通向多维空间的解，目标光子计数值的求解过程中，需要最大化的函数可局部近似为多维二次函数m(f)，则有，m(f)＝-(f-fs)ta(f-fs)(21)▽m(f)＝-2a(f-fs)(22)式(22)中，()t表示转置，a表示平方、对称、正定矩阵，根据式(21)和(22)，采用差分形式给出第k次迭代引入的变化gk，如式(23)所示，gk＝ε▽m(fk)＝-2εa(fk-fs)(23)式(19)中，根据矩阵对角化原理，通过矩阵a的特征值构成的对角矩阵λa和特征向量s给出矩阵a，即，a＝sλas-1(24)设i为单位矩阵，λb是矩阵b的特征值构成的对角矩阵，b＝i-2εa,λb＝i-2ελa，f0表示估计值的初值，λbk为λb的k次方，ε表示一个用来防止λb出现负值的向量，根据式(20)、式(23)和式(24)，fk为，根据式(25)得到第k-1次迭代引入的变化gk-1，k-δ-1次迭代引入的变化gk-1-δ，k-δ次迭代的估计值fk-δ，gk-1＝sλbk(i-λb-1)s-1(f0-fs) fs(26)fk-δ＝sλbk-δs-1(f0-fs) fs(27)gk-δ-1＝sλbk-δ(i-λb-1)s-1(f0-fs)(28)式(28)中，λbk-δ为λb的k-δ次方，且有i-λb-δ≈δ(i-λb-1)；s4.4、计算估计值之差估计值fk和估计值fk-δ之差hk,-δ表示为，hk,-δ＝fk-fk-δ＝sλbk(i-λb-δ)s-1(f0-fs)(29)对于第k δ次迭代的估计值fk δ的预测为，fk δ＝sλbk δs-1(f0-fs) fs(30)第k δ-1次迭代引入的变化gk δ-1为，gk δ-1＝sλbk δ(i-λb-1)s-1(f0-fs)(31)估计值fk δ和估计值fk之差hk,δ为，hk,δ＝fk δ-fk＝sλbk δ(i-λb-δ)s-1(f0-fs)(32)式(32)中，λbk δ为λb的k δ次方，λb-δ为λb的-δ次方，λb-1为λb的-1次方，s-1为s的-1次方；s4.5、估计加速参数基于泰勒级数的一阶矢量外推法，通过最小二乘投影αkhk,-δ预测hk,δ，得，hk,δ＝αkhk,-δ(33)(hk,δ-αkhk,-δ)thk,-δ＝0(34)式(34)中，αk为加速参数，求解式(34)，得到，又因为i-λb-δ≈δ(i-λb-1),将i-λb-δ≈δ(i-λb-1)代入式(32)可得，hk,δ≈δsλbk(i-λb-1)s-1(f0-fs)＝δgk-1(36)以此类推得到，hk,-δ≈δgk-δ-1(37)将式(35)和式(37)代入式(35)，可得加速参数αk为，s4.6、估计携带目标特征的目标光子计数值设第k次待估计的目标光子计数值为lk，根据第k次估计中已经得到的fk、fk-1、αk，代入式lk＝fk αk(fk-fk-1)，得到的lk，k 1后，将lk代入得到fk 1，同理通过式(20)得到gk 1，通过式(38)得到αk 1，再代入lk＝fk αk(fk-fk-1)的变换式lk 1＝fk 1 αk 1(fk 1-fk)，得到lk 1，以此类推，对lk进行循环更新，直到达到最大迭代次数kmax，获得最终携带目标特征的目标光子计数值lkmax。s5、根据最终携带目标特征的目标光子计数值分布，通过残差去除﹑峰值搜索和质心拟合方法提取距离数值，具体为，s5.1、通过最终携带目标特征的目标光子计数值lkmax，最高幅度值选取阈值β，提取出满足lkmax＜β对应的所有k值，对应k的lkmax幅度值置为0，得到的数据记为lkmax_opt，即第k个深度单元对应的目标光子计数值；s5.2、搜索lkmax_opt中的峰值，并且提取峰值对应的k值；s5.3、在搜索到的峰值附近提取l个非零数值点，采用质心拟合算法计算距离数值d，式(39)中，τ为最小的距离分辨单元。s6、逐点重构距离数值，得到重构的二维深度图像，具体为，将准直器固定在二维导轨平台上，扫描获得每个像素点的光子到达时间点，采用步骤s3～步骤s6所述的方法，如图3所示，逐点重构距离数值，得到重构的二维深度图像。为了描述单点测距的性能，定义m次测量的单点距离均方误差为：其中，为第j次估计的距离值，dactual为实际距离数值，tc为码元宽度，c为光速。定义m次测量的单点距离标准差为：其中，为估计值的均值。为了初始化最大迭代次数，需要定义纵向距离间隔均方误差(rmse)为：为算法估计的两个目标的纵向距离间隔，dactualseparation为两个目标实际的纵向距离间隔。以估计的间隔和实际间隔的均方误差作为评价算法性能的指标，并确定距离间隔均方误差最低的迭代次数作为最大迭代次数。设为算法估计的两个目标的平均纵向距离间隔，则纵向距离间隔标准差(std)表示为:实验中，首先，对7米处的白墙成像，每个像素点获取的信号光子计数值为5400个，噪声信号光子计数值为600个，总共为6000个。分别采用50％、30％、10％三种压缩比，则每个像素点获取的光子计数值分别为3000个、1800个和600个，采用图3中带噪声高效重构方法得到扩频光子计数分布。测量50次，得到该分布的峰值对应距离数值的rmsesingle和stdsingle如表1所示。同等条件下，采用最佳检测的方法重构的扩频光子计数分布的峰值对应距离数值的rmsesingle和stdsingle分别为3.1cm和2.7cm。也即采用50％的光子计数值即可达到最佳检测的方法的距离精确度，如表1所示。其次，对放置在白墙前面的积木成一维距离像，积木和白墙间隔8厘米，成像目标如图1(a)所示。将激光打在积木边缘，根据图3高效率高分辨力成像方法流程，采用步骤3-6，先通过带噪声高效重构方法得到扩频光子计数值分布，以最小的rmse得到最大迭代次数kmax为5次，如图5所示，根据步骤s2计算得到拟合的冲激响应，再采用快速迭代重构目标光子计数值方法得到高分辨力距离像。表2给出8厘米间隔，50次测量的间隔均方误差。表3为在同等条件下，采用传统的最佳检测获得扩频光子计数值分布，再采用ancombe变换和快速反卷积算法提高分辨能力。比较表2和3，本发明提出的方法只需要获取60％的光子计数值，即可以达到原方法估计的距离间隔均方误差，省去了单独的去噪声算法，降低了数据存取复杂度，提高运算速度和运算效率。表1带噪高效重建方法距离估计值统计量压缩比rmsesingle[cm]stdsingle[cm]50％32.130％98.610％14.412.7表2高效率高分辨力间隔估计值统计量压缩比rmse[cm]std[cm]60％3.92.950％8.56.330％11.510.910％20.818.9表3传统最佳检测高分辨力间隔估计值统计量rmse[cm]std[cm]3.82.7以上显示和描述了本发明的基本原理、主要特征及优点。本行业的技术人员应该了解，本发明不受上述实施例的限制，上述实施例和说明书中描述的只是说明本发明的原理，在不脱离本发明精神和范围的前提下，本发明还会有各种变化和改进，这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。当前第1页1 2 3
技术特征：1.一种光子计数目标高效率高分辨力深度重建方法，基于一体化光纤式伪随机码幅度调制深度矫正装置的成像系统，其特征在于，包括以下步骤：
s1、将成像系统的准直器固定在二维导轨平台上，去掉a单光子探测器、固定衰减器和分路器，伪随机序列发生器产生一路同步高电平脉冲，该脉冲产生时间和开始发送伪随机序列的时间一致，将同步高电平脉冲连接光子时间到达记录仪，作为开始计时的触发信号，伪随机序列发生器驱动激光器发送激光脉冲，激光脉冲信号通过成像系统后得到单点扩频光子计数值分布，表达为系统冲激响应和携带目标特征的目标光子计数值分布的卷积积分；
s2、改变s1中的成像系统结构，两个准直器面对面发射并接收光子，获取平均的系统冲激响应并拟合；
s3、恢复成像系统结构到初始状态，基于最佳检测的运算机制，通过构建合适的稀疏基，高效反演目标的扩频光子计数值分布；
s4、采用基于贝叶斯理论和最大似然估计的迭代重建算法，引入矢量外推加速算法，估计最终携带目标特征的目标光子计数值分布；
s5、根据最终携带目标特征的目标光子计数值分布，通过残差去除﹑峰值搜索和质心拟合方法提取距离数值；
s6、逐点重构距离数值，得到重构的二维深度图像。
2.根据权利要求1所述的一种光子计数目标高效率高分辨力深度重建方法，其特征在于：所述步骤s1将准直器固定在二维导轨平台上，去掉a单光子探测器、固定衰减器和分路器，伪随机序列发生器产生一路同步高电平脉冲，该脉冲产生时间和开始发送伪随机序列的时间一致，将同步高电平脉冲连接光子时间到达记录仪，作为开始计时的触发信号，伪随机序列发生器驱动激光器发送激光脉冲，激光脉冲信号通过成像系统后得到的单点扩频光子计数值分布，表达为系统冲激响应和携带目标特征的目标光子计数值分布的卷积积分，具体为，
将准直器固定在二维导轨平台上，去掉a单光子探测器、固定衰减器和分路器，伪随机序列发生器产生一路同步高电平脉冲，该脉冲产生时间和开始发送伪随机序列的时间一致，将同步高电平脉冲连接光子时间到达记录仪，作为开始计时的触发信号，伪随机序列发生器驱动激光器发送激光脉冲，激光脉冲发射后，先后通过大气介质、目标和接收系统，在短距离下在到达接收系统前，脉冲回波信号s(n)由目标冲激响应决定；
通过接收系统后，目标冲激响应ho(n)和探测器的冲激响应hd(n)是造成接收的目标信号产生瞬时失真的主要因素，则伪随机光脉冲碰到目标后，往返一次，回波到达接收系统，探测到的单点扩频光子计数值分布r(n)为，
式中，r(n)为携带目标信息并且被系统冲激响应模糊的单点扩频光子计数值分布，c(n)为互相关函数，为系统冲激响应，表征系统自身产生的瞬时干扰，f(n)＝ho(n)为携带目标特征的目标光子计数值分布，b为噪声光子。
3.根据权利要求2所述的一种光子计数目标高效率高分辨力深度重建方法，其特征在于：所述步骤s2改变s1中的成像系统结构，两个准直器面对面发射并接收光子，获取平均的系统冲激响应并拟合，具体为，
s2.1、在成像系统上，断开可调衰减器和光环行器的连接，将准直器b和可调衰减器连接，作为收准直器端，将连接光环行器的准直器a作为发准直器端，将准直器a和准直器b面对面放置，通过微调和粗调准直器的位置，使得从准直器a中发出的光信号导入接收的准直器b中，再由b单光子探测器接收后，得到接收路光子时间到达点；
s2.2、将matlab的序列重构方法嵌入labview，基于labview，获取接收路的光子时间到达点，并重构0﹑1序列，表示为y(n)，伪随机序列发生器中存储的模板序列p(n)与y(n)进行快速傅里叶变换，得到扩频光子计数值分布，即系统冲激响应h，其代表系统的瞬时干扰，
式中，表示傅里叶变换，为傅里叶逆变换；
s2.3、测量若干次的系统冲激响应并作平均，得到平均的系统冲激响应；
s2.4、采用高斯拟合法拟合平均的系统冲激响应，
设平均的系统冲激响应数据(oi,hi)i＝1,2,3,…,n，用高斯函数描述为，
式中，hmax、omax和s为待估参数，hmax代表高斯曲线的峰高，omax代表高斯曲线峰位置，s代表半高全宽；
将(3)式两边取自然对数，化为，
令
则(4)式化为二次多项式拟合函数，
考虑全部数据和测量误差εi，并以矩阵形式表示如下，
简记为zn×1＝on×3b3×1 en×1，
根据最小二乘原理，求得拟合常数b0﹑b1﹑b2构成的矩阵b的广义最小二乘解为，
b＝(oto)-1otz(8)
进而根据(5)式，求出待估参数hmax﹑omax和s，代入(3)式，从而得到拟合的高斯函数。
4.根据权利要求3所述的一种光子计数目标高效率高分辨力深度重建方法，其特征在于：所述步骤s3恢复成像系统结构并成像，基于最佳检测的运算机制，通过稀疏基表示方法，高效反演目标的扩频光子计数值分布，具体为，
恢复成像系统结构，对目标成像，获得接收路光子到达点时间，
s3.1、根据步骤s2.2的序列重构方法，重构接收路光子到达时间点，一个周期的序列长度为na，接收的经过延时后的光子到达时间信号记为
s3.2、生成数据压缩比α，令m＝αna，以此为参量，随机生成测量矩阵
s3.3、基于最佳检测的运算机制，生成稀疏基ψ，
基于最佳检测的运算机制，采用一路伪随机序列p(n)作为发送的参考码型，另一路为接收的经过延时后的光子到达时间点，重构后得到序列x(n)，两个序列最佳检测的运算表示为，
其中，na为一个周期的序列长度，m为时间单元，
扩频光子计数值分布r(m)的峰值对应横坐标即为时间延时，其代表目标所处的距离信息；基于最佳检测的运算方法，以参考序列为模板序列，拟采用包含时间位移的hankel矩阵结构，构造稀疏基ψ，该稀疏基的列元素包含所有可能的参考序列的时间位移，分以下两种情况建立，
1)当积分的时间ta和序列的周期tp相等，即模板序列和接收序列的长度相等，则构建的稀疏基表示为，
其中，ψ(:,i)代表稀疏基的第i列向量，zeros代表零向量，[]t代表矩阵的转置，p(i:j)为从第i个元素到第j个元素的参考向量；
2)当积分时间ta大于序列的周期tp，向量p右侧补零，直至p和x长度相等；
s3.4、计算稀疏基的相干性δ(ψ)
采用δ(ψ)表示稀疏基的相干性，最劣的相干性即为两个不同稀疏基元素之间内积的绝对值的最大值，δ(ψ)的取值范围为[0，1]，计算式(11)以防止出现相干的稀疏基元素对，
式中，ψi为第i列元素构成的向量，ψj为第j列元素构成的向量；
s3.5、重构扩频光子计数值分布
设测量信号为为接收的经过延时后的光子到达时间信号，为变换系数，即待恢复的扩频光子计数值分布，其描述光子在时间轴上的分布特征，其具有尖锐的峰值，在时间分布上具有稀疏特性，即||r||0≤k，是k稀疏的，则有在稀疏基的表示下是k稀疏的，基于压缩感知的测量过程理解为对经过延时后的光子到达时间信号的m次独立测量，每次以线性函数的形式测量接收数据的线性投影m个测量向量构成的测量矩阵为其将高维信号x映射到低维，则测量过程表述为
y＝φx＝φ(ψr)(12)
设v表示噪声，带噪观测情况下，式(12)表示为，
y＝ar v(13)
其中a＝φψ；
s3.5.1式(13)转化为以下lasso标准问题的求解，λ为代价参数，
s3.5.2计算rk＝(ata ρi)-1[aty ρ(zk-1-uk-1)]，
其中rk为第k次迭代计算得到的扩频光子计数值分布，即式(1)中的r(n)，ρ为惩罚系数，i为单位矩阵，()-1为矩阵的逆，()t为矩阵的转置，zk-1和uk-1均为第k-1次迭代的长度为2na-1的向量；
s3.5.3计算
其中prox为关于变量z的l1范数近端算子，具有显式解；
s3.5.4计算uk＝uk-1 rk-zk；
s3.6、判断δ(ψ)＜ε是否成立，如果小于设定数值ε，则k＝k 1，继续采用s3.5的步骤迭代更新r，如果大于设定数值ε，则停止迭代，输出r。
5.根据权利要求4所述的一种光子计数目标高效率高分辨力深度重建方法，其特征在于：所述步骤s4采用基于贝叶斯理论和最大似然估计的迭代重建算法，引入矢量外推加速算法，估计最终携带目标特征的目标光子计数值分布，具体为，
s4.1、初始化拟合的冲激响应﹑最大迭代次数kmax和加速参数α；
s4.2、设携带目标特征的目标光子计数值分布为f，简称为目标光子计数值，由贝叶斯模型可得，
式(14)中，p(f/r)为已知扩频光子计数值分布情况下，目标光子计数值的概率分布情况，p(r/f)为已知目标光子计数值的情况下，扩频光子计数值概率分布情况，p(f)为目标光子计数概率分布；当p(f/r)取得最大数值，则恢复的目标信号接近理想的波形，
式(14)的分母部分为常量，通过最大似然估计法求解其最大值，
ml(f)＝max[p(r/f)](15)
式(15)中，ml(f)表示采用最大似然估计的方法估计f，max为最大值，
设待恢复的波形，即目标光子计数值f，其与系统的冲激响应h的卷积表示为t表示每个像素的积分时间，则联合概率密度分布l(r/f)为，
式(16)等号两边求对数得到，
对等式(17)两边求导，并令其为0，根据退化函数归一化的性质，可得，
引入迭代运算，则有，
式(19)是迭代计算公式的一般形式，其中，ψ(fk)表示线性搜索的迭代算法，对于一维扩频光子计数值分布，fk表示第k次迭代计算时产生的目标光子计数值的估计值，fk 1表示第k 1次迭代计算时产生的目标光子计数值的估计值；
s4.3、基于差分表示的估计值变化引入
根据式(19)从现有的估计值fk-1产生估计值fk的迭代出发，给出一阶泰勒矢量，
fk＝ψ(fk-1)＝fk-1 gk(20)
设fk逐渐收敛到第s次迭代计算时产生的估计值fs，在第k次迭代中引入的变化gk朝着解的梯度方向生成一条通向多维空间的解，目标光子计数值的求解过程中，需要最大化的函数可局部近似为多维二次函数m(f)，则有，
m(f)＝-(f-fs)ta(f-fs)(21)
式(22)中，()t表示转置，a表示平方、对称、正定矩阵，
根据式(21)和(22)，采用差分形式给出第k次迭代引入的变化gk，如式(23)所示，
式(19)中，根据矩阵对角化原理，通过矩阵a的特征值构成的对角矩阵λa和特征向量s给出矩阵a，即，
a＝sλas-1(24)
设i为单位矩阵，λb是矩阵b的特征值构成的对角矩阵，b＝i-2εa,λb＝i-2ελa，f0表示估计值的初值，λbk为λb的k次方，ε表示一个用来防止λb出现负值的向量，根据式(20)、式(23)和式(24)，fk为，
根据式(25)得到第k-1次迭代引入的变化gk-1，k-δ-1次迭代引入的变化gk-1-δ，k-δ次迭代的估计值fk-δ，
gk-1＝sλbk(i-λb-1)s-1(f0-fs) fs(26)
fk-δ＝sλbk-δs-1(f0-fs) fs(27)
gk-δ-1＝sλbk-δ(i-λb-1)s-1(f0-fs)(28)
式(28)中，λbk-δ为λb的k-δ次方，且有i-λb-δ≈δ(i-λb-1)；
s4.4、计算估计值之差
估计值fk和估计值fk-δ之差hk,-δ表示为，
hk,-δ＝fk-fk-δ＝sλbk(i-λb-δ)s-1(f0-fs)(29)
对于第k δ次迭代的估计值fk δ的预测为，
fk δ＝sλbk δs-1(f0-fs) fs(30)
第k δ-1次迭代引入的变化gk δ-1为，
gk δ-1＝sλbk δ(i-λb-1)s-1(f0-fs)(31)
估计值fk δ和估计值fk之差hk,δ为，
hk,δ＝fk δ-fk＝sλbk δ(i-λb-δ)s-1(f0-fs)(32)
式(32)中，λbk δ为λb的k δ次方，λb-δ为λb的-δ次方，λb-1为λb的-1次方，s-1为s的-1次方；
s4.5、估计加速参数
基于泰勒级数的一阶矢量外推法，通过最小二乘投影αkhk,-δ预测hk,δ，得，
hk,δ＝αkhk,-δ(33)
(hk,δ-αkhk,-δ)thk,-δ＝0(34)
式(34)中，αk为加速参数，
求解式(34)，得到，
又因为i-λb-δ≈δ(i-λb-1),将i-λb-δ≈δ(i-λb-1)代入式(32)可得，
hk,δ≈δsλbk(i-λb-1)s-1(f0-fs)＝δgk-1(36)
以此类推得到，
hk,-δ≈δgk-δ-1(37)
将式(35)和式(37)代入式(35)，可得加速参数αk为，
s4.6、估计携带目标特征的目标光子计数值
设第k次待估计的目标光子计数值为lk，根据第k次估计中已经得到的fk、fk-1、αk，代入式lk＝fk αk(fk-fk-1)，得到的lk，k 1后，将lk代入得到fk 1，同理通过式(20)得到gk 1，通过式(38)得到αk 1，再代入lk＝fk αk(fk-fk-1)的变换式lk 1＝fk 1 αk 1(fk 1-fk)，得到lk 1，以此类推，对lk进行循环更新，直到达到最大迭代次数kmax，获得最终携带目标特征的目标光子计数值lkmax。
6.根据权利要求5所述的一种光子计数目标高效率高分辨力深度重建方法，其特征在于：所述步骤s5根据最终携带目标特征的目标光子计数值分布，通过残差去除﹑峰值搜索和质心拟合方法提取距离数值，具体为，
s5.1、通过最终携带目标特征的目标光子计数值lkmax，最高幅度值选取阈值β，提取出满足lkmax＜β对应的所有k值，对应k的lkmax幅度值置为0，得到的数据记为lkmax_opt，即第k个深度单元对应的目标光子计数值；
s5.2、搜索lkmax_opt中的峰值，并且提取峰值对应的k值；
s5.3、在搜索到的峰值附近提取l个非零数值点，采用质心拟合算法计算距离数值d，
式(39)中，τ为最小的距离分辨单元。
7.根据权利要求1所述的一种光子计数目标高效率高分辨力深度重建方法，其特征在于：所述步骤s6逐点重构距离数值，得到重构的二维深度图像，具体为，
将准直器固定在二维导轨平台上，扫描获得每个像素点的光子到达时间点，采用步骤s3～步骤s5所述的方法，逐点重构距离数值，得到重构的二维深度图像。
技术总结本发明公开了一种光子计数目标高效率高分辨力深度重建方法，包括以下步骤，S1、将成像系统的准直器固定在二维导轨平台上，去掉A单光子探测器、固定衰减器和分路器，伪随机序列发生器产生一路同步高电平脉冲，该脉冲产生时间和开始发送伪随机序列的时间一致，将同步高电平脉冲连接光子时间到达记录仪，作为开始计时的触发信号，伪随机序列发生器驱动激光器发送激光脉冲，激光脉冲信号通过成像系统后得到单点扩频光子计数值分布，表达为系统的冲激响应和携带目标特征的目标光子计数值分布的卷积积分。省去单独的噪声去除步骤，以最快的速度方向，加速估计下一次迭代的数值，运算速度快于线性搜索的迭代算法。
技术研发人员：沈姗姗;孙肖琳;苏适;何瑞清;冒添翼
受保护的技术使用者：南京工业职业技术大学
技术研发日：2021.04.30
技术公布日：2021.08.03
