[发明专利]一种快速分组光子强度关联方法

权利要求书

1.一种快速分组光子强度关联方法，其特征在于，包括如下步骤:

步骤a：记录同一观测源辐射出的光子序列分别到达空间中两个探测器的时间，得到两组光子到达时间序列信息；将得到的两组光子到达时间序列组成一组数据，然后按照时间顺序排序；

步骤b：将排序后的数据[a₁₁,a₁₂，…，a_1U,a₂₁,a₂₂，…，a_2U]分割成2U个相对应的数据单元，其中a_1i与a_2i(i＝1,2...U)互相对应；将时间序列数据[a₁₁,a₁₂，…，a_1U]标记为第一通道数据单元，时间序列数据[a₂₁,a₂₂，…，a_2U]标记为第二通道数据单元；

步骤c：确定光程差τ_g的搜索范围[t_a,t_b]及搜索步长t_min，设定初始搜索值t_j的大小为t_a，令j＝1，进行迭代操作；

步骤d：将第二通道数据单元[a₂₁,a₂₂，…，a_2U]中的每个数据都加上t_j，当达到迭代终止条件t_j＝t_b时跳转至步骤i，否则进入步骤e；

步骤e：将数据[a₁₁,a₂₁],[a₁₂,a₂₂]，…，[a_1U,a_2U]输入强度关联方法队列进行并行计算，得到计算结果G_e；

步骤f：在第一通道数据单元和第二通道数据单元交界处交叉对应提取部分数据，组成新数据单元A；

步骤g：将得到的新数据单元A输入强度关联方法队列进行并行计算，得到计算结果G_g；

步骤h：令G_i＝G_e+G_g，t_j+1＝t_j+t_min，j＝j+1,跳转至步骤d；

步骤i：输出强度关联曲线，纵坐标为G_i，横坐标为t_i，曲线峰值对应的横坐标即为光程差τ_g的最终迭代值。

2.根据权利要求1所述的快速分组光子强度关联方法，其特征在于，步骤c中根据观测源距离第一个探测器和第二个探测器的距离差来对所述光程差τ_g的搜索范围[t_a,t_b]进行初始化。

3.根据权利要求1所述的快速分组光子强度关联方法，其特征在于，步骤f中通过将每个通道中一段时间的累计探测到的光子数据分成U个数据块，每个数据块进行独立计算，再对相邻数据块中漏掉的关联光子进行补充，最后获得整体的关联光子数；补充计算考虑交叉的情况时a_1,i的尾部数据与a_2,i+1的头部数据关联，并且a_1,i+1的头部数据与a_2,i的尾部数据关联,补充计算选取头部和尾部数据块的大小为最后2倍光子符合门限时长的数据。

4.根据权利要求3所述的快速分组光子强度关联方法，其特征在于，当n个探测器进行两两强度关联时，则每两个通道之间进行独立关联，获得个关联数据值，从而通过分布式并行计算获得多通道的光子关联数据。

5.根据权利要求1所述的快速分组光子强度关联方法，其特征在于，步骤e和步骤g中强度关联方法的内容为：

1)创建两个数组a[chal]和b[chbl]，数组a表示第一个探测器探测到的chal个光子的到达时间序列，数组b表示第二个探测器探测到的chbl个光子的到达时间序列；

2)令变量i_a初始值为0,j_b初始值为0，la初始值为chal，lb初始值为chbl，ans初始值为0；

3)如果a[i_a]大于b[j_b]，令j_b＝j_b+1，如果a[i_a]小于b[j_b]，令i_a＝i_a+1；两种情况都不满足时令i_a＝i_a+1，j_b＝j_b+1且ans值自增1；

4)判断是否i_a＜la且j_b＜lb，如果结果为是，则进入步骤3)，如果结果为否则输出ans值，结束运算。