CN110687585B - 获取晶体效率的方法、装置、计算机设备和存储介质 - Google Patents

获取晶体效率的方法、装置、计算机设备和存储介质 Download PDF

Info

Publication number
CN110687585B
CN110687585B CN201910899826.9A CN201910899826A CN110687585B CN 110687585 B CN110687585 B CN 110687585B CN 201910899826 A CN201910899826 A CN 201910899826A CN 110687585 B CN110687585 B CN 110687585B
Authority
CN
China
Prior art keywords
event data
crystal
detector
coincidence
estimated
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
Application number
CN201910899826.9A
Other languages
English (en)
Other versions
CN110687585A (zh
Inventor
冯涛
何鎏春
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai United Imaging Healthcare Co Ltd
Original Assignee
Shanghai United Imaging Healthcare Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shanghai United Imaging Healthcare Co Ltd filed Critical Shanghai United Imaging Healthcare Co Ltd
Priority to CN201910899826.9A priority Critical patent/CN110687585B/zh
Publication of CN110687585A publication Critical patent/CN110687585A/zh
Application granted granted Critical
Publication of CN110687585B publication Critical patent/CN110687585B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T7/00—Details of radiation-measuring instruments
    • G01T7/005—Details of radiation-measuring instruments calibration techniques
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/02—Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03—Computerised tomographs
    • A61B6/037—Emission tomography
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/29—Measurement performed on radiation beams, e.g. position or section of the beam; Measurement of spatial distribution of radiation
    • G01T1/2914—Measurement of spatial distribution of radiation
    • G01T1/2985—In depth localisation, e.g. using positron emitters; Tomographic imaging (longitudinal and transverse section imaging; apparatus for radiation diagnosis sequentially in different planes, steroscopic radiation diagnosis)
    • 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

Abstract

本申请涉及一种获取晶体效率的方法、装置、计算机设备和存储介质。所述方法包括：获取探测器扫描得到的衰减图像、符合事件数据以及每个探测器晶体的测量单事件数据；根据衰减图像以及符合事件数据，得到每个探测器晶体的估计单事件数据；根据所述测量单事件数据和估计单事件数据，计算得到相应探测器晶体的晶体效率。通过临床数据能够精确的计算探测器晶体的晶体效率，并且能够减少放射源的使用次数，减少成本。

Description

获取晶体效率的方法、装置、计算机设备和存储介质
技术领域
本申请涉及医学图像校正技术领域，特别是涉及一种获取晶体效率的方法、装置、计算机设备和存储介质。
背景技术
随着医疗技术的发展，正电子发射型计算机断层显像技术(PET)是核医学领域比较先进的临床检查影像技术，为了能获取准确的PET重建图像，需要提前获取每个探测器晶体的晶体效率，获取错误的晶体效率可能会带来重建图像的伪影。
然而，目前的现有技术通常是使用固定形状的标准模体进行质量控制(QC，quality control)，利用标准模体进行质量控制，由于测量用的标准模体结构已知，再加上放射源的存在，能充分的模拟系统正常工作环境下的反应。但由于需要额外的放射源，并且在QC的过程中不能使用仪器，因此通过对标准模体进行QC获取晶体效率的方法成本较大，往往不能经常使用。
发明内容
基于此，有必要针对上述技术问题，提供一种能够节省成本的获取晶体效率的方法、装置、计算机设备和存储介质。
一种获取晶体效率的方法，所述方法包括：获取探测器扫描得到的衰减图像、符合事件数据以及每个探测器晶体的测量单事件数据；根据衰减图像以及符合事件数据，得到每个探测器晶体的估计单事件数据；根据所述测量单事件数据和估计单事件数据，计算得到相应探测器晶体的晶体效率。
在其中一个实施例中，根据衰减图像以及符合事件数据，得到每个探测器晶体的估计单事件数据包括：根据所述衰减图像以及符合事件数据进行蒙特卡罗模拟，得到每个探测器晶体的估计单事件数据。
在其中一个实施例中，所述根据所述衰减图像以及符合事件数据进行蒙特卡罗模拟，得到每个探测器晶体的估计单事件数据包括：根据所述衰减图像以及符合事件数据确定蒙特卡罗模拟边界以及初始的光子信息；根据所述蒙特卡罗模拟边界以及初始的光子信息进行光子运动模拟，得到每个探测器晶体的估计单事件数据。
在其中一个实施例中，根据所述测量单事件数据和估计单事件数据，计算得到相应探测器晶体的晶体效率包括：构建所述测量单事件数据、所述估计单事件数据与所述晶体效率的第一预设关系；根据所述第一预设关系、测量单事件数据以及估计单事件数据，计算得到相应探测器晶体的晶体效率。
在其中一个实施例中，所述第一预设关系基于最大似然法建立。
在其中一个实施例中，所述根据所述晶体效率方程、测量单事件数据以及估计单事件数据，得到相应探测器晶体的晶体效率包括：获取探测器扫描得到的符合事件数据，构建所述符合事件数据与活度图像的第二预设关系；所述第二预设关系与所述第一预设关系交替运算，得到相应探测器晶体的晶体效率。
在其中一个实施例中，所述获取探测器扫描得到的衰减图像、符合事件数据以及每个探测器晶体的测量单事件数据包括：对注射氟代脱氧葡萄糖的头部进行扫描，获取探测器扫描得到的衰减图像、符合事件数据以及每个探测器晶体的测量单事件数据。
一种获取晶体效率的装置，所述装置包括：获取模块，用于获取探测器扫描得到的衰减图像、符合事件数据以及每个探测器晶体的测量单事件数据；单事件数据计算模块，用于根据衰减图像以及符合事件数据，得到每个探测器晶体的估计单事件数据；晶体效率计算模块，用于根据所述测量单事件数据和估计单事件数据，计算得到相应探测器晶体的晶体效率。
一种计算机设备，包括存储器和处理器，所述存储器存储有计算机程序，所述处理器执行所述计算机程序时实现上述任一种所述方法的步骤。
一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现上述任一种所述的方法的步骤。
上述获取晶体效率的方法、装置、计算机设备和存储介质，通过临床数据的探测器扫描得到的衰减图像、符合事件数据以及每个探测器晶体的测量单事件数据。在根据衰减图像以及符合事件数据，得到每个探测器晶体的估计单事件数据。最后根据测量单事件数据以及估计单事件数据，计算得到相应探测器晶体的晶体效率。通过临床数据能够精确的计算探测器晶体的晶体效率，并且能够减少放射源的使用次数，减少成本。
附图说明
图1为一个实施例中获取晶体效率的方法流程示意图；
图2为一个实施例中根据单事件数据计算晶体效率的方法流程示意图；
图3为一个实施例中获取晶体效率的装置的结构框图；
图4为一个实施例中计算机设备的内部结构图。
具体实施方式
为了便于理解本申请，为使本申请的上述目的、特征和优点能够更加明显易懂，下面结合附图对本申请的具体实施方式做详细的说明。在下面的描述中阐述了很多具体细节以便于充分理解本申请，附图中给出了本申请的较佳实施方式。但是，本申请可以以许多不同的形式来实现，并不限于本文所描述的实施方式。相反地，提供这些实施方式的目的是使对本申请的公开内容理解的更加透彻全面。本申请能够以很多不同于在此描述的其它方式来实施，本领域技术人员可以在不违背本申请内涵的情况下做类似改进，因此本申请不受下面公开的具体实施例的限制。
此外，术语“第一”、“第二”仅用于描述目的，而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。在本申请的描述中，“多个”的含义是至少两个，例如两个，三个等，除非另有明确具体的限定。在本申请的描述中，“若干”的含义是至少一个，例如一个，两个等，除非另有明确具体的限定。
除非另有定义，本文所使用的所有的技术和科学术语与属于本申请的技术领域的技术人员通常理解的含义相同。本文中所使用的术语只是为了描述具体的实施方式的目的，不是旨在于限制本申请。本文所使用的术语“及/或”包括一个或多个相关的所列项目的任意的和所有的组合。
正电子发射型计算机断层显像(Positron Emission Computed Tomography，PET)，是核医学领域比较先进的临床检查影像技术。是将某种物质，一般是生物生命代谢中必须的物质，如：葡萄糖、蛋白质、核酸、脂肪酸，标记上短寿命的放射性核素(如18F，11C等)，注入人体后，放射性核素在衰变过程中释放出正电子，一个正电子在行进十分之几毫米到几毫米后遇到一个电子后发生湮灭，从而产生方向相反的一对能量为511KeV的光子。这对光子，通过高度灵敏的照相机捕捉，并经计算机进行散射和随机信息的校正。经过对不同的正电子进行相同的分析处理，我们可以得到在生物体内聚集情况的三维图像，从而达到诊断的目的。
在一个实施例中，如图1所示，提供了一种获取晶体效率的方法，包括以下步骤：
步骤S102，获取探测器扫描得到的衰减图像、符合事件数据以及每个探测器晶体的测量单事件数据。
具体地，在进行PET扫描时放射性核素在人体内衰变，并释放出正电子，正电子在运动过程中与电子碰撞发生湮灭，湮灭后产生两个运动方向相反的光子被探测器接收，探测器接收到一对光子的事件称为符合事件，探测器接收到一对光子的数据称为符合事件数据。探测器由许多的晶体构成，每个晶体接收到一个光子的事件称为单事件，每个晶体接收到一个光子的数据称为单事件数据。衰减图像为表示身体各部分组织衰减效率的图像，其中，衰减图可以由电子计算机断层扫描设备(CT)对待检测物体进行扫描得到；衰减图还可以由磁共振成像设备(MRI)扫描待检测物体，再对生成的图像进行图像分割后得到；衰减图还可以由PET系统对待扫描物体进行扫描，通过MLAA算法得到。在PET的临床应用中，通过PET对待扫描人体进行扫描，得到符合事件数据以及每个探测器晶体的测量单事件数据。为了获取更为准确的单事件数据，优选的，对注射氟代脱氧葡萄糖的头部进行扫描，获取探测器扫描得到的衰减图像、符合事件数据以及每个探测器晶体的测量单事件数据。由于对其他部位进行扫描时，在系统的视场范围(FOV)外可能存在较低的活度，并且有较低的概率出现衰减-活度不匹配的问题。因此，最优选的是对注射氟代脱氧葡萄糖的头部进行扫描，获取相应数据，得到的单事件数据更加的精确。
步骤S104，根据衰减图像以及符合事件数据，得到每个探测器晶体的估计单事件数据；
具体地，估计单事件数据可以有两种获取方式。其中一种得到估计单事件数据的方式为：根据所述衰减图像以及符合事件数据进行图像重建，获取估计的活度图像。根据估计的活度图像和衰减图像进行蒙特卡罗模拟，得到每个探测器晶体的估计单事件数据。更具体地，根据所述衰减图像以及符合事件数据进行图像重建，得到估计的活度图像，再根据估计的活度图像以及衰减图像确定蒙特卡罗模拟边界以及初始的光子信息。其中，初始光子信息包括光子的初始位置、光子的初始方向以及光子的初始能量。根据所述蒙特卡罗模拟边界以及初始的光子信息进行光子运动模拟，从而得到光子到达探测器的位置、光子到达探测器的方向、光子到达探测器的能量，也就是得到每个探测器晶体的估计单事件数据。另一种得到估计单事件数据的方式为：根据所述衰减图像以及符合事件数据进行图像重建，获取估计的活度图像。根据活度图像以及衰减图像进行蒙特卡罗模拟得到估计散射单事件数据；或者根据活度图像以及衰减图像进行单散射模拟得到估计散射单事件数据。根据活度图像以及投影矩阵得到估计真实单事件数据。将估计散射单事件数据和估计真实单事件数据求和，得到估计单事件数据。
步骤S106，根据所述测量单事件数据和估计单事件数据，计算得到相应探测器晶体的晶体效率。
具体地，构建所述测量单事件数据、所述估计单事件数据与所述晶体效率的第一预设关系；其中，第一预设关系基于最大似然法建立。最大似然法，也叫最大似然估计，是一种重要而普遍的求估计量的方法。最大似然法明确地使用概率模型，其目标是寻找能够以较高概率产生观察数据的系统发生树。根据所述第一预设关系、测量单事件数据以及估计单事件数据，计算得到相应探测器晶体的晶体效率。也就是将测量单事件数据以及估计单事件数据代入第一预设关系，求解得到相应探测器晶体的晶体效率。
上述获取晶体效率的方法、装置、计算机设备和存储介质，通过临床数据的探测器扫描得到的衰减图像、符合事件数据以及每个探测器晶体的测量单事件数据。在根据衰减图像以及符合事件数据，得到每个探测器晶体的估计单事件数据。最后根据测量单事件数据以及估计单事件数据，计算得到相应探测器晶体的晶体效率。通过临床数据能够精确的计算探测器晶体的晶体效率，并且能够减少放射源的使用次数，减少成本。
在一个实施例中，如图2所示，提供了一种根据单事件数据计算晶体效率的方法，包括以下步骤：
步骤S202，构建所述测量单事件数据、所述估计单事件数据与所述晶体效率的第一预设关系。
具体地，测量单事件数据以及估计单事件数据满足如下方程，也就是第一预设关系：
s(i)＝E(i)(∑j∑kHi，j，kAi，j，kx(k)+R(i))
其中，s(i)为每个探测器晶体的测量单事件数据，Ai，j，k为衰减图像从k点到探测器的投影，Hi，j，k为投影矩阵，E(i)为第i个探测器晶体的晶体效率，x(k)为是未知的像素点k的活度图像，R(i)为每个探测器晶体的估计单事件数据。
根据第一预设关系，计算每个探测器对应单事件数据的似然值，具体计算公式如下：
Figure BDA0002211466530000061
其中，Likelihood为单事件数据的似然值，s(i)为每个探测器晶体的测量单事件数据，Ai，j，k为衰减图像从k点到探测器的投影，Hi，j，k为投影矩阵，E(i)为第i个探测器晶体的晶体效率，x(k)为是未知的像素点k的活度图像，fi(μ，x)为每个探测器晶体的估计单事件数据。
将最大的似然值对应的探测器晶体的晶体效率作为第一参数。
再计算最大的似然值对应的似然值对数，具体计算公式如下：
Figure BDA0002211466530000071
其中，L为似然值对数，s(i)为每个探测器晶体的测量单事件数据，Ai，j，k为衰减图像从k点到探测器的投影，Hi，j，k为投影矩阵，E(i)为第i个探测器晶体的晶体效率，x(k)为是未知的像素点k的活度图像，fi(μ，x)为每个探测器晶体的估计单事件数据。
将最大的似然值的对数作为第二参数。若第一参数与第二参数的偏导等于零，则根据第一预设关系计算相应探测器晶体的晶体效率。也即当似然值最大的晶体效率满足
Figure BDA0002211466530000072
则：
Figure BDA0002211466530000073
由此可以得到探测器晶体计算公式：
Figure BDA0002211466530000074
其中，s(i)为每个探测器晶体的测量单事件数据，Ai，j，k为衰减图像从k点到探测器的投影，Hi，j，k为投影矩阵，E(i)为第i个探测器晶体的晶体效率，x(k)为是未知的像素点k的活度图像，fi(μ，x)为每个探测器晶体的估计单事件数据。
步骤S204，根据所述第一预设关系、测量单事件数据以及估计单事件数据，计算得到相应探测器晶体的晶体效率。
具体地，获取探测器扫描得到的符合事件数据，构建所述符合事件数据与活度图像的第二预设关系。
获取探测器扫描得到的符合事件数据满足一下方程：
y(i，j)＝E(i)*E(j)(Ai，j∑kHi，j，kx(k)+R(i，j))
其中，y(i，j)为探测器扫描得到的符合事件数据，E(i)为第i个探测器晶体的晶体效率，E(j)为第i个探测器晶体的晶体效率，x(k)为是未知的像素点k的活度图像，Hi，j，k为投影矩阵，Ai，j是衰减图的投影，R(i，j)是由于散射和随机事件带来的额外影响。R(i，j)可以通过延迟时间估计随机效应得到，也可以通过蒙特卡罗估计散射效应得到。
在假定一直探测器效率的前提下，使用最大似然法可得到活度图像的迭代方程，也即第二预设关系：
Figure BDA0002211466530000081
其中，y(i，j)为探测器扫描得到的符合事件数据，E(i)为第i个探测器晶体的晶体效率，E(j)为第j个探测器晶体的晶体效率，x(k)为是未知的像素点k的活度图像，Hi，j，k为投影矩阵，Ai，j是衰减图的投影，R(i，j)是由于散射和随机事件带来的额外影响。
所述第二预设关系与所述第一预设关系交替运算，得到相应探测器晶体的晶体效率。其中，第二预设关系中E(j)是未知项，在第一预设关系中x(k)是未知项，为了能同事获取正确的E(j)以及x(k)，第二预设关系与第一预设关系需要交替运算。
更具体地，探测器晶体计算公式中的分母可以通过多种方法计算得到。第一种，通过分析法得到，使用投影矩阵计算∑kHi，j，kAi，j，kx(k)部分，使用单散射模型计算fi(μ，x)部分。第二种，使用蒙特卡罗模拟的方法，直接计算∑j∑kHi，j，kAi，j，kx(k)+fi(μ，x)。
在其中一个实施例中，当探测器晶体的晶体效率接近1的前提下，为获取准确的晶体效率，需要进行确定绝对值的步骤，具体公式如下：
Figure BDA0002211466530000082
其中，E(i)为第i个探测器晶体的晶体效率，N为探测器晶体的数量。
在其中一个实施例中，对上述晶体效率进行正确性理论验证。在忽略散射效应、随机效应、衰减效应以及距离效应的时候，获取的符合事件可以通过第一公式简单表示为：
y(i，j)＝E(i)*E(j)∑kHi，j，kx(k)
其中，y(i，j)为探测器扫描得到的符合事件数据，E(i)为第i个探测器晶体的晶体效率，E(j)为第j个探测器晶体的晶体效率，x(k)为是未知的像素点k的活度图像，Hi，j，k为投影矩阵。
获取的单事件数据可以通过第二公式表示为：
s(i)＝E(i)∑j∑kHi，j，kx(k)
其中，s(i)为每个探测器晶体的测量单事件数据，Hi，j，k为投影矩阵，E(i)为第i个探测器晶体的晶体效率，x(k)为是未知的像素点k的活度图像。
将第一公式代入第二公式可以得到：
Figure BDA0002211466530000091
由上述算式可知，只要矩阵y(i，j)可逆，就可以通过单事件数据和符合事件数据获取晶体效率，虽然不能保证所有扫描获取的y(i，j)可逆，但由于在点源扫描时为y(i，j)对角阵，因此可以认为在扫描物体尺度有限的时候(例如头部扫描)，矩阵y(i，j)可逆。在实际情况下，散射，衰减，距离等效应并不能忽略，故以上推导只能近似的证明此问题可解，不能用于实际计算。
上述根据单事件数据计算晶体效率的方法，能够通过测量单事件数据以及估计单事件数据精确的计算得到晶体探测器效率。并且计算的到的晶体探测器效率能够用于PET图像的重建。
应该理解的是，虽然图1-2的流程图中的各个步骤按照箭头的指示依次显示，但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明，这些步骤的执行并没有严格的顺序限制，这些步骤可以以其它的顺序执行。而且，图1-2中的至少一部分步骤可以包括多个子步骤或者多个阶段，这些子步骤或者阶段并不必然是在同一时刻执行完成，而是可以在不同的时刻执行，这些子步骤或者阶段的执行顺序也不必然是依次进行，而是可以与其它步骤或者其它步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。
在一个实施例中，如图3所示，提供了一种获取晶体效率的装置，包括：获取模块100、单事件数据计算模块200和晶体效率计算模块300，其中：
获取模块100，用于获取探测器扫描得到的衰减图像、符合事件数据以及每个探测器晶体的测量单事件数据。
单事件数据计算模块200，用于根据衰减图像以及符合事件数据，得到每个探测器晶体的估计单事件数据。
晶体效率计算模块300，用于根据所述测量单事件数据和估计单事件数据，计算得到相应探测器晶体的晶体效率。
单事件数据计算模块200，还用于根据所述衰减图像以及符合事件数据进行蒙特卡罗模拟，得到每个探测器晶体的估计单事件数据。
单事件数据计算模块200，还用于根据所述衰减图像以及符合事件数据确定蒙特卡罗模拟边界以及初始的光子信息；根据所述蒙特卡罗模拟边界以及初始的光子信息进行光子运动模拟，得到每个探测器晶体的估计单事件数据。
晶体效率计算模块300包括：第一构件单元以及晶体效率计算单元。
关系构件单元，用于构建所述测量单事件数据、所述估计单事件数据与所述晶体效率的第一预设关系。
晶体效率计算单元，用于根据所述第一预设关系、测量单事件数据以及估计单事件数据，计算得到相应探测器晶体的晶体效率。
晶体效率计算单元包括：关系构件子单元以及晶体效率计算子单元。
关系构件子单元，用于获取探测器扫描得到的符合事件数据，构建所述符合事件数据与活度图像的第二预设关系。
晶体效率计算子单元，用于所述第二预设关系与所述第一预设关系交替运算，得到相应探测器晶体的晶体效率。
获取模块100，还用于对注射氟代脱氧葡萄糖的头部进行扫描，获取探测器扫描得到的衰减图像、符合事件数据以及每个探测器晶体的测量单事件数据。
关于获取晶体效率的装置的具体限定可以参见上文中对于获取晶体效率的方法的限定，在此不再赘述。上述获取晶体效率的装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中，也可以以软件形式存储于计算机设备中的存储器中，以便于处理器调用执行以上各个模块对应的操作。
在一个实施例中，提供了一种计算机设备，该计算机设备可以是终端，其内部结构图可以如图4所示。该计算机设备包括通过系统总线连接的处理器、存储器、网络接口、显示屏和输入装置。其中，该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作系统和计算机程序。该内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该计算机设备的网络接口用于与外部的终端通过网络连接通信。该计算机程序被处理器执行时以实现一种获取晶体效率的方法。该计算机设备的显示屏可以是液晶显示屏或者电子墨水显示屏，该计算机设备的输入装置可以是显示屏上覆盖的触摸层，也可以是计算机设备外壳上设置的按键、轨迹球或触控板，还可以是外接的键盘、触控板或鼠标等。
本领域技术人员可以理解，图4中示出的结构，仅仅是与本申请方案相关的部分结构的框图，并不构成对本申请方案所应用于其上的计算机设备的限定，具体的计算机设备可以包括比图中所示更多或更少的部件，或者组合某些部件，或者具有不同的部件布置。
在一个实施例中，提供了一种计算机设备，包括存储器和处理器，存储器中存储有计算机程序，该处理器执行计算机程序时实现以下步骤：
获取探测器扫描得到的衰减图像、符合事件数据以及每个探测器晶体的测量单事件数据；根据衰减图像以及符合事件数据，得到每个探测器晶体的估计单事件数据；根据所述测量单事件数据和估计单事件数据，计算得到相应探测器晶体的晶体效率。
在一个实施例中，处理器执行计算机程序时还实现以下步骤：
根据所述衰减图像以及符合事件数据进行蒙特卡罗模拟，得到每个探测器晶体的估计单事件数据。
在一个实施例中，处理器执行计算机程序时还实现以下步骤：
根据所述衰减图像以及符合事件数据确定蒙特卡罗模拟边界以及初始的光子信息；根据所述蒙特卡罗模拟边界以及初始的光子信息进行光子运动模拟，得到每个探测器晶体的估计单事件数据。
在一个实施例中，处理器执行计算机程序时还实现以下步骤：
构建所述测量单事件数据、所述估计单事件数据与所述晶体效率的第一预设关系；根据所述第一预设关系、测量单事件数据以及估计单事件数据，计算得到相应探测器晶体的晶体效率。
在一个实施例中，处理器执行计算机程序时还实现以下步骤：
获取探测器扫描得到的符合事件数据，构建所述符合事件数据与活度图像的第二预设关系；所述第二预设关系与所述第一预设关系交替运算，得到相应探测器晶体的晶体效率。
在一个实施例中，处理器执行计算机程序时还实现以下步骤：
对注射氟代脱氧葡萄糖的头部进行扫描，获取探测器扫描得到的衰减图像、符合事件数据以及每个探测器晶体的测量单事件数据。
在一个实施例中，提供了一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现以下步骤：
获取探测器扫描得到的衰减图像、符合事件数据以及每个探测器晶体的测量单事件数据；根据衰减图像以及符合事件数据，得到每个探测器晶体的估计单事件数据；根据所述测量单事件数据和估计单事件数据，计算得到相应探测器晶体的晶体效率。
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：
根据所述衰减图像以及符合事件数据进行蒙特卡罗模拟，得到每个探测器晶体的估计单事件数据。
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：
根据所述衰减图像以及符合事件数据确定蒙特卡罗模拟边界以及初始的光子信息；根据所述蒙特卡罗模拟边界以及初始的光子信息进行光子运动模拟，得到每个探测器晶体的估计单事件数据。
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：
构建所述测量单事件数据、所述估计单事件数据与所述晶体效率的第一预设关系；根据所述第一预设关系、测量单事件数据以及估计单事件数据，计算得到相应探测器晶体的晶体效率。
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：
获取探测器扫描得到的符合事件数据，构建所述符合事件数据与活度图像的第二预设关系；所述第二预设关系与所述第一预设关系交替运算，得到相应探测器晶体的晶体效率。
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：
对注射氟代脱氧葡萄糖的头部进行扫描，获取探测器扫描得到的衰减图像、符合事件数据以及每个探测器晶体的测量单事件数据。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一非易失性计算机可读取存储介质中，该计算机程序在执行时，可包括如上述各方法的实施例的流程。其中，本申请所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用，均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器(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 (10)

1.一种获取晶体效率的方法，其特征在于，所述方法包括：
获取探测器扫描得到的衰减图像、符合事件数据以及每个探测器晶体的测量单事件数据；
根据衰减图像以及符合事件数据，得到每个探测器晶体的估计单事件数据；
根据所述测量单事件数据和估计单事件数据，计算得到相应探测器晶体的晶体效率。
2.根据权利要求1所述的方法，其特征在于，根据衰减图像以及符合事件数据，得到每个探测器晶体的估计单事件数据包括：
根据所述衰减图像以及符合事件数据进行蒙特卡罗模拟，得到每个探测器晶体的估计单事件数据。
3.根据权利要求2所述的方法，其特征在于，所述根据所述衰减图像以及符合事件数据进行蒙特卡罗模拟，得到每个探测器晶体的估计单事件数据包括：
根据所述衰减图像以及符合事件数据确定蒙特卡罗模拟边界以及初始的光子信息；
根据所述蒙特卡罗模拟边界以及初始的光子信息进行光子运动模拟，得到每个探测器晶体的估计单事件数据。
4.根据权利要求1所述的方法，其特征在于，根据所述测量单事件数据和估计单事件数据，计算得到相应探测器晶体的晶体效率包括：
构建所述测量单事件数据、所述估计单事件数据与所述晶体效率的第一预设关系；
根据所述第一预设关系、测量单事件数据以及估计单事件数据，计算得到相应探测器晶体的晶体效率。
5.根据权利要求4所述的方法，其特征在于，所述第一预设关系基于最大似然法建立。
6.根据权利要求4所述的方法，其特征在于，根据所述第一预设关系、测量单事件数据以及估计单事件数据，计算得到相应探测器晶体的晶体效率包括：
获取探测器扫描得到的符合事件数据，构建所述符合事件数据与活度图像的第二预设关系；
所述第二预设关系与所述第一预设关系交替运算，得到相应探测器晶体的晶体效率。
7.根据权利要求1所述的方法，其特征在于，所述获取探测器扫描得到的衰减图像、符合事件数据以及每个探测器晶体的测量单事件数据包括：
对注射氟代脱氧葡萄糖的头部进行扫描，获取探测器扫描得到的衰减图像、符合事件数据以及每个探测器晶体的测量单事件数据。
8.一种获取晶体效率的装置，其特征在于，所述装置包括：
获取模块，用于获取探测器扫描得到的衰减图像、符合事件数据以及每个探测器晶体的测量单事件数据；
单事件数据计算模块，用于根据衰减图像以及符合事件数据，得到每个探测器晶体的估计单事件数据；
晶体效率计算模块，用于根据所述测量单事件数据和估计单事件数据，计算得到相应探测器晶体的晶体效率。
9.一种计算机设备，包括存储器和处理器，所述存储器存储有计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至7中任一项所述方法的步骤。
10.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现权利要求1至7中任一项所述的方法的步骤。
CN201910899826.9A 2019-09-23 2019-09-23 获取晶体效率的方法、装置、计算机设备和存储介质 Active CN110687585B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910899826.9A CN110687585B (zh) 2019-09-23 2019-09-23 获取晶体效率的方法、装置、计算机设备和存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910899826.9A CN110687585B (zh) 2019-09-23 2019-09-23 获取晶体效率的方法、装置、计算机设备和存储介质

Publications (2)

Publication Number Publication Date
CN110687585A CN110687585A (zh) 2020-01-14
CN110687585B true CN110687585B (zh) 2021-05-11

Family

ID=69109909

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910899826.9A Active CN110687585B (zh) 2019-09-23 2019-09-23 获取晶体效率的方法、装置、计算机设备和存储介质

Country Status (1)

Country Link
CN (1) CN110687585B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112269205B (zh) * 2020-10-13 2023-05-23 中国原子能科学研究院 用于确定放射性探测器参数的方法
CN117203555A (zh) * 2021-06-08 2023-12-08 中加健康工程研究院(合肥)有限公司 用于展示pet设备的晶体效率的方法、装置、电子设备及介质

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103315763A (zh) * 2013-07-04 2013-09-25 沈阳东软医疗系统有限公司 对成像设备中的扫描数据进行正规化校正的方法和装置
CN106725573A (zh) * 2017-01-16 2017-05-31 沈阳东软医疗系统有限公司 一种连续进床pet系统的计数丢失校正方法和装置
CN108520542A (zh) * 2018-03-29 2018-09-11 明峰医疗系统股份有限公司 一种pet/ct数据时相匹配的重建方法
CN108831546A (zh) * 2018-06-22 2018-11-16 上海联影医疗科技有限公司 一种数据处理方法、装置和非暂态计算机可读存储介质
WO2019057572A1 (en) * 2017-09-22 2019-03-28 Koninklijke Philips N.V. MANAGING DETECTOR PIXEL PERFORMANCE VARIATION IN DIGITAL POSITRON EMISSION TOMOGRAPHY
CN110009707A (zh) * 2019-04-09 2019-07-12 上海联影医疗科技有限公司 图像散射校正方法、装置、计算机设备和存储介质
CN110025329A (zh) * 2019-04-24 2019-07-19 上海联影医疗科技有限公司 符合计数弦图生成方法、装置、计算机设备和存储介质
CN110215223A (zh) * 2019-05-16 2019-09-10 上海联影医疗科技有限公司 散射校正方法、系统、可读存储介质和设备

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7890282B2 (en) * 2007-09-17 2011-02-15 Siemens Medical Solutions Usa, Inc. Estimation of crystal efficiency with axially compressed sinogram
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
CN106251313B (zh) * 2016-08-15 2020-06-26 上海联影医疗科技有限公司 医学成像方法及系统
CN106388841B (zh) * 2016-09-23 2019-05-21 东软医疗系统股份有限公司 一种图像重建方法、装置及设备
CN107411768B (zh) * 2017-07-31 2020-11-10 东软医疗系统股份有限公司 一种设备校准方法和装置

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103315763A (zh) * 2013-07-04 2013-09-25 沈阳东软医疗系统有限公司 对成像设备中的扫描数据进行正规化校正的方法和装置
CN106725573A (zh) * 2017-01-16 2017-05-31 沈阳东软医疗系统有限公司 一种连续进床pet系统的计数丢失校正方法和装置
WO2019057572A1 (en) * 2017-09-22 2019-03-28 Koninklijke Philips N.V. MANAGING DETECTOR PIXEL PERFORMANCE VARIATION IN DIGITAL POSITRON EMISSION TOMOGRAPHY
CN108520542A (zh) * 2018-03-29 2018-09-11 明峰医疗系统股份有限公司 一种pet/ct数据时相匹配的重建方法
CN108831546A (zh) * 2018-06-22 2018-11-16 上海联影医疗科技有限公司 一种数据处理方法、装置和非暂态计算机可读存储介质
CN110009707A (zh) * 2019-04-09 2019-07-12 上海联影医疗科技有限公司 图像散射校正方法、装置、计算机设备和存储介质
CN110025329A (zh) * 2019-04-24 2019-07-19 上海联影医疗科技有限公司 符合计数弦图生成方法、装置、计算机设备和存储介质
CN110215223A (zh) * 2019-05-16 2019-09-10 上海联影医疗科技有限公司 散射校正方法、系统、可读存储介质和设备

Also Published As

Publication number Publication date
CN110687585A (zh) 2020-01-14

Similar Documents

Publication Publication Date Title
CN106491151B (zh) Pet图像获取方法及系统
US9155514B2 (en) Reconstruction with partially known attenuation information in time of flight positron emission tomography
CN110584698B (zh) 探测器质量控制效验方法、装置、计算机设备和存储介质
CN108932741B (zh) 动态pet参数成像方法、装置、系统和计算机可读存储介质
CN108615250B (zh) 图像重建方法、装置、系统和计算机可读存储介质
US10492737B2 (en) Cross-calibration for quantitative functional imaging
US20090088620A1 (en) Quantitative clinical and pre-clinical imaging
CN111080584A (zh) 医学图像的质控方法、计算机设备和可读存储介质
WO2020146032A1 (en) Determination of metabolic rate from static pet scan
CN107392876B (zh) 获取归一化校正因子的方法、装置及医学成像方法
US10126439B2 (en) Reconstruction with multiple photopeaks in quantitative single photon emission computed tomography
CN110749916A (zh) 获取pet探测器晶体时间延迟量的方法、装置和计算机设备
CN111588399B (zh) 医学成像设备状态监控的方法、设备和计算机设备
CN110687585B (zh) 获取晶体效率的方法、装置、计算机设备和存储介质
CN110047116B (zh) Pet图像校正方法、装置、计算机设备和存储介质
WO2021004157A1 (zh) 医学扫描成像方法、装置、存储介质及计算机设备
CN110136076B (zh) 医学扫描成像方法、装置、存储介质及计算机设备
CN110327065B (zh) 单事件计数率校正方法、装置、计算机设备和存储介质
CN110223247B (zh) 图像衰减校正方法、装置、计算机设备和存储介质
CN112365479B (zh) Pet参数图像处理方法、装置、计算机设备及存储介质
CN108873047B (zh) 检测放射源活度的方法、系统、计算机设备和存储介质
CN110400361B (zh) 子集划分及图像重建的方法、装置和计算机设备
US11354831B2 (en) Image reconstruction using tracer and various parameters
CN113712583A (zh) 死时间校正方法、装置、计算机设备和存储介质
CN108703769B (zh) Tof数据的校正方法、装置、系统和计算机可读存储介质

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