1.一种重建探测器所探测的能谱的方法,包括步骤:用射线照射被检查物体;利用探测器接收穿透所述被检查物体的射线,将接收到的射线转换为表示探测到的能谱的能谱数据;以及利用事先创建的探测器响应模型通过统计迭代算法基于所述能谱数据重建所述探测器的能谱,其中通过在多个能量段中采集不同靶材料产生的X射线荧光能谱对所述探测器响应模型的参数进行标定,得到所述探测器响应模型的显式表示,其中所述探测器响应模型是基于单能峰、探测器的材料的特征峰、逃逸峰以及基线模型而生成的。
2.如权利要求1所述的方法,其中所述单能峰、探测器的材料的特征峰以及逃逸峰中的每一个基于高斯峰来表示,所述高斯峰表示为:其中,E D 为实际探测到的能谱能量,E是峰位的能量,σ(E)是标准差。
3.如权利要求2所述的方法,其中所述单能峰、探测器的材料的特征峰以及逃逸峰的分布表示为:R peak (E D ,E)=G(E D ,E)+c 11 R DXRF (E D )+c 12 R escape (E D ,E)其中,R peak (E D ,E)是分布表示,G(E D ,E)是所述单能峰,R DXRF (E D )是特征峰以及R escape (E D ,E)是逃逸峰,所述特征峰和所述逃逸峰分别由具有两个或多个不同参数的高斯峰相加得到,所述高斯峰的数量由探测器的材料的元素数量决定,其中,c 11 ,c 12 是常数参数。
4.如权利要求3所述的方法,其中所述探测器是碲化镉探测器。
5.如权利要求4所述的方法,其中所述分布表示中的特征峰和逃逸峰分别由以下等式获得:其中,E Cd 和E Te 分别为碲化镉晶体中碲、镉两种元素对应的荧光光子能量,c 6 为常数参数。
6.如权利要求1所述的方法,其中所述基线模型部分表示为:其中,w(E D ,E)是权重参数,使用二次函数进行拟合,w(E D ,E)表示为:w(E D ,E)=c 3 E D 2 +c 4 E D +(c 5 E+1)其中,c 3 ,c 4 ,c 5 是常数参数,E D 为实际探测到的能谱能量,E是峰位的能量,σ(E)是标准差。
7.如权利要求1所述的方法,还包括:在建立了所述探测器响应模型后,通过采集探测器的多种材料的单能荧光光子信号进行常数参数的标定以得到最终的探测器响应函数。
8.如权利要求1所述的方法,其中基于下式进行统计迭代重建:其中D i 表示探测到的能谱数据,R ij 表示所述探测器响应模型的矩阵表示,S j 表示重建的探测器能谱,n为迭代次数。
9.一种重建探测器所探测的能谱的设备,包括:射线源,产生射线,照射被检查物体;探测器,接收穿透所述被检查物体的射线,并且将接收到的射线转换为表示探测到的能谱的能谱数据;以及数据处理设备,利用事先创建的探测器响应模型通过统计迭代算法基于所述能谱数据重建所述探测器的能谱,其中通过在多个能量段中采集不同靶材料产生的X射线荧光能谱对所述探测器响应模型的参数进行标定,得到所述探测器响应模型的显式表示,其中所述探测器响应模型是基于单能峰、探测器的材料的特征峰、逃逸峰以及基线模型而生成的。
10.如权利要求9所述的设备,其中所述单能峰、探测器的材料的特征峰以及逃逸峰中的每一个基于高斯峰来表示,所述高斯峰表示为:其中,E D 为实际探测到的能谱能量,E是峰位的能量,σ(E)是标准差。
11.如权利要求10所述的设备,其中所述单能峰、探测器的材料的特征峰以及逃逸峰的分布表示为:R peak (E D ,E)=G(E D ,E)+c 11 R DXRF (E D )+c 12 R escape (E D ,E)其中,R peak (E D ,E)是分布表示,G(E D ,E)是所述单能峰,R DXRF (E D )是特征峰以及R escape (E D ,E)是逃逸峰,所述特征峰和所述逃逸峰分别由具有两个或多个不同参数的高斯峰相加得到,所述高斯峰的数量由探测器的材料的元素数量决定,其中,c 11 ,c 12 是常数参数。
12.如权利要求11所述的设备,其中所述探测器是碲化镉探测器。
13.如权利要求12所述的设备,其中所述分布表示中的特征峰和逃逸峰分别由以下等式获得:其中,E Cd 和E Te 分别为碲化镉晶体中碲、镉两种元素对应的荧光光子能量,c 6 为常数参数。
14.如权利要求9所述的设备,其中所述基线模型部分表示为:其中,w(E D ,E)是权重参数,使用二次函数进行拟合,w(E D ,E)表示为:w(E D ,E)=c 3 E D 2 +c 4 E D +(c 5 E+1)其中,c 3 ,c 4 ,c 5 是常数参数,E D 为实际探测到的能谱能量,E是峰位的能量,σ(E)是标准差。
15.如权利要求9所述的设备,其中所述数据处理设备基于下式进行统计迭代重建:其中D i 表示探测到的能谱数据,R ij 表示所述探测器响应模型的矩阵表示,S j 表示重建的探测器能谱,n为迭代次数。