一种构造虚拟DOI及相应系统矩阵提高PET图像重建质量的方法技术方案

本发明专利技术公开一种构造虚拟DOI及相应系统矩阵提高PET图像重建质量的方法，包括如下步骤：1)虚拟DOI的划分：将PET探测器环中的闪烁晶体划分成若干个区域并对应不同的作用深度2)计算虚拟DOI的探测概率：将虚拟划分出的DOI区域进行配对，计算每一对在不同信号源位置被检测到的概率3)划分虚拟DOI中的虚拟LOR事件：将原先一对晶体中的一条LOR事件在虚拟DOI中划分为多个子事件4)对常规系统矩阵的优化：考虑了多个影响图像重建质量的因素5)图像重建：利用统计迭代算法重建图像6)利用GPU平台进行分布运算：可大大降低数据处理的时间。通过本方法可极大提高PET类设备图像重建后的分辨率，精确定位肿瘤位置，具有非常重要的临床价值与意义。

【技术实现步骤摘要】
一种构造虚拟DOI及相应系统矩阵提高PET图像重建质量的方法
本专利技术涉及核医学PET图像重建领域，尤其涉及一种基于构造虚拟DOI及相应系统矩阵的PET图像重建与优化方法。
技术介绍
PET(positronemissiontomography,正电子发射计算机断层扫描)是核医学领域重要的临床检查影像技术。PET的探测原理是将能发射正电子的核素作为特异性示踪剂注入病人体内，示踪剂发射的正电子与人体内的电子湮没产生一对方向相反的γ光子，这对光子被探测装置记录并储存下来，通过图像重建就可以得到示踪剂在病人体内的分布情况，从而实现肿瘤的早期发现和定位。PET设备中的探测装置一般由若干个闪烁晶体与光子探测器组成的近似环形结构组成，可以探测各种角度发射的光子对。具体的工作原理是：γ光子在闪烁晶体中衰减一段距离，被转换为可见光信号，然后被光子探测器捕获，由光电倍增管转换为电信号记录，探测器每记录下的一次信号被称为一次LOR事件(Lineofresponse)。对信号传统的处理方式是：忽略探测器的尺寸，将一对探测器(互成180°)捕捉到的LOR事件归纳为同一类事件，即光子的发射位置处于这一对探测器的中心连线某点上。这种方法可以大大提高探测器的计数效率，并且计算简单方便，但同时也忽略了光子在闪烁晶体内被探测到时的位置深度，会造成计算光子发射位置时存在很大误差，造成图像的空间分辨率下降。而PET图像的分辨率过低也正是PET技术发展中的一大诟病。如图1所示：传统PET成像设备将某一次由源A所产生的LOR事件，归纳为由这对探测器P1P2的中心连线上某点产生的，计算出源处于P1P2连线上某点，而实际上，该次LOR事件的源处于A点，故经过图像重建后会造成肿瘤位置的偏差。为了解决这一问题，可以在物理硬件上重新设计制造出能检测光子在晶体中的作用深度DOI(depthofinteraction)的探测器，但物理方法有几个致命的缺点：一是目前探测器的制作工艺困难，太薄难以捕获光子，太厚又无意义；二是会使探测器更复杂导致成本更高；三是很多已有的PET设备不具备升级硬件的条件，不仅需要重新匹配PET系列中的其他模块，还需要重新测试系统性能。
技术实现思路
本专利技术所要解决的技术问题是提供一种构造虚拟DOI及相应系统矩阵来进行PET图像重建的新方法。首先对晶体虚拟创建出DOI，再对每一对DOI检测到γ光子的概率进行计算，增大了检测到PET数据(LOR事件)的几何范围；并根据虚拟DOI构造出相应的系统矩阵用以后续的PET图像重建迭代算法，这个系统矩阵将不同于以往的简便计算方法，考虑了更多的影响图像质量的因素，计算精度将大大提高；最后在云计算和图形处理单元(GPU)平台上，通过图像重建迭代算法对数据进行分布式图像重建计算，降低了计算时间。具体的技术细节如下：步骤1，虚拟DOI的划分，将PET探测器环中的闪烁晶体划分成若干个区域并对应不同的作用深度；步骤2，计算虚拟DOI的探测概率，将虚拟划分出的DOI区域进行配对，计算每一对在不同信号源位置被检测到的概率；步骤3，划分虚拟DOI中的虚拟LOR事件，将原先一对晶体中的一条LOR事件在虚拟DOI中划分为多个子事件，并计算每个虚拟子事件的发生概率，将其作为后续重建图像时每条LOR子事件的权重系数；步骤4，对PET重建中的常规系统矩阵进行优化，将步骤2计算获得的虚拟DOI对应的探测概率加入到系统矩阵；步骤5，利用步骤4获得的系统矩阵对图像进行重建，获得优化后的图像。进一步的，步骤1中根据闪烁晶体的几何排列和尺寸大小进行虚拟DOI的划分。进一步的，步骤2中计算虚拟DOI的探测概率的实现方式如下，通过连接信号源与某对虚拟DOI中心的连线，计算覆盖在晶体中的长度，确定检测概率D，对于一个光子，它在晶体内被检测概率D与其在之前相邻晶体中穿透过的距离A和最终某个晶体内部穿越的距离L相关：D＝(1-e-μL)e-μA，其中μ为γ光子在某种晶体中的衰减系数，为常数；假设源M处发生的一次LOR事件，由一对探测器P1和P2探测到，γ光子在探测器P1中先穿透距离A1，最后才在某个晶体内穿越距离L1后被捕获，将A1和L1带入公式可以计算出光子在晶体P1内被检测到的概率D1，同样的，计算出光子在晶体P2内被检测到的概率D2，则光子在这对探测器中被检测到的概率D＝D1D2。进一步的，步骤4中优化后的系统矩阵的计算公式如下，Pij＝Pdet·blu，P，eom其中Pdet.blur代表光子在晶体不同区域检测到的概率，由步骤2获得；Pgeom代表光子通过像素点的几何概率，由像素j到光子投影线i的几何距离的远近决定。进一步的，步骤5中利用图像重建迭代算法OSEM对图像进行重建。进一步的，步骤4和步骤5在GPU平台上进行分布运算。与现有技术相比，本专利技术的优点和有益效果如下：1)能有效的提高重建后的PET图像的分辨率2)极大降低了制造成本，本专利技术通过数学方法构造虚拟DOI，在提供相当性能的同时，能极大降低设备的制造成本3)本专利技术所示的方法流程简单方便，适用范围广，可以轻松应用在各种PET设备上，均能提高相应的图像质量。适用范围：构造虚拟DOI所适用的PET设备，不局限于设备的结构，可以是全身，局部器官或小动物的设备；不局限于晶体与探测器的几何分布，可以是环形，半环形，平板或多边形分布。附图说明图1为传统PET探测器处理LOR事件忽略DOI产生误差示意图；图2为本专利技术实施例中PET探测器的虚拟DOI划分示意图；图3为本专利技术实施例中光子在闪烁晶体不同位置(DOI)被检测到的概率计算示意图；图4为本专利技术具体实施的流程图；具体实施方式下面结合附图对本
技术实现思路
作进一步介绍，本专利技术实施例的具体实现步骤如下：步骤1，虚拟DOI的划分：将PET探测器环中的闪烁晶体划分成若干个区域并对应不同的作用深度(DOI)，以增大探测器探测到的几何范围。区别于可以检测到实际DOI的物理探测器，这些区域是虚拟的。根据探测器的角度和覆盖的几何范围，虚拟DOI的个数可以自由调整，可由PET探测器的几何分布与晶体大小决定划分的数量。例如探测器本身所覆盖的几何范围很小，可以划分很少的虚拟DOI。而对于边缘部分的晶体(需要更精确的探测更大范围内的γ光子)，可以大量划分虚拟DOI。如图2所示，将一对互成180°的探测器P1,P2虚拟划分为P1,1；P1,2；P2,1；P2,2这4个子探测器，对原先探测器仅能计算出源分布在P1P2连线上，通过划分虚拟DOI后计算出源分布可以在P1,1P2,1；P1,1P2,2；P1,2P2,1；P1,2P2,2这4条连线上，分辨率较原先有了极大提高。优化虚拟DOI的划分：不同于物理硬件上根据晶体深度方向划分DOI，虚拟DOI的划分不局限在晶体的深度方向上，可以取决于不同的几何形态，为了模拟所有γ光子可能被检测到的几何位置，虚拟DOI可以划分在晶体内的其他方向，例如细分在晶体长度与宽度方向上，如果晶体非长方体，也可根据不同几何形态划分，这些是物理探测器达不到的。步骤2，计算虚拟DOI的探测概率：将虚拟划分出的DOI区域进行配对，计算每一对在不同信号源位置被检测到的概率。以图2所示为例，即将一对探测器P1P2实际所采集到的所有LOR事件，分解到P1本文档来自技高网...

【技术保护点】
1.一种构造虚拟DOI及相应系统矩阵提高PET图像重建质量的方法，其特征在于，包括如下步骤：步骤1，虚拟DOI的划分，将PET探测器环中的闪烁晶体划分成若干个区域并对应不同的作用深度；步骤2，计算虚拟DOI的探测概率，将虚拟划分出的DOI区域进行配对，计算每一对在不同信号源位置被检测到的概率；步骤3，划分虚拟DOI中的虚拟LOR事件，将原先一对晶体中的一条LOR事件在虚拟DOI中划分为多个子事件，并计算每个虚拟子事件的发生概率，将其作为后续重建图像时每条LOR子事件的权重系数；步骤4，对PET重建中的常规系统矩阵进行优化，将步骤2计算获得的虚拟DOI对应的探测概率加入到系统矩阵；步骤5，利用步骤4获得的系统矩阵对图像进行重建，获得优化后的图像。

【技术特征摘要】
1.一种构造虚拟DOI及相应系统矩阵提高PET图像重建质量的方法，其特征在于，包括如下步骤：步骤1，虚拟DOI的划分，将PET探测器环中的闪烁晶体划分成若干个区域并对应不同的作用深度；步骤2，计算虚拟DOI的探测概率，将虚拟划分出的DOI区域进行配对，计算每一对在不同信号源位置被检测到的概率；步骤3，划分虚拟DOI中的虚拟LOR事件，将原先一对晶体中的一条LOR事件在虚拟DOI中划分为多个子事件，并计算每个虚拟子事件的发生概率，将其作为后续重建图像时每条LOR子事件的权重系数；步骤4，对PET重建中的常规系统矩阵进行优化，将步骤2计算获得的虚拟DOI对应的探测概率加入到系统矩阵；步骤5，利用步骤4获得的系统矩阵对图像进行重建，获得优化后的图像。2.如权利要求1所述的一种构造虚拟DOI及相应系统矩阵提高PET图像重建质量的方法，其特征在于：步骤1中根据闪烁晶体的几何排列和尺寸大小进行虚拟DOI的划分。3.如权利要求1所述的一种构造虚拟DOI及相应系统矩阵提高PET图像重建质量的方法，其特征在于：步骤2中计算虚拟DOI的探测概率的实现方式如下，通过连接信号源与某对虚拟DOI中心的连线，计算覆盖在晶体中的长度，确定检测概率D，对于一个光子，它在晶体内被检测概率...

【专利技术属性】
技术研发人员：彭浩，王成，
申请(专利权)人：武汉大学，
类型：发明
国别省市：湖北,42