一种晶体条与伽马光子作用位置映射表生成方法及系统技术方案

本发明专利技术提供了一种晶体条与伽马光子作用位置映射表生成方法及系统，以点源或泛源持续映射晶体阵列，对PMT接收信号进行计算得到反映符合事件作用位置的散点图，经数字化处理去除二维位置谱中的奇异点后，对二维位置谱采用局部最大值算法估算每个晶体条的中心位置，及对数据进行二值化处理，得到二值化位置图谱；对二值化位置图谱计算每一个像素点距离最近亮点的欧几里得距离，得到距离分布图；沿X方向和Y方向分别对距离分布图求局部最大值并划分各个晶体条的边界位置，这样，每个区域与晶体条阵列一一对应，进而确定出系统的位置映射表，该方法可覆盖所有的位置，不会出现划分死区。区。区。

【技术实现步骤摘要】
一种晶体条与伽马光子作用位置映射表生成方法及系统


[0001]本专利技术属于正电子发射断层扫描
，特别涉及一种晶体条与伽马光子作用位置映射表生成方法及系统。

技术介绍

[0002]正电子发射断层扫描技术是一种新颖的三维核医学成像技术，属于发射型断层扫描技术的一种，它利用正电子核素标记葡萄糖等人体代谢物作为显像剂，通过病灶对显像剂的摄取来反映其代谢变化，提供疾病的生物代谢信息。它可以更早期、灵敏、准确的准断和指导治疗多种疾病，在生物学、神经科学、肿瘤学、药物动力学、临床诊断和治疗评价等方面具有重要的应用价值。其基本原理是将微量的正电子核素示踪剂注入人体中，正电子与人体内的负电子湮灭产生一对背对背发射的γ光子，用体外探测装置探测这一对光子来探知放射性药物在体内的分布。目前常用的探测器系统是由晶体条组成的晶体阵列与光电倍增管耦合的方式实现的。当正负电子湮灭产生的γ光子进入晶体条，与晶体发生反应，晶体原子退激时会产生荧光，荧光沿着晶体阵列传播然后被光电倍增管转换为电信号。该电信号携带了入射γ光子的位置信息，利用定位算法得到入射γ光子发生反应的晶体条的位置信息，从而确定湮灭事件符合线的位置。位置信息的精确度将直接影响探测器的空间分辨率，进而影响PET的成像质量和图像定量分析精度。所以，γ光子发生能量沉积的晶体条的准确定位至关重要。
[0003]但是实际上，因闪烁晶体排列方式、分光光导的切割方法、光电倍增管信号处理能力的差异，及前端电子信号采集线路等多方面原因，经定位算法计算之后的位置信息容易发生畸变，因此必须通过相应的位置映射算法来校正。

技术实现思路

[0004]针对现有技术的不足，本专利技术旨在提供一种晶体条与伽马光子作用位置映射表生成方法，包括以下步骤：步骤一：对入射光子事件的位置信息进行二维统计，得到二维位置谱，所述二维位置谱中每个像素点的亮度代表该位置的事例数；步骤二：利用3*3的二维模板对获取的二维位置谱进行中值滤波处理；步骤三：利用1*3和3*1模板对二维位置谱分别沿X方向和Y方向做均值滤波处理；步骤四：对二维位置谱采用局部最大值算法估算每个晶体条的中心位置，及对数据进行二值化处理，得到二值化位置图谱；步骤五：对二值化位置图谱计算每一个像素点距离最近亮点的欧几里得距离，得到距离分布图；步骤六：沿X方向和Y方向分别对距离分布图求局部最大值并划分各个晶体条的边界位置；步骤七：采用区域标识算法，将每个区域与晶体条编号一一对应，完成散点图与晶体条的映射关系；及步骤八：采用边界提取位置映射算法，生成位置映射表。
[0005]根据本专利技术上述的一种晶体条与伽马光子作用位置映射表生成方法，进一步地，上述步骤四包括：对二维位置谱分别在X方向和Y方向的进行一维寻峰；将X方向和Y方向的寻峰结果在同一像素点作“与运算”，如果该像素点同时是X方向和Y方向的峰值点，那么将
该像素点值设为1,否则设为0，完成了二值化处理。
[0006]根据本专利技术上述的一种晶体条与伽马光子作用位置映射表生成方法，进一步地，上述步骤六包括：采用微商寻峰法对距离分布图沿X方向和Y方向分别进行一维寻峰；将X方向和Y方向的寻峰结果，在同一像素点作“或运算”，数值最大对应各个晶体条边界位置，得到边界位置分布图。
[0007]根据本专利技术上述的一种晶体条与伽马光子作用位置映射表生成方法，进一步地，上述步骤七包括：区域标识、边界归属和顺序编码。
[0008]根据本专利技术上述的一种晶体条与伽马光子作用位置映射表生成方法，进一步地，区域标识包括：将边界位置分布图进行区域标识，形成的区域标识并不是按照顺序而是随机编码的，只是计算出哪些像素点属于同一个区域，属于同一区域的像素点采用相同编码，并且每个区域间编码保证不重复。
[0009]根据本专利技术上述的一种晶体条与伽马光子作用位置映射表生成方法，进一步地，边界归属包括：判断边界位置对应的像素点的区域归属时，根据对距离分布图求局部最大值中采用一维寻峰算法特点，按照以下规则进行边界归属判断：沿X方向的边界点属于边界右边的区域，沿Y方向的边界点属于边界下方的区域。
[0010]根据本专利技术上述的一种晶体条与伽马光子作用位置映射表生成方法，进一步地，顺序编码包括采用顺序编码算法，使得每个区域与晶体表编码顺序对应，完成散点图与晶体条的映射关系。
[0011]本专利技术旨还提供了一种晶体条与伽马光子作用位置映射表生成系统，包括阵列晶体、数据统计探测器、中值滤波处理器、均值滤波处理器、晶体条的中心位置估算模块、距离计算模块、晶体条边界确定模块、区域标识模块、边界归属模块、区域顺序编码器、位置映射表生成模块和位置映射表输出模块，阵列晶体具有多个晶体条，数据统计探测器与阵列晶体耦合，数据统计探测器、中值滤波处理器、均值滤波处理器、晶体条的中心位置估算模块、距离计算模块、晶体条边界确定模块、区域标识模块、边界归属模块、区域顺序编码器、位置映射表生成模块和位置映射表输出模块依次连接。
[0012]根据本专利技术上述的一种晶体条与伽马光子作用位置映射表生成系统，进一步地，数据统计探测器为PET探测器。
[0013]根据本专利技术上述的一种晶体条与伽马光子作用位置映射表生成系统，进一步地，中值滤波处理器和均值滤波处理器均为高斯滤波处理器。
[0014]根据本专利技术上述的一种晶体条与伽马光子作用位置映射表生成系统，进一步地，中值滤波处理器为3*3的二维模板的中值滤波处理器。
[0015]根据本专利技术上述的一种晶体条与伽马光子作用位置映射表生成系统，进一步地，均值滤波处理器包括X方向均值滤波处理器和Y方向均值滤波处理器，X方向均值滤波处理器和Y方向均值滤波处理器均与中值滤波处理器和晶体条的中心位置估算模块连接，X方向均值滤波处理器为1*3模板的均值滤波处理器，Y方向均值滤波处理器为3*1模板的均值滤波处理器。
[0016]根据本专利技术上述的一种晶体条与伽马光子作用位置映射表生成系统，进一步地，晶体条的中心位置估算模块具有“与运算”处理器和二值化处理器，所述“与运算”处理器和所述二值化处理器串联集成至所述晶体条的中心位置估算模块。
[0017]根据本专利技术上述的一种晶体条与伽马光子作用位置映射表生成系统，进一步地，距离计算模块为欧几里得距离计算器。
[0018]根据本专利技术上述的一种晶体条与伽马光子作用位置映射表生成系统，进一步地，晶体条边界确定模块包括X方向一维寻峰单元、Y方向一维寻峰单元和“或运算”处理器，X方向一维寻峰单元、所述Y方向一维寻峰单元和“或运算”处理器集成至晶体条边界确定模块，X方向一维寻峰单元和Y方向一维寻峰单元两端均分别连接至距离计算模块和“或运算”处理器，“或运算”处理器另一端与区域标识模块连接。
[0019]根据本专利技术上述的一种晶体条与伽马光子作用位置映射表生成系统，进一步地，区域标识模块的编码器为随机编码器。
附图说明
[0020]上述说明仅是本专利技术技术方案的概述，为了能够更清楚了解本专利技术的技术手段，而可依照说明书的内容予以实施，并且本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种晶体条与伽马光子作用位置映射表生成方法，其特征在于：包括以下步骤：步骤一：对入射光子事件的位置信息进行二维统计，得到二维位置谱，所述二维位置谱中每个像素点的亮度代表该位置的事例数；步骤二：利用3*3的二维模板对获取的二维位置谱进行中值滤波处理；步骤三：利用1*3和3*1模板对二维位置谱分别沿X方向和Y方向做均值滤波处理；步骤四：对二维位置谱采用局部最大值算法估算每个晶体条的中心位置，及对数据进行二值化处理，得到二值化位置图谱；步骤五：对二值化位置图谱计算每一个像素点距离最近亮点的欧几里得距离，得到距离分布图；步骤六：沿X方向和Y方向分别对距离分布图求局部最大值并划分各个晶体条的边界位置；步骤七：采用区域标识算法，将每个区域与晶体条编号一一对应，完成散点图与晶体条的映射关系；及步骤八：采用边界提取位置映射算法，生成位置映射表。2.如权利要求1所述的一种晶体条与伽马光子作用位置映射表生成方法，其特征在于：所述步骤四包括：对二维位置谱分别在X方向和Y方向的进行一维寻峰；将X方向和Y方向的寻峰结果在同一像素点作“与运算”，如果该像素点同时是X方向和Y方向的峰值点，那么将该像素点值设为1,否则设为0，完成了二值化处理。3.如权利要求1所述的一种晶体条与伽马光子作用位置映射表生成方法，其特征在于：所述步骤六包括：采用微商寻峰法对距离分布图沿X方向和Y方向分别进行一维寻峰；将X方向和Y方向的寻峰结果，在同一像素点作“或运算”，数值最大对应各个晶体条边界位置，得到边界位置分布图。4.如权利要求1所述的一种晶体条与伽马光子作用位置映射表生成方法，其特征在于：所述步骤七包括：区域标识、边界归属和顺序编码。5.如权利要求4所述的一种晶体条与伽马光子作用位置映射表生成方法，其特征在于：所述区域标识包括：将边界位置分布图进行区域标识，形成的区域标识并不是按照顺序而是随机编码的，只是计算出哪些像素点属于同一个区域，属于同一区域的像素点采用相同编码，并且每个区域间编码保证不重复。6.如权利要求4所述的一种晶体条与伽马光子作用位置映射表生成方法，其特征在于：所述边界归属包括：判断边界位置对应的像素点的区域归属时，根据对距离分布图求局部最大值中采用一维寻峰算法特点，按照以下规则进行边界归属判断：沿X方向的边界点属于边界右边的区域，沿Y方向的边界点属于边界下方的区域。7.如权利要求4所述的一种晶体条与伽马光子作用位置映射表生成方法，其特征在于：所述顺序编码包括采用顺序编码算法，使得每个区域与晶体表编码顺序对应，完成散点图与晶体条的映射关系。8.一种晶体条与伽马光子作用位...

【专利技术属性】
技术研发人员：范胜男，刘小平，李景涛，严李李，
申请(专利权)人：北京大基康明医疗设备有限公司，
类型：发明
国别省市：