1.一种放射性物质辐射剂量率的二维角度分布的测量方法，其特征在于，包括：
伽马相机探测目标区域内从放射性物质定义的目标角平面(α,β)各个方向入射的伽马光子，并根据所述入射的伽马光子及其相应的探测作用位置重组生成矢量投影数据p，以及根据所述入射的伽马光子在所述伽马相机的探测器中沉积能量重组生成矢量能谱数据E；
根据所述矢量能谱数据E获取入射光子能量Epeak，并根据所述矢量投影数据p和所述入射光子能量Epeak重建得到伽马辐射图像；
根据所述伽马辐射图像获取实际入射光子计数率的二维角度分布；
根据所述实际入射光子计数率的二维角度分布获取所述目标角平面(α,β)各个方向伽马光子的注量率的二维角度分布φ(Epeak)(α,β)；
根据所述入射光子能量Epeak查找相应的转换系数κ(Epeak)，并根据所述转换系数κ(Epeak)和所述伽马光子的注量率的二维角度分布φ(Epeak)(α,β)通过以下公式获取目标区域内所述放射性物质对伽马相机位置辐射剂量率的二维角度分布D(α,β)：
D(α,β)＝κ(Epeak)×φ(Epeak)(α,β)；
其中，所述根据伽马辐射图像获取实际入射光子计数率的二维角度分布具体包括：
根据伽马辐射图像 获取所述目标角平面(α,β)各个方向实际入射光子数量的二维角度分布N(α,β)，并根据所述实际入射光子数量的二维角度分布N(α,β)获取所述实际入射光子计数率的二维角度分布N(α,β)/T，其中，T为探测从所述目标角平面(α,β)各个方向入射的伽马光子的探测时间；
所述根据所述实际入射光子计数率的二维角度分布获取所述目标角平面各个方向伽马光子的注量率的二维角度分布具体包括：
通过以下公式获取所述伽马光子的注量率的二维角度分布φ(Epeak)(α,β)：
其中，S为所述探测器的有效探测面积。
2.如权利要求1所述的方法，其特征在于，所述根据所述矢量能谱数据E获取入射光子能量Epeak具体包括：
根据所述矢量能谱数据E获取每个所述伽马光子的光子能量；
通过以下公式获取伽马光子平均能量Eave：
其中，Ei为所述探测到的伽马光子中第i个伽马光子的光子能量；
根据所述伽马光子平均能量Eave获取所述入射光子能量Epeak，其中，
Epeak＝Eave×ae+be，
其中，ae和be为拟合系数。
3.如权利要求1所述的方法，其特征在于，所述根据所述矢量投影数据和所述入射光子能量Epeak重建得到伽马辐射图像具体包括：
根据所述入射光子能量Epeak查找相应的传输矩阵M(α,β)以及相应的从所述目标角平面各个方向入射的伽马光子的探测效率的二维角度分布ζ(Epeak)(α,β)；
根据所述传输矩阵M(α,β)和所述探测效率的二维角度分布ζ(Epeak)(α,β)获取考虑实际传输效率的传输矩阵M；
根据所述考虑实际传输效率的传输矩阵M对所述放射性物质的矢量投影数据p进行重建以得到所述伽马辐射图像 其中，
其中，i＝1,2,…,I，j＝1,2,…,J，pj为所述探测器的第j个像素在所述目标角平面(α,β)各个方向探测到的伽马光子数量组成的向量，M的第j行、i列的元素，Mji表示从所述目标角平面(α,β)的(αi,βi)方向入射的伽马光子在所述伽马相机的所述探测器的第j个像素上被探测到的概率；
通过以下期望最大化的统计迭代算法对所述伽马辐射图像 求解：
其中， 为第n次迭代得到的伽马辐射图像中第i个像素的像素估计值， i＝1,
2,…,I。
4.如权利要求3所述的方法，其特征在于，通过以下步骤预先建立入射光子能量Epeak与伽马光子的探测效率的二维角度分布ζ(Epeak)(α,β)的对应关系：
获取所述目标角平面(α,β)各个方向入射到达所述探测器前表面的能量为Epeak的伽马光子被探测到的二维角度分布的第一概率εdet(Epeak)(α,β)；
获取所述目标角平面(α,β)各个方向入射的伽马光子穿透所述探测器的准直器到达所述探测器表面的二维角度分布的第二概率εcol(Epeak)(α,β)；
根据所述第一概率εdet(Epeak)(α,β)和所述第二概率εcol(Epeak)(α,β)获取与所述入射光子能量Epeak对应的伽马光子的探测效率的二维角度分布ζ(Epeak)(α,β)，其中，ζ(Epeak)(α,β)＝εcol(Epeak)(α,β)×εdet(Epeak)(α,β)。
5.如权利要求1所述的方法，其特征在于，所述根据所述矢量投影数据和所述入射光子能量Epeak重建得到伽马辐射图像具体包括：
根据所述入射光子能量Epeak查找相应的传输矩阵M(α,β)以及相应的从所述目标角平面(α,β)的预设角度(αm,βn)入射的伽马光子的探测效率ζ(Epeak)(m,n)；
根据所述传输矩阵M(α,β)和所述探测效率ζ(Epeak)(m,n)获取考虑实际传输效率的传输矩阵M'；
根据所述考虑实际传输效率的传输矩阵M'对所述放射性物质的矢量投影数据p进行重建以得到伽马辐射图像 其中，
其中，i＝1,2,…,I，j＝1,2,…,J，pj为所述探测器的第j个像素在所述目标角平面(α,β)各个方向探测到的伽马光子数量组成的向量，M'的第j行、i列的元素，M'ji表示从所述目标角平面(α,β)的(αi,βi)方向入射的伽马光子在所述伽马相机的所述探测器的第j个像素上被探测到的概率；
通过以下期望最大化的统计迭代算法对所述伽马辐射图像 求解：
其中， 为第n次迭代得到的伽马辐射图像中第i个像素的像素估计值， i＝1,
2,…,I。
6.如权利要求5所述的方法，其特征在于，
所述根据伽马辐射图像获取实际入射光子计数率的二维角度分布具体包括：
根据所述伽马辐射图像 获取所述目标角平面(α,β)各个方向实际入射光子数量的二维角度分布N'(α,β)，并根据所述实际入射光子数量的二维角度分布N'(α,β)获取所述实际入射光子计数率的二维角度分布N'(α,β)/T，其中，T为探测从所述目标角平面(α,β)各个方向入射的伽马光子的探测时间；
所述根据所述实际入射光子计数率的二维角度分布获取所述目标角平面各个方向伽马光子的注量率的二维角度分布具体包括：
通过以下公式获取所述伽马光子的注量率的二维角度分布φ(Epeak)(α,β)：
其中，S为所述探测器的有效探测面积。
7.如权利要求5或6所述的方法，其特征在于，通过以下步骤预先建立入射光子能量Epeak与所述目标角平面(α,β)的预设角度(αm,βn)入射的伽马光子的探测效率ζ(Epeak)(m,n)的对应关系：
获取从所述目标角平面的预设角度(αm,βn)入射到达所述探测器前表面的能量为Epeak的伽马光子被探测到的第一概率εdet(Epeak)(m,n)；
获取从所述目标角平面的预设角度(αm,βn)的入射光子穿透所述探测器的准直器到达所述探测器表面的第二概率εcol(Epeak)(m,n)；
根据所述第一概率εdet(Epeak)(m,n)和所述第二概率εcol(Epeak)(m,n)获取与所述入射光子能量Epeak对应的伽马光子的探测效率ζ(Epeak)(m,n)，其中，
ζ(Epeak)(m,n)＝εcol(Epeak)(m,n)×εdet(Epeak)(m,n)。
8.一种放射性物质辐射剂量率的二维角度分布的测量设备，其特征在于，包括：
伽马相机，用于探测目标区域内从放射性物质定义的目标角平面(α,β)各个方向入射的伽马光子，并根据所述入射的伽马光子及其相应的探测作用位置重组生成矢量投影数据p，以及根据所述入射的伽马光子在所述伽马相机的探测器中沉积能量重组生成矢量能谱数据E；
投影数据分析模块，用于根据所述矢量能谱数据E获取入射光子能量Epeak，并根据所述矢量投影数据p和所述入射光子能量Epeak重建得到伽马辐射图像；
注量率分布获取模块，用于根据所述伽马辐射图像获取实际入射光子计数率的二维角度分布，并根据所述实际入射光子计数率的二维角度分布获取所述目标角平面(α,β)各个方向伽马光子的注量率的二维角度分布φ(Epeak)(α,β)；
辐射剂量率分布获取模块，用于根据所述入射光子能量Epeak查找相应的转换系数κ(Epeak)，并根据所述转换系数κ(Epeak)和所述伽马光子的注量率的二维角度分布φ(Epeak)(α,β)通过以下公式获取目标区域内所述放射性物质对伽马相机位置辐射剂量率的二维角度分布D(α,β)：
D(α,β)＝κ(Epeak)×φ(Epeak)(α,β)；
其中，所述注量率分布获取模块具体用于：
根据伽马辐射图像 获取所述目标角平面(α,β)各个方向实际入射光子数量的二维角度分布N(α,β)，并根据所述实际入射光子数量的二维角度分布N(α,β)获取所述实际入射光子计数率的二维角度分布N(α,β)/T，其中，T为探测从所述目标角平面(α,β)各个方向入射的伽马光子的探测时间；
通过以下公式获取所述伽马光子的注量率的二维角度分布φ(Epeak)(α,β)：
其中，S为所述探测器的有效探测面积。
9.如权利要求8所述的设备，其特征在于，所述投影数据分析模块具体用于：
根据所述矢量能谱数据E获取每个所述伽马光子的光子能量；
通过以下公式获取伽马光子平均能量Eave：
其中，Ei为所述探测到的伽马光子中第i个伽马光子的光子能量；
根据所述伽马光子平均能量Eave获取所述入射光子能量Epeak，其中，
Epeak＝Eave×ae+be，
其中，ae和be为拟合系数。
10.如权利要求8所述的设备，其特征在于，所述投影数据分析模块具体用于：
根据所述入射光子能量Epeak查找相应的传输矩阵M(α,β)以及相应的从所述目标角平面各个方向入射的伽马光子的探测效率的二维角度分布ζ(Epeak)(α,β)；
根据所述传输矩阵M(α,β)和所述探测效率的二维角度分布ζ(Epeak)(α,β)获取考虑实际传输效率的传输矩阵M；
根据所述考虑实际传输效率的传输矩阵M对所述放射性物质的矢量投影数据p进行重建以得到所述伽马辐射图像 其中，
其中，i＝1,2,…,I，j＝1,2,…,J，pj为所述探测器的第j个像素在所述目标角平面(α,β)各个方向探测到的伽马光子数量组成的向量，M的第j行、i列的元素，Mji表示从所述目标角平面(α,β)的(αi,βi)方向入射的伽马光子在所述伽马相机的所述探测器的第j个像素上被探测到的概率；
通过以下期望最大化的统计迭代算法对所述伽马辐射图像 求解：
其中， 为第n次迭代得到的伽马辐射图像中第i个像素的像素估计值， i＝1,
2,…,I。
11.如权利要求10所述的设备，其特征在于，通过以下步骤预先建立入射光子能量Epeak与伽马光子的探测效率的二维角度分布ζ(Epeak)(α,β)的对应关系：
获取所述目标角平面(α,β)各个方向入射到达所述探测器前表面的能量为Epeak的伽马光子被探测到的二维角度分布的第一概率εdet(Epeak)(α,β)；
获取所述目标角平面(α,β)各个方向入射的伽马光子穿透所述探测器的准直器到达所述探测器表面的二维角度分布的第二概率εcol(Epeak)(α,β)；
根据所述第一概率εdet(Epeak)(α,β)和所述第二概率εcol(Epeak)(α,β)获取与所述入射光子能量Epeak对应的伽马光子的探测效率的二维角度分布ζ(Epeak)(α,β)，其中，ζ(Epeak)(α,β)＝εcol(Epeak)(α,β)×εdet(Epeak)(α,β)。
12.如权利要求8所述的设备，其特征在于，所述投影数据分析模块具体用于：
根据所述入射光子能量Epeak查找相应的传输矩阵M(α,β)以及相应的从所述目标角平面(α,β)的预设角度(αm,βn)入射的伽马光子的探测效率ζ(Epeak)(m,n)；
根据所述传输矩阵M(α,β)和所述探测效率ζ(Epeak)(m,n)获取考虑实际传输效率的传输矩阵M'；
根据所述考虑实际传输效率的传输矩阵M'对所述放射性物质的矢量投影数据p进行重建以得到伽马辐射图像 其中，
其中，i＝1,2,…,I，j＝1,2,…,J，pj为所述探测器的第j个像素在所述目标角平面(α,β)各个方向探测到的伽马光子数量组成的向量，M'的第j行、i列的元素，M'ji表示从所述目标角平面(α,β)的(αi,βi)方向入射的伽马光子在所述伽马相机的所述探测器的第j个像素上被探测到的概率；
通过以下期望最大化的统计迭代算法对所述伽马辐射图像 求解：
其中， 为第n次迭代得到的伽马辐射图像中第i个像素的像素估计值， i＝1,
2,…,I。
13.如权利要求12所述的设备，其特征在于，所述注量率分布获取模块具体用于：
根据所述伽马辐射图像 获取所述目标角平面(α,β)各个方向实际入射光子数量的二维角度分布N'(α,β)；
根据所述实际入射光子数量的二维角度分布N'(α,β)获取所述实际入射光子计数率的二维角度分布N'(α,β)/T，其中，T为探测从所述目标角平面(α,β)各个方向入射的伽马光子的探测时间；
通过以下公式获取所述伽马光子的注量率的二维角度分布φ(Epeak)(α,β)：
其中，S为所述探测器的有效探测面积。
14.如权利要求12或13所述的设备，其特征在于，通过以下步骤预先建立入射光子能量Epeak与所述目标角平面(α,β)的预设角度(αm,βn)入射的伽马光子的探测效率ζ(Epeak)(m,n)的对应关系：
获取从所述目标角平面的预设角度(αm,βn)入射到达所述探测器前表面的能量为Epeak的伽马光子被探测到的第一概率εdet(Epeak)(m,n)；
获取从所述目标角平面的预设角度(αm,βn)入射的伽马光子穿透所述探测器的准直器到达所述探测器表面的第二概率εcol(Epeak)(m,n)；
根据所述第一概率εdet(Epeak)(m,n)和所述第二概率εcol(Epeak)(m,n)获取与所述入射光子能量Epeak对应的伽马光子的探测效率ζ(Epeak)(m,n)，其中，
ζ(Epeak)(m,n)＝εcol(Epeak)(m,n)×εdet(Epeak)(m,n)。