CN102640016B - Pet数据处理方法、pet数据处理装置以及数据处理方法 - Google Patents
Pet数据处理方法、pet数据处理装置以及数据处理方法 Download PDFInfo
- Publication number
- CN102640016B CN102640016B CN201180004335.5A CN201180004335A CN102640016B CN 102640016 B CN102640016 B CN 102640016B CN 201180004335 A CN201180004335 A CN 201180004335A CN 102640016 B CN102640016 B CN 102640016B
- Authority
- CN
- China
- Prior art keywords
- mentioned
- sinogram
- triangle
- pet
- former
- 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
- 238000003672 processing method Methods 0.000 title claims abstract description 15
- 238000012545 processing Methods 0.000 title claims abstract description 5
- 238000000034 method Methods 0.000 claims abstract description 63
- 238000005070 sampling Methods 0.000 claims abstract description 27
- 238000002600 positron emission tomography Methods 0.000 description 30
- 230000008569 process Effects 0.000 description 16
- 230000005251 gamma ray Effects 0.000 description 13
- 239000013078 crystal Substances 0.000 description 11
- 238000013480 data collection Methods 0.000 description 9
- 238000001514 detection method Methods 0.000 description 8
- 230000011218 segmentation Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 238000013461 design Methods 0.000 description 3
- 229940121896 radiopharmaceutical Drugs 0.000 description 3
- 239000012217 radiopharmaceutical Substances 0.000 description 3
- 230000002799 radiopharmaceutical effect Effects 0.000 description 3
- 238000009825 accumulation Methods 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000002596 correlated effect Effects 0.000 description 2
- 230000000875 corresponding effect Effects 0.000 description 2
- 239000003814 drug Substances 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000000704 physical effect Effects 0.000 description 2
- 238000003325 tomography Methods 0.000 description 2
- 238000012879 PET imaging Methods 0.000 description 1
- 230000004071 biological effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 238000007599 discharging Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 230000000155 isotopic effect Effects 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000002285 radioactive effect Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- -1 suck or oral uptake Substances 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
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/40—Scaling of whole images or parts thereof, e.g. expanding or contracting
- G06T3/4007—Scaling of whole images or parts thereof, e.g. expanding or contracting based on interpolation, e.g. bilinear interpolation
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Nuclear Medicine (AREA)
Abstract
提供一种能够快速地生成均匀地采样的正弦图的PET数据处理方法、PET数据处理装置以及数据处理方法。为了重建PET图像而将PET数据进行插补的方法包括：收集通过使用PET扫描仪来对被检体进行扫描而得到的PET现象数据的步骤；根据收集到的PET现象数据，生成二维的LOR正弦图的步骤；为了确定连接LOR正弦图内的采样点的多个三角形，确定在LOR正弦图内相邻的采样点在水平方向与每列连接的多个三角形的步骤；使用这些被确定的多个三角形，根据LOR正弦图确定均匀地采样的正弦图的步骤。
Description
技术领域
本说明书所记载的实施方式涉及使用三角形分割插补法来重建将收集到的检测器信息作为材料的图像的伽马射线扫描仪***(scannersystem)。
背景技术
正电子发射断层摄影(Positron Emission Tomography：PET)的用途在图像诊断领域正在扩大。在PET成像中，通过注射、吸入、或经口摄取，将放射性医药品导入应该形成图像的被检体。投放放射性医药品后，根据其物理性质以及分子生物学性质，集中于被检体内的特定部位。药剂的实际的空间分布、积聚区域中的药品的浓度、以及从投放到最终排出的过程的动态均为临床上具有重要意义的因素。在该过程中，附着于放射性医药品上的一个正电子放出体根据半衰期、分支比等同位素的物理性质而放出正电子。
放射性核素放出正电子。当一个放正电子与一个电子碰撞时，发生一个湮没现象，该正电子以及电子被破坏。大多数情况下，通过一个湮没现象，实质上发生向大致180度相反方向行进的两条511keV伽马射线。
通过检测这两条伽马射线，并引出连结检测部位彼此的线，即引出同时计数线(Line-Of-Response：LOR)，从而能够发现概然性高的原来的分解部位。该过程仅仅是单纯地识别可能发生相互作用的线的过程，如果大量蓄积这些线，则能够通过使用重建断层的过程，来推定原来的分布。如果在除了上述2个闪烁现象的部位之外的部位，还能够准确利用定时(几百皮秒以内)，则能够通过计算飞行时间 。
发明内容
然而，作为插补原LOR正弦图时发现的难点的例子，可以列举出(1)由于它们没有存在于矩形的格子上，因此难以确定所希望的采样点的最近的点，以及(2)处理不规则地分布的采样点，不存在***及方法。
进而，实现线性三角形分割插补是创造性的工作。例如，形成三角形时，哪一原采样点都不能存在于其他的三角形的内部。由于不规则，离散的数据的三角形分割是含糊的，因此，希望找到生成最大的具有最小角度的三角形的最优的三角形分割方法。另外，采样数据点的三角形分割易受到各变量中所使用的单元的影响。甚至以最优的三角形分割法的德洛涅的方法，根据情况，有时生成不同的三角形分割的组。
一种为了重建实施方式涉及的PET图像而将PET数据进行插补的PET数据处理方法，包括：收集通过使用PET扫描仪将被检体进行扫描而得到的PET现象数据，根据上述收集到的PET现象数据生成原正弦图，确定连接上述原正弦图内的采样点的多个三角形，使用上述确定的多个三角形，根据上述原正弦图确定均匀地采样后的正弦图。
能够提供在对原正弦图进行插补时为找到所希望的采样点的最近的点的***及方法。
与附图一起考察，由于通过参照以下的详细的说明加深对本发明的理解，因此，大概易于更完全地理解本发明以及附属的许多诸优点。
图1是表示具有不是存在于矩形的格子上，而存在于菱形的采样格子上的采样点的原LOR正弦图的图。
图2是表示用于得到插补正弦图的均匀地采样后的格子的图。
图3是表示针对在中央径向槽中均匀地分布的、被配置在摄像区域的中心的点源(小型磁盘)的2D正弦图，将线性插补与三角形分割插补进行比较的图。
图4是表示按照本说明书所记载的一实施方式形成的三角形的实例的图。
图5是表示按照本说明书所记载的一实施方式，根据图1所示的正弦图进行三角形分割的结果的图。
图6是表示使用本说明书所记载的一实施方式的某一组件形成的三角形的图。
图7是表示使用本说明书所记载的一实施方式的另一组件形成的三角形的图。
图8是表示按照本说明书所记载的代替实施方式，根据图1所示的正弦图进行三角形分割的结果的图。
图9是表示按照本说明书所记载的另一代替实施方式，根据图1所示的正弦图进行三角形分割的结果的图。
图10A是表示本说明书所记载的一实施方式涉及的方法的步骤的图。
图10B是表示本说明书所记载的一实施方式涉及的方法的步骤的图。
图10C是表示本说明书所记载的一实施方式涉及的方法的步骤的图。
图10D是表示本说明书所记载的一实施方式涉及的方法的步骤的图。
图11是表示为了按照本说明书所记载的实施方式得到PET现象信息而能够使用的伽马射线检测***的图。
具体实施方式
在一实施方式中，为了重建PET图像而插补PET数据的方法包含：(1)收集通过使用PET扫描仪来将被检体进行扫描而得到的PET现象数据的步骤；(2)根据收集到的PET现象数据生成原正弦图的步骤；(3)确定连接该原正弦图内的采样点的多个三角形的步骤，(4)使用这些被确定的多个三角形，根据原正弦图确定均匀地采样后的正 弦图的步骤。
在另一实施方式中，生成原正弦图的步骤包括：确定存储对于原正弦图内的各采样点的矢径值的第1列表的步骤；确定存储对于原正弦图内的各采样点的角度值的第2列表的步骤，确定多个三角形的上述步骤包括，为了确定这些多个三角形，在原正弦图内，将相邻的采样点在水平方向与每列进行连接的步骤。
在另一实施方式中，生成原正弦图的步骤包括：(1)将原正弦图内的第1采样点，即，使用一对指数(irad，iphi)而从各第1以及第2列表得到矢径值以及角度值的第1采样点，选择为第1三角形点的步骤；(2)将第2三角形点确定为使用指数(irad+1，iphi)而从各第1以及第2列表得到的原采样点的步骤；(3)将第3三角形点确定为使用指数(irad+2，iphi)而从各第1以及第2列表得到的原采样点的步骤；(4)将这些各第1、第2、以及第3三角形点作为上述多个三角形中的某个三角形的各顶点进行存储的步骤。
另外，根据另一实施方式，确定多个三角形的步骤包括，在irad指数值的某个范围内，以固定iphi指数的状态，对于通过使指数irad增加1而选择的多个第1原采样点，重复上述的选择步骤、第1确定步骤、第2确定步骤、以及存储步骤的步骤，与在iphi指数值的某个范围内重复上述的重复步骤的步骤。
在另一实施方式中，确定均匀地采样后的正弦图的步骤包括：(1)对均匀地采样后的正弦图的每个均匀的采样点，确定上述多个三角形中的对应的三角形的步骤；(2)针对上述多个三角形的各三角形，使用该三角形的各顶点的数据值(例如，现象的度数)来实施使用了重心坐标的插补方法，并对于与该三角形对应的均匀采样点生成插补数据值的步骤。
另外，在另一实施方式中，上述的方法包含使用均匀地采样后的正弦图以及插补数据值来实施重建，生成PET图像的步骤。
在另一实施方式中，为了重建图像而将断层法数据进行插补的方法包括：(1)收集通过使用扫描仪将被检体进行扫描而得到的现象数 据的步骤；(2)根据收集到的现象数据生成原正弦图的步骤；(3)为了确定连接该原正弦图内的采样点的多个三角形，确定在该原正弦图内相邻的采样点在水平方向与每列连接的多个三角形的步骤；(4)使用这些被确定的多个三角形，根据原正弦图确定均匀地采样后的正弦图的步骤。
上述方法之后还详细地进行说明，也能够扩展到三维空间。
在具备有模块样式的多边形的环体的PET***中，原正弦图的采样点不能完全地自由设定。没有使用如德洛涅三角形分割那样的方法，而使用根据晶体指数的直接三角形分割法。详细而言，该方法是连接相邻的采样点，在水平方向在每列构成三角形的方法，其结果构成图5所示的三角形。该方法大幅度地简化了三角形分割过程及其实施方式。另外，以这样的方式直接形成三角形时，为了使三角形的最小角度最大化，使用采样图案的信息。
接着，将三角形分割插补法的一实施方式与伪码相关联地进行说明。在具有多个环体的所提供的PET检测器***中，设2D原正弦图坐标在每个环体为irad以及iphi。这些是整数，irad＝1，2，...，nrad、iphi＝1，2，...，nphi。nrad以及nphi是irad以及iphi中的指数的总数。设S[irad，iphi]是具备了具有指数irad以及iphi的所提供的原采样点的s值的列表。同样地，设P[irad，iphi]是具备了具有指数irad以及iphi的所提供的原采样点的p值的列表。进而，设一个三角形具有3个顶点(x(1)，y(1))，(x(2)，y(2))，以及(x(3)，y(3))。
在一实施方式中，为了生成三角形，使用具有与Rad相关的两个内环以及与Phi相关的一个外环的下述的伪码。
在iphi＝1(角度为0°附近)或者iphi＝nphi(角度为180°附近)的Phi的边界，三角形如iphi＝1与iphi＝nphi相互彼此相邻地形成。应用了这样的结束条件，但为了简化，省略图示。
图4表示使用之前所示的伪码形成的三角形。具体而言，上面的两个三角形通过第1内环生成，下面的两个三角形通过第2内环生成。另外，图4的底部所示的两个三角形分别通过之前所示的第2内环的各第1以及第2框(block)而生成。对于指数iphi的某个固定值，伴随着irad指数增加，要注意P(irad，iphi)值从左向右在两个值相互间切换。从而，如图4的三角形所示，P(irad，iphi)＝P(irad+2，iphi)。
例如，图4的左下方的三角形根据下述的伪码(之前所示的第2内框)生成。
x(1)＝S(irad，iphi)；
x(2)＝S(irad+2，iphi)；
x(3)＝S(irad+1，iphi+1)；
y(1)＝P(irad，iphi)；
y(2)＝P(irad+2，iphi)；
y(3)＝P(irad+1，iphi+1)；
具体而言，在上述的伪码中，当irad＝ir且iphi＝jp时，通过添加了图4的(ir，jp)的标记的点，来给出三角形的第1顶点(x(1)，y(1))。另外，根据伪码，第2顶点(x(2)，y(2))通过添加了图4的(ir+2，jp)的标记的点而给出，第3的顶点(x(3)，y(3))通过添加了图4的(ir+1，jp+1)的标记的点而给出。
另外，图4的右下方的三角形根据下述的伪码(之前所示的第3内框)生成。
x(1)＝S(irad+2，iphi)；
x(2)＝S(irad+1，iphi+1)；
x(3)＝S(irad+3，iphi+1)；
y(1)＝P(irad+2，iphi)；
y(2)＝P(irad+1，iphi+1)；
y(3)＝P(irad+3，iphi+1)；
具体而言，在上述的伪码中，当irad＝ir且iphi＝jp时，通过添 加了图4的(ir+2，jp)的标记的点，来给出三角形的第1顶点(x(1)，y(1))。进而，根据伪码，第2顶点(x(2)，y(2))通过添加了图4的(ir+1，jp+1)的标记的点而给出，第3顶点(x(3)，y(3))通过添加了图4的(ir+3，jp+1)的标记的点而给出。
图1表示作为对之前所示的伪码的输入而提供的原采样空间的仅仅的一部分。图5表示对于相同的采样空间区域，进行上述的方法时的三角形分割的结果。
图6表示之前所示的伪码的、只执行第1内环时的结果，图7表示只执行第2内环时的结果。
在使用下述的伪码生成一部分的三角形的代替实施方式中，生成图8所示的三角形。
在使用下述的伪码生成一部分的三角形的另一代替实施方式中，生成图9所示的三角形。
另外，这样得到的几个图8以及图9的三角形变为小的三角形，这是不怎么希望的结果。
在3D收集中，3D原正弦图具有(Rad，Phi，ringCombination) 的坐标。关联的插补正弦图具有四个维度(s，p，z，tilt)。此时，z是LOR的轴方向的位置的平均，tilt是LOR对于横断面的轴倾斜角度。3D LOR正弦图的ringCombination在插补正弦图中，与(z，tilt)相关。
在本说明书中所公开的三角形分割插补法的诸实施方式中，具有若干优点。
第1、它们为找到所希望的采样点的最近的点提供了***及方法，能够将它们应用于离散的数据分布。
第2、上述的方法是局部的方法，只使用最近的点，来确定所希望的采样点中的插补值。对此，由于附近的槽内部的原采样点全部对于所希望的采样点有益，因此，线性插补法是广域的(相对而言)的方法。称为局部的方法的上述的方法的性质对提高活动密度的边界附近的插补结果起作用。例如，一个极端的事例是：配置在摄像区域(field-of-view)的中心的点源(小型磁盘)的2D正弦图。如图3所示，这是在径向槽的中央均匀地分布的线。图3的上侧所示的线性插补法归结于非常不均匀地分布，但如图3的下侧所示，三角形分割插补提供了极为完美的结果。
第3、实现公开的三角形分割法的实施方式的方法简单明了，易于实现，带来希望的成果。
图10A～10D表示为了按照一实施方式重建PET图像，而对PET数据进行插补的方法的诸步骤。
图10A综合地表示本实施方式的方法。
在步骤201中，获得通过使用PET扫描仪对被检体进行扫描而得到的PET现象数据。PET现象数据以前接着PET扫描而存储在存储器中，并从该存储器中取出。如相关申请第12/571,562号所说明的那样，现象数据能够编成适合清单模式的重建的清单(list)。该申请的内容通过在本说明书中参照而进行援引。
在步骤202中，根据收集到的PET现象数据，生成原同时计数线(Line-Of-Response：LOR)正弦图。如之前所叙述的那样，二维的 情况下，生成原LOR正弦图的步骤包含，确定存储对于原LOR正弦图内的各原采样点的矢径值的第1列表S[irad，iphi]，和确定存储对于原LOR正弦图内的各原采样点的角度值的第2列表P[irad，iphi]。原采样点的矢径值以及角度值使用一对指数(irad，iphi)，从各第1以及第2列表得到。
在步骤203中，使用如图10B～10D所示的那样的各种方法，从原采样点生成多个三角形。
在步骤208中，使用这些被确定的多个三角形，根据原LOR正弦图确定均匀地采样的正弦图。具体而言，步骤208包含：对均匀地采样的正弦图的每个均匀采样点，选择上述多个三角形中的对应的三角形的步骤，和对上述多个三角形中的每个该三角形，使用该三角形的各顶点的数据值实施使用了重心坐标的插补方法，对于与该三角形对应的均匀采样点生成插补数据值的步骤。
在步骤209中，为了生成PET图像，使用均匀地采样后的正弦图以及插补数据值，实施重建。
图10B表示生成上述多个三角形的方法。
在步骤301中，根据各第1以及第2列表的大小，将参数nphi以及nrad进行初始化。指数iphi也初始化为1。
在步骤302中，使用第1方法，即、使用上述的第1内环来生成三角形组，并进行存储。参照图10C的步骤。
在步骤303中，使用第2方法，即、使用上述的第2内环来生成另一三角形组，并进行存储。参照图10D的步骤。
在步骤304中，当iphi＝nphi-1时结束生成三角形的过程，通常，本方法进入步骤305。因此，iphi指数增加1，本过程再次进入步骤302。
图10C表示按照之前所示的伪码的第1内环生成三角形的第1方法。
在步骤401中，irad指数初始化为1。
在步骤402中，将原LOR正弦图内的第1原采样点，即，使用 一对指数(irad，iphi)而从各第1以及第2列表得到矢径值以及角度值的第1原采样点选择为第1三角形点。
在步骤403中，第2三角形点被确定为使用指数(irad+1，iphi)而从各第1以及第2列表得到的原采样点。
在步骤404中，第3三角形点被确定为使用指数(irad+2，iphi)而从各第1以及第2列表得到的原采样点。
在步骤405中，将这些各第1、第2、以及第3三角形点作为上述多个三角形中的某个三角形的各顶点存储在存储器中。
在步骤406中，如果irad与nrad-2相等，则结束本过程。否则，在步骤407中将irad指数增加1，本过程再次进入步骤402。
图10D是表示按照之前所示的伪码的第2内环生成三角形的第2方法。
在步骤501中，将irad指数初始化为1。
在步骤502中，将原LOR正弦图内的第1原采样点，即，使用一对指数(irad，iphi)从各第1以及第2列表得到矢径值以及角度值的第1原采样点选择为第1三角形点。
在步骤503中，将第2三角形点确定为使用指数(irad+2，iphi)而从各第1以及第2列表得到的原采样点。
在步骤504中，将第3三角形点确定为使用指数(irad+1，iphi+1)而从各第1以及第2列表得到的原采样点。
在步骤505中，将这些各第1、第2、以及第3三角形点作为上述多个三角形中的某个三角形的各顶点存储在存储器中。
在步骤506中，将原LOR正弦图内的第1原采样点，即，使用一对指数(irad+2，iphi)而从各第1以及第2列表得到矢径值以及角度值的第1原采样点选择为第1三角形点。
在步骤507中，将第2三角形点确定为使用指数(irad+1，iphi+1)而从各第1以及第2列表得到的原采样点。
在步骤508中，将第3三角形点确定为使用指数(irad+3，iphi+1)而从各第1以及第2列表得到的原采样点。
在步骤509中，将这些各第1、第2、以及第3三角形点作为上述多个三角形中的某个三角形的各顶点存储在存储器中。
在步骤510中，如果irad为nrad-3以上，则结束本过程。否则，在步骤511将irad指数增加2，本过程再次进入步骤502。
图11是按照本实施方式得到伽马射线或者PET现象信息而能够使用的伽马射线检测***的示意图。在图11中，在光导130上配置有光电倍增管135以及140，在光导130下配置有第1闪烁晶体的阵列105。另外，第2闪烁晶体的阵列125与闪烁晶体105对置地配设，在之上配置有光导115和光电倍增管195以及110。由光电倍增管、光导、以及闪烁晶体能够形成检测器模块，该图的伽马射线检测***具备有内设于环体的多个检测器模块。
在图11中，从被检体(省略图示)放射出伽马射线时，这些伽马射线在相差约180°，向反方向上移动。由闪烁晶体100以及120同时进行伽马射线检测，当在规定的限制时间内由闪烁晶体100以及120检测这些伽马射线时，确定闪烁现象。从而，检测伽马射线的定时的***由闪烁晶体100以及120同时检测伽马射线。其中，为了全面简化，伽马射线检测与闪烁晶体100相关联地进行了说明。然而，针对在本说明书中形成的闪烁晶体100的说明，同样适用于在闪烁晶体120的伽马射线检测，如果是本领域的技术人员，会注意到这一点。
各光电倍增管110、135、140、以及195分别与数据收集单元150连接。数据收集单元具备有为了处理来自光电倍增管的信号而构成的硬件。数据收集单元150测定伽马射线的到达时间。将对于***时钟(省略图示)的判别脉冲的时间进行编码化的数据收集单元150生成两个输出(一个组合光电倍增管135/140，一个组合光电倍增管110/195)。在飞行时间PET***用中，数据收集单元150典型地以15～25ps的精度生成时间戳。数据收集单元测定各PMT的信号(来自数据收集单元150的4个输出)的振幅。
数据收集单元150的输出被提供给CPU(Central ProcessingUnit)170，在此进行处理。该处理有时包含由根据数据收集单元输出 来推定能量以及位置，同时根据对每个现象由时间戳输出推定到达时间而成，并根据以前的校正来应用大量的校正步骤，提高能量、位置、以及时间的推定精度。如是本领域的技术人员将注意到那样，CPU170能够作为个别逻辑门、面向特定用途的集成电路(Application SpecificIntegrated Circuit：ASIC)、现场可编程逻辑门阵列(FieldProgrammable Gate Array：FPGA)、或者其他的复杂可编程逻辑器件(Complex Programmable Logic Device：CPLD)来实现。FPGA或者CPLD的实施方式能够通过VHDL、Verilog、或者其他的任意的硬件描述语言来来进行编码化，该编码能够存储在直接搭载于该FPGA或者CPLD内部的电子存储器或者作为独立的电子存储器的电子存储器中。另外，电子存储器能够使ROM(Read Only Memory)、EPROM(Erasable Programmable Read Only Memory)、EEPROM(Electronically Erasable and Programmable Read Only Memory)、或者闪存等非易失性的存储器。另外，电子存储器也能够是静态RAM(Random Access Memory)或者动态RAM等易失性的存储器，为了进行电子存储器的管理以及FPGA或者CPLD与电子存储器之间的相互作用的管理，也可以设置微控制器或者微处理器那样的处理器。
或者，CPU170也可以执行包含实施本说明书所记载的功能的一组计算机可读的命令的计算机程序，即，也可以执行存储于上述的非暂时的电子存储器及/或硬盘驱动器、CD(Compact Disc)、DVD(Digital Versatile Disc)、闪盘、或者其他的任意的公知的存储媒体中的任一个的程序。另外，计算机可读的命令能够提供基于美国因特尔公司的Xeon(注册商标)处理器或者基于美国AMD公司的Opteron(注册商标)处理器等处理器、以及MicrosoftVISTA(注册商标)、UNIX(注册商标)、Solaris(注册商标)、LINUX(注册商标)、Apple(注册商标)、MAC-OS(注册商标)、以及本领域的技术人员公知的其他的操作***等操作***一起工作的、应用程序、背景程式、或者操作***的构成要素、或者它们的组合。
一旦通过CPU170来进行处理，则该处理完成的信号被存储在电 子存储部180，且/或显示在显示器145。如是本领域的技术人员则会注意到那样，电子存储部180能够作为硬盘驱动器、CD-ROM驱动器、DVD驱动器、闪盘、RAM、ROM、或者在该技术领域公知的其他的任意的电子存储部。显示器145能够作为LCD(Liquid CrystalDisplay)显示器、CRT(Cathode-Ray Tube)显示器、等离子显示器、OLED(Organic Light Emitted Diode)、LED(Light Emitted Diode)、或者在该技术领域公知的其他的任意的显示器来实现。从而，在本说明书中形成的电子存储部180以及显示器145的说明仅仅是示例，绝不是限定本进步的技术的范围。
另外，图11也具备接口175，伽马射线检测***经由其进行其他的外部装置以及/或用户的互动。例如，接口175可以作为USB(Universal Serial Bus)接口、PCMCIA(Personal Computer MemoryCard International Association)接口、以太网(注册商标)接口、或者本技术领域公知的其他的任意的接口。接口175也能够是有线式或者无线式，也可以包含用于与键盘以及/或者鼠标或者用户进行对话的该技术领域公知的其他的人机接口。
虽然说明了本发明的几个实施方式，但这些实施方式是作为例子而提示的，并不意图限定本发明的范围。实际上，本说明书记载的新颖的方法及***的能够以其他的各种形态进行实施。并且，在不脱离本发明的精神的范围内，本说明书记载的方法及***的形式能够进行各种的省略、置换、变更。后附的权利要求书及其等同内容，意味着作为对象而包含属于本发明的范围及精神的方式或变形例。
Claims (10)
1.一种为了重建PET图像而对PET数据进行插补的PET数据处理方法，其特征在于，包括：
收集通过使用PET扫描仪对被检体进行扫描而得到的PET现象数据，
根据上述收集到的PET现象数据生成原正弦图，
确定连接上述原正弦图内的采样点的多个三角形，
使用上述被确定的多个三角形，根据上述原正弦图来确定均匀地采样后的正弦图，
上述原正弦图的生成包括：
确定存储对于上述原正弦图内的各采样点的矢径值的第1列表；以及
确定存储对于上述原正弦图内的各采样点的角度值的第2列表，
上述多个三角形的确定包含：在上述原正弦图内，将相邻的采样点在水平方向针对每列连接，确定上述多个三角形。
2.根据权利要求1所述的PET数据处理方法，其特征在于，
上述多个三角形的确定包括：
将上述原正弦图内的第1采样点，即，使用一对指数(irad，iphi)从各上述第1以及第2列表得到矢径值以及角度值的第1采样点，选择为第1三角形点；
将第2三角形点确定为使用指数(irad+1，iphi)而从各上述第1以及第2列表得到的采样点；
将第3三角形点确定为使用指数(irad+2，iphi)而从各上述第1以及第2列表得到的采样点；以及
将各上述第1、第2、以及第3三角形点作为上述多个三角形中的某个三角形的各顶点而进行存储。
3.根据权利要求1所述的PET数据处理方法，其特征在于，
上述多个三角形的确定包括：
将上述原正弦图内的第1采样点，即，使用一对指数(irad，iphi)而从各上述第1以及第2列表得到矢径值以及角度值的第1采样点，选择为第1三角形点；
将第2三角形点确定为使用指数(irad+2，iphi)而从各上述第1以及第2列表得到的采样点；
将第3三角形点确定为使用指数(irad+1，iphi+1)而从各上述第1以及第2列表得到的采样点；以及
将各上述第1、第2、以及第3三角形点作为上述多个三角形中的某个三角形的各顶点而进行存储。
4.根据权利要求1所述的PET数据处理方法，其特征在于，
上述多个三角形的确定包括：
将上述原正弦图内的第1采样点，即，使用一对指数(irad+2，iphi)而从各上述第1以及第2列表得到矢径值以及角度值的第1采样点，选择为第1三角形点；
将第2三角形点确定为使用指数(irad+1，iphi+1)而从各上述第1以及第2列表得到的采样点；
将第3三角形点确定为使用指数(irad+3，iphi+1)而从各上述第1以及第2列表得到的采样点；以及
将各上述第1、第2、以及第3三角形点作为上述多个三角形中的某个三角形的各顶点进行存储。
5.根据权利要求2所述的PET数据处理方法，其特征在于，
上述多个三角形的确定包括：
以在irad指数值的某个范围内固定iphi指数的状态，对于通过使irad指数增加1而选择的多个第1采样点，重复上述第1三角形点的选择、上述第2三角形点的确定、上述第3三角形点的确定、以及各三角形点的存储，
在iphi指数值的某个范围内，重复上述重复。
6.根据权利要求1所述的PET数据处理方法，其特征在于，
均匀地采样后的正弦图的确定包括：
对上述均匀地采样后的正弦图的每个均匀采样点，选择上述多个三角形中的对应的三角形；以及
针对上述多个三角形的各三角形，使用上述三角形的各顶点的数据值来实施使用了重心坐标的插补方法，对于与上述三角形对应的上述均匀采样点生成插补数据值。
7.根据权利要求1所述的PET数据处理方法，其特征在于，
为了生成上述PET图像，也包含使用上述均匀地采样后的正弦图以及上述插补数据值来实施重建。
8.一种为了重建PET图像而将PET数据进行插补的PET数据处理装置，其特征在于，具备：
存储部，存储通过使用PET扫描仪将被检体进行扫描而得到的PET现象数据；
生成部，根据上述存储的PET现象数据生成原正弦图；
第1确定部，确定连接上述原正弦图内的采样点的多个三角形；
第2确定部，使用上述确定的多个三角形，根据上述原正弦图确定均匀地采样后的正弦图；以及
插补部，对上述均匀地采样后的正弦图的每个均匀采样点，选择上述多个三角形中的对应的三角形，并针对上述多个三角形的各三角形，使用上述三角形的各顶点的数据值来实施使用了重心坐标的插补方法，并对于与上述三角形对应的上述均匀采样点生成插补数据值。
9.根据权利要求8所述的PET数据处理装置，其特征在于
还具备重建部，为了生成上述PET图像，使用上述均匀地采样后的正弦图以及上述插补数据值来实施重建。
10.一种插补用于图像重建的断层数据的数据处理方法，其特征在于，包括：
收集通过使用扫描仪将被检体进行扫描而得到的现象数据；
根据上述收集到的现象数据生成原正弦图；
为了确定连接上述原正弦图内的采样点的多个三角形，确定在上述原正弦图内将相邻的采样点在水平方向针对每列连接的多个三角形；以及
使用上述被确定的多个三角形，根据上述原正弦图确定均匀地采样后的正弦图。
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US12/958,762 | 2010-12-02 | ||
| US12/958,762 US8503753B2 (en) | 2010-12-02 | 2010-12-02 | System and method for triangular interpolation in image reconstruction for PET |
| PCT/JP2011/077722 WO2012074028A1 (ja) | 2010-12-02 | 2011-11-30 | Ｐｅｔデータ処理方法、ｐｅｔデータ処理装置、コンピュータ読み取り可能な記録媒体、及びデータ処理方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102640016A CN102640016A (zh) | 2012-08-15 |
| CN102640016B true CN102640016B (zh) | 2015-08-19 |
Family
ID=46162283
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201180004335.5A Active CN102640016B (zh) | 2010-12-02 | 2011-11-30 | Pet数据处理方法、pet数据处理装置以及数据处理方法 |
Country Status (5)
| Country | Link |
|---|---|
| US (1) | US8503753B2 (zh) |
| EP (1) | EP2629117B1 (zh) |
| JP (1) | JP6009761B2 (zh) |
| CN (1) | CN102640016B (zh) |
| WO (1) | WO2012074028A1 (zh) |
Families Citing this family (28)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| SE533704C2 (sv) | 2008-12-05 | 2010-12-07 | Flatfrog Lab Ab | Pekkänslig apparat och förfarande för drivning av densamma |
| US8718973B2 (en) * | 2011-09-09 | 2014-05-06 | Kabushiki Kaisha Toshiba | Method, device, and system for calculating a geometric system model using an area-simulating-volume algorithm in three dimensional reconstruction |
| US10168835B2 (en) | 2012-05-23 | 2019-01-01 | Flatfrog Laboratories Ab | Spatial resolution in touch displays |
| KR20140025922A (ko) * | 2012-08-23 | 2014-03-05 | 삼성전기주식회사 | 터치패널 |
| WO2014168567A1 (en) * | 2013-04-11 | 2014-10-16 | Flatfrog Laboratories Ab | Tomographic processing for touch detection |
| CN103565459B (zh) * | 2013-09-25 | 2015-07-01 | 沈阳东软医疗***有限公司 | 一种判断成像质量方法和装置 |
| US10126882B2 (en) | 2014-01-16 | 2018-11-13 | Flatfrog Laboratories Ab | TIR-based optical touch systems of projection-type |
| US10146376B2 (en) | 2014-01-16 | 2018-12-04 | Flatfrog Laboratories Ab | Light coupling in TIR-based optical touch systems |
| JP2017512997A (ja) | 2014-03-28 | 2017-05-25 | コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． | 陽電子放射断層撮影（ｐｅｔ）における欠落画素の補償 |
| CN106233336B (zh) * | 2014-04-23 | 2020-11-03 | 皇家飞利浦有限公司 | Pet随机符合校正 |
| WO2015199602A1 (en) | 2014-06-27 | 2015-12-30 | Flatfrog Laboratories Ab | Detection of surface contamination |
| WO2016122385A1 (en) | 2015-01-28 | 2016-08-04 | Flatfrog Laboratories Ab | Dynamic touch quarantine frames |
| US10318074B2 (en) | 2015-01-30 | 2019-06-11 | Flatfrog Laboratories Ab | Touch-sensing OLED display with tilted emitters |
| EP3537269A1 (en) | 2015-02-09 | 2019-09-11 | FlatFrog Laboratories AB | Optical touch system |
| US10401546B2 (en) | 2015-03-02 | 2019-09-03 | Flatfrog Laboratories Ab | Optical component for light coupling |
| US10775937B2 (en) | 2015-12-09 | 2020-09-15 | Flatfrog Laboratories Ab | Stylus identification |
| US10761657B2 (en) | 2016-11-24 | 2020-09-01 | Flatfrog Laboratories Ab | Automatic optimisation of touch signal |
| EP3667475B1 (en) | 2016-12-07 | 2022-09-07 | FlatFrog Laboratories AB | A curved touch device |
| CN110300950B (zh) | 2017-02-06 | 2023-06-16 | 平蛙实验室股份公司 | 触摸感测***中的光学耦合 |
| WO2018174788A1 (en) | 2017-03-22 | 2018-09-27 | Flatfrog Laboratories | Object characterisation for touch displays |
| EP4036697A1 (en) | 2017-03-28 | 2022-08-03 | FlatFrog Laboratories AB | Optical touch sensing apparatus |
| WO2019045629A1 (en) | 2017-09-01 | 2019-03-07 | Flatfrog Laboratories Ab | IMPROVED OPTICAL COMPONENT |
| WO2019172826A1 (en) | 2018-03-05 | 2019-09-12 | Flatfrog Laboratories Ab | Improved touch-sensing apparatus |
| CN109636873B (zh) * | 2018-12-12 | 2023-07-18 | 上海联影医疗科技股份有限公司 | 用于医学图像重建的数据处理方法和医学图像重建方法 |
| WO2020153890A1 (en) | 2019-01-25 | 2020-07-30 | Flatfrog Laboratories Ab | A videoconferencing terminal and method of operating the same |
| CN109829869B (zh) * | 2019-03-04 | 2020-11-13 | 广东明峰医疗科技有限公司 | 一种降低pet图像噪声的数值域插值方法 |
| WO2021162602A1 (en) | 2020-02-10 | 2021-08-19 | Flatfrog Laboratories Ab | Improved touch-sensing apparatus |
| CN112907664B (zh) * | 2021-02-08 | 2023-04-25 | 明峰医疗***股份有限公司 | 图像像素位置查找方法、***及计算机可读存储介质 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101034479A (zh) * | 2006-03-10 | 2007-09-12 | Ge医疗***环球技术有限公司 | 图像重建方法和x射线ct设备 |
| CN101267768A (zh) * | 2005-07-22 | 2008-09-17 | 断层放疗公司 | 对接受放射疗法的患者的呼吸时相进行检测的***和方法 |
| CN101467180A (zh) * | 2006-04-27 | 2009-06-24 | 马科外科公司 | 轮廓三角测量***和方法 |
Family Cites Families (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6002738A (en) * | 1995-07-07 | 1999-12-14 | Silicon Graphics, Inc. | System and method of performing tomographic reconstruction and volume rendering using texture mapping |
| WO2000011603A2 (en) * | 1998-08-20 | 2000-03-02 | Apple Computer, Inc. | Graphics processor with pipeline state storage and retrieval |
| US6625249B1 (en) * | 2002-09-20 | 2003-09-23 | General Electric Company | Method, system and program product for tomographic backprojection using multi-color rendering engine |
| ATE439653T1 (de) * | 2003-09-09 | 2009-08-15 | Univ Illinois | Schnelle hierarchische tomographieverfahren und - vorrichtungen |
| US7444011B2 (en) * | 2004-02-10 | 2008-10-28 | University Of Chicago | Imaging system performing substantially exact reconstruction and using non-traditional trajectories |
| GB0512252D0 (en) * | 2005-06-16 | 2005-07-27 | Hammersmith Imanet Ltd | Method of and software for conducting motion correction in tomographic scanning and system for tomographic scanning using the method |
| US7848559B2 (en) * | 2006-05-17 | 2010-12-07 | Siemens Medical Solutions Usa, Inc. | Discrete axial re-binning of time-of-flight positron emission tomography data |
| GB2450073B (en) * | 2006-08-25 | 2009-11-04 | Siemens Molecular Imaging Ltd | Regional reconstruction of spatially distributed functions |
| US20080095465A1 (en) * | 2006-10-18 | 2008-04-24 | General Electric Company | Image registration system and method |
| US20110142316A1 (en) * | 2009-10-29 | 2011-06-16 | Ge Wang | Tomography-Based and MRI-Based Imaging Systems |
| TW201203052A (en) * | 2010-05-03 | 2012-01-16 | Flatfrog Lab Ab | Touch determination by tomographic reconstruction |
-
2010
- 2010-12-02 US US12/958,762 patent/US8503753B2/en active Active
-
2011
- 2011-11-30 WO PCT/JP2011/077722 patent/WO2012074028A1/ja active Application Filing
- 2011-11-30 CN CN201180004335.5A patent/CN102640016B/zh active Active
- 2011-11-30 JP JP2011262973A patent/JP6009761B2/ja active Active
- 2011-11-30 EP EP11844146.8A patent/EP2629117B1/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101267768A (zh) * | 2005-07-22 | 2008-09-17 | 断层放疗公司 | 对接受放射疗法的患者的呼吸时相进行检测的***和方法 |
| CN101034479A (zh) * | 2006-03-10 | 2007-09-12 | Ge医疗***环球技术有限公司 | 图像重建方法和x射线ct设备 |
| CN101467180A (zh) * | 2006-04-27 | 2009-06-24 | 马科外科公司 | 轮廓三角测量***和方法 |
Non-Patent Citations (3)
| Title |
|---|
| Gap Filling Strategies for 3-D-FBP Reconstructions of High-Resolution Research Tomograph Scans;Floris H. P. van Velden等;《IEEE TRANSACTIONS ON MEDICAL IMAGING》;20080731;第27卷(第7期);第934页 * |
| On interpolation and approximation by bivariate quintic splines with boundary conditions;Huanwen Liu等;《Computer-Aided Design and Computer Graphics》;20090830;第274页 * |
| Scattered data interpolation methods for electronic imaging systems: a survey;Isaac Amidror等;《Journal of Electronic Imaging》;20020430;第11卷(第2期);第157、159页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| US8503753B2 (en) | 2013-08-06 |
| EP2629117A1 (en) | 2013-08-21 |
| JP2012118079A (ja) | 2012-06-21 |
| EP2629117B1 (en) | 2019-10-23 |
| WO2012074028A1 (ja) | 2012-06-07 |
| JP6009761B2 (ja) | 2016-10-19 |
| US20120141001A1 (en) | 2012-06-07 |
| CN102640016A (zh) | 2012-08-15 |
| EP2629117A4 (en) | 2017-01-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102640016B (zh) | Pet数据处理方法、pet数据处理装置以及数据处理方法 | |
| US6674083B2 (en) | Positron emission tomography apparatus | |
| CN1272638C (zh) | Pet装置及pet装置的图像生成方法 | |
| US10216701B2 (en) | Image-based point-spread-function modelling in time-of-flight positron-emission-tomography iterative list-mode reconstruction | |
| CN102755172B (zh) | 核医学成像方法以及核医学成像装置 | |
| US9291725B2 (en) | Random coincidence reduction in positron emission tomography using tangential time-of-flight mask | |
| US8558181B2 (en) | Positron emission tomography system with hybrid detection geometries and sampling | |
| CN114666495B (zh) | 高分辨率康普顿相机成像方法、装置、电子设备及介质 | |
| CN102540239B (zh) | γ射线检测***以及正电子发射断层摄影*** | |
| US20140209804A1 (en) | Apparatuses for and methods of generating images | |
| US8483992B2 (en) | Analytical computation of geometric normalization factor for TOF-PET | |
| US9241678B2 (en) | Random estimation in positron emission tomography with tangential time-of-flight mask | |
| JP2005043104A (ja) | 放射線位置検出器の校正方法 | |
| US8866086B2 (en) | PET scanner with emission and transmission structures in a checkerboard configuration | |
| Jin et al. | Experimental evaluation of a 3-D CZT imaging spectrometer for potential use in Compton-Enhanced PET imaging | |
| WO1997008569A1 (en) | An imaging apparatus | |
| US8787644B2 (en) | Method and device for calculating voxels defining a tube-of-response using a central-ray-filling algorithm | |
| CN111080737A (zh) | 图像重建方法、装置及pet扫描*** | |
| Stefaan et al. | State of the art in total body PET | |
| Schug et al. | Fast and unbiased 3D calibration method of arbitrary scintillator based PET detectors | |
| Freire et al. | Calibration of PET Detectors Based on Monolithic Blocks Using Voronoi Diagrams | |
| Samanta | Enhancing PET Functionalities with Novel Geometries and Multimodal Imaging Techniques | |
| Johnson | Development of a Small-Animal SPECT System with a High-Purity Germanium Detector | |
| CN110368011A (zh) | 用于计算针对晶间事件的响应线的pet成像***和方法 | |
| Chiesa et al. | Advances in nuclear medicine imaging in oncology |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | 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 | ||
| C41 | Transfer of patent application or patent right or utility model | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20160803 Address after: Japan Tochigi Patentee after: Toshiba Medical System Co., Ltd. Address before: Tokyo, Japan, Japan Patentee before: Toshiba Corp Patentee before: Toshiba Medical System Co., Ltd. |