CN103645491B - 放射源定位方法、装置及*** - Google Patents
放射源定位方法、装置及*** Download PDFInfo
- Publication number
- CN103645491B CN103645491B CN201310602906.6A CN201310602906A CN103645491B CN 103645491 B CN103645491 B CN 103645491B CN 201310602906 A CN201310602906 A CN 201310602906A CN 103645491 B CN103645491 B CN 103645491B
- Authority
- CN
- China
- Prior art keywords
- crystal
- radioactive source
- count
- counting
- angle
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Measurement Of Radiation (AREA)
Abstract
本发明公开了一种放射源定位方法、装置及***。该方法包括:根据模拟数据获取探测器的晶体阵列中晶体计数比与放射源角度的关系;测量待测放射源放射射线时所述晶体阵列的实测晶体计数;在计数噪声的变化范围内,根据所述实测晶体计数以及各所述计数噪声获取多个去噪晶体计数比,根据所述晶体计数比与放射源角度的关系获取与各所述去噪晶体计数比相对应的各实测放射源角度,并选取计数噪声最优时的实测放射源角度为所述待测放射源的放射源角度。该装置包括:处理模块、测量模块和计算模块。该***包括:放射源定位装置、待测放射源和探测器。本发明采用晶体计数比作为计算参数,采用有效的扣除噪声技术,提高了放射源定位的灵敏度和准确度。
Description
技术领域
本发明涉及放射源探测技术,特别涉及一种放射源定位方法、装置及***。
背景技术
放射源探测可广泛应用于各种放射性环境,如实验室、核电厂、核燃料厂、大型重离子加速器、放射源库等。用于探测伽马射线的探测器多种多样,目前较先进的技术是利用伽马射线中的粒子能够被碘化铯(CsI)或锗酸铋(BGO)等晶体吸收的原理,采用四块晶体组成2×2晶体阵列,对同一放射源进行探测并计数,根据各块晶体的计数值来反推出放射源的角度,以确定放射源的方位。
但是,当放射源的位置与晶体阵列不在同一平面上时,计算出的放射源角度不够准确,使得定位准确度低;同时,由于晶体计数值随角度的变化不明显,使得现有的放射源定位方法灵敏度不高,分辨能力偏低。
发明内容
在下文中给出关于本发明的简要概述,以便提供关于本发明的某些方面的基本理解。应当理解,这个概述并不是关于本发明的穷举性概述。它并不是意图确定本发明的关键或重要部分,也不是意图限定本发明的范围。其目的仅仅是以简化的形式给出某些概念,以此作为稍后论述的更详细描述的前序。
本发明提供一种放射源定位方法、装置及***,用以解决现有技术中放射源定位准确度低及灵敏度低的缺陷,实现高灵敏度、高准确度的放射源定位。
本发明提供了一种放射源定位方法,包括:
步骤S1、根据理拟数据获取探测器的晶体阵列中晶体计数比与放射源角度的关系;
步骤S2、测量待测放射源放射射线时所述晶体阵列的实测晶体计数;
步骤S3、在计数噪声的变化范围内,根据所述实测晶体计数以及各所述计数噪声获取多个去噪晶体计数比,根据所述晶体计数比与放射源角度的关系获取与各所述去噪晶体计数比相对应的各实测放射源角度,并选取所述计数噪声最优时的实测放射源角度为所述待测放射源的放射源角度。
本发明还提供了一种放射源定位装置,包括:
处理模块,用于根据模拟数据获取探测器的晶体阵列中晶体计数比与放射源角度的关系;
测量模块,用于测量待测放射源放射射线时所述晶体阵列的实测晶体计数;
计算模块,用于在计数噪声的变化范围内,根据各所述计数噪声以及所述测量模块测量得到的所述实测晶体计数获取多个去噪晶体计数比,根据所述处理模块获取的所述晶体计数比与放射源角度的关系获取与各所述去噪晶体计数比相对应的各实测放射源角度,并选取所述计数噪声最优时的实测放射源角度为所述待测放射源的放射源角度。
本发明还提供了一种放射源定位***,包括:上述放射源定位装置、待测放射源和探测器。
本发明提供的放射源定位方法、装置及***,采用晶体计数比作为计算参数,使参数随放射源角度的变化明显,采用有效的扣除噪声技术,对于放射源不在晶体阵列平面上的情况,能够有效提高参数中与放射源角度有关的变化因子,避免了现有技术的缺陷,提高了放射源定位的灵敏度和准确度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明放射源定位方法第一实施例的流程图;
图2为本发明放射源定位方法第二实施例的流程图;
图3为本发明实施例采用的2×2晶体阵列的示意图;
图4为本发明实施例放射源角度与晶体计数比的关系图;
图5为本发明放射源定位装置第一实施例的结构示意图;
图6为本发明放射源定位装置第二实施例的结构示意图;
图7为本发明放射源定位***实施例的***框图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。在本发明的一个附图或一种实施方式中描述的元素和特征可以与一个或更多个其它附图或实施方式中示出的元素和特征相结合。应当注意,为了清楚的目的,附图和说明中省略了与本发明无关的、本领域普通技术人员已知的部件和处理的表示和描述。基于本发明中的实施例,本领域普通技术人员在没有付出创造性劳动的前提下所获得的所有其他实施例,都属于本发明保护的范围。
图1为本发明放射源定位方法第一实施例的流程图。如图1所示,本发明实施例提供了一种放射源定位方法,包括:
步骤S10、根据模拟数据获取探测器的晶体阵列中晶体计数比与放射源角度的关系,其中,晶体计数比为晶体阵列中相邻两块晶体的计数和与晶体阵列的计数和之比;
步骤S20、测量待测放射源放射射线时晶体阵列的实测晶体计数,其中,实测晶体计数即为晶体阵列实际测量得到的晶体计数;
步骤S30、在计数噪声的变化范围内,根据实测晶体计数以及各计数噪声获取多个去噪晶体计数比,根据晶体计数比与放射源角度的关系获取与各去噪晶体计数比相对应的各实测放射源角度,并选取计数噪声最优时的实测放射源角度为待测放射源的放射源角度,其中,去噪晶体计数比为利用去噪后的晶体计数计算出的晶体计数比。
在本发明实施例中,首先根据模拟数据获取晶体计数比与放射源角度的关系,该模拟数据是通过模拟软件模拟放射源与晶体阵列置于同一平面上的情况而获得,由此获得的关系中计数噪声可视为零。然后,针对某待测环境,利用上述探测器的晶体阵列测量实测晶体计数。当待测放射源与晶体阵列不在同一平面上时,所有晶体上的计数均将高于在同一平面上的计数,这个计数的增加(即为,计数噪声)使得晶体计数比偏离上述晶体计数比与放射源角度的关系。因此,可以根据具体实测数据确定计数噪声的变化范围。对实测晶体计数进行去噪处理,得到去噪晶体计数比,并获取与其对应的实测放射源角度。最后,在计数噪声的变化范围内,选取计数噪声最优时的实测放射源角度,作为待测放射源的放射源角度。
本发明实施例提供的放射源定位方法,采用晶体计数比作为计算参数,使参数随放射源角度的变化明显,采用有效的扣除噪声技术,对于放射源不在晶体阵列平面上的情况,能够有效提高参数中与放射源角度有关的变化因子,避免了现有技术的缺陷,提高了放射源定位的灵敏度和准确度。
图2为本发明放射源定位方法第二实施例的流程图。如图2所示,在上述方法第一实施例的基础上,上述晶体阵列可以为2×2晶体阵列,上述步骤S10可以进一步具体包括:
步骤S101、在2×2晶体阵列中任选一块晶体作为参照晶体;
步骤S102、将参照晶体及与其相邻的一块晶体组合为第一晶体组,将参照晶体及与其相邻的另一块晶体组合为第二晶体组;
步骤S103、利用模拟软件测量第一晶体组的第一晶体计数比AB和第二晶体组的第二晶体计数比BC,其中,AB=(A+B)/(A+B+C+D),BC=(B+C)/(A+B+C+D),A、B、C和D依次为该2×2晶体阵列中四块晶体的晶体计数;
步骤S104、分别获取第一晶体计数比AB与放射源角度θ的离散关系以及第二晶体计数比BC与放射源角度θ的离散关系;
步骤S105、将上述离散关系拟合为线性关系。
图3为本发明实施例采用的2×2晶体阵列的示意图。如图3所示,a、b、c和d为2×2晶体阵列中的四块晶体。假如,选晶体b为参照晶体,将晶体a与晶体b组合为第一晶体组合,将晶体b和晶体c组合为第二晶体组合。A、B、C和D依次为晶体a、b、c和d的计数,利用模拟软件测量第一晶体计数比AB和第二晶体计数比BC。分别获取AB与θ的离散关系以及BC与θ的离散关系。图4为本发明实施例放射源角度与晶体计数比的关系图。如图4所示,AB与θ的关系可以分为两个部分:θ∈(-90°,90°)与θ∈(90°,270°),在两个区间内,AB与θ关系均接近线性关系。同样,BC在θ∈(0°,180°)与θ∈(180°,360°)两个区间内也是接近线性关系。因此,对其进行拟合,得到线性关系θ=(AB-y1)/x1,θ=(BC-y2)/x2,其中,x1、x2、y1和y2为直线拟合参数。
本发明实施例在探测器工作时,可以通过测量直接计算出AB和BC,因此可以通过上述公式计算出θ,计算极其简单,同时从图4中可以看出,晶体计数比与放射源角度的变化明显,因此定位灵敏度很高。
更进一步地,待测放射源放射射线时,晶体a、b、c和d的实测晶体计数可以依次计为A’、B’、C’和D’,计数噪声N的变化范围可以设为Nmin至Nmax。一般Nmin可以设为0,Nmax可以设为A’、B’、C’和D’中最小值的60%。上述步骤S30可以具体包括:
步骤S301、将计数噪声N设置为Nmin;
步骤S302、根据AB’=(A’+B’-2×N)/(A’+B’+C’+D’-4×N)获取第一去噪晶体计数比,根据BC’=(B’+C’-2×N)/(A’+B’+C’+D’-4×N)获取第二去噪晶体计数比;
步骤S303、根据线性关系获取AB’对应的放射源角度θ1和BC’对应的放射源角度θ2;
步骤S304、计数噪声N自加1;
步骤S305、比较|θ1-θ2|和的大小,和为预设参数,若则设置参数否则,参数和不变;
步骤S306、判断计数噪声N是否大于Nmax,若是,则执行步骤S307,否则,执行步骤S302;
步骤S307、根据获取待测放射源的放射源角度。
当待测放射源与晶体阵列不在同一平面上时,所有晶体上的计数均将高于在同一平面上的计数,这个计数的增加(即为,计数噪声)使得晶体计数比AB和BC偏离图4中所示的线性曲线。本发明实施例中,假设晶体a、b、c和d上的计数噪声均为N,则得到去噪后的晶体计数比:
AB’=(A’-N+B’-N)/(A’-N+B’-N+C’-N+D’-N)
=(A’+B’-2×N)/(A’+B’+C’+D’-4×N),
BC’=(B’-N+C’-N)/(A’-N+B’-N+C’-N+D’-N)
=(B’+C’-2×N)/(A’+B’+C’+D’-4×N)。
根据实际情况由于对应同一放射源,角度应该是相同的。因此,计数噪声的确定方法为,在一定范围内连续改变噪声大小,使得通过曲线AB与曲线BC计算出的角度最接近,此时的计数噪声,即为最接近真实的最优计数噪声。此时计算出的角度即为最接近真实的放射源实际角度。
本发明实施例中的待测放射源可以为伽马射线放射源。
本发明实施例提供的放射源定位方法,采用晶体计数比作为计算参数,使参数随放射源角度的变化明显,采用有效的扣除噪声技术,对于放射源不在晶体阵列平面上的情况,有效地去除无益的计数噪声,能够有效提高参数中与放射源角度有关的变化因子,避免了现有技术的缺陷,提高了放射源定位的灵敏度和准确度。
图5为本发明放射源定位装置第一实施例的结构示意图。如图5所示,本发明实施例提供了一种放射源定位装置,包括:处理模块51、测量模块52和计算模块53。其中,处理模块51用于根据模拟数据获取探测器的晶体阵列中晶体计数比与放射源角度的关系;测量模块52用于测量待测放射源放射射线时晶体阵列的实测晶体计数;计算模块53用于在计数噪声的变化范围内,根据各计数噪声以及测量模块52测量得到的实测晶体计数获取多个去噪晶体计数比,根据处理模块51获取的晶体计数比与放射源角度的关系获取与各去噪晶体计数比相对应的各实测放射源角度,并选取计数噪声最优时的实测放射源角度为待测放射源的放射源角度。
在本发明实施例中,首先处理模块51根据模拟数据获取晶体计数比与放射源角度的关系,其中,晶体计数比为晶体阵列中相邻两块晶体的计数和与晶体阵列的计数和之比。该模拟数据是通过模拟软件模拟放射源与晶体阵列置于同一平面上的情况而获得,由此获得的晶体计数比与放射源角度的关系中计数噪声可视为零。然后,针对某待测环境,测量模块52利用上述探测器的晶体阵列测量实测晶体计数。当待测放射源与晶体阵列不在同一平面上时,所有晶体上的计数均将高于在同一平面上的计数,这个计数的增加(即为,计数噪声)使得晶体计数比偏离上述线性关系。因此,计算模块53可以根据具体实测数据确定计数噪声的变化范围。对实测晶体计数进行去噪处理,得到去噪晶体计数比,并获取与其对应的实测放射源角度。最后,在计数噪声的变化范围内,选取计数噪声最优时的实测放射源角度,作为待测放射源的放射源角度。
本发明实施例提供的放射源定位装置,采用晶体计数比作为计算参数,使参数随放射源角度的变化明显,采用有效的扣除噪声技术,对于放射源不在晶体阵列平面上的情况,能够有效提高参数中与放射源角度有关的变化因子,避免了现有技术的缺陷,提高了放射源定位的灵敏度和准确度。
图6为本发明放射源定位装置第二实施例的结构示意图。如图6所示,在上述装置第一实施例的基础上,上述晶体阵列可以为2×2晶体阵列,上述处理模块51可以进一步包括:选择子模块61、测量子模块62、第一处理子模块63和拟合子模块64。其中,选择子模块61用于在2×2晶体阵列中任选一块晶体作为参照晶体,将参照晶体及与其相邻的一块晶体组合为第一晶体组,将参照晶体及与其相邻的另一块晶体组合为第二晶体组;测量子模块62用于利用模拟软件测量第一晶体组的第一晶体计数比AB和第二晶体组的第二晶体计数比BC,其中,AB=(A+B)/(A+B+C+D),BC=(B+C)/(A+B+C+D),A、B、C和D依次为该2×2晶体阵列中四块晶体的晶体计数;第一处理子模块63用于分别获取第一晶体计数比AB与放射源角度θ的离散关系以及第二晶体计数比BC与放射源角度θ的离散关系;拟合子模块64用于将上述第一处理子模块63获取到的离散关系拟合为线性关系。
在图3中,a、b、c和d为2×2晶体阵列中的四块晶体。假如,选择子模块61选晶体b为参照晶体,将晶体a与晶体b组合为第一晶体组合,将晶体b和晶体c组合为第二晶体组合。A、B、C和D依次为晶体a、b、c和d的计数,测量子模块62利用模拟软件测量第一晶体计数比AB和第二晶体计数比BC。第一处理子模块63分别获取AB与θ的离散关系以及BC与θ的离散关系。如图4所示,AB与θ的关系可以分为两个部分:θ∈(-90°,90°)与θ∈(90°,270°),在两个区间内,AB与θ关系均接近线性关系。同样,BC在θ∈(0°,180°)与θ∈(180°,360°)两个区间内也是接近线性关系。因此,拟合子模块64对其进行拟合,得到线性关系θ=(AB-y1)/x1,θ=(BC-y2)/x2,其中,x1、x2、y1和y2为直线拟合参数。
本发明实施例在探测器工作时,可以通过测量直接计算出AB和BC,因此可以通过上述公式计算出θ,计算极其简单,同时从图4中可以看出,晶体计数比与放射源角度的变化明显,因此定位灵敏度很高。
更进一步地,待测放射源放射射线时,晶体a、b、c和d的实测晶体计数可以依次计为A’、B’、C’和D’,计数噪声N的变化范围可以设为Nmin至Nmax。一般Nmin可以设为0,Nmax可以设为A’、B’、C’和D’中最小值的60%。
上述计算模块53可以具体包括:设置子模块65、第一计算子模块66、第二处理子模块67、自加子模块68、比较子模块69、判断子模块71和第二计算子模块72。其中,设置子模块65用于将计数噪声N设置为Nmin;第一计算子模块66用于根据AB’=(A’+B’-2×N)/(A’+B’+C’+D’-4×N)获取第一去噪晶体计数比,根据BC’=(B’+C’-2×N)/(A’+B’+C’+D’-4×N)获取第二去噪晶体计数比;第二处理子模块67用于根据拟合子模块64拟合的线性关系获取AB’对应的放射源角度θ1和BC’对应的放射源角度θ2;自加子模块68用于将计数噪声N自加1;比较子模块69用于比较|θ1-θ2|和的大小,和为预设参数,若则设置参数否则,参数和不变;判断子模块71用于判断自加子模块68得出的计数噪声N是否大于Nmax;第二计算子模块72用于当计数噪声N大于Nmax时,根据获取待测放射源的放射源角度。
当待测放射源与晶体阵列不在同一平面上时,所有晶体上的计数将均高于在同一平面上的计数,这个计数的增加(即为,计数噪声)使得晶体计数比AB和BC偏离图4中所示的线性曲线。本发明实施例中,假设晶体a、b、c和d上的计数噪声均为N,则得到去噪后的晶体计数比:
AB’=(A’-N+B’-N)/(A’-N+B’-N+C’-N+D’-N)
=(A’+B’-2×N)/(A’+B’+C’+D’-4×N),
BC’=(B’-N+C’-N)/(A’-N+B’-N+C’-N+D’-N)
=(B’+C’-2×N)/(A’+B’+C’+D’-4×N)。
根据实际情况由于对应同一放射源,角度应该是相同的。因此,在一定范围内连续改变噪声大小,使得通过曲线AB与曲线BC计算出的角度最接近,此时的计数噪声,即为最接近真实的最优计数噪声。此时计算出的角度即为最接近真实的放射源实际角度。
本发明实施例中的待测放射源可以为伽马射线放射源。
本发明实施例提供的放射源定位装置,采用晶体计数比作为计算参数,使参数随放射源角度的变化明显,采用有效的扣除噪声技术,对于放射源不在晶体阵列平面上的情况,有效地去除无益的计数噪声,能够有效提高参数中与放射源角度有关的变化因子,避免了现有技术的缺陷,提高了放射源定位的灵敏度和准确度。
图7为本发明放射源定位***实施例的***框图。如图7所示,本发明实施例提供一种放射源地位***,包括:放射源定位装置73、待测放射源74和探测器75。放射源定位装置73根据模拟数据获取探测器75中晶体计数比与放射源角度的关系,然后根据该关系计算待测放射源74的放射源角度。
本***实施例中放射源定位装置73的功能如上述图5和图6所示装置实施例中的具体描述,在此不再赘述。
本发明实施例中的待测放射源74可以为伽马射线放射源。
本发明实施例提供的放射源定位***,采用晶体计数比作为计算参数,使参数随放射源角度的变化明显,采用有效的扣除噪声技术,对于放射源不在晶体阵列平面上的情况,有效地去除无益的计数噪声,能够有效提高参数中与放射源角度有关的变化因子,避免了现有技术的缺陷,提高了放射源定位的灵敏度和准确度。
在本发明上述各实施例中,实施例的序号仅仅便于描述,不代表实施例的优劣。对各个实施例的描述都各有侧重,某个实施例中没有详述的部分,可以参见其他实施例的相关描述。
本领域普通技术人员可以理解:实现上述方法实施例的全部或部分步骤可以通过程序指令相关的硬件来完成,前述的程序可以存储于一计算机可读取存储介质中,该程序在执行时,执行包括上述方法实施例的步骤;而前述的存储介质包括:只读存储器(Read-Only Memory,简称ROM)、随机存取存储器(Random Access Memory,简称RAM)、磁碟或者光盘等各种可以存储程序代码的介质。
在本发明的装置和方法等实施例中,显然,各部件或各步骤是可以分解、组合和/或分解后重新组合的。这些分解和/或重新组合应视为本发明的等效方案。同时,在上面对本发明具体实施例的描述中,针对一种实施方式描述和/或示出的特征可以以相同或类似的方式在一个或更多个其它实施方式中使用,与其它实施方式中的特征相组合,或替代其它实施方式中的特征。
应该强调,术语“包括/包含”在本文使用时指特征、要素、步骤或组件的存在,但并不排除一个或更多个其它特征、要素、步骤或组件的存在或附加。
最后应说明的是:虽然以上已经详细说明了本发明及其优点,但是应当理解在不超出由所附的权利要求所限定的本发明的精神和范围的情况下可以进行各种改变、替代和变换。而且,本发明的范围不仅限于说明书所描述的过程、设备、手段、方法和步骤的具体实施例。本领域内的普通技术人员从本发明的公开内容将容易理解,根据本发明可以使用执行与在此所述的相应实施例基本相同的功能或者获得与其基本相同的结果的、现有和将来要被开发的过程、设备、手段、方法或者步骤。因此,所附的权利要求旨在在它们的范围内包括这样的过程、设备、手段、方法或者步骤。
Claims (8)
1.一种放射源定位方法,其特征在于,包括:
步骤S1、根据模拟数据获取探测器的晶体阵列中晶体计数比与放射源角度的关系;
步骤S2、测量待测放射源放射射线时所述晶体阵列的实测晶体计数;
步骤S3、在计数噪声的变化范围内,根据所述实测晶体计数以及各所述计数噪声获取多个去噪晶体计数比,根据所述晶体计数比与放射源角度的关系获取与各所述去噪晶体计数比相对应的各实测放射源角度,并选取所述计数噪声最优时的实测放射源角度为所述待测放射源的放射源角度;其中,
所述晶体阵列为2×2晶体阵列,所述步骤S1包括:
步骤S11、在所述2×2晶体阵列中任选一块晶体作为参照晶体;
步骤S12、将所述参照晶体及与其相邻的一块晶体组合为第一晶体组,将所述参照晶体及与其相邻的另一块晶体组合为第二晶体组;
步骤S13、利用模拟软件测量所述第一晶体组的第一晶体计数比AB和所述第二晶体组的第二晶体计数比BC,其中,AB=(A+B)/(A+B+C+D),BC=(B+C)/(A+B+C+D),A、B、C和D依次为所述2×2晶体阵列中四块晶体的晶体计数;
步骤S14、分别获取所述第一晶体计数比AB与所述放射源角度θ的离散关系以及所述第二晶体计数比BC与所述放射源角度θ的离散关系;
步骤S15、将所述离散关系拟合为线性关系。
2.根据权利要求1所述的放射源定位方法,其特征在于,所述待测放射源放射射线时,所述2×2晶体阵列中四块晶体的实测晶体计数依次为A’、B’、C’和D’,所述计数噪声的变化范围为Nmin至Nmax,所述步骤S3包括:
步骤S31、将所述计数噪声N设置为Nmin;
步骤S32、根据AB’=(A’+B’-2×N)/(A’+B’+C’+D’-4×N)获取第一去噪晶体计数比,根据BC’=(B’+C’-2×N)/(A’+B’+C’+D’-4×N)获取第二去噪晶体计数比;
步骤S33、根据所述线性关系获取AB’对应的放射源角度θ1和BC’对应的放射源角度θ2;
步骤S34、所述计数噪声N自加1;
步骤S35、比较|θ1-θ2|和的大小,和为预设参数,若 则设置参数否则,参数和不变;
步骤S36、判断所述计数噪声N是否大于Nmax,若是,则执行步骤S37,否则,执行步骤S32;
步骤S37、根据获取所述待测放射源的放射源角度。
3.根据权利要求1或2所述的放射源定位方法,其特征在于,所述线性关系为θ=(AB-y1)/x1,θ=(BC-y2)/x2,其中,x1、x2、y1和y2为直线拟合参数。
4.一种放射源定位装置,其特征在于,包括:
处理模块,用于根据模拟数据获取探测器的晶体阵列中晶体计数比与放射源角度的关系;
测量模块,用于测量待测放射源放射射线时所述晶体阵列的实测晶体计数;
计算模块,用于在计数噪声的变化范围内,根据各所述计数噪声以及所述测量模块测量得到的所述实测晶体计数获取多个去噪晶体计数比,根据所述处理模块获取的所述晶体计数比与放射源角度的关系获取与各所述去噪晶体计数比相对应的各实测放射源角度,并选取所述计数噪声最优时的实测放射源角度为所述待测放射源的放射源角度;其中,
所述晶体阵列为2×2晶体阵列,所述处理模块包括:
选择子模块,用于在所述2×2晶体阵列中任选一块晶体作为参照晶体,将所述参照晶体及与其相邻的一块晶体组合为第一晶体组,将所述参照晶体及与其相邻的另一块晶体组合为第二晶体组;
测量子模块,用于利用模拟软件测量所述第一晶体组的第一晶体计数比AB和所述第二晶体组的第二晶体计数比BC,其中,AB=(A+B)/(A+B+C+D),BC=(B+C)/(A+B+C+D),A、B、C和D依次为所述2×2晶体阵列中四块晶体的晶体计数;
第一处理子模块,用于分别获取所述第一晶体计数比AB与所述放射源角度θ的离散关系以及所述第二晶体计数比BC与所述放射源角度θ的离散关系;
拟合子模块,用于将所述第一处理子模块获取到的所述离散关系拟合为线性关系。
5.根据权利要求4所述的放射源定位装置,其特征在于,所述待测放射源放射射线时,所述2×2晶体阵列中四块晶体的实测晶体计数依次为A’、B’、C’和D’,所述计数噪声的变化范围为Nmin至Nmax,所述计算模块包括:
设置子模块,用于将所述计数噪声N设置为Nmin;
第一计算子模块,用于根据AB’=(A’+B’-2×N)/(A’+B’+C’+D’-4×N)获取第一去噪晶体计数比,根据BC’=(B’+C’-2×N)/(A’+B’+C’+D’-4×N)获取第二去噪晶体计数比;
第二处理子模块,用于根据所述拟合子模块拟合的所述线性关系获取AB’对应的放射源角度θ1和BC’对应的放射源角度θ2;
自加子模块,用于将所述计数噪声N自加1;
比较子模块,用于比较|θ1-θ2|和的大小,和为预设参数,若则设置参数否则,参数和不变;
判断子模块,用于判断所述自加子模块得出的所述计数噪声N是否大于Nmax;
第二计算子模块,用于当所述计数噪声N大于Nmax时,根据 获取所述待测放射源的放射源角度。
6.根据权利要求4或5所述的放射源定位装置,其特征在于,所述拟合子模块拟合的所述线性关系为θ=(AB-y1)/x1,θ=(BC-y2)/x2,其中,x1、x2、y1和y2为直线拟合参数。
7.一种放射源定位***,其特征在于,包括:如权利要求4-6中任一所述的放射源定位装置、待测放射源和探测器。
8.根据权利要求7所述的放射源定位***,其特征在于,所述待测放射源为伽马射线放射源。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310602906.6A CN103645491B (zh) | 2013-11-25 | 2013-11-25 | 放射源定位方法、装置及*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310602906.6A CN103645491B (zh) | 2013-11-25 | 2013-11-25 | 放射源定位方法、装置及*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103645491A CN103645491A (zh) | 2014-03-19 |
CN103645491B true CN103645491B (zh) | 2016-12-07 |
Family
ID=50250739
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310602906.6A Active CN103645491B (zh) | 2013-11-25 | 2013-11-25 | 放射源定位方法、装置及*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103645491B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104155681A (zh) * | 2014-07-21 | 2014-11-19 | 北京辛耕普华医疗科技有限公司 | 用于监测放射源的定位仪及监测放射源的方法 |
CN104155675A (zh) * | 2014-08-27 | 2014-11-19 | 中国科学院高能物理研究所 | 一种射线源定位成像装置 |
CN104330814B (zh) * | 2014-10-29 | 2017-06-20 | 西南科技大学 | 一种放射源定位方法及*** |
CN104597068A (zh) * | 2015-01-08 | 2015-05-06 | 中国科学院高能物理研究所 | 一种空间主动激发x荧光谱仪的准直结构 |
CN104965216B (zh) * | 2015-06-11 | 2017-07-07 | 山东航天电子技术研究所 | 基于射线源在轨标定***的探测器探测效率标定的方法 |
CN105068109B (zh) * | 2015-07-28 | 2017-08-01 | 济南中威仪器有限公司 | 一种用于搜寻遗失放射源的定位方法 |
CN106483545B (zh) * | 2016-09-14 | 2021-12-24 | 中国科学院高能物理研究所 | 放射源监控方法及*** |
CN107064987B (zh) * | 2017-01-16 | 2019-10-22 | 北京科技大学 | 一种放射源定位***及定位方法 |
CN107436446B (zh) * | 2017-08-07 | 2020-06-19 | 西南科技大学 | 一种放射源定位方法及定位*** |
CN107884805B (zh) * | 2017-09-28 | 2019-08-23 | 苏州瑞派宁科技有限公司 | 一种射源定位的方法及装置 |
CN108445527B (zh) * | 2018-01-30 | 2020-04-14 | 张岚 | 一种射线源位置坐标的获取方法 |
CN114371494B (zh) * | 2022-03-22 | 2022-06-24 | 西南科技大学 | 面向自主寻源机器人的放射源场景模拟方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN86106216A (zh) * | 1986-08-23 | 1988-03-02 | 株式会社岛津制作所 | 光电信增管放大倍数的稳定装置 |
CN1338048A (zh) * | 1999-01-29 | 2002-02-27 | 特罗克斯勒电子实验有限公司 | 校正核辐射计用的设备和方法 |
WO2003007019A1 (en) * | 2001-07-09 | 2003-01-23 | National Institute Of Radiological Sciences Office Of Public Relations | Radiation three-dimensional position detector |
JP3795041B2 (ja) * | 2003-12-16 | 2006-07-12 | 三菱重工業株式会社 | 放射性物質の含有量測定方法及び測定装置 |
CN201237645Y (zh) * | 2008-04-30 | 2009-05-13 | 宋跃进 | 放射源在线定位监控装置 |
CN101460864A (zh) * | 2006-09-19 | 2009-06-17 | 株式会社岛津制作所 | 核医学诊断装置 |
CN101688918A (zh) * | 2007-04-23 | 2010-03-31 | 独立行政法人放射线医学综合研究所 | 利用能量和位置信息的辐射检测方法及装置 |
CN102759408A (zh) * | 2011-04-25 | 2012-10-31 | 中国科学院空间科学与应用研究中心 | 一种单光子计数成像***及其方法 |
CN102944890A (zh) * | 2012-11-06 | 2013-02-27 | 中国科学院高能物理研究所 | 基于ps-pmt的探测器信号读出方法和*** |
CN202815229U (zh) * | 2012-08-08 | 2013-03-20 | 北京辛耕普华医疗科技有限公司 | 用于伽马放射源定位仪的伽马相机及其伽马放射源定位仪 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4568818B2 (ja) * | 2005-09-29 | 2010-10-27 | 独立行政法人 日本原子力研究開発機構 | ガンマ線源を利用した可視化装置 |
JP5588190B2 (ja) * | 2010-02-08 | 2014-09-10 | 日立アロカメディカル株式会社 | 放射線検出モジュール |
-
2013
- 2013-11-25 CN CN201310602906.6A patent/CN103645491B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN86106216A (zh) * | 1986-08-23 | 1988-03-02 | 株式会社岛津制作所 | 光电信增管放大倍数的稳定装置 |
CN1338048A (zh) * | 1999-01-29 | 2002-02-27 | 特罗克斯勒电子实验有限公司 | 校正核辐射计用的设备和方法 |
WO2003007019A1 (en) * | 2001-07-09 | 2003-01-23 | National Institute Of Radiological Sciences Office Of Public Relations | Radiation three-dimensional position detector |
JP3795041B2 (ja) * | 2003-12-16 | 2006-07-12 | 三菱重工業株式会社 | 放射性物質の含有量測定方法及び測定装置 |
CN101460864A (zh) * | 2006-09-19 | 2009-06-17 | 株式会社岛津制作所 | 核医学诊断装置 |
CN101688918A (zh) * | 2007-04-23 | 2010-03-31 | 独立行政法人放射线医学综合研究所 | 利用能量和位置信息的辐射检测方法及装置 |
CN201237645Y (zh) * | 2008-04-30 | 2009-05-13 | 宋跃进 | 放射源在线定位监控装置 |
CN102759408A (zh) * | 2011-04-25 | 2012-10-31 | 中国科学院空间科学与应用研究中心 | 一种单光子计数成像***及其方法 |
CN202815229U (zh) * | 2012-08-08 | 2013-03-20 | 北京辛耕普华医疗科技有限公司 | 用于伽马放射源定位仪的伽马相机及其伽马放射源定位仪 |
CN102944890A (zh) * | 2012-11-06 | 2013-02-27 | 中国科学院高能物理研究所 | 基于ps-pmt的探测器信号读出方法和*** |
Also Published As
Publication number | Publication date |
---|---|
CN103645491A (zh) | 2014-03-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103645491B (zh) | 放射源定位方法、装置及*** | |
CN103163548B (zh) | 基于伽马相机的放射性物质探测方法及其装置和*** | |
CN104536029B (zh) | 基于多个NaI(T1)探测器的放射性物质二维定位方法 | |
US20210173096A1 (en) | Method and device for multi-dimensional direction measurement of gamma radiation in the far field | |
CN103901052B (zh) | 一种sgs与tgs联合测量装置及准直器优化方法 | |
Canadas et al. | Validation of a small-animal PET simulation using GAMOS: a GEANT4-based framework | |
CN105496436A (zh) | 用于pet装置的时间校正方法和装置 | |
Jan et al. | PSF reconstruction for Compton-based prompt gamma imaging | |
Surti et al. | Impact of event positioning algorithm on performance of a whole-body PET scanner using one-to-one coupled detectors | |
CN100399051C (zh) | 辐射探测器无源效率刻度的方法 | |
CN107195342A (zh) | 一种测量核电厂核素活度分布的方法 | |
JP5523407B2 (ja) | 放射線検出装置及び検出方法 | |
CN109581472B (zh) | Czt谱仪特征能量峰道址范围选择和确定方法 | |
Xiao et al. | A study on scattering correction for γ-photon 3D imaging test method | |
Li et al. | A SPECT system simulator built on the SolidWorksTM 3D-Design package | |
CN108445029A (zh) | 提高康普顿散射成像性能的方法 | |
Kowash | A rotating modulation imager for the orphan source search problem | |
Dong et al. | Development and preliminary results of a large-pixel two-layer labr3 compton camera prototype | |
Hetzel et al. | Near-field coded-mask technique and its potential for proton therapy monitoring | |
Shirakawa et al. | Remote sensing of nuclear accidents using a direction finding detector | |
Heath et al. | Gantryless Associated-Particle Neutron Radiography | |
Niknami et al. | Design and simulation of a semiconductor detector-based Compton imaging system with efficiency analysis | |
El Bitar et al. | The performance of a hybrid analytical-Monte Carlo system response matrix in pinhole SPECT reconstruction | |
Toivonen et al. | On the symmetry of photon detection arrays: A directionally sensitive 3D model | |
Sjölin et al. | A method for geometric calibration of edge‐on detectors in a CT‐gantry |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |