CN111815734B - 一种pet扫描数据修复方法、装置和计算机设备 - Google Patents
一种pet扫描数据修复方法、装置和计算机设备 Download PDFInfo
- Publication number
- CN111815734B CN111815734B CN202010847753.1A CN202010847753A CN111815734B CN 111815734 B CN111815734 B CN 111815734B CN 202010847753 A CN202010847753 A CN 202010847753A CN 111815734 B CN111815734 B CN 111815734B
- Authority
- CN
- China
- Prior art keywords
- chord
- detector
- data
- chord graph
- graph
- 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
- 238000000034 method Methods 0.000 title claims abstract description 73
- 230000004044 response Effects 0.000 claims abstract description 70
- 239000013078 crystal Substances 0.000 claims abstract description 64
- 230000010354 integration Effects 0.000 claims abstract description 17
- 238000012545 processing Methods 0.000 claims abstract description 14
- 230000008439 repair process Effects 0.000 claims description 44
- 238000004590 computer program Methods 0.000 claims description 14
- 238000009826 distribution Methods 0.000 claims description 14
- 238000012937 correction Methods 0.000 claims description 12
- 230000006870 function Effects 0.000 claims description 11
- 238000012567 pattern recognition method Methods 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 6
- 230000001502 supplementing effect Effects 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 abstract description 7
- 229920000139 polyethylene terephthalate Polymers 0.000 description 81
- 239000005020 polyethylene terephthalate Substances 0.000 description 81
- 238000003860 storage Methods 0.000 description 9
- 238000010586 diagram Methods 0.000 description 8
- 230000002159 abnormal effect Effects 0.000 description 6
- 230000008901 benefit Effects 0.000 description 3
- 238000005094 computer simulation Methods 0.000 description 3
- 238000003759 clinical diagnosis Methods 0.000 description 2
- 150000001875 compounds Chemical class 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000009499 grossing Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 108010001267 Protein Subunits Proteins 0.000 description 1
- 230000005856 abnormality Effects 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 210000000988 bone and bone Anatomy 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 238000013135 deep learning Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000000338 in vitro Methods 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 230000004060 metabolic process Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 210000003205 muscle Anatomy 0.000 description 1
- 238000011056 performance test Methods 0.000 description 1
- -1 polyethylene terephthalate Polymers 0.000 description 1
- 238000002600 positron emission tomography Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000012549 training Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/008—Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/416—Exact reconstruction
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Nuclear Medicine (AREA)
- Measurement Of Radiation (AREA)
Abstract
本申请涉及一种PET扫描数据修复方法、装置和计算机设备，其中，该PET扫描数据修复方法包括：通过获取PET探测器接收到的符合事件信息，根据符合事件信息获取基于探测器晶体响应线的第一弦图以及基于探测器晶体响应线整合处理后形成的第二弦图，根据各探测器模块接收到的符合事件信息，确定第一探测器模块，第一探测器模块的符合事件计数率低于第一阈值，修复第二弦图中与所述第一探测器模块相关的缺失数据，获得第二弦图修复数据，根据所述第二弦图修复数据修复所述第一弦图，得到修复后的第一弦图，解决了相关技术中失效探测器模块对重建图像质量的影响，消除了图像伪影，提高了PET图像的定量准确度。
Description
技术领域
本申请涉及医疗影像技术领域，特别是涉及一种PET扫描数据修复方法、装置和计算机设备。
背景技术
正电子发射计算机断层扫描(Positron Emission Tomography，简称PET)是一种利用向生物体内部注入正电子放射性同位素标记的化合物，而在体外测量它们的空间分布和时间特性的三维成像无损检测技术，轴向视野的不断增长是PET的主要趋势，一方面，更长的轴向视野带来了探测灵敏度的巨大提升，有助于进一步降低被扫描人受到的辐射剂量，或是进行更精细的动态代谢研究；但在另一方面，更多的探测器模块也大大增加了PET系统出现失效探测器模块(也就是坏通道)的概率，失效探测器模块会导致与该模块相关的响应线(Lines Of Resolution，简称LOR)上的计数异常，在弦图中的表现为一条或若干条计数为零的计数异常的黑线，使用这样的弦图数据直接进行PET图像重建，会造成图像伪影和定量异常，严重影响临床诊断。
在相关技术中，针对失效探测器模块会导致与该模块相关的响应线(Lines OfResolution，简称LOR)上的计数异常，在弦图中的表现为一条或若干条计数为零的计数异常的黑线，一种解决方式是在重建系统建模中去除与该失效探测器模块相关的响应线，该方式可以在一定程度上缓解图像伪影，但是无法修复定量偏差，进而存在重建的PET扫描图像精度低的问题。
目前针对相关技术中在重建系统建模中去除与失效探测器模块相关响应线的方式来消除失效探测器模块的影响，存在重建的PET扫描图像精度低的问题，尚未提出有效的解决方案。
发明内容
本申请实施例提供了一种PET扫描数据修复方法、装置和计算机设备，解决了相关技术中失效探测器模块对重建图像质量的影响，消除了图像伪影，提高了PET图像的定量准确度。
第一方面，本申请实施例提供了一种PET扫描数据修复方法，所述方法包括：获取PET探测器接收到的符合事件信息，所述PET探测器由多个探测器模块组成；
根据所述符合事件信息获取基于探测器晶体响应线的第一弦图以及基于探测器晶体响应线整合处理后形成的第二弦图；
根据所述探测器模块接收到的符合事件信息，确定第一探测器模块，其中，所述第一探测器模块的符合事件计数率低于第一阈值；
修复所述第二弦图中与所述第一探测器模块相关的缺失数据，获得第二弦图修复数据；
根据所述第二弦图修复数据修复所述第一弦图，得到修复后的第一弦图。
在其中一些实施例中，所述方法还包括：所述基于探测器晶体响应线整合处理后形成的第二弦图包括：
以探测器模块为单元整合处理探测器晶体响应线。
在其中一些实施例中，所述修复所述第二弦图中与所述第一探测器模块相关的缺失数据，得到第二弦图修复数据包括：
通过弦图空间内与第一探测器模块弦图数据相邻的无损坏弦图数据的数据分布情况补充所述第二弦图中与所述第一探测器模块相关的缺失数据，将所述补充的缺失数据作为第二弦图修复数据。
在其中一些实施例中，所述根据符合事件信息获取基于探测器晶体响应线整合处理后形成的第二弦图之后，所述方法还包括：
对所述第二弦图进行随机校正和/或散射校正。
在其中一些实施例中，所述修复所述第二弦图中与所述第一探测器模块相关的缺失数据，得到第二弦图修复数据包括：
通过插值法或模式识别法修复所述第二弦图中所述第一探测器模块的缺失数据，得到第二弦图修复数据。
在其中一些实施例中，根据所述第二弦图修复数据修复所述第一弦图，得到修复后的第一弦图之前，所述方法还包括：
获取所述第二弦图修复数据的平滑函数值，若平滑函数值小于第二阈值，则根据所述第二弦图修复数据修复所述第一弦图，得到修复后的第一弦图。
在其中一些实施例中，所述根据所述第二弦图修复数据修复所述第一弦图，得到修复后的第一弦图包括：
将所述第二弦图修复数据随机分配给第一弦图中相应的晶体响应线，得到修复后的第一弦图。
第二方面，本申请实施例提供了一种PET扫描数据的重建方法，所述方法对通过第一方面所述的PET扫描数据修复方法获得的PET扫描数据进行图像重建，获得PET扫描图像。
第三方面，本申请实施例提供了一种PET扫描数据修复装置，所述装置包括获取模块、第一修复模块和第二修复模块；
所述获取模块，用于获取PET探测器接收到的符合事件信息，所述PET探测器由多个探测器模块组成，所述获取模块，还用于根据所述符合事件信息获取基于探测器晶体响应线的第一弦图以及基于探测器晶体响应线整合处理后形成的第二弦图；
所述第一修复模块，用于根据所述探测器模块接收到的符合事件信息，确定第一探测器模块，其中，所述第一探测器模块的符合事件计数率低于第一阈值；所述第一修复模块，还用于修复所述第二弦图中与所述第一探测器模块相关的缺失数据，获得第二弦图修复数据；
所述第二修复模块，用于根据所述第二弦图修复数据修复所述第一弦图，得到修复后的第一弦图。
第四方面，本申请实施例提供了一种计算机设备，包括存储器、处理器以及存储在所述存储器上并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现如上述第一方面所述的PET扫描数据修复方法。
相比于相关技术，本申请实施例提供的一种PET扫描数据修复方法、装置和计算机设备，通过获取PET探测器接收到的符合事件信息，所述PET探测器由多个探测器模块组成，根据所述符合事件信息获取基于探测器晶体响应线的第一弦图以及基于探测器晶体响应线整合处理后形成的第二弦图，根据所述探测器模块接收到的符合事件信息，确定第一探测器模块，所述第一探测器模块的符合事件计数率低于第一阈值，修复所述第二弦图中与所述第一探测器模块相关的缺失数据，获得第二弦图修复数据，根据所述第二弦图修复数据修复所述第一弦图，得到修复后的第一弦图，解决了相关技术中失效探测器模块对重建图像质量的影响，消除了图像伪影，提高了PET图像的定量准确度。
附图说明
此处所说明的附图用来提供对本申请的进一步理解，构成本申请的一部分，本申请的示意性实施例及其说明用于解释本申请，并不构成对本申请的不当限定。在附图中：
图1是根据本申请实施例的PET扫描数据修复方法的流程图一；
图2是根据本申请实施例的PET扫描数据修复方法的流程图二；
图3是根据本申请实施例的PET扫描数据修复方法的流程图三；
图4是根据本申请实施例的PET扫描数据修复方法的流程图四；
图5是根据本申请实施例的PET扫描数据修复方法的流程图五；
图6是根据本申请实施例的PET扫描数据修复方法的流程图六；
图7是根据本申请实施例的PET扫描数据修复装置的结构框图；
图8是根据本申请实施例的计算机设备的内部结构示意图。
具体实施方式
为了使本申请的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本申请进行描述和说明。应当理解，此处所描述的具体实施例仅仅用以解释本申请，并不用于限定本申请。基于本申请提供的实施例，本领域普通技术人员在没有作出创造性劳动的前提下所获得的所有其他实施例，都属于本申请保护的范围。
显而易见地，下面描述中的附图仅仅是本申请的一些示例或实施例，对于本领域的普通技术人员而言，在不付出创造性劳动的前提下，还可以根据这些附图将本申请应用于其他类似情景。此外，还可以理解的是，虽然这种开发过程中所作出的努力可能是复杂并且冗长的，然而对于与本申请公开的内容相关的本领域的普通技术人员而言，在本申请揭露的技术内容的基础上进行的一些设计，制造或者生产等变更只是常规的技术手段，不应当理解为本申请公开的内容不充分。
在本申请中提及“实施例”意味着，结合实施例描述的特定特征、结构或特性可以包含在本申请的至少一个实施例中。在说明书中的各个位置出现该短语并不一定均是指相同的实施例，也不是与其它实施例互斥的独立的或备选的实施例。本领域普通技术人员显式地和隐式地理解的是，本申请所描述的实施例在不冲突的情况下，可以与其它实施例相结合。
除非另作定义，本申请所涉及的技术术语或者科学术语应当为本申请所属技术领域内具有一般技能的人士所理解的通常意义。本申请所涉及的“一”、“一个”、“一种”、“该”等类似词语并不表示数量限制，可表示单数或复数。本申请所涉及的术语“包括”、“包含”、“具有”以及它们任何变形，意图在于覆盖不排他的包含；例如包含了一系列步骤或模块(单元)的过程、方法、系统、产品或设备没有限定于已列出的步骤或单元，而是可以还包括没有列出的步骤或单元，或可以还包括对于这些过程、方法、产品或设备固有的其它步骤或单元。本申请所涉及的“连接”、“相连”、“耦接”等类似的词语并非限定于物理的或者机械的连接，而是可以包括电气的连接，不管是直接的还是间接的。本申请所涉及的“多个”是指两个或两个以上。“和/或”描述关联对象的关联关系，表示可以存在三种关系，例如，“A和/或B”可以表示：单独存在A，同时存在A和B，单独存在B这三种情况。字符“/”一般表示前后关联对象是一种“或”的关系。本申请所涉及的术语“第一”、“第二”、“第三”等仅仅是区别类似的对象，不代表针对对象的特定排序。
在本申请的实施例中，PET系统中的探测器包括若干个探测器模块，探测器模块包括晶体阵列，一对符合事件的探测器模块的晶体之间形成响应线；随着更长的轴向视野，更多的探测器模块也大大增加了PET系统出现失效探测器模块的概率，失效探测器模块会导致与该模块相关的响应线(Lines Of Resolution，简称LOR)上的计数异常，在弦图中的表现为一条或若干条计数为零的计数异常的黑线，使用这样的弦图数据直接进行PET图像重建，会造成图像伪影和定量异常，严重影响临床诊断，在相关技术中，一种解决方式是在重建系统建模中去除与该失效探测器模块相关的响应线，该方式可以在一定程度上缓解图像伪影，但是无法修复定量偏差，进而存在重建的PET扫描图像精度低的问题；本申请提供的PET扫描数据修复方法，通过获取PET探测器接收到的符合事件信息，PET探测器由多个探测器模块组成，根据符合事件信息获取基于探测器晶体响应线的第一弦图以及基于探测器晶体响应线整合处理后形成的第二弦图，并根据探测器模块接收到的符合事件信息，确定符合事件计数率低于第一阈值的第一探测器模块，并修复第二弦图中与第一探测器模块相关的缺失数据以获得第二弦图修复数据，并根据第二弦图修复数据修复第一弦图，得到修复后的第一弦图，利用修复后的第一弦图的弦图数据重建PET扫描图像，解决了相关技术中失效探测器模块对重建图像质量的影响，消除了图像伪影，提高了PET图像的定量准确度。
本实施例提供了一种PET扫描数据修复方法，图1是根据本申请实施例的PET扫描数据修复方法的流程图一，如图1所示，该方法包括如下步骤：
步骤S101，获取PET探测器接收到的符合事件信息，PET探测器由多个探测器模块组成，根据符合事件信息获取基于探测器晶体响应线的第一弦图以及基于探测器晶体响应线整合处理后形成的第二弦图；
需要说明的是，第一弦图是基于探测器中单个晶体响应线为单位的弦图，由于探测器模块包括晶体阵列，因此第二弦图可以是基于探测器中一定角度的晶体阵列响应线为单位的弦图，也可以是基于探测器中整个晶体阵列响应线为单位的弦图。
步骤S102，根据探测器模块接收到的符合事件信息，确定第一探测器模块，其中，第一探测器模块的符合事件计数率低于第一阈值；
需要说明的是，PET系统中若某探测器模块失效，则该失效探测器模块对应的全部晶体响应线均存在异常，因此PET系统中符合事件计算率低于第一阈值的探测器模块便是失效探测器模块，失效探测器模块也就是本实施例中的第一探测器模块。
步骤S103，修复第二弦图中与第一探测器模块相关的缺失数据，获得第二弦图修复数据；
其中，弦图等价于各个角度观察物体后得到的物体信息的线积分分布，在已知失效探测器模块也就是第一探测器模块的情况下，可根据第二弦图中的数据分布，预测第二弦图中与第一探测器模块相关的缺失数据，该缺失数据便是第二弦图修复数据。
步骤S104，根据第二弦图修复数据修复第一弦图，得到修复后的第一弦图；
需要说明的是，在已知失效探测器模块的缺失数据的情况下，将该缺失数据分配给失效探测器模块中相关的晶体响应线，并将失效探测器模块中相关的晶体响应线数据加入至基于探测器中单个晶体响应线为单位的第一弦图中，以得到基于完整探测器晶体响应线的第一弦图。
通过上述步骤S101至步骤S104，获取基于探测器晶体响应线的第一弦图以及基于探测器晶体响应线整合处理后形成的第二弦图，并预测失效探测器模块在第二弦图中的缺失数据，并根据缺失数据修复第一弦图，以得到基于完整探测器晶体响应线的第一弦图，利用修复后的第一弦图的弦图数据重建PET扫描图像，解决了相关技术中失效探测器模块对重建图像质量的影响，消除了图像伪影，提高了PET图像的定量准确度。
在其中一些实施例中，基于探测器晶体响应线整合处理后形成的第二弦图包括：以探测器模块为单元整合处理探测器晶体响应线；
其中，由于探测器模块包括晶体阵列，且若探测器模块失效则失效的探测器模块对应的全部晶体响应线均存在异常，以探测器模块为单元整合处理探测器晶体响应线得到第二弦图，也就是将一对失效探测器模块中的多根晶体响应线合并为一根响应线，以增大该响应线上的计数统计性。
在其中一些实施例中，图2是根据本申请实施例的PET扫描数据修复方法的流程图二，如图2所示，修复第二弦图中与第一探测器模块相关的缺失数据，得到第二弦图修复数据的方法包括如下步骤：
步骤S201，通过弦图空间内与第一探测器模块弦图数据相邻的无损坏弦图数据的数据分布情况补充第二弦图中与第一探测器模块相关的缺失数据，将补充的缺失数据作为第二弦图修复数据；
其中，弦图等价于各个角度观察物体后得到的物体信息的线积分分布，进而它的排列分布决定了在弦图域上这些信息一定是连续平滑的，进而可通过弦图空间内与第一探测器模块弦图数据相邻的无损坏弦图数据的数据分布情况补充第二弦图中与第一探测器模块相关的缺失数据。
通过上述步骤S201，根据弦图中物体信息的连续平滑的排列分布，进而可通过弦图空间内与第一探测器模块弦图数据相邻的无损坏弦图数据的数据分布情况补充第二弦图中与第一探测器模块相关的缺失数据，进而方便获取与第一探测器模块相关的且具有一定可靠性的缺失数据。
在其中一些实施例中，图3是根据本申请实施例的PET扫描数据修复方法的流程图三，如图3所示，根据符合事件信息获取基于探测器晶体响应线整合处理后形成的第二弦图之后，该方法还包括如下步骤：
步骤S301，对第二弦图进行随机校正和/或散射校正；
需要说明的是，PET系统是复杂精密的大型医用诊断设备，在使用过程中，需定期做性能测试和标定，不断更新校正参数以确保PET系统处于良好的工作状态，对第二弦图进行随机校正可以是采用对应的算法进行同位素的衰变校正、体内肌肉和骨骼等对光子的吸收衰减校正等。
通过上述步骤S301，对第二弦图进行随机校正和/或散射校正，以提高第二弦图中弦图数据的真实度。
在其中一些实施例中，图4是根据本申请实施例的PET扫描数据修复方法的流程图四，如图4所示，修复第二弦图中与第一探测器模块相关的缺失数据，得到第二弦图修复数据的方法包括如下步骤：
步骤S401，通过插值法或模式识别法修复第二弦图中第一探测器模块的缺失数据，得到第二弦图修复数据；
需要说明的是，弦图等价于各个角度观察物体后得到的物体信息的线积分分布，进而它的排列分布决定了在弦图域上这些信息一定是连续平滑的，使用插值法可以计算缺失数据的估计值，所谓的插值法就是通过两点(x0，y0)，(x1，y1)估计中间点的值，假设y＝f(x)是一条直线，通过已知的两点来计算函数f(x)，然后只要知道x就能求出y，以此方法来估计缺失值，也可以假设f(x)不是直线而是其他函数；需要进一步说明的是，模式识别法是指利用已有的正常数据获得的无损弦图，人为指定无损弦图中部分模块损坏，从而获得无损弦图对应的损坏弦图，将大量的类似组合构成数据集进行训练，获得从损坏弦图到无损弦图的映射关系，进而根据已训练的映射关系，在已知第二弦图中第一探测器模块的缺失数据的情况下，得到第二弦图修复数据。
通过步骤S401，基于插值法或模式识别法修复第二弦图中第一探测器模块的缺失数据，得到第二弦图修复数据，进而方便获取与第一探测器模块相关的且具有一定可靠性的缺失数据。
在其中一些实施例中，图5是根据本申请实施例的PET扫描数据修复方法的流程图五，如图5所示，根据第二弦图修复数据修复第一弦图，得到修复后的第一弦图之前，该方法还包括如下步骤：
步骤S501，获取第二弦图修复数据的平滑函数值，若平滑函数值小于第二阈值，则根据第二弦图修复数据修复第一弦图，得到修复后的第一弦图。
需要说明的是，平滑函数值，即每次迭代中更新第二弦图修复数据前后，第二弦图修复数据中每个像素值变化量的均方差，因此若均方差小于预设的第二阈值，则第二弦图修复数据更加符合真实性。
通过步骤S501，在将第二弦图修复数据修复第一弦图之前，先对第二弦图修复数据进行判断，在第二弦图修复数据的平滑函数值小于第二阈值的情况下，将第二弦图修复数据修复第一弦图，以确保用来修复第一弦图的第二弦图修复数据的真实性。
在其中一些实施例中，图6是根据本申请实施例的PET扫描数据修复方法的流程图六，如图6所示，根据第二弦图修复数据修复第一弦图，得到修复后的第一弦图的方法包括如下步骤：
步骤S601，将第二弦图修复数据随机分配给第一弦图中相应的晶体响应线，得到修复后的第一弦图；
其中，第二弦图修复数据是整个失效探测器模块的响应线数据，因此在已知失效探测器模块的缺失数据的情况下，将该缺失数据分配给失效探测器模块中相关的晶体响应线，且失效探测器模块中相关的晶体响应线数据加入至基于探测器中单个晶体响应线为单位的第一弦图中，以得到基于完整探测器晶体响应线的第一弦图。
通过步骤S601，将第二弦图修复数据分配给失效探测器模块中相关的晶体响应线，且失效探测器模块中相关的晶体响应线数据加入至基于探测器中单个晶体响应线为单位的第一弦图中，以得到基于完整探测器晶体响应线的第一弦图，进而使得利用修复后的第一弦图的弦图数据重建PET扫描图像，提高了重建的PET扫描图像的精度。
需要说明的是，在上述流程中或者附图的流程图中示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行，并且，虽然在流程图中示出了逻辑顺序，但是在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤。
本实施例提供了PET扫描数据的重建方法，对通过上述PET扫描数据修复方法中获得的PET扫描数据进行图像重建，获得PET扫描图像，需要说明的是，在本实施例中，可以将通过上述PET扫描数据修复方法中修复后的第一弦图，输入至经训练的深度学习图像重建模型，以获取PET扫描图像，相较于相关技术中在重建系统建模中去除与失效探测器模块相关响应线的方式来消除失效探测器模块的影响，提高了重建的PET扫描图像的精度。
本实施例还提供了一种PET扫描数据修复装置，该装置用于实现上述实施例及优选实施方式，已经进行过说明的不再赘述。如以下所使用的，术语“模块”、“单元”、“子单元”等可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。
图7是根据本申请实施例的PET扫描数据修复装置的结构框图，如图7所示，装置包括获取模块71、第一修复模块72和第二修复模块73；
获取模块71，用于获取PET探测器接收到的符合事件信息，PET探测器由多个探测器模块组成；获取模块71，还用于根据符合事件信息获取基于探测器晶体响应线的第一弦图以及基于探测器晶体响应线整合处理后形成的第二弦图；
第一修复模块72，用于根据探测器模块接收到的符合事件信息，确定第一探测器模块，其中，第一探测器模块的符合事件计数率低于第一阈值；第一修复模块72，还用于修复第二弦图中与第一探测器模块相关的缺失数据，获得第二弦图修复数据；
第二修复模块73，用于根据第二弦图修复数据修复第一弦图，得到修复后的第一弦图。
通过上述的PET扫描数据修复装置，获取基于探测器晶体响应线的第一弦图以及基于探测器晶体响应线整合处理后形成的第二弦图，并预测失效探测器模块在第二弦图中的缺失数据，并根据缺失数据修复第一弦图，以得到基于完整探测器晶体响应线的第一弦图，利用修复后的第一弦图的弦图数据重建PET扫描图像，解决了相关技术中失效探测器模块对重建图像质量的影响，消除了图像伪影，提高了PET图像的定量准确度。
在其中一些实施例中，获取模块71、第一修复模块72和第二修复模块73还用于实现上述各实施例提供的PET扫描数据修复方法中的步骤，在这里不再赘述。
需要说明的是，上述各个模块可以是功能模块也可以是程序模块，既可以通过软件来实现，也可以通过硬件来实现。对于通过硬件来实现的模块而言，上述各个模块可以位于同一处理器中；或者上述各个模块还可以按照任意组合的形式分别位于不同的处理器中。
在一个实施例中，提供了一种计算机设备，该计算机设备可以是终端。该计算机设备包括通过系统总线连接的处理器、存储器、网络接口、显示屏和输入装置。其中，该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作系统和计算机程序。该内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该计算机设备的网络接口用于与外部的终端通过网络连接通信。该计算机程序被处理器执行时以实现一种PET扫描数据修复方法。该计算机设备的显示屏可以是液晶显示屏或者电子墨水显示屏，该计算机设备的输入装置可以是显示屏上覆盖的触摸层，也可以是计算机设备外壳上设置的按键、轨迹球或触控板，还可以是外接的键盘、触控板或鼠标等。
在一个实施例中，图8是根据本申请实施例的计算机设备的内部结构示意图，如图8所示，提供了一种计算机设备，该计算机设备可以是服务器，其内部结构图可以如图8所示。该计算机设备包括通过系统总线连接的处理器、存储器、网络接口和数据库。其中，该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作系统、计算机程序和数据库。该内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该计算机设备的数据库用于存储数据。该计算机设备的网络接口用于与外部的终端通过网络连接通信。该计算机程序被处理器执行时以实现一种PET扫描数据修复方法。
本领域技术人员可以理解，图8中示出的结构，仅仅是与本申请方案相关的部分结构的框图，并不构成对本申请方案所应用于其上的计算机设备的限定，具体的计算机设备可以包括比图中所示更多或更少的部件，或者组合某些部件，或者具有不同的部件布置。
在一个实施例中，提供了一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，处理器执行计算机程序时实现上述各实施例提供的PET扫描数据修复方法中的步骤。
在一个实施例中，提供了一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现上述各个实施例提供的PET扫描数据修复方法中的步骤。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，该计算机程序可存储于一非易失性计算机可读取存储介质中，该计算机程序在执行时，可包括如上述各方法的实施例的流程。其中，本申请所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用，均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器(ROM)、可编程ROM(PROM)、电可编程ROM(EPROM)、电可擦除可编程ROM(EEPROM)或闪存。易失性存储器可包括随机存取存储器(RAM)或者外部高速缓冲存储器。作为说明而非局限，RAM以多种形式可得，诸如静态RAM(SRAM)、动态RAM(DRAM)、同步DRAM(SDRAM)、双数据率SDRAM(DDRSDRAM)、增强型SDRAM(ESDRAM)、同步链路(Synchlink)DRAM(SLDRAM)、存储器总线(Rambus)直接RAM(RDRAM)、直接存储器总线动态RAM(DRDRAM)、以及存储器总线动态RAM(RDRAM)等。
以上所述实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本申请的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本申请构思的前提下，还可以做出若干变形和改进，这些都属于本申请的保护范围。因此，本申请专利的保护范围应以所附权利要求为准。
Claims (9)
1.一种PET扫描数据修复方法，其特征在于，所述方法包括：
获取PET探测器接收到的符合事件信息，所述PET探测器由多个探测器模块组成；
根据所述符合事件信息获取基于探测器晶体响应线的第一弦图以及基于探测器晶体响应线整合处理后形成的第二弦图；所述第一弦图为基于探测器中单个晶体响应线为单位的弦图；所述第二弦图为基于探测器中一定角度晶体阵列或整个晶体阵列响应线为单位的弦图；
根据所述探测器模块接收到的符合事件信息，确定第一探测器模块，其中，所述第一探测器模块的符合事件计数率低于第一阈值；
修复所述第二弦图中与所述第一探测器模块相关的缺失数据，获得第二弦图修复数据；所述修复所述第二弦图中与所述第一探测器模块相关的缺失数据，获得第二弦图修复数据包括：通过插值法或模式识别法修复所述第二弦图中所述第一探测器模块的缺失数据，获得第二弦图修复数据；
根据所述第二弦图修复数据修复所述第一弦图，得到修复后的第一弦图。
2.根据权利要求1所述的方法，其特征在于，所述基于探测器晶体响应线整合处理后形成的第二弦图包括：
以探测器模块为单元整合处理探测器晶体响应线。
3.根据权利要求1所述的方法，其特征在于，所述修复所述第二弦图中与所述第一探测器模块相关的缺失数据，获得第二弦图修复数据包括：
通过弦图空间内与第一探测器模块弦图数据相邻的无损坏弦图数据的数据分布情况补充所述第二弦图中与所述第一探测器模块相关的缺失数据，将所述补充的缺失数据作为第二弦图修复数据。
4.根据权利要求1所述的方法，其特征在于，所述根据符合事件信息获取基于探测器晶体响应线整合处理后形成的第二弦图之后，所述方法还包括：
对所述第二弦图进行随机校正和/或散射校正。
5.根据权利要求1所述的方法，其特征在于，根据所述第二弦图修复数据修复所述第一弦图，得到修复后的第一弦图之前，所述方法还包括：
获取所述第二弦图修复数据的平滑函数值，若平滑函数值小于第二阈值，则根据所述第二弦图修复数据修复所述第一弦图，得到修复后的第一弦图。
6.根据权利要求1所述的方法，其特征在于，所述根据所述第二弦图修复数据修复所述第一弦图，得到修复后的第一弦图包括：
将所述第二弦图修复数据随机分配给第一弦图中相应的晶体响应线，得到修复后的第一弦图。
7.一种PET扫描数据的重建方法，其特征在于，对权利要求1至6任一项所述的方法获得的PET扫描数据进行图像重建，获得PET扫描图像。
8.一种PET扫描数据修复装置，其特征在于，所述装置包括获取模块、第一修复模块和第二修复模块；
所述获取模块，用于获取PET探测器接收到的符合事件信息，所述PET探测器由多个探测器模块组成，所述获取模块，还用于根据所述符合事件信息获取基于探测器晶体响应线的第一弦图以及基于探测器晶体响应线整合处理后形成的第二弦图；所述第一弦图为基于探测器中单个晶体响应线为单位的弦图；所述第二弦图为基于探测器中一定角度晶体阵列或整个晶体阵列响应线为单位的弦图；
所述第一修复模块，用于根据所述探测器模块接收到的符合事件信息，确定第一探测器模块，其中，所述第一探测器模块的符合事件计数率低于第一阈值；所述第一修复模块，还用于修复所述第二弦图中与所述第一探测器模块相关的缺失数据，获得第二弦图修复数据；所述修复所述第二弦图中与所述第一探测器模块相关的缺失数据，得到第二弦图修复数据包括：通过插值法或模式识别法修复所述第二弦图中所述第一探测器模块的缺失数据，获得第二弦图修复数据；
所述第二修复模块，用于根据所述第二弦图修复数据修复所述第一弦图，得到修复后的第一弦图。
9.一种计算机设备，包括存储器、处理器以及存储在所述存储器上并可在所述处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现如权利要求1至6中任一项所述的PET扫描数据修复方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010847753.1A CN111815734B (zh) | 2020-08-21 | 2020-08-21 | 一种pet扫描数据修复方法、装置和计算机设备 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010847753.1A CN111815734B (zh) | 2020-08-21 | 2020-08-21 | 一种pet扫描数据修复方法、装置和计算机设备 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111815734A CN111815734A (zh) | 2020-10-23 |
| CN111815734B true CN111815734B (zh) | 2024-04-16 |
Family
ID=72860460
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010847753.1A Active CN111815734B (zh) | 2020-08-21 | 2020-08-21 | 一种pet扫描数据修复方法、装置和计算机设备 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111815734B (zh) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102509322A (zh) * | 2011-11-11 | 2012-06-20 | 刘华锋 | 一种基于卡尔曼滤波的pet图像重建方法 |
| CN102968762A (zh) * | 2012-10-24 | 2013-03-13 | 浙江理工大学 | 一种基于稀疏化和泊松模型的pet重建方法 |
| WO2018033003A1 (en) * | 2016-08-15 | 2018-02-22 | Shenzhen United Imaging Healthcare Co., Ltd. | System and method for pet data correction |
| WO2018108641A1 (en) * | 2016-12-14 | 2018-06-21 | Koninklijke Philips N.V. | Dead pixel correction for digital pet reconstruction |
| CN110025329A (zh) * | 2019-04-24 | 2019-07-19 | 上海联影医疗科技有限公司 | 符合计数弦图生成方法、装置、计算机设备和存储介质 |
| CN110584698A (zh) * | 2019-09-23 | 2019-12-20 | 上海联影医疗科技有限公司 | 探测器质量控制效验方法、装置、计算机设备和存储介质 |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7381959B2 (en) * | 2005-08-17 | 2008-06-03 | General Electric Company | Technique for reconstructing PET scan images |
| US7718954B2 (en) * | 2006-01-31 | 2010-05-18 | Koninklijke Philips Electronics N.V. | Component method and system for PET detector efficiency normalization |
| US7858943B2 (en) * | 2007-09-18 | 2010-12-28 | Siemens Medical Solutions Usa, Inc. | Methods of performing PET sinogram quality check |
| US9507033B2 (en) * | 2013-02-05 | 2016-11-29 | Siemens Medical Solutions Usa, Inc. | Method and apparatus for compensating for scattering of emission gamma photons for PET imaging |
| US9804275B2 (en) * | 2013-11-19 | 2017-10-31 | Toshiba Medical Systems Corporation | Method and system of random-event-based count-rate adaptive normalization for PET detector efficiency |
| US9483850B2 (en) * | 2013-11-19 | 2016-11-01 | Toshiba Medical Systems Corporation | Cylinder source software-positioning method for PET calibration and image quality assurance |
-
2020
- 2020-08-21 CN CN202010847753.1A patent/CN111815734B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102509322A (zh) * | 2011-11-11 | 2012-06-20 | 刘华锋 | 一种基于卡尔曼滤波的pet图像重建方法 |
| CN102968762A (zh) * | 2012-10-24 | 2013-03-13 | 浙江理工大学 | 一种基于稀疏化和泊松模型的pet重建方法 |
| WO2018033003A1 (en) * | 2016-08-15 | 2018-02-22 | Shenzhen United Imaging Healthcare Co., Ltd. | System and method for pet data correction |
| WO2018108641A1 (en) * | 2016-12-14 | 2018-06-21 | Koninklijke Philips N.V. | Dead pixel correction for digital pet reconstruction |
| CN110025329A (zh) * | 2019-04-24 | 2019-07-19 | 上海联影医疗科技有限公司 | 符合计数弦图生成方法、装置、计算机设备和存储介质 |
| CN110584698A (zh) * | 2019-09-23 | 2019-12-20 | 上海联影医疗科技有限公司 | 探测器质量控制效验方法、装置、计算机设备和存储介质 |
Non-Patent Citations (1)
| Title |
|---|
| 基于FMECA的PET-CT可靠性评价指标研究;陈颖等;《中国医疗设备》;第34卷(第08期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111815734A (zh) | 2020-10-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111080584B (zh) | 医学图像的质控方法、计算机设备和可读存储介质 | |
| CN106491151B (zh) | Pet图像获取方法及系统 | |
| CN110063742B (zh) | 散射校正方法、装置、计算机设备和存储介质 | |
| CN110584698A (zh) | 探测器质量控制效验方法、装置、计算机设备和存储介质 | |
| CN110415310A (zh) | 医学扫描成像方法、装置、存储介质及计算机设备 | |
| CN110136076A (zh) | 医学扫描成像方法、装置、存储介质及计算机设备 | |
| CN111161182B (zh) | Mr结构信息约束的非局部均值引导的pet图像部分容积校正方法 | |
| US7186984B2 (en) | Method for quantifying the radioactivity of living structures of small dimensions by employing emission tomography | |
| CN111243052A (zh) | 图像重建方法、装置、计算机设备和存储介质 | |
| Pan et al. | Iterative residual optimization network for limited-angle tomographic reconstruction | |
| CN110687585B (zh) | 获取晶体效率的方法、装置、计算机设备和存储介质 | |
| CN111588399B (zh) | 医学成像设备状态监控的方法、设备和计算机设备 | |
| CN111815734B (zh) | 一种pet扫描数据修复方法、装置和计算机设备 | |
| CN116071506B (zh) | 四维血管造影的重建方法、装置、计算机设备和存储介质 | |
| CN110473241B (zh) | 图像配准方法、存储介质及计算机设备 | |
| Pretorius et al. | Observer studies of image quality of denoising reduced-count cardiac single photon emission computed tomography myocardial perfusion imaging by three-dimensional Gaussian post-reconstruction filtering and deep learning | |
| CN112488949A (zh) | 一种低剂量pet图像还原方法、系统、设备和介质 | |
| CN116721045A (zh) | 一种多ct图像融合的方法及装置 | |
| Van Eyndhoven et al. | Combined motion estimation and reconstruction in tomography | |
| Mesquita et al. | Choosing the ART relaxation parameter for Clear-PEM 2D image reconstruction | |
| CN114943782A (zh) | 一种正电子发射断层成像方法、系统和装置 | |
| JP7459243B2 (ja) | １以上のニューラルネットワークとしての画像形成のモデル化による画像再構成 | |
| CN112365479A (zh) | Pet参数图像处理方法、装置、计算机设备及存储介质 | |
| US20240029324A1 (en) | Method for image reconstruction, computer device and storage medium | |
| CN111914393B (zh) | 死时间校正方法、装置、计算机设备和存储介质 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB02 | Change of applicant information |
Address after: 201807 Shanghai City, north of the city of Jiading District Road No. 2258 Applicant after: Shanghai Lianying Medical Technology Co.,Ltd. Address before: 201807 Shanghai City, north of the city of Jiading District Road No. 2258 Applicant before: SHANGHAI UNITED IMAGING HEALTHCARE Co.,Ltd. |
|
| CB02 | Change of applicant information | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |