专利类型 | 发明授权 | 法律事件 | 公开; 实质审查; 授权; |
专利有效性 | 有效专利 | 当前状态 | 授权 |
申请号 | CN202111056477.8 | 申请日 | 2021-09-09 |
公开(公告)号 | CN113777648B | 公开(公告)日 | 2024-04-12 |
申请人 | 南京航空航天大学; | 申请人类型 | 学校 |
发明人 | 王泽宇; 汤晓斌; 龚频; 王鹏; 胡志猛; 梁大戬; 沈小磊; | 第一发明人 | 王泽宇 |
权利人 | 南京航空航天大学 | 权利人类型 | 学校 |
当前权利人 | 南京航空航天大学 | 当前权利人类型 | 学校 |
省份 | 当前专利权人所在省份:江苏省 | 城市 | 当前专利权人所在城市:江苏省南京市 |
具体地址 | 当前专利权人所在详细地址:江苏省南京市御道街29号 | 邮编 | 当前专利权人邮编:210016 |
主IPC国际分类 | G01T1/29 | 所有IPC国际分类 | G01T1/29 ; G03B42/00 ; G06N3/084 |
专利引用数量 | 0 | 专利被引用数量 | 0 |
专利权利要求数量 | 8 | 专利文献类型 | B |
专利代理机构 | 北京邦中知识产权代理有限公司 | 专利代理人 | 方岩; 张君; |
摘要 | 本 申请 公开了一种基于随机编码与神经网络探测器成像的方法,包括:S1、获得多个穿过编码 准直 器 的编码图像信息;S2、将所述编码图像信息划分为训练集和测试集,利用所述训练集构建反向神经网络模型;S3、利用所述训练集对所述反向神经网络模型进行训练,获得训练后的神经网络模型,再利用所述测试集对所述训练后的神经网络模型进行测试;S4、若测试结果中 放射源 位置 识别准确率、放射源数量识别准确率、解码图像 信噪比 满足预设值时,则神经网络模型训练完成后即为探测器成像方法,若测试结果达不到要求,则重复步骤S2、S3、S4。使用本 发明 方法得到的图像结果可以准确反应放射源的真实位置信息,同时相比于传统方法具有更高的信噪比和图像清晰度。 | ||
权利要求 | 1.一种基于随机编码与神经网络探测器成像的方法,其特征在于,包括: |
||
说明书全文 | 一种基于随机编码与神经网络探测器成像的方法及伽马相机技术领域背景技术[0002] 随着我国国民经济和核能利用的大力发展,核能开发、国土安全与核反恐备受关注,这就使得对核设施辐射安全、放射源的检测与监控提出了更高的要求。在一些放射事故中,若发生放射源丢失或放射源弥散的情况,即需要迅速判断放射性核素所在。在核探测技术领域,相对于传统的数值、谱线等辐射探测手段,射线成像技术可在远距离给出反映辐射热点分布情况的二维图像,并可以与光学图像相融合,具有效率高、效果直观等突出优点。 [0003] 编码孔径成像技术(Coded Aperture Imaging)是在单针孔成像的基础上发展起来的一种高灵敏度的射线成像方法。传统的单孔成像技术虽然可以获得具有较高信噪比的图像,但是成像时间长,效率低下。其通过采用预设图样的多开孔准直器,大大提高了系统的探测效率和图像信噪比,显著降低了系统的成像时间,并逐渐在放射源检测与监控、核能开发和生物医学等领域广泛应用。 [0004] 编码孔径成像过程分两步:投影和重建,光源透过编码准直器后将编码阵列图案投影至探测器平面,将探测数据与编码阵列做相关运算即可获得光源的空间位置。编码准直器和位置灵敏探测器分别对应两个步骤的关键设备,他们的参数影响着相机的成像性能,如图1所示。 [0005] (1)重建图像的信噪比,取决于编码阵列的自相关性,图像的信噪比比越高,则图像越清晰,分辨细节的能力受噪声影响越小; [0007] (3)编码视场的几何半角宽: (全编码), (部分编码),此指标指相机可成像的面积的多少,即对多大的区域进行成像; [0008] 其中,各个字母含义如下,D:码版到探测器的距离;d:码版元素的大小;S:码版的大小;x:探测器的大小。 [0009] 但现有的编码孔成像技术要求探测器在任意时刻必须探测到完整的编码阵列才可以重建最好的图像,例如像素阵列为n×n的探测器,则必须要求编码阵列N1×N2(N1≤n、N2≤n)。在应用过程中就很难与常见探测器探测器进行完全匹配,这就造成了探测器像素的浪费,同时无法最大化编码阵列,造成嵌套后成像系统视场角(FOV)受限。 [0010] 编码准直器上常用两个编码阵列,修正的均匀冗余阵列(Modified Uniformly Redundant Arrays,MURA)与均匀冗余阵列(Uniformly Redundant Arrays,URA)。他们具有开孔率大(约 50%)、在有限大小内的具有理论的相关特性等优点从而被广泛应用。URA可以构建N1×N2的阵列,N1、N2互质且N1‑N2=2。MURA可以构建N×N的阵列。在实际应用中为了扩大全编码视场,一般在原有码版排列的基础上向四周循环复制得到嵌套码板,嵌套后码版阵列为(2×N1‑1)×(2×N2‑1),其上任意连续N1×N2区域保持原有编码信息。 [0011] 上述原理表明,编码阵列最大只能约为探测器阵列的四倍。常见商用探测器阵列造价高,且阵列固定,由于MURA和URA阵列的特殊性在应用过程中就很难与常见探测器探测器进行完全匹配,这就造成了探测器像素的浪费,同时无法最大化编码阵列,造成嵌套后成像系统视场角FOV受限。 [0012] 随机阵列(Random Array,RA)可以不受行列的限制获得足够大的阵列,从而不受设计限制最大化系统FOV,同时能最大化利用探测器像素。但是传统的相关运算重建应用到基于随机编码的成像中时,会在重建过程中引入固有的大量的噪声,使得图像信噪比低下。如利用迭代法解码,可以获得优异信噪比的图像,但同时会消耗大量运算时间,一定程度上限制了其在实时、快速成像方面应用。 [0013] 由上可知,现有技术下编码孔径成像装置受限于编码阵列编码规则、探测器阵列大小等因素,无法兼顾相机成像质量和系统FOV。基于上述,需要一种能拓展相机成像视场角FOV性能并且能快速获得最优化重建图像信噪比的重建技术发明内容 [0014] 针对现有技术的不足,本发明提供一种基于随机编码与神经网络探测器成像的方法和一种伽马相机,通过采用基于反向神经网络的编码阵列规则,可以兼顾相机成像质量和系统视场角。 [0015] 本发明提供了一种基于随机编码与神经网络探测器成像的方法,包括:S1、获得多个穿过编码准直器的编码图像信息;S2、将所述编码图像信息划分为训练集和测试集,利用所述训练集构建反向神经网络模型;S3、利用所述训练集对所述反向神经网络模型进行训练,获得训练后的神经网络模型,再利用所述测试集对所述训练后的神经网络模型进行测试;S4、若测试结果中放射源位置识别准确率、放射源数量识别准确率、解码图像信噪比满足预设值时,则神经网络模型训练完成后即为探测器成像方法,若测试结果达不到要求,则重复步骤S2、S3、S4。 [0016] 进一步地,获得所述多个穿过所述编码准直器的编码图像信息的方法包括:使用蒙特卡罗软件构建包括随机编码阵列的编码孔准直器模型和探测器模型;模拟不同放射源发射射线穿过所述编码孔准直器模型;所述探测器模型接收所述不同位置放射源的编码图像信息。 [0017] 进一步地,模拟的所述不同位置的放射源位于所述编码孔准直器模型和所述探测器模型的视场角内。 [0018] 进一步地,将同一时间的两个或多个编码图像叠加并归一化处理,直到达到训练集和调试集所需的编码图像数据数量。 [0019] 进一步地,所述蒙特卡罗方法采用MORSE、MCNP、EGS、GEANT4、FLUKA、SuperMC、Phits和GADRAS中的一种或几种。 [0020] 进一步地,所述反向神经网络包括:输入层,所述输入层节点数为探测器像素个数;输出层,所述输出层节点数为接收的所述编码图像信息的像素个数;利用所述输入层、输出层计算所述隐含层节点数,并利用所述激活函数构建所述反向神经网络。 [0021] 进一步地,所述反向传播神经网络框架为DeepLearnToolbox、Caffe、CNTK、TensorFlow、Theano、Torch、Keras、Lasagne、DSSTNE、MXNet、DeepLearning4J、ConvNetJS、Chainer和ScikitLearn中的一种或多种。 [0022] 进一步地,放射源位置识别准确率预设值为大于或等于95%,放射源数量识别准确率预设值为大于或等于95%,解码图像信噪比预设值为大于等于40。 [0023] 本申请还公开了一种利用上述的方法的伽马相机,用于发射射线的放射源识别,包括:编码准直器,所述射线穿过编码准直器;探测组件,包括接收器、处理器及存储器,所述接收器接收穿过编码准直器的伽马射线投影的编码阵列图案,所述存储器存储有计算机指令,当所述计算机指令被所述处理器执行时,使得所述处理器执行对所述投影的编码阵列图案进行成像的方法。 [0025] 图1为编码孔成像示意图。 [0026] 图2为基于随机编码与神经网络探测器成像的方法流程图。 [0027] 图3为编码孔伽马相机的模型及探测到地数据模型。 [0028] 图4为迅速获得大量训练数据的方法示意图。 [0029] 图5展示了利用本申请中的方法对2‑5个放射源同时存在的解码图像结果,和相关运算重建方法的解码图像结果。 具体实施方式[0030] 以下结合附图和实施例,对本发明的具体实施方式进行更加详细的说明,以便能够更好地理解本发明的方案以及其各个方面的优点。然而,以下描述的具体实施方式和实施例仅是说明的目的,而不是对本发明的限制。 [0031] 下文的公开提供了许多不同的实施方式或例子用来实现本发明的不同结构。为了简化本发明的公开,下文中对特定例子的部件和设置进行描述。当然,它们仅仅为示例,并且目的不在于限制本发明。此外,本发明可以在不同例子中重复参考数字和/或参考字母,这种重复是为了简化和清楚的目的,其本身不指示所讨论各种实施方式和/或设置之间的关系。此外,本发明提供了的各种特定的工艺和材料的例子,但是本领域普通技术人员可以意识到其他工艺的应用和/或其他材料的使用。 [0032] 首先,需解释下BPNN(反向传播神经网络)的由来。人工神经网络(Artificial Neural Network,即ANN ),是20世纪80 年代以来人工智能领域兴起的研究热点。最近十多年来,人工神经网络的研究工作不断深入,已经取得了很大的进展,而反向传播算法(Back‑Propagation)作为最早且最普遍使用的多层感知机网络训练方法,早在70、80 年代就被提出来,并一直沿用至今。以误差反向传播算法为核心的反向传播人工神经网络(BPNN)也同样应用广泛。 [0033] 反向神经网络(BP)模型拓扑结构包括输入层、隐含层以及输出层,输入层的神经元由样本属性的维度决定,输出层的神经元个数由样本分类个数决定,隐含层则可以依据实际情况自己决定。从本质上讲,BP算法就是以网络误差平方目标函数、采用梯度下降法来计算目标函数的最小值。基本BP算法包括信号的前向传播和误差的反向传播两个过程。BPNN(反向传播神经网络)具有诸多优点,比如非线性映射能力强、高度的自学习和自适应能力、容错能力强等。使用BPNN,可以较为有效的提高整个剂量率高度修正的效率,实现核素剂量率高度修正的一体式修正。 [0034] 此外,还需简单介绍下本申请中利用的MCNP模拟软件。在MCNP软件是一种利用蒙特卡洛模拟程方法的软件。核技术的学科特点导致很多情况下无法进行实际的测量或实验工作,此时计算机模拟计算就显示出了独特的优势。同时对于结构和反应机制复杂的核反应问题,一般的数值方法难以求解,而MC(Monte Carlo)方法能够准确地模拟实际中的物理过程,解决以往数值方法难以解决的问题,因此该方法广泛地应用于涉核领域的研究中。上世纪 40 年代中期随着原子能事业的发展,MC方法逐渐成长起来,其基本思想是基于随机数选择的统计抽样方法,由于传统的经验方法不能逼近真实的物理过程, 难以得到满意的结果,在解决粒子输运问题方面,蒙特卡罗方法有着显著的优势,随后人们逐渐开发了众多的 MC 模拟程序,目前主流的有EGS、MCNP、GEANT4 等,其中MCNP是由美国洛斯阿拉莫斯国家实验室开发,主要用于三维几何结构中光子、电子、中子等在物质中输运问题的蒙特卡罗模拟程序,其适用的光子能量范围在1E‑3MeV—1E5MeV,电子能量范围在1E‑3MeV—1E3MeV,中子能量范围1E‑11MeV—20MeV。MCNP的程序功能较全,各种物质反应截面数据丰富,减小方差的方法多样,通用性强且使用简单。 [0035] 图2示出了基于随机编码与神经网络探测器成像的方法流程图。 [0036] 如图2所示,首先获取多个穿过编码准直器的编码图像信息(104步),这些图像信息可以划分为训练集和测试集,利用训练集中的训练数据构建反向神经网络模型,并利用训练数据进行神经网络训练(105步)。获得训练后的神经网络模型,再利用所述测试集对所述训练后的神经网络模型进行测试。若测试结果中射源位置识别准确率大于或等于95%,放射源数量识别准确率大于或等于95%,解码图像信噪比大于等于40时(106步),则训练后的神经网络模型即为探测器成像方法(107步)。若测试结果不满足上述要求,则需要重复105‑107步过程。 [0037] 需要注意的是,当为单个放射源时,位置识别准确率需大于或等于99%,当为多个放射源时,位置识别准确率大于或等于95%, [0038] 如图2所示,流程图中还示出了模型搭建的过程。首先需确定阵列探测器类型以及参数和规格,编码孔径阵列大小及孔径、厚度设计等参数(101、102步)。参考以上参数,使用蒙特卡罗软件构建包括随机编码阵列的编码孔准直器模型和探测器模型(103步),和模拟编码孔径的成像过程,获得训练数据。具体过程为模拟不同放射源发射射线穿过所述编码孔准直器模型;所述探测器模型接收所述不同位置放射源的编码图像信息。编码图像信息中包括有放射源地位置信息,此时编码图像信息即为训练数据。 [0039] 进一步地,在步骤101,需确定选用的阵列探测器的参数。然后根据步骤102对编码孔径参数进行初始设计,所需考虑的设计参数包括准直器的材料、尺寸、孔径的大小、孔的形状、孔的排列方式等。需要说明的是,101‑102步可以根据使用者所具备的客观条件特殊设计。 [0041] 对于步骤104,本实施例采用Co‑57(活度为798E4 Bq)点源,放置于相机视野内距离相机2m处,分别模拟测量时间为120s、60 s、30s、10s、的情况,统计接收到的编码图像信息。再将模拟获得的编码图像和归一化处理至(‑1,1),结合对应空间位置作为输出标签输入到网络模型中,作为训练数据。 [0042] 步骤105采用的深度学习网络框架是DeepLearnToolbox。利用上述的训练数据构建神经网络。 [0043] 进一步地,根据步骤106对初步训练的深度学习模型进行测试,如所选取的测试样本为模拟测量时间为10 s,1s的已经归一化的样本,且测试样本随机包含1‑9多点源成像数据。需要检验的性能参数包括放射源位置识别准确率、放射源数量识别准确率、解码图像信噪比。如果测试参数满足要求,则获得训练完成的深度学习网络。本实施例中设置单放射源位置识别准确率为0.99、多放射源(5放射源)位置识别准确率为0.95、放射源数量(1‑5源)识别准确率为0.99、信噪比为40条件来判断网络是否完成训练并是否可进行应用,如不满足,则重复步骤105‑107。 [0044] 下面结合实际模型来详细描述成像过程。图3示出了编码孔伽马相机的模型及探测到地数据模型。如图3所示,编码孔伽马相机包括阵列探测器201和编码孔准直器202,此外还包括屏蔽或者相机外壳203,阵列探测器201和编码孔准直器202置于屏蔽或者相机外壳203中。其中编码孔准直器202中具有编码板。 [0045] 进一步地,204为伽马相机全编码探测目标区域,205为伽马相机部分编码探测目标区域,本实施例中只考虑全编码探测区域即204区域。206、207所示为放射源所在位置;206处于全编码视场角FOV内,探测器探测到的数据如图中208所示。例如,此时探测数据为 16×16矩阵。207处于部分编码视场角FOV内,探测器数据不采用。图示209中,示出一个实施例,如只考虑全编码视场角FOV,则全编码探测区域204被离散成17×17的像素化平面,即探测器可以探测到17×17个连续的且不同的16×16编码阵列,放射源所处位置为1(如图中白色点),其余为0(如图中黑色区域),并构建对应位置标签。获得探测数据及对应位置标签作为训练数据的输入数据,即可进行上述步骤105,构建反向传播神经网络模型,训练网络。 [0046] 上述表示了一个训练数据的获得过程,但是一个完整的具有高识别率的反向传播人工神经网络需要大量的训练数据进行训练。本申请还示出了如何迅速获得大量训练数据。图4为迅速获得大量训练数据的方法示意图。图4即对上述的图2中104步(蒙特卡罗模拟编码孔径成像过程,获取训练数据)做了进一步说明。 [0047] 如图4所示,通过不断改变放射源206在204全编码探测目标区域内放射源位置,即可以在阵列探测器201上获得大量编码图像,但是还远远不够网络训练的需求量。步骤301本实施例中模拟了不同时间下单个放射源在离散探测目标平面内的所有17×17种可能成像结果。步骤302在301的基础上,随机抽取同时间下301中2种探测数据进行叠加并归一化,对应位置标签叠加并归一化,本实施例重复以上操作1000000次获得2点源数据。步骤303在301的基础上,随机抽取同时间下301中3种探测数据进行叠加并归一化,对应位置标签叠加并归一化,本实施例重复以上操作1000000次获得3点源数据。步骤304在301的基础上,随机抽取同时间下301中4种探测数据进行叠加并归一化,对应位置标签叠加并归一化,本实施例重复以上操作1000000次获得4点源数据。步骤305在301的基础上,随机抽取同时间下301中5种探测数据进行叠加并归一化,对应位置标签叠加并归一化,本实施例重复以上操作 1000000次获得5点源数据。将上述原始单点源数据301,以及叠加所得302、303、304、305多点源数据打包输入数据,对构建的反向神经网络进行训练,即图2中的步骤105。 [0048] 还需特别指出的是在本申请实施例中,全编码成像构建的反向传播神经网络参数如下: [0049] 1层输入层,输入层节点数需要为探测器的像素个数,如节点数为16×16。1层输出层,输出层节点数等于输出图像的像素个数,如本实施例中节点数为17×17。为了保证训练速度,设置1层隐含层,隐含层节点依照经验公式,设置为300。 [0050] 本实施例中还设置了所需的其他参数,具体如下: [0051] 全连接层激活函数为Sigmoid函数。激活函数学习速率0.01。平均下降梯度步长200。训练数据全扫描次数500。 [0052] 本申请在不同的使用场景中,需调整对应的构建的反向网络参数: [0053] 场景一:探测器阵列固定、编码准直器至探测器平面距离固定、码孔尺寸固定的情况下。假定设计编码阵列为N1×N2、探测器阵列为N3×N4。N1、N2的取值决定了全编码FOV的大小。对应全编码FOV被离散成(N1+1‑N3)×(N2+1‑N4)的像素化平面。在进行步骤104、105时,需要调整对应的输入层节点以及输出层节点。 [0054] 场景二:编码准直器阵列固定、编码准直器至探测器平面距离固定、码孔尺寸固定的情况下,假定设计编码阵列为N1×N2、探测器阵列为N3×N4。N3、N4的取值决定了全编码FOV以及部分编码FOV的大小。对应全编码FOV需被离散成(N1+1‑N3)×(N2+1‑N4)的像素化平面。在进行步骤104、105时,需要调整对应的输入层节点以及输出层节点。 [0055] 本申请还涉及一种伽马相机,用于发射射线的放射源识别,包括:编码准直器,射线穿过编码准直器。探测组件,包括接收器和处理器及存储器,接收器接收穿过编码准直器的伽马射线投影的编码阵列图案,存储器存储有计算机指令,当计算机指令被所述处理器执行时,使得处理器执行对投影的编码阵列图案进行成像的方法。 [0056] 本申请具有以下优点: [0057] (1)用采用反向传播神经网络作为重建算法,可以有效抑制重建背景噪声,获得高信噪比图像,且重建效果不受编码阵列编码规则的限制。相比于最大似然估计法,重建迅速;相比于相关运算重建,重建图像信噪比优异。 [0058] (2)基于固定的探测器阵列,使用随机编码作为相机码板的编码规则,可以灵活设计编码准直器的阵列大小、形状以及开孔率。突出的优点是获得较大的编码阵列,从而获得较大的理论成像FOV视场角。 [0059] (3)基于固定的编码准直器,可以通过相对小阵列的阵列探测器重建较大视野的射线热点图像,即且获得更大的成像视场角FOV的同时减少探测器成本。 [0060] 图5展示了利用本申请中的方法对2‑5个放射源同时存在的解码图像结果,和相关运算重建方法的解码图像结果。 [0061] 如图5所示,可以见得,使用本发明所述的反向传播神经网络重建方式得到的图像结果可以准确反应放射源的真实位置信息,同时相比于传统方法具有更高的信噪比和图像清晰度。 [0062] 显然,上述实施例仅仅是为清楚地说明本发明所作的举例,而并非对实施方式的限定。对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。而由此所引申出的显而易见的变化或变动仍处于本发明的保护范围之中。 |