一种基于双编码板的伽马相机及其用于放射性物质定位的
方法
技术领域
背景技术
[0002] 随着核工业与核技术应用深入到国家经济发展的各个领域,放射性物质的安全监管和核事故应急能
力成为核安全和核安保行业特别关注的问题。对放射性物质空间分布的精确描述对于提升核安全监管能力至关重要,在工业及医用
放射源安全管理、核事故应急处置、环境辐射监测、公共安全等领域具有广泛应用。
[0003] 编码孔径伽马射线成像系统又称伽马相机,是一种采用编码板与阵列式
辐射探测器为核心器件的放射性物质定位装置。伽马
光子通过按照特定方式排列的编码板后,在阵列探测器上形成编码图像,编码图像经过特殊
算法进行解码可以还原待测区域的二维辐射分布情况。传统的伽马相机的空间
位置分辨能力取决于阵列探测器的
像素个数,想要提高放射性图像的空间位置分辨能力,只有提高阵列探测器的像素数目。
[0004] 在实际应用场景中,对放射性物质进行精确定位是十分重要的,但是限于目前伽马辐射阵列探测器的像素个数较少、编码板的孔径不能太小等问题,尚无法在物理测量层面实际突破探测器像素个数限制实现更加高分辨的放射性成像。
发明内容
[0005] 本发明提供了一种基于双编码板的伽马相机和放射性物质定位方法;可以获得高于阵列探测器像素个数的高
分辨率辐射图像。
[0006] 为实现以上目的,本发明采用以下技术方案:一种基于双编码板的伽马相机,包括:
数据采集处理系统1、阵列探测器2、探测器位置变换装置3、第一编码板4、屏蔽结构5、第二编码板6、
深度相机7、编码板切换装置24;数据采集处理系统1、阵列探测器2、探测器位置变换装置3都位于屏蔽结构5围成的六面体空间内,屏蔽结构5六个面中的前端一面用编码板代替,探测器位置变换装置3是指承载安装探测器的以及实现探测器位置变换的整体机构,可以看做是一
块平板通过卡扣或螺丝固定到屏蔽结构5的内侧以及数据采集处理系统1的前端,阵列探测器2通过引脚插槽固定到探测器位置变换装置3上的插针上,深度相机7利用卡扣和螺丝固定到屏蔽结构5的下侧,编码板切换装置24通过螺丝或者粘接的方式分别与屏蔽结构5的前端两侧边连接,编码板切换装置24为内置微型
电机的可以实现自动开合的
铰链结构,第一编码板4和第二编码板6分别通过编码板切换装置24与屏蔽结构5连接,数据采集处理系统1置于阵列探测器2后端,阵列探测器
2通过利用数据线连接到数据采集处理系统1,数据采集处理系统1包括数据获取模块、数据传输模块和
数据处理模块,数据采集处理系统1通过所述数据传输模块与数据处理装置与阵列探测器2相连接。
[0007] 以上所述结构中,第一编码板4为基于16阶编码阵列的31阶嵌套编码板,第二编码板6为基于17阶修正均匀冗余阵列的33阶嵌套阵列,阵列探测器2探测单元的大小与编码板的孔径单元大小相同,第一编码板4、第二编码板6和屏蔽结构5所使用的材料为铅、钨、
铜、
铁或其
合金;阵列探测器2为阵列式辐射探测器,阵列探测器2的探测材料包括有机闪烁晶体、无机闪烁晶体、
半导体或气体电离室,阵列探测器2的阵列结构为线阵式、像素式或多层阻性板式;所述的数据获取、传输、处理模块的设计和搭建采用的核心部件为ASIC芯片、FPGA芯片、PCB板、ARM板、Windows
主板、无线射频模块、GPRS数据传输模块中的一种或几种组合。
[0008] 以上所述基于双编码板的伽马相机用于放射性物质定位的方法,包括以下步骤:步骤1:首先将第一编码板和处于位置一处的阵列探测器结合成第一伽马相机,对待测区域放射性物质进行成像,所述位置一为阵列探测器所在平面的中心位置,该位置与第一编码板平行且二者中心点连线垂直于探测器;
步骤2:将所述阵列探测器移动到位置二处,将第二编码板与所述移动后的阵列探测器结合成第二伽马相机,对待测区域放射性物质进行成像,所述位置二为所述阵列探测器所在平面的左上、左下、右上、右下四个
角落中的某一个方向平移半个探测器像素大小;
步骤3:保持伽马相机方向和位置不变,利用深度相机对待测区域进行光学成像并测量待测区域与成像系统的距离;
步骤4:利用数据采集处理系统将第一伽马相机和第二伽马相机获得的放射性成像结果进行
叠加,获得高空间分辨放射性成像结果;
步骤5:利用数据采集处理系统将步骤4获得的高空间分辨放射性成像结果与步骤3获得的光学成像结果进行配准,获得能够直接指示出放射性物质空间位置的复合图像;
步骤6:利用阵列探测器测到的伽马能谱数据和深度相机测得的距离信息,计算出放射性物质的核素种类、放射性活度、
剂量率。
[0009] 以上所述步骤中,步骤1和步骤2中所述第一编码板、所述第二编码板分别与所述阵列探测器结合构成第一伽马相机和第二伽马相机时,编码板的中心点所处位置相同,均位于处于位置一处的阵列探测器的正前方,第一编码板和第二编码板分别与阵列探测器结合组成伽马相机的方式利用铰链或合页装置自动地将编码板转动至阵列探测器前方相同位置,驱动编码板转动的机构为可遥控操作的电动马达,编码板和阵列探测器的大小以及二者之间的距离还要依据伽马相机的使用场景进行设计,所需考虑参数包括成像距离范围、视场角大小、图像的空间分辨率和角分辨率要求,所述编码板的各项参数需要根据阵列探测器的参数进行设计并优化,所需考虑的设计参数包括
准直器的材料、
准直器的大小和厚度、孔的形状、孔的大小;步骤3中所述深度相机在使用前通过校准使得深度相机与第一伽马相机的
视野中心重合,确保伽马相机和深度相机对同一个区域进行成像,同时深度相机的光学成像视野角大于等于伽马相机的视野角,确保光学可成像区域大于等于伽马相机的可成像区域,深度相机为光学相机、双目或多目视觉系统、结构光视觉系统、TOF(Time of Flight)视觉系统中的一种或几种;
步骤4中所述将第一伽马相机和第二伽马相机获得的放射性成像结果进行叠加时,将第二伽马相机获得的17阶放射性图像与第一伽马相机获得的16阶图像的每一像素分割成四份并按照中心点对齐,直接进行数值相加,获得一幅新的放射性图像;
步骤5中的
图像配准采用的方法为基于灰度信息的方法、基于空间变换域的方法或基于图像特征的方法,图像配准所采用的计算平台为OpenCV
计算机视觉库、ArcGIS或MATLAB中;
步骤6中利用阵列探测器测得的伽马能谱数据分析
放射性核素种类的方法包括寻峰匹配法、对称零面积法、
人工神经网络方法或
机器学习方法;通过阵列探测器测得的剂量率和距离信息反推放射性核素活度和种类的方法为利用核素识别结果判断伽马射线
能量,利用特征伽马射线计数和距离平方反比关系反推放射性物质活度和一米远处空气剂量率;
以上基于双编码板的伽马相机用于放射性物质定位方法中还利用了用于远程控制伽马相机和深度相机的上位机下位机系统、
无线通信模块、GPS定位系统;上位机系统指的是操作人员远程操控时使用的软
硬件系统,用于远程控制伽马相机和深度相机进行测量,控制探测器位置变换装置和编码板切换装置。下位机系统指的是安装在伽马相机上的用于接收上位机系统发出的控制指令的软硬件系统。无线通信模块用于
信号传输和指令发送,GPS定位系统用于实时定位伽马相机地理位置,安装在伽马相机下位机系统上。
[0010] 有益效果:本发明提供了一种基于双编码板的伽马相机和放射性物质定位方法,采用双编码板设计,利用一个阵列探测器和两次成像方法,可以更加精确地判断放射性物质的方向和位置;采用双编码板自动切换方法可以在短时间内遥控操作,提高效率、方便使用,减少了工作人员的受辐照
风险;而且本发明采用光学成像和伽马射线成像技术共同使用的方法,在获得放射性成像结果之后还可以通过与光学图像的叠加准确判断放射性物质所处场景位置;本发明将深度检测方法引入伽马相机系统中,利用测距信息、伽马能谱测量结果可以在成像的同时识别放射性位置的核素种类、放射性活度和一米远处空气剂量率,为处置工作提供参考,本发明的伽马相机既可以人工操作一台单独的伽马相机,也可以搭载在无人机、无人
汽车、
机器人、无人艇和便携式移动平台,
载荷设备的控制系统与上下位机系统通过
信号传输模块连接,用以为载荷设备进行移动轨迹规划,实现一种多功能多用途的放射性位置定位方法。
附图说明
[0011] 图1是本发明伽马相机的结构示意图;图2是本发明实施中放射性物质定位方法的
流程图;
图3是本发明中第一编码板和第二编码板的孔径排列示意图,其中,301为第一编码板孔径排列示意图,302为第二编码板孔径排列示意图;
图4是本发明
实施例中将阵列探测器从位置一移动到位置二的结构装置示意图;
图5是本发明实施例中第一编码板和阵列探测器构成第一伽马相机并对中心位置处放射源进行成像的示意图;
图6是本发明实施例中第二编码板和阵列探测器构成第二伽马相机并对中心位置处放射源进行成像的示意图;
图7是本发明实施例中对第一伽马相机、第二伽马相机的成像结果进行叠加获得可以更高
精度显示放射源位置的辐射图像;
图中,1为数据采集系统,2为阵列探测器,3为探测器位置变换装置,4为第一编码板,5为屏蔽结构,6为第二编码板,7为深度相机,8为搭载阵列探测器平台,9为电
磁铁,10 为金属轴,11为阵列探测器中间4×4个像素a,12为第一编码板中间3×3个单元,13为放射源位置a,14为阵列探测器中间4×4个像素b,15为第二编码板中间3×3个单元,16为17阶阵列探测器,17为放射源位置b,18为第二伽马相机获得的17阶辐射热点图像,19为第二伽马相机获得的17阶辐射热点图像右上角16×16个像素,20为第一伽马相机获得的16阶辐射热点图像,21为第二伽马相机获得的17阶辐射热点图像右上角16×16个像素中间4×4个像素,22为第一伽马相机获得的16阶辐射热点图像中间4×4个像素,23为放射源在辐射热点图象上的对应位置,24为编码板切换装置。
具体实施方式
[0012] 下面结合附图和具体实施例对本发明进行详细说明:如图1所示,包括:数据采集处理系统1、阵列探测器2、探测器位置变换装置3、第一编码板4、屏蔽结构5、第二编码板6、深度相机7、编码板切换装置24;数据采集处理系统1、阵列探测器2、探测器位置变换装置3都位于屏蔽结构5围成的六面体空间内,屏蔽结构5六个面中的前端一面用编码板代替,探测器位置变换装置3是指承载安装探测器的以及实现探测器位置变换的整体机构,可以看做是一块平板通过卡扣或螺丝固定到屏蔽结构5的内侧以及数据采集处理系统1的前端,阵列探测器2通过引脚插槽固定到探测器位置变换装置3上的插针上,深度相机7利用卡扣和螺丝固定到屏蔽结构5的下侧,编码板切换装置24通过螺丝或者粘接的方式分别与屏蔽结构5的前端两侧边连接,编码板切换装置24为内置微型电机的可以实现自动开合的铰链结构,第一编码板4和第二编码板6分别通过编码板切换装置24与屏蔽结构5连接,数据采集处理系统1置于阵列探测器2后端,阵列探测器2通过利用数据线连接到数据采集处理系统1,数据采集处理系统1包括数据获取模块、数据传输模块和数据处理模块,数据采集处理系统1通过所述数据传输模块与数据处理装置与阵列探测器2相连接。
[0013] 以上所述结构中,第一编码板4为基于16阶编码阵列的31阶嵌套编码板,第二编码板6为基于17阶修正均匀冗余阵列的33阶嵌套阵列,阵列探测器2探测单元的大小与编码板的孔径单元大小相同,第一编码板4、第二编码板6和屏蔽结构5所使用的材料为铅、钨、铜、铁或其合金;阵列探测器2为阵列式辐射探测器,阵列探测器2的探测材料包括有机闪烁晶体、无机闪烁晶体、半导体或气体电离室,阵列探测器2的阵列结构为线阵式、像素式或多层阻性板式;所述的数据获取、传输、处理模块的设计和搭建采用的核心部件为ASIC芯片、FPGA芯片、PCB板、ARM板、Windows主板、无线射频模块、GPRS数据传输模块中的一种或几种组合。
[0014] 如图2所示,以上所述基于双编码板的伽马相机用于放射性物质定位的方法,包括以下步骤:步骤201:首先将第一编码板和处于位置一处的阵列探测器结合成第一伽马相机,对待测区域放射性物质进行成像,所述位置一为阵列探测器所在平面的中心位置,该位置与第一编码板平行且二者中心点连线垂直于探测器;
步骤202:将所述阵列探测器移动到位置二处,将第二编码板与所述移动后的阵列探测器结合成第二伽马相机,对待测区域放射性物质进行成像,所述位置二为所述阵列探测器所在平面的左上、左下、右上、右下四个角落中的某一个方向平移半个探测器像素大小;
步骤203:保持伽马相机方向和位置不变,利用深度相机对待测区域进行光学成像并测量待测区域与成像系统的距离;
步骤204:利用数据采集处理系统将第一伽马相机和第二伽马相机获得的放射性成像结果进行叠加,获得高空间分辨放射性成像结果;
步骤205:利用数据采集处理系统将步骤4获得的高空间分辨放射性成像结果与步骤3获得的光学成像结果进行配准,获得能够直接指示出放射性物质空间位置的复合图像;
步骤206:利用阵列探测器测到的伽马能谱数据和深度相机测得的距离信息,计算出放射性物质的核素种类、放射性活度、剂量率。
[0015] 以上所述步骤中,步骤1和步骤2中所述第一编码板、所述第二编码板分别与所述阵列探测器结合构成第一伽马相机和第二伽马相机时,编码板的中心点所处位置相同,均位于处于位置一处的阵列探测器的正前方,第一编码板和第二编码板分别与阵列探测器结合组成伽马相机的方式利用铰链或合页装置自动地将编码板转动至阵列探测器前方相同位置,驱动编码板转动的机构为可遥控操作的电动马达,编码板和阵列探测器的大小以及二者之间的距离还要依据伽马相机的使用场景进行设计,所需考虑参数包括成像距离范围、视场角大小、图像的空间分辨率和角分辨率要求,所述编码板的各项参数需要根据阵列探测器的参数进行设计并优化,所需考虑的设计参数包括准直器的材料、准直器的大小和厚度、孔的形状、孔的大小;步骤3中所述深度相机在使用前通过校准使得深度相机与第一伽马相机的视野中心重合,确保伽马相机和深度相机对同一个区域进行成像,同时深度相机的光学成像视野角大于等于伽马相机的视野角,确保光学可成像区域大于等于伽马相机的可成像区域,深度相机为光学相机、双目或多目视觉系统、结构光视觉系统、TOF(Time of Flight)视觉系统中的一种或几种;
步骤4中所述将第一伽马相机和第二伽马相机获得的放射性成像结果进行叠加时,将第二伽马相机获得的17阶放射性图像与第一伽马相机获得的16阶图像的每一像素分割成四份并按照中心点对齐,直接进行数值相加,获得一幅新的放射性图像;
步骤5中的图像配准采用的方法为基于灰度信息的方法、基于空间变换域的方法或基于图像特征的方法,图像配准所采用的计算平台为OpenCV计算机视觉库、ArcGIS或MATLAB中;
步骤6中利用阵列探测器测得的伽马能谱数据分析放射性核素种类的方法包括寻峰匹配法、对称零面积法、人工神经网络方法或机器学习方法;通过阵列探测器测得的剂量率和距离信息反推放射性核素活度和种类的方法为利用核素识别结果判断伽马射线能量,利用特征伽马射线计数和距离平方反比关系反推放射性物质活度和一米远处空气剂量率。
[0016] 以上基于双编码板的伽马相机用于放射性物质定位方法中还利用了用于远程控制伽马相机和深度相机的上位机下位机系统、无线通信模块、GPS定位系统;上位机系统指的是操作人员远程操控时使用的软硬件系统,用于远程控制伽马相机和深度相机进行测量,控制探测器位置变换装置和编码板切换装置。下位机系统指的是安装在伽马相机上的用于接收上位机系统发出的控制指令的软硬件系统。无线通信模块用于信号传输和指令发送,GPS定位系统用于实时定位伽马相机地理位置,安装在伽马相机下位机系统上。
[0017] 所述第一编码板为基于16阶编码阵列的31阶嵌套编码板,其孔径排列图案如图3中的301所示;所述第二编码板为基于17阶修正均匀冗余阵列的33阶嵌套阵列,其孔径排列图案如图3中的302所示,其中黑色部分表示射线无法直接通过的屏蔽材料,白色部分表示射线可以直接通过的小孔,所述阵列探测器的行列数为16,阵列探测器探测单元的大小与编码板的孔径单元大小相同,均为1 cm。
[0018] 图4是本发明实施例中将阵列探测器从位置一处移动到位置二处的结构装置示意图,其中,搭载阵列探测器平台8前部的开槽为16.5 cm,阵列探测器2长宽为16 cm;9为平台上的开孔并安装在孔周围的电磁铁,10为与探测器相连接的金属轴,当阵列探测器2位于搭载阵列探测器平台8前部开槽的左下方时,金属轴10被磁铁9的右下方吸引过去,此时阵列探测器位于步骤201中所述位置一;当金属轴10被磁铁9的左上方吸引过去时,阵列探测器2随之移动到阵列探测器平台8前部开槽的右上方,即阵列探测器向右上方移动了半个像素大小,此时阵列探测器位于步骤202所述的位置二处。
[0019] 所述第一、二编码板与屏蔽结构通过电动合页连接,通过外部上位机
软件控制编码板的移动,使之分别与阵列探测器组成第一、二伽马相机。
[0020] 使用深度相机获得探测目标区域的
深度图和光学图像,在本实施例中,使用商业化深度检测视觉产品Microsoft Kinect,可以通过配套软件系统直接获得探测目标区域的深度图和光学图像,步骤203中所述的光学相机在使用前已经通过校准使得光学成像视野中心与第一伽马相机的视野中心重合,确保对同一个区域进行成像,同时通过结构设计已使得光学成像的视野角大于等于伽马相机的视野角,确保光学可成像区域大于等于伽马相机的可成像区域。
[0021] 图5是本实施例中第一编码板和阵列探测器构成第一伽马相机并对中心位置处放射源进行成像的示意图,其中4是第一编码板,2阵列探测器,11指示的实线方格代表阵列探测器2的中间4×4个像素,12指示的虚线方格代表第一编码板4的中间3×3个单元,13指示的两个黑色圆点代表放射源的位置,该图为从放射源所在平面正视第一伽马相机的结构示意图,可以看出编码板的单元与阵列探测器的像素不是一一对齐的,图中的浅灰色区域代表左下角的放射源在实际成像结果中反映出来的辐射热点区域,深灰色区域代表右上角的放射源在实际成像结果中反映出来的辐射热点区域。
[0022] 图6是本实施例中第二编码板和阵列探测器构成第二伽马相机并对中心位置处放射源进行成像的示意图,其中6是第二编码板,第二编码板位置与图5中所示第一编码板位置相同,2是阵列探测器(16阶),16指示的续断线方框是16阶阵列探测器加上左侧一列及下侧一行虚构出的17阶阵列探测器,在实际操作中左侧一列和下侧一行数据被赋值为其他所有像素值的平均值,14指示的实线方格代表阵列探测器2的中间4×4个像素,15指示的虚线方格代表第二编码板的中间3×3个单元,17指示的两个黑色圆点代表放射源的位置,放射源位置与图5中所示放射源位置13相同,该图为从放射源所在平面正视第二伽马相机的结构示意图,可以看出编码板的单元与阵列探测器的像素是一一对齐的,即14指示的实线方格与15指示的虚线方格是重合的,图中的浅灰色区域代表右上角的放射源在实际成像结果中反映出来的辐射热点区域,深灰色区域代表左下角角的放射源在实际成像结果中反映出来的辐射热点区域。
[0023] 图7是本实施例中对第一、二伽马相机的成像结果进行叠加获得可以更高精度显示放射源位置的辐射图像,18指示的续断线方框为第二伽马相机获得的17阶辐射热点图像,19指示的细实线方框为第二伽马相机获得的17阶辐射热点图像的右上角16×16像素,20指示的粗实线方框代表第一伽马相机获得的16阶辐射热点图像,21指示的虚线方格是第二伽马相机获得的17阶辐射热点图像的右上角16×16像素中间4×4个像素,22指示的实线方格是第一伽马相机获得的16阶辐射热点图像中间4×4个像素,23指示的是放射源在辐射热点图像上的对应位置;图中灰色部分为图5和图6中辐射热点成像结果的叠加,在进行像素值叠加前,每个像素都被分割成了四份,然后按照图7所示位置将对应像素依次求和,获得叠加图像结果,可以看出通过对第一伽马相机和第二伽马相机成像结果的叠加处理,可以更加清晰地指示出放射源实际所在位置,将热点范围缩小在放射源所在像素,图像分辨率也增大了一倍。
[0024] 阵列探测器和编码板的周围使用伽马屏蔽材料作为
外壳,采用刚性连接,数据采集处理系统置于阵列探测器后端,通过信号传输模块与
图像处理装置相连接。
[0025] 在本实施例中,使用深度检测视觉系统获得探测目标区域的深度图和光学图像。在本实施例中,使用商业化深度检测视觉产品Microsoft Kinect,可以通过配套软件系统直接获得探测目标区域的深度图和光学图像。
[0026] 将辐射热点分布图像、深度图和光学图像进行图像配准,得到可以精确反映放射性空间分布和周围环境特征的复合图像。
[0027] 步骤5采用基于空间变换域的方法,结合OpenCV计算机视觉库和MATLAB软件,通过一系列实验设置确定配准模型,然后将辐射热点分布图像、深度图和光学图像在配准模型中进行配准。
[0028] 以上仅仅是本发明优选实施方式,本发明所属技术领域的技术人员可以对所描述的具体实施例做各种
修改或补充或采用类似的方式替代,只要不偏离本发明
说明书的内容或者超越本
权利要求书所定义的范围,均应属于本发明的保护范围。