CN109998582A - 符合判选方法、装置、设备和介质 - Google Patents
符合判选方法、装置、设备和介质 Download PDFInfo
- Publication number
- CN109998582A CN109998582A CN201910299086.5A CN201910299086A CN109998582A CN 109998582 A CN109998582 A CN 109998582A CN 201910299086 A CN201910299086 A CN 201910299086A CN 109998582 A CN109998582 A CN 109998582A
- Authority
- CN
- China
- Prior art keywords
- event
- tdc
- time
- value
- period
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 39
- 230000000737 periodic effect Effects 0.000 claims description 40
- 238000012163 sequencing technique Methods 0.000 claims description 14
- 238000012216 screening Methods 0.000 claims description 10
- 230000001360 synchronised effect Effects 0.000 claims description 9
- 238000004590 computer program Methods 0.000 claims description 2
- 230000015654 memory Effects 0.000 description 15
- 230000004044 response Effects 0.000 description 10
- 238000010586 diagram Methods 0.000 description 9
- 230000006870 function Effects 0.000 description 5
- 239000000523 sample Substances 0.000 description 5
- 238000001514 detection method Methods 0.000 description 2
- 238000002600 positron emission tomography Methods 0.000 description 2
- 230000000630 rising effect Effects 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000003851 biochemical process Effects 0.000 description 1
- 239000011248 coating agent Substances 0.000 description 1
- 238000000576 coating method Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000010006 flight Effects 0.000 description 1
- 238000000338 in vitro Methods 0.000 description 1
- 239000003068 molecular probe Substances 0.000 description 1
- 230000008450 motivation Effects 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 230000035479 physiological effects, processes and functions Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/037—Emission tomography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/42—Arrangements for detecting radiation specially adapted for radiation diagnosis
- A61B6/4266—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a plurality of detector units
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5205—Devices using data or image processing specially adapted for radiation diagnosis involving processing of raw data to produce diagnostic data
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Radiology & Medical Imaging (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Physics & Mathematics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Nuclear Medicine (AREA)
- Measurement Of Radiation (AREA)
Abstract
本发明实施例公开了一种符合判选方法、装置、设备和介质，涉及医学图像技术领域。该方法包括：依据探测器接收光子的事件发生周期，确定第一事件发生时刻所在周期和第二事件发生时刻所在周期的周期差；依据所述周期差、第一事件的TDC值和第二事件的TDC值，确定第一事件和第二事件发生时间差；若所述时间差小于设定时间阈值，则确定第一事件和第二事件具有符合关系。本发明实施例提供的一种符合判选方法、装置、设备和介质，实现了对不同探测器单元间或相同探测器单元内的任意两个探测器的符合判断。
Description
技术领域
本发明实施例涉及医学图像技术领域，尤其涉及一种符合判选方法、装置、设备和介质。
背景技术
正电子发射断层成像(Positron Emission Tomography，以下均简称PET)是一种非侵入式的造影方法。它成像的基本原理是将正电子放射性核素标记于分子探针上，当放射性核素衰变产生的正电子与生物体内的负电子碰撞湮灭后，发出一对能量为511KeV，运动方向近似相反的γ光子。PET采用环绕生物体的探测器环将入射的γ光子转换为电信号，从而获得其能量、位置和时间信息。通过湮灭符合技术，得到湮灭事件所在响应线的位置，并通过二维或三维断层重建算法获得正电子核素在生物体中的分布，从而在体外观测生物体内的生理和生化过程。如图1所示，一个正电子湮灭后产生互为180度的两个伽马光子，被探测器对M0和M10接收到，M0和M10间这一条光子飞行的直线叫做响应线，光子到达时间差就是两个光子飞行的路程差。
传统的PET设备通常只有一个PET探测器单元，该单元内部由二十个以上的探测器组成探测器环。因此，传统的符合装置只能处理一个单元内部各探测器数据间的小窗宽符合判选。
然而，传统只有一个PET探测器单元的PET设备的几何立体角覆盖比例较小，有大量γ光子从环状探测器两端空隙逃逸，没有被探测器收集到，从而符合事件的探测效率很低。为了提高检测灵敏度，在环状探测器两端对称添加多个PET探测器单元。
针对由多个PET探测器单元构成的探测器设备，需要符合处理装置能完成不同探测器单元间的任意两个探测器的符合处理。由于两个单元上的探测器组成的响应线较同一单元内部两个探测器间响应线更长，所以传统的小窗宽符合判选装置已经不能实现对不同探测器单元间的任意两个探测器的符合处理。
发明内容
本发明实施例提供一种符合判选方法、装置、设备和介质，以实现对不同探测器单元间或相同探测器单元内的任意两个探测器的符合判断。
第一方面，本发明实施例提供了一种符合判选方法，该方法包括：
依据探测器接收光子的事件发生周期，确定第一事件发生时刻所在周期和第二事件发生时刻所在周期的周期差；
依据所述周期差、第一事件的TDC(time-to-digital converter,时间数字转换)值和第二事件的TDC值，确定第一事件和第二事件发生时间差；
若所述时间差小于设定时间阈值，则确定第一事件和第二事件具有符合关系。
第二方面，本发明实施例还提供了一种符合判选装置，该装置包括：
周期差确定模块，用于依据探测器接收光子的事件发生周期，确定第一事件发生时刻所在周期和第二事件发生时刻所在周期的周期差；
时间差确定模块，用于依据所述周期差、第一事件的TDC时间数字转换值和第二事件的TDC值，确定第一事件和第二事件发生时间差；
符合确定模块，用于若所述时间差小于设定时间阈值，则确定第一事件和第二事件具有符合关系。
第三方面，本发明实施例还提供了一种设备，所述设备包括：
一个或多个处理器；
存储装置，用于存储一个或多个程序，
当所述一个或多个程序被所述一个或多个处理器执行，使得所述一个或多个处理器实现如本发明实施例中任一所述的一种符合判选方法。
第四方面，本发明实施例还提供了一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现如本发明实施例中任一所述的一种符合判选方法。
本发明实施例通过依据第一事件发生时刻所在周期和第二事件发生时刻所在周期的周期差、第一事件的TDC值和第二事件的TDC值，确定第一事件和第二事件发生时间差；依据发生时间差确定两事件的是否具有符合关系。从而不仅可以实现对相同探测器单元内的任意两个探测器的符合判断，还可以实现对不同探测器单元间的任意两个探测器的符合判断。
附图说明
图1为PET设备中单个探测器单元内部响应线示意图；
图2是PET设备中两个探测器单元间跨单元响应线示意图；
图3是符合事件出现在同一个时钟周期且A事件先发生的示意图；
图4是符合事件出现在同一个时钟周期且B事件先发生的示意图；
图5是符合事件出现在不同时钟周期且B事件先A事件发生一个时钟周期的示意图；
图6是符合事件出现在不同时钟周期且A事件先B事件发生一个时钟周期的示意图；
图7是符合事件出现在不同时钟周期且B事件先A事件发生两个时钟周期的示意图；
图8是符合事件出现在不同时钟周期且A事件先B事件发生两个时钟周期的示意图；
图9是本发明实施例一提供的一种符合判选方法的流程图；
图10是本发明实施例二提供的一种符合判选方法的流程图；
图11是本发明实施例三提供的一种符合判选方法的流程图；
图12是本发明实施例三提供的一种符合判选装置的结构示意图；
图13是本发明实施例四提供的一种符合判选装置的结构示意图；
图14为本发明实施例五提供的一种设备的结构示意图。
具体实施方式
下面结合附图和实施例对本发明作进一步的详细说明。可以理解的是，此处所描述的具体实施例仅仅用于解释本发明，而非对本发明的限定。另外还需要说明的是，为了便于描述，附图中仅示出了与本发明相关的部分而非全部结构。
参见图2，发明人在实现本发明的过程中发现，当PET设备存在多个探测器单元时，存在跨探测器单元的两个探测器间的响应线，即探测器单元U0中探测器M0与探测器单元U3中探测器M10之间的响应线。由于该响应线相比同一探测器单元内的两个探测器间的响应线会更长，所以两个光子可能出现的飞行时间差会更大。为了把这样的符合事件甄别出来，需要符合处理装置具备处理大窗宽的符合判断的能力。
为此发明人遍历所有可能出现的情况，针对每种情况分析得出TDC_A、TDC_B和Twindow间的关系如下。其中TDC_A是探测器A事件(以下简称A事件)发生时距离下一个时钟周期上升沿的时间的TDC量化值，也即A事件的TDC值。TDC_B是探测器B事件(以下简称B事件)发生时距离下一个时钟周期上升沿的时间的TDC量化值，也即B事件的TDC值。Twindow是符合窗宽的大小，也即设定时间阈值，此处以超过一个时钟周期为例。
图3示出了符合事件出现在同一个时钟周期，且A事件先B事件发生，两事件的发生时间和TDC值。从图3可以看出若A事件和B事件为符合事件，则需要满足Tab＝TDC_A-TDC_B<Twindow。其中Tab是探测器A和探测器B接收到的两个事件的飞行时间差。
图4示出了符合事件出现在同一个时钟周期，且B事件先A事件发生，两事件的发生时间和TDC值。从图4可以看出若A事件和B事件为符合事件，则需满足Tab＝TDC_B-TDC_A<Twindow。
综合图3和图4，若两个符合事件出现在同一个时钟周期内，则需满足|TDC_A-TDC_B|<Twindow，该公式也变形为：-Twindow<TDC_A–TDC_B<Twindow。
图5示出了符合事件出现在不同时钟周期，且B事件先A事件发生一个时钟周期，两事件的发生时间和TDC值。若两个事件为符合事件，则需满足Tab＝TDC_B+T-TDC_A<Twindow，公式变形得TDC_A-TDC_B>-Twindow+T。其中T为***同步时钟周期，以10ns为例。
图6示出了符合事件出现在不同时钟周期，且A事件先B事件发生一个时钟周期，两事件的发生时间和TDC值。若两个事件为符合事件，则需满足Tab＝TDC_A+T-TDC_B<Twindow，公式变形为TDC_A-TDC_B<-T+Twindow。
图7示出了符合事件出现在不同时钟周期，且B事件先A事件发生两个时钟周期，两事件的发生时间和TDC值。若两个事件为符合事件，则需满足Tab＝TDC_B+2T-TDC_A<Twindow，公式变形得TDC_A-TDC_B>2T-Twindow。
图8示出了符合事件出现在不同时钟周期，且A事件先B事件发生两个时钟周期，两事件的发生时间和TDC值。若两个事件为符合事件，则需满足Tab＝TDC_A+2T-TDC_B<Twindow，公式变形为TDC_A-TDC_B<Twindow-2T。
综上所述：
1.如果相较于B事件，A事件发生在下N个时钟周期，N为整数，那么若两个事件为符合事件，则需满足Tab＝TDC_B+N×T-TDC_A<Twindow即TDC_A-TDC_B>N×T-Twindow；
2.如果相较于A事件，B事件发生在下N个时钟周期，N为整数，那么若两个事件为符合事件，则需满足Tab＝TDC_A+N×T-TDC_B<Twindow即TDC_A-TDC_B<-N×T+Twindow；
3.N＝0时候，表明A事件、B事件发生在同一个时钟周期内。
对于Twindow的取值没有限制，可以根据需要取任意大。
实施例一
图9是本发明实施例一提供的一种符合判选方法的流程图。本实施例可适用于对不同探测器单元间或相同探测器单元内的任意两个探测器的符合判断的情况。该方法可以由一种符合判选装置来执行，该装置可以由软件和/或硬件的方式实现。参见图9，本实施例提供的符合判选方法包括：
S110、依据探测器接收光子的事件发生周期，确定第一事件发生时刻所在周期和第二事件发生时刻所在周期的周期差。
其中，第一事件和第二事件是待符合判断的两个探测器接收光子的事件。且第一事件的所属探测器与第二事件的所属探测器可以位于同一探测器单元，也可以位于不同探测器单元。
具体地，探测器接收光子的事件发生周期可以通过事件发生时刻***的时钟脉冲所属周期确定。
S120、依据所述周期差、第一事件的TDC值和第二事件的TDC值，确定第一事件和第二事件发生时间差。
具体地，若所述周期差为0，说明第一事件和第二事件位于同一周期，则第一事件和第二事件发生时间差等于第一事件的TDC值和第二事件的TDC值的绝对值。
若所述周期差为非零，说明第一事件和第二事件是跨周期事件，且第一事件先于第二事件发生，则第一事件和第二事件发生时间差等于第一事件的TDC值+所述周期差×时钟周期-第二事件的TDC值。
S130、若所述时间差小于设定时间阈值，则确定第一事件和第二事件具有符合关系。
其中，设定时间阈值是指两个具有符合关系的事件的发生时间差的最大值。若将设定时间阈值看作是时间窗口，则本实施例可适应任意大窗宽的符合判选，具体可以根据实际需要设定。
本发明实施例的技术方案，通过依据第一事件发生时刻所在周期和第二事件发生时刻所在周期的周期差、第一事件的TDC值和第二事件的TDC值，确定第一事件和第二事件发生时间差；依据发生时间差确定两事件的是否具有符合关系。从而不仅可以实现对相同探测器单元内的任意两个探测器的符合判断，还可以实现对不同探测器单元间的任意两个探测器的符合判断。对于未来的PET设备，无论探测器单元数有多少，任意两个单元间的数据符合判断，均适用。
为实现对后续符合事件的处理，所述若所述时间差小于设定时间阈值，则确定第一事件和第二事件具有符合关系之后，还包括：
将第一事件和第二事件各自的时间、能量和位置信息合并成一个符合数据并输出。
需要强调的是，本实施例可适用于后端处理，与探测器连接，由硬件电路实现，电路印刷在符合电路板上。相比将第一事件和第二事件输出至前端利用滑窗进行符合判断，本实施例节省了第一事件和第二事件的输出，从而提高了符合事件判断的实时性和处理速度。
实施例二
图10是本发明实施例二提供的一种符合判选方法的流程图。本实施例是在上述实施例的基础上，提出的一种可选方案。参见图10，本实施例提供的符合判选方法包括：
S210、依据探测器接收光子的事件发生周期，确定第一事件发生时刻所在周期和第二事件发生时刻所在周期的周期差。
S220、依据所述周期差或探测器接收光子的事件发生周期，确定第一事件和第二事件发生时间的先后顺序。
具体地，第一事件和第二事件发生时间的先后顺序可以依据所述周期差的正负确定，也可以根据第一事件和第二事件的发生周期确定。
例如，第一事件发生在第一时钟周期，第二事件发生在第二时钟周期，那么就可以确定第一事件先于第二事件发生。
S230、依据确定的先后顺序确定目标时间差公式。
具体地，所述依据确定的先后顺序确定目标时间差公式，包括：
若第一事件先于第二事件发生，则确定目标时间周期差公式为：第一事件的TDC值+周期差的绝对值×***同步时钟周期-第二事件的TDC值；
若第二事件先于第一事件发生，则确定目标时间周期差公式为：第二事件的TDC值+周期差的绝对值×***同步时钟周期-第一事件的TDC值。
S240、将所述周期差、第一事件的TDC值和第二事件的TDC值带入所述目标时间差公式，得出第一事件和第二事件发生时间差。
S250、若所述时间差小于设定时间阈值，则确定第一事件和第二事件具有符合关系。
本发明实施例的技术方案，通过根据第一事件和第二事件发生时间的先后顺序确定目标时间差公式，将所述周期差、第一事件的TDC值和第二事件的TDC值带入所述目标时间差公式，得出第一事件和第二事件发生时间差。若所述时间差小于设定时间阈值，则确定第一事件和第二事件具有符合关系。从而实现对两事件的符合判断。
实施例三
图11是本发明实施例三提供的一种符合判选方法的流程图。本实施例是在上述实施例的基础上，以待符合判断的两事件分别为探测器A和探测器B采集的A事件与B事件为例，提出的一钟可选方案。参见图11，本实施例提供的符合判选方法包括：
S310、确定A事件与B事件发生时刻的周期差。
S320、确定A事件与B事件发生的先后顺序。
S330、计算A事件的TDC值与B事件的TDC值之间的差值，并赋值给一个有符号寄存器。
S340、若A事件先于B事件发生，选用不等式TDC_A-TDC_B<-N×T+Twindow，若B事件先于A事件发生，选用不等式TDC_A-TDC_B>N×T-Twindow，进行符合判断，其中TDC_A是A事件的TDC值，TDC_B是B事件的TDC值，N是A事件与B事件发生时刻的周期差，T是***时钟周期，Twindow是设定时间阈值。
S350、从有符号寄存器获取A事件的TDC值与B事件的TDC值之间的差值，根据A事件的TDC值与B事件的TDC值之间的差值以及确定的不等式，确定A事件与B事件是否具有符合关系；把两个具有符合关系的单事件合并成一个数据，输出具有符合关系的两个事件各自的时间、能量、位置信息。
本发明实施例对各步骤的执行顺序不做限定。可选地，S340可以先于S330执行。
参见图12，基于本实施例描述的符合判选方法，符合装置的可以包括顺序确定模块、作差模块、符合不等式确定模块和符合判选及符合输出模块。
其中，顺序模块用于确定A事件与B事件发生的先后顺序；作差模块用于计算A事件的TDC值与B事件的TDC值之间的差值，并赋值给一个有符号寄存器；符合不等式确定模块用于若A事件先于B事件发生，选用不等式TDC_A-TDC_B<-N×T+Twindow，若B事件先于A事件发生，选用不等式TDC_A-TDC_B>N×T-Twindow，进行符合判断；符合判选及符合输出模块用于从有符号寄存器获取A事件的TDC值与B事件的TDC值之间的差值，根据A事件的TDC值与B事件的TDC值之间的差值以及确定的不等式，确定A事件与B事件是否具有符合关系；把两个具有符合关系的单事件合并成一个数据，输出具有符合关系的两个事件各自的时间、能量、位置信息。
本实施例的技术方案可以实现如下效果，本实施例能适应任意大窗宽的符合判选，覆盖了所有的窗宽情形，窗宽取值的不同就是一个变化例。
另外，本实施例能够处理跨任意单元数的符合判选，对于未来的PET设备，无论探测器单元数有多少，任意两个单元间的数据符合判断，均适用。
需要说明的是，基于上述实施例的技术教导，本领域技术人员有动机将上述实施方式进行组合，以实现对不同探测器单元间或相同探测器单元内的任意两个探测器的符合判断。
实施例四
图13是本发明实施例四提供的一种符合判选装置的结构示意图。参见图13，本实施例提供的符合判选装置包括：周期差确定模块10、时间差确定模块20、符合确定模块30。
其中，周期差确定模块10，用于依据探测器接收光子的事件发生周期，确定第一事件发生时刻所在周期和第二事件发生时刻所在周期的周期差；
时间差确定模块20，用于依据所述周期差、第一事件的TDC时间数字转换值和第二事件的TDC值，确定第一事件和第二事件发生时间差；
符合确定模块30，用于若所述时间差小于设定时间阈值，则确定第一事件和第二事件具有符合关系。
本发明实施例通过依据第一事件发生时刻所在周期和第二事件发生时刻所在周期的周期差、第一事件的TDC值和第二事件的TDC值，确定第一事件和第二事件发生时间差；依据发生时间差确定两事件的是否具有符合关系。从而不仅可以实现对相同探测器单元内的任意两个探测器的符合判断，还可以实现对不同探测器单元间的任意两个探测器的符合判断。
进一步地，所述时间差确定模块，包括：发生顺序确定单元、目标公式确定单元和时间差确定单元。
其中，发生顺序确定单元，用于依据所述周期差或探测器接收光子的事件发生周期，确定第一事件和第二事件发生时间的先后顺序；
目标公式确定单元，用于依据确定的先后顺序确定目标时间差公式；
时间差确定单元，用于将所述周期差、第一事件的TDC值和第二事件的TDC值带入所述目标时间差公式，得出第一事件和第二事件发生时间差。
进一步地，所述目标公式确定单元具体用于：
若第一事件先于第二事件发生，则确定目标时间周期差公式为：第一事件的TDC值+周期差的绝对值×***同步时钟周期-第二事件的TDC值；
若第二事件先于第一事件发生，则确定目标时间周期差公式为：第二事件的TDC值+周期差的绝对值×***同步时钟周期-第一事件的TDC值。
进一步地，所述的装置，还包括：数据合并模块和符合数据输出模块。
数据合并模块，用于所述若所述时间差小于设定时间阈值，则确定第一事件和第二事件具有符合关系之后，将第一事件和第二事件各自的时间、能量和位置信息合并成一个符合数据；
符合数据输出模块，用于将符合数据输出。
本发明实施例所提供的符合判选装置可执行本发明任意实施例所提供的符合判选方法，具备执行方法相应的功能模块和有益效果。
实施例五
图14为本发明实施例五提供的一种设备的结构示意图。如图14所示，该设备包括处理器70、存储器71、输入装置72和输出装置73；设备中处理器70的数量可以是一个或多个，图14中以一个处理器70为例；设备中的处理器70、存储器71、输入装置72和输出装置73可以通过总线或其他方式连接，图14中以通过总线连接为例。
存储器71作为一种计算机可读存储介质，可用于存储软件程序、计算机可执行程序以及模块，如本发明实施例中的符合判选方法对应的程序指令/模块(例如，符合判选装置中的周期差确定模块10、时间差确定模块20、符合确定模块30)。处理器70通过运行存储在存储器71中的软件程序、指令以及模块，从而执行设备的各种功能应用以及数据处理，即实现上述的符合判选方法。
存储器71可主要包括存储程序区和存储数据区，其中，存储程序区可存储操作***、至少一个功能所需的应用程序；存储数据区可存储根据终端的使用所创建的数据等。此外，存储器71可以包括高速随机存取存储器，还可以包括非易失性存储器，例如至少一个磁盘存储器件、闪存器件、或其他非易失性固态存储器件。在一些实例中，存储器71可进一步包括相对于处理器70远程设置的存储器，这些远程存储器可以通过网络连接至设备。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。
输入装置72可用于接收输入的数字或字符信息，以及产生与设备的用户设置以及功能控制有关的键信号输入。输出装置73可包括显示屏等显示设备。
实施例六
本发明实施例六还提供一种包含计算机可执行指令的存储介质，所述计算机可执行指令在由计算机处理器执行时用于执行一种符合判选方法，该方法包括：
依据探测器接收光子的事件发生周期，确定第一事件发生时刻所在周期和第二事件发生时刻所在周期的周期差；
依据所述周期差、第一事件的TDC时间数字转换值和第二事件的TDC值，确定第一事件和第二事件发生时间差；
若所述时间差小于设定时间阈值，则确定第一事件和第二事件具有符合关系
当然,本发明实施例所提供的一种包含计算机可执行指令的存储介质,其计算机可执行指令不限于如上所述的方法操作,还可以执行本发明任意实施例所提供的符合判选方法中的相关操作.
通过以上关于实施方式的描述，所属领域的技术人员可以清楚地了解到，本发明可借助软件及必需的通用硬件来实现，当然也可以通过硬件实现，但很多情况下前者是更佳的实施方式。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在计算机可读存储介质中，如计算机的软盘、只读存储器(Read-Only Memory,ROM)、随机存取存储器(RandomAccess Memory,RAM)、闪存(FLASH)、硬盘或光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述的方法。
值得注意的是，上述符合判选装置的实施例中，所包括的各个单元和模块只是按照功能逻辑进行划分的，但并不局限于上述的划分，只要能够实现相应的功能即可；另外，各功能单元的具体名称也只是为了便于相互区分，并不用于限制本发明的保护范围。
注意，上述仅为本发明的较佳实施例及所运用技术原理。本领域技术人员会理解，本发明不限于这里所述的特定实施例，对本领域技术人员来说能够进行各种明显的变化、重新调整和替代而不会脱离本发明的保护范围。因此，虽然通过以上实施例对本发明进行了较为详细的说明，但是本发明不仅仅限于以上实施例，在不脱离本发明构思的情况下，还可以包括更多其他等效实施例，而本发明的范围由所附的权利要求范围决定。
Claims (10)
1.一种符合判选方法，其特征在于，包括：
依据探测器接收光子的事件发生周期，确定第一事件发生时刻所在周期和第二事件发生时刻所在周期的周期差；
依据所述周期差、第一事件的TDC时间数字转换值和第二事件的TDC值，确定第一事件和第二事件发生时间差；
若所述时间差小于设定时间阈值，则确定第一事件和第二事件具有符合关系。
2.根据权利要求1所述方法，其特征在于，所述依据所述周期差、第一事件的TDC值和第二事件的TDC值，确定第一事件和第二事件发生时间差，包括：
依据所述周期差或探测器接收光子的事件发生周期，确定第一事件和第二事件发生时间的先后顺序；
依据确定的先后顺序确定目标时间差公式；
将所述周期差、第一事件的TDC值和第二事件的TDC值带入所述目标时间差公式，得出第一事件和第二事件发生时间差。
3.根据权利要求2所述的方法，其特征在于，所述依据确定的先后顺序确定目标时间差公式，包括：
若第一事件先于第二事件发生，则确定目标时间周期差公式为：第一事件的TDC值+周期差的绝对值×***同步时钟周期-第二事件的TDC值；
若第二事件先于第一事件发生，则确定目标时间周期差公式为：第二事件的TDC值+周期差的绝对值×***同步时钟周期-第一事件的TDC值。
4.根据权利要求1-3中任一所述的方法，其特征在于，所述若所述时间差小于设定时间阈值，则确定第一事件和第二事件具有符合关系之后，还包括：
将第一事件和第二事件各自的时间、能量和位置信息合并成一个符合数据并输出。
5.一种符合判选装置，其特征在于，包括：
周期差确定模块，用于依据探测器接收光子的事件发生周期，确定第一事件发生时刻所在周期和第二事件发生时刻所在周期的周期差；
时间差确定模块，用于依据所述周期差、第一事件的TDC时间数字转换值和第二事件的TDC值，确定第一事件和第二事件发生时间差；
符合确定模块，用于若所述时间差小于设定时间阈值，则确定第一事件和第二事件具有符合关系。
6.根据权利要求5所述装置，其特征在于，所述时间差确定模块，包括：
发生顺序确定单元，用于依据所述周期差或探测器接收光子的事件发生周期，确定第一事件和第二事件发生时间的先后顺序；
目标公式确定单元，用于依据确定的先后顺序确定目标时间差公式；
时间差确定单元，用于将所述周期差、第一事件的TDC值和第二事件的TDC值带入所述目标时间差公式，得出第一事件和第二事件发生时间差。
7.根据权利要求6所述的装置，其特征在于，所述目标公式确定单元具体用于：
若第一事件先于第二事件发生，则确定目标时间周期差公式为：第一事件的TDC值+周期差的绝对值×***同步时钟周期-第二事件的TDC值；
若第二事件先于第一事件发生，则确定目标时间周期差公式为：第二事件的TDC值+周期差的绝对值×***同步时钟周期-第一事件的TDC值。
8.根据权利要求5-7中任一所述的装置，其特征在于，还包括：
数据合并模块，用于所述若所述时间差小于设定时间阈值，则确定第一事件和第二事件具有符合关系之后，将第一事件和第二事件各自的时间、能量和位置信息合并成一个符合数据；
符合数据输出模块，用于将符合数据输出。
9.一种设备，其特征在于，所述设备包括：
一个或多个处理器；
存储装置，用于存储一个或多个程序，
当所述一个或多个程序被所述一个或多个处理器执行，使得所述一个或多个处理器实现如权利要求1-4中任一所述的一种符合判选方法。
10.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该程序被处理器执行时实现如权利要求1-4中任一所述的一种符合判选方法。
Priority Applications (4)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910299086.5A CN109998582B (zh) | 2019-04-15 | 2019-04-15 | 符合判选方法、装置、设备和介质 |
| EP20791156.1A EP3937785A4 (en) | 2019-04-15 | 2020-04-15 | SYSTEMS AND METHODS FOR DETERMINING TRUE COINCIDENTAL EVENTS |
| PCT/CN2020/084861 WO2020211764A1 (en) | 2019-04-15 | 2020-04-15 | Systems and methods for determining true coincidence events |
| US17/450,837 US20220031272A1 (en) | 2019-04-15 | 2021-10-14 | Systems and methods for determining true coincidence events |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910299086.5A CN109998582B (zh) | 2019-04-15 | 2019-04-15 | 符合判选方法、装置、设备和介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109998582A true CN109998582A (zh) | 2019-07-12 |
| CN109998582B CN109998582B (zh) | 2022-05-13 |
Family
ID=67171764
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910299086.5A Active CN109998582B (zh) | 2019-04-15 | 2019-04-15 | 符合判选方法、装置、设备和介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109998582B (zh) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110432923A (zh) * | 2019-08-14 | 2019-11-12 | 上海联影医疗科技有限公司 | 图像重建方法、装置、计算机设备和存储介质 |
| CN110507344A (zh) * | 2019-09-04 | 2019-11-29 | 上海联影医疗科技有限公司 | 一种符合判选方法、装置、设备及存储介质 |
| WO2020211764A1 (en) * | 2019-04-15 | 2020-10-22 | Shanghai United Imaging Healthcare Co., Ltd. | Systems and methods for determining true coincidence events |
| CN112102426A (zh) * | 2020-08-28 | 2020-12-18 | 上海联影医疗科技股份有限公司 | 本底符合事件判选方法、装置、设备及可读存储介质 |
| CN114780168A (zh) * | 2022-03-30 | 2022-07-22 | 全球能源互联网研究院有限公司南京分公司 | 智能终端容器安全策略动态变更的方法、装置及电子设备 |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070206721A1 (en) * | 2004-09-13 | 2007-09-06 | Tkaczyk John E | Photon counting x-ray detector with overrange logic control |
| CN101750621A (zh) * | 2008-12-08 | 2010-06-23 | 皇家飞利浦电子股份有限公司 | 对来自像素化探测器的测量的评估 |
| CN102631212A (zh) * | 2012-04-28 | 2012-08-15 | 中国科学院高能物理研究所 | 正电子发射断层扫描仪及其中的符合判选方法 |
| CN102933150A (zh) * | 2011-06-08 | 2013-02-13 | 株式会社东芝 | 正电子放射断层摄影方法、正电子放射断层摄影装置及非暂时性计算机可读存储介质 |
| US20140008542A1 (en) * | 2011-03-30 | 2014-01-09 | The Board Of Trustees Of The Stanford Junior University | Methods and Systems for increasing the sensitivity of simultaneous multi-isotope positron emission tomography |
| CN105450215A (zh) * | 2015-11-09 | 2016-03-30 | 科大国盾量子技术股份有限公司 | 一种符合测量***及方法 |
| CN106108934A (zh) * | 2016-08-31 | 2016-11-16 | 清华大学 | 多伽马光子同时发射药物时间符合核医学成像***及方法 |
| CN106974671A (zh) * | 2016-06-12 | 2017-07-25 | 上海联影医疗科技有限公司 | 正电子发射断层成像***及其图像重建方法 |
| CN107464270A (zh) * | 2017-07-17 | 2017-12-12 | 沈阳东软医疗***有限公司 | 一种图像重建方法和装置 |
| US20180123773A1 (en) * | 2016-11-01 | 2018-05-03 | Shenyang Neusoft Medical Systems Co., Ltd. | Clock synchronizing |
| CN108013895A (zh) * | 2017-11-28 | 2018-05-11 | 上海联影医疗科技有限公司 | 一种符合判断方法、装置、设备和介质 |
| CN109350099A (zh) * | 2018-09-13 | 2019-02-19 | 中山市明峰医疗器械有限公司 | 一种应用于临床pet***的随机事件去除处理方法 |
-
2019
- 2019-04-15 CN CN201910299086.5A patent/CN109998582B/zh active Active
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070206721A1 (en) * | 2004-09-13 | 2007-09-06 | Tkaczyk John E | Photon counting x-ray detector with overrange logic control |
| CN101750621A (zh) * | 2008-12-08 | 2010-06-23 | 皇家飞利浦电子股份有限公司 | 对来自像素化探测器的测量的评估 |
| US20140008542A1 (en) * | 2011-03-30 | 2014-01-09 | The Board Of Trustees Of The Stanford Junior University | Methods and Systems for increasing the sensitivity of simultaneous multi-isotope positron emission tomography |
| CN102933150A (zh) * | 2011-06-08 | 2013-02-13 | 株式会社东芝 | 正电子放射断层摄影方法、正电子放射断层摄影装置及非暂时性计算机可读存储介质 |
| CN102631212A (zh) * | 2012-04-28 | 2012-08-15 | 中国科学院高能物理研究所 | 正电子发射断层扫描仪及其中的符合判选方法 |
| CN105450215A (zh) * | 2015-11-09 | 2016-03-30 | 科大国盾量子技术股份有限公司 | 一种符合测量***及方法 |
| CN106974671A (zh) * | 2016-06-12 | 2017-07-25 | 上海联影医疗科技有限公司 | 正电子发射断层成像***及其图像重建方法 |
| CN106108934A (zh) * | 2016-08-31 | 2016-11-16 | 清华大学 | 多伽马光子同时发射药物时间符合核医学成像***及方法 |
| US20180123773A1 (en) * | 2016-11-01 | 2018-05-03 | Shenyang Neusoft Medical Systems Co., Ltd. | Clock synchronizing |
| CN107464270A (zh) * | 2017-07-17 | 2017-12-12 | 沈阳东软医疗***有限公司 | 一种图像重建方法和装置 |
| CN108013895A (zh) * | 2017-11-28 | 2018-05-11 | 上海联影医疗科技有限公司 | 一种符合判断方法、装置、设备和介质 |
| CN109350099A (zh) * | 2018-09-13 | 2019-02-19 | 中山市明峰医疗器械有限公司 | 一种应用于临床pet***的随机事件去除处理方法 |
Non-Patent Citations (2)
| Title |
|---|
| CHAO WANG等: "A Real Time Coincidence System for High Count-Rate TOF or Non-TOF PET Cameras Using Hybrid Method Combining ANDLogic and Time-Mark Technology", 《IEEE TRANS NUCL SCI》 * |
| 边兆英: "基于区域时空先验的动态PET重建及PET图像恢复算法研究", 《中国优秀硕士学位论文全文数据库信息科技辑》 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2020211764A1 (en) * | 2019-04-15 | 2020-10-22 | Shanghai United Imaging Healthcare Co., Ltd. | Systems and methods for determining true coincidence events |
| CN110432923A (zh) * | 2019-08-14 | 2019-11-12 | 上海联影医疗科技有限公司 | 图像重建方法、装置、计算机设备和存储介质 |
| CN110507344A (zh) * | 2019-09-04 | 2019-11-29 | 上海联影医疗科技有限公司 | 一种符合判选方法、装置、设备及存储介质 |
| CN110507344B (zh) * | 2019-09-04 | 2022-05-13 | 上海联影医疗科技股份有限公司 | 一种符合判选方法、装置、设备及存储介质 |
| CN112102426A (zh) * | 2020-08-28 | 2020-12-18 | 上海联影医疗科技股份有限公司 | 本底符合事件判选方法、装置、设备及可读存储介质 |
| CN112102426B (zh) * | 2020-08-28 | 2024-03-26 | 上海联影医疗科技股份有限公司 | 本底符合事件判选方法、装置、设备及可读存储介质 |
| CN114780168A (zh) * | 2022-03-30 | 2022-07-22 | 全球能源互联网研究院有限公司南京分公司 | 智能终端容器安全策略动态变更的方法、装置及电子设备 |
| CN114780168B (zh) * | 2022-03-30 | 2023-04-28 | 全球能源互联网研究院有限公司南京分公司 | 智能终端容器安全策略动态变更的方法、装置及电子设备 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109998582B (zh) | 2022-05-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109998582A (zh) | 符合判选方法、装置、设备和介质 | |
| CN102631212B (zh) | 正电子发射断层扫描仪及其中的符合判选方法 | |
| Sander et al. | Towards increased trustworthiness of deep learning segmentation methods on cardiac MRI | |
| CN102648856B (zh) | 正电子发射断层扫描仪及其中的符合判选方法 | |
| CN102783964B (zh) | Pet成像中单事件列表式数据的同步方法及*** | |
| Markiewicz et al. | Rapid processing of PET list-mode data for efficient uncertainty estimation and data analysis | |
| Yuan et al. | A two-stage hog feature extraction processor embedded with SVM for pedestrian detection | |
| CN110327067A (zh) | 图像重建方法、装置、终端设备及pet*** | |
| Fysikopoulos et al. | A Spartan 6 FPGA-based data acquisition system for dedicated imagers in nuclear medicine | |
| Luciw et al. | Automated generation of cerebral blood flow and arterial transit time maps from multiple delay arterial spin‐labeled MRI | |
| Kirkland et al. | Imaging from temporal data via spiking convolutional neural networks | |
| US8117142B2 (en) | Method of real-time crystal peak tracking for positron emission tomography (PET) avalanche-photodiodes (APD) detector | |
| CN110507344A (zh) | 一种符合判选方法、装置、设备及存储介质 | |
| CN108013895B (zh) | 一种符合判断方法、装置、设备和介质 | |
| US10367653B2 (en) | Transmitting data in PET system | |
| CN104146726B (zh) | Pet***符合探测响应的生成方法 | |
| US8704190B2 (en) | Radiation detection signal processing method and system | |
| Zhao et al. | An fpga-based hardware accelerator of ransac algorithm for matching of images feature points | |
| Esteve Bosch et al. | The Event Detection System in the NEXT-White Detector | |
| Zhang et al. | A hardware-oriented histogram of oriented gradients algorithm and its VLSI implementation | |
| Gu et al. | HDNet: a hierarchically decoupled network for crowd counting | |
| RU146670U1 (ru) | Система мониторинга данных изображений объектов динамической компьютерной томографии | |
| Raath et al. | A permutation entropy analysis of voyager interplanetary magnetic field observations | |
| Chen et al. | PET timing calibration using low rank constraint | |
| Chen et al. | Dynamic supervisor for cross-dataset object detection |
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 | ||
| 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. |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |