CN105212954A - 一种脉冲堆积事件实时处理方法与系统 - Google Patents
一种脉冲堆积事件实时处理方法与系统 Download PDFInfo
- Publication number
- CN105212954A CN105212954A CN201510751817.7A CN201510751817A CN105212954A CN 105212954 A CN105212954 A CN 105212954A CN 201510751817 A CN201510751817 A CN 201510751817A CN 105212954 A CN105212954 A CN 105212954A
- Authority
- CN
- China
- Prior art keywords
- pulse
- pile
- mvt
- value
- sampling
- 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
- 238000003672 processing method Methods 0.000 title claims abstract description 28
- 238000000034 method Methods 0.000 claims abstract description 101
- 238000005070 sampling Methods 0.000 claims abstract description 64
- 238000012545 processing Methods 0.000 claims abstract description 19
- 238000004519 manufacturing process Methods 0.000 claims abstract description 5
- 230000010354 integration Effects 0.000 claims description 58
- 230000008569 process Effects 0.000 claims description 25
- 230000000630 rising effect Effects 0.000 claims description 13
- 238000009825 accumulation Methods 0.000 claims description 12
- 230000008859 change Effects 0.000 claims description 8
- 238000012937 correction Methods 0.000 claims description 8
- 230000002085 persistent effect Effects 0.000 claims description 5
- 230000000694 effects Effects 0.000 abstract description 7
- 238000011084 recovery Methods 0.000 abstract description 6
- 238000010586 diagram Methods 0.000 description 7
- 239000013078 crystal Substances 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 238000013016 damping Methods 0.000 description 3
- 238000001514 detection method Methods 0.000 description 3
- 230000006866 deterioration Effects 0.000 description 3
- 239000002245 particle Substances 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 238000009206 nuclear medicine Methods 0.000 description 2
- 238000002600 positron emission tomography Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000005316 response function Methods 0.000 description 2
- 230000011664 signaling Effects 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- 229910052688 Gadolinium Inorganic materials 0.000 description 1
- 206010019133 Hangover Diseases 0.000 description 1
- BPQQTUXANYXVAA-UHFFFAOYSA-N Orthosilicate Chemical compound [O-][Si]([O-])([O-])[O-] BPQQTUXANYXVAA-UHFFFAOYSA-N 0.000 description 1
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical compound [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- ORCSMBGZHYTXOV-UHFFFAOYSA-N bismuth;germanium;dodecahydrate Chemical compound O.O.O.O.O.O.O.O.O.O.O.O.[Ge].[Ge].[Ge].[Bi].[Bi].[Bi].[Bi] ORCSMBGZHYTXOV-UHFFFAOYSA-N 0.000 description 1
- 235000000332 black box Nutrition 0.000 description 1
- JJWKPURADFRFRB-UHFFFAOYSA-N carbonyl sulfide Chemical compound O=C=S JJWKPURADFRFRB-UHFFFAOYSA-N 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000005713 exacerbation Effects 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- UIWYJDYFSGRHKR-UHFFFAOYSA-N gadolinium atom Chemical compound [Gd] UIWYJDYFSGRHKR-UHFFFAOYSA-N 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000010365 information processing Effects 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 238000012856 packing Methods 0.000 description 1
- 230000002285 radioactive effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 229910052710 silicon Inorganic materials 0.000 description 1
- 239000010703 silicon Substances 0.000 description 1
- 230000001960 triggered effect Effects 0.000 description 1
- 229910052727 yttrium Inorganic materials 0.000 description 1
- VWQVUPCCIRVNHF-UHFFFAOYSA-N yttrium atom Chemical compound [Y] VWQVUPCCIRVNHF-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- 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/16—Measuring radiation intensity
- G01T1/17—Circuit arrangements not adapted to a particular type of detector
- G01T1/171—Compensation of dead-time counting losses
-
- 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/16—Measuring radiation intensity
- G01T1/161—Applications in the field of nuclear medicine, e.g. in vivo counting
- G01T1/164—Scintigraphy
- G01T1/1641—Static instruments for imaging the distribution of radioactivity in one or two dimensions using one or several scintillating elements; Radio-isotope cameras
- G01T1/1642—Static instruments for imaging the distribution of radioactivity in one or two dimensions using one or several scintillating elements; Radio-isotope cameras using a scintillation crystal and position sensing photodetector arrays, e.g. ANGER cameras
-
- 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/16—Measuring radiation intensity
- G01T1/161—Applications in the field of nuclear medicine, e.g. in vivo counting
-
- 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/16—Measuring radiation intensity
- G01T1/161—Applications in the field of nuclear medicine, e.g. in vivo counting
- G01T1/164—Scintigraphy
-
- 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/16—Measuring radiation intensity
- G01T1/161—Applications in the field of nuclear medicine, e.g. in vivo counting
- G01T1/164—Scintigraphy
- G01T1/1641—Static instruments for imaging the distribution of radioactivity in one or two dimensions using one or several scintillating elements; Radio-isotope cameras
- G01T1/1648—Ancillary equipment for scintillation cameras, e.g. reference markers, devices for removing motion artifacts, calibration devices
-
- 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/16—Measuring radiation intensity
- G01T1/20—Measuring radiation intensity with scintillation detectors
- G01T1/2006—Measuring radiation intensity with scintillation detectors using a combination of a scintillator and photodetector which measures the means radiation intensity
-
- 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/16—Measuring radiation intensity
- G01T1/24—Measuring radiation intensity with semiconductor detectors
- G01T1/248—Silicon photomultipliers [SiPM], e.g. an avalanche photodiode [APD] array on a common Si substrate
-
- 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
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D10/00—Energy efficient computing, e.g. low power processors, power management or thermal management
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Molecular Biology (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Measurement Of Radiation (AREA)
- Nuclear Medicine (AREA)
Abstract
一种脉冲堆积事件实时处理方法，其包括步骤：生成脉冲下降沿的拟合基线值查询表和拟合能量值查询表，利用多电压阈值采样方法完成堆积脉冲的识别和脉冲信号处理过程的触发；利用脉冲先验信息和已采集脉冲信息以查找表的方式获得因脉冲堆积未正确采样部分的信息，完成了堆积脉冲信息的实时恢复。本发明首先利用多电压阈值采样方法完成了高计数率下堆积脉冲的识别和脉冲信号处理过程的触发，接下来利用脉冲先验信息和已采集脉冲信息以查找表的方式获得因脉冲堆积未正确采样部分的信息，完成了堆积脉冲信息的实时恢复。本发明的方法简洁高效，易于在探测器级别的实时数据采集系统上实现，在低采样率的情况下仍然能够较优的堆积脉冲信息恢复效果。
Description
技术领域
本发明涉及数字信号处理领域，尤其涉及一种基于多电压阈值采样技术的脉冲堆积事件实时处理方法与系统。
背景技术
高能粒子探测领域以及正电子发射断层成像(PositronEmissionTomography，简称PET)等核医学成像领域需要完成伽马光子信息的探测，通常利用闪烁晶体将伽马光子转变成可见光光子，再用光电转换器件转换成闪烁脉冲电信号，最终经数据采集电路提取闪烁脉冲电信号所包含的对应伽马光子的信息。
计数性能是高能粒子探测器的基础性能指标之一，反映了系统在单位时间内能采集到的事件数量极限。在诸如高辐射通量宇宙射线探测以及短半衰期放射性示踪剂、动态扫描等高活度、高灵敏度要求的核医学成像应用中，将对高能粒子探测器的计数性能极限提出甚高要求。
而随着探测器的计数率升高，闪烁脉冲的堆积现象将越发严重，并成为影响探测器性能和计数极限的主要原因。脉冲堆积定义为：同一个探测通道在一个闪烁脉冲持续时间内探测到两个以上的伽马光子事件。脉冲堆积使得多个闪烁脉冲波形重叠，影响数据采集系统对每个闪烁脉冲信息的独立求取，将造成计数丢失、能量信息恶化和位置信息恶化等一系列问题，严重恶化高计数率下探测器的性能表现。
为解决高计数率下脉冲堆积引发的探测器性能恶化问题，研究人员提出了多种方法，如动态积分方法:对两个脉冲触发的时间间隔进行检测，若在一个积分时间内有第二个脉冲到来，则结束对第一个脉冲的积分，并通过脉冲模型对第一个脉冲未完成积分的部分进行估算，在第一个脉冲积分值上加上该估算值，并从第二个脉冲的积分值中减去该估算值，最终得到堆积两个脉冲的实际积分值，参见文献[1]。
文献[2]中提出的数字化单事件信息恢复(DigitalSingleEventReconstruction，简称DiSER)方法可看作对动态积分方法的改进，DiSER方法基于闪烁脉冲的模型进行堆积事件中单个脉冲的信息恢复，该脉冲模型由一个快速上升的上升沿和一个指数衰减的下降沿构成。DiSER方法通过对采样得到的一帧数字化脉冲波形求导，确定闪烁脉冲的个数和发生时刻的间隔，两个闪烁脉冲上升沿之间的采样点，被认为是前一个脉冲的下降沿并被用于重建其全部下降沿，在求取后一个闪烁脉冲的波形时将减去该下降沿，以恢复堆积事件中各个脉冲的真实波形。
另一种思路是将打入探测器的伽马光子看作一个冲激函数，而将探测器视作一个黑盒系统，并根据其构成具有特定的系统响应函数，将探测器输出到电子学的闪烁脉冲电信号视作冲激函数与系统响应函数进行卷积后的结果，考虑到噪声的存在，无法通过解析的方式求解出原脉冲，可引入迭代的方法从闪烁脉冲电信号中恢复出伽马光子的信息。由于将探测器视作一个整体系统考虑，该方法从工作原理上不受限于特定的脉冲形状，且由于恢复后的脉冲持续时间相当短，适用于超高计数率的探测器应用，参见文献[3]。
另一些研究工作针对脉冲堆积引起的基线漂移进行处理，如[4]中对堆积引起的基线漂移问题，采用了统计基线值的分布的方法:由于电路白噪声符合高斯分布而闪烁脉冲的拖尾符合泊松分布，同时考虑到即使在高活度下，无脉冲拖尾影响的点仍然为整个模数转换器(Analog-to-DigitalConverter，简称ADC)采样波形中的绝大部分，即从基线值得统计分布上可以分辨出无脉冲基线构成的高斯峰，该高斯峰的均值即为当前的基线值，因此通过实时更新一个ADC采样波形中基线部分的基线库，并统计其幅值分布，可追踪当前基线值并反馈到脉冲采样端进行基线校正。
动态积分方法和DiSER方法在应用于位置敏感型的探测器时，将面临一个伽马光子事件所产生的多个同源闪烁脉冲信号的处理同步问题，如位置敏感型光电倍增管(PhotomultiplierTube，简称PMT)接收到伽马光子后将产生一组包含位置信息的角脉冲信号，且同组的角脉冲信号波形存在差异，采用上述方法很可能造成只有部分脉冲能够得到触发和后续处理，且由于各角脉冲信号独立触发，造成触发时刻不同步，恶化能量及位置信息的求取精度。此外，以上两种方法需要采集完成一定时长的数据帧后再进行数据回放来确定堆积脉冲的发生时刻，不能在探测器上进行堆积脉冲的实时处理。
利用探测器系统响应函数和迭代算法进行脉冲恢复的方法对数字化闪烁脉冲过程的采样率有较高要求，且迭代算法耗时较多，使得该方法在需要实时处理的探测器及系统上的实现存在困难。
基于脉冲基线值统计的脉冲堆积基线漂移处理方法的基线更新过程需要一定响应时间，只能减轻脉冲拖尾造成的基线漂移影响，无法在发生脉冲堆积时对脉冲进行实时的恢复处理。
因此，针对上述技术问题，有必要提供一种改良结构的脉冲堆积实时处理方法，以克服上述缺陷。
[1]LewellenTK；BiceAN；PollardKR,etal."EvaluationofaClinicalScintillationCamerawithPulseTailExtrapolationElectronics,"TheJournalofNuclearMedicine,vol.30,no.9,pp.1554-1558,1989.
[2]X.Wang；Q.G.Xie；Y.B.Chen；M.Niu；andP.Xiao,"AdvantagesofDigitallySamplingScintillationPulsesinPileupProcessinginPET,"IEEETransactionsonNuclearScience,vol.59,pp.498-506,Jun2012.
[3]Z.Z.Deng；Q.G.Xie；Z.W.DuanandP.Xiao,"Scintillationeventenergymeasurementviaapulsemodelbasediterativedeconvolutionmethod,"Phys.Med.Biol,vol.58,pp.7815,2013.
[4]H.D.Li；C.Wang；BaghaeiH.；Y.X.Zhang；RamirezR.；S.T.Liu；S.H.An；Wai-HoiWong,"ANewStatistics-BasedOnlineBaselineRestorerforaHighCount-RateFullyDigitalSystem,"IEEETransactionsonNuclearScience,vol.57,no.2,pp.550-555,April2010.
发明内容
有鉴于此，本发明的目的在于提供一种基于多电压阈值采样(Multi-voltageThreshold，简称MVT)技术的脉冲堆积事件实时处理方法，该方法简洁高效，易于在探测器级别的实时数据采集系统上实现，在低采样率的情况下仍然能够较优的堆积脉冲信息恢复效果。
为实现上述目的，本发明提供如下技术方案：
一种脉冲堆积事件实时处理方法，其包括步骤：
S1：获取非堆积脉冲数字波形数据库，生成脉冲下降沿的数字波形查询表，并存储在数据采集系统中；
S2：在进行有堆积脉冲的数据采集时，利用多电压阈值采样方法完成堆积脉冲的识别和脉冲信号处理过程的触发；
S3：在脉冲信号处理过程中，利用已采集的未堆积部分脉冲信息查找数字波形表以获得因脉冲堆积未正确采样部分的信息，完成了堆积脉冲的实时处理。
上述的脉冲堆积事件实时处理方法，优选地，所述步骤S1包括步骤：
P1：在无堆积脉冲的情况下，用多电压阈值采样方法作为采集正常脉冲的触发信号；
P2：触发后用固定时间间隔采样(Regular-timesampling，均简称RTS)方法进行脉冲波形的采样，得到闪烁脉冲数字波形；
P3：提取闪烁脉冲数字波形库中各个闪烁脉冲的下降沿的幅值变化情况，并生成对应的拟合基线值查询表和拟合能量值查询表；
P4：将上述表存储在数据采集系统中。
上述的脉冲堆积事件实时处理方法，优选地，所述步骤S2包括步骤：
P1：数据采集系统对高计数下的闪烁脉冲电信号进行采集，由于闪烁脉冲具有一定持续时间，且出现时间不固定，因此相距过近的两个或多个闪烁脉冲将发生脉冲堆积；
P2：闪烁脉冲电信号分别被MVT方法采样，MVT方法利用多个电压阈值与闪烁脉冲波形快速上升沿进行电压比较，无论该脉冲是否为堆积脉冲都能在其上升沿得到电平翻转信号；
P3：将MVT方法得到的电平翻转信号作为脉冲的识别信号，触发RTS方法对脉冲进行采样和数字积分求取脉冲能量信息。
上述的脉冲堆积事件实时处理方法，优选地，所述步骤S3具体包括步骤：
a.在数据采集系统进行高计数下的数据采样过程中，当有堆积事件发生时，前一个MVT触发信号作为前一个脉冲的数字积分开始信号；
b.当第二个MVT触发信号发生时，利用堆积脉冲中前一脉冲未堆积部分脉冲下降沿作为输入参数查找拟合能量值查询表以得到前一脉冲堆积部分能量值，并与已采集的未堆积部分能量值加和作为前一脉冲的能量值；
c.同时利用前一脉冲未堆积部分脉冲下降沿作为输入参数查询拟合基线值查询表以得到前一脉冲在堆积发生后的脉冲幅值，并在每个采样时刻将该脉冲幅值按照脉冲下降沿衰减规律进行更新；
d.在后一脉冲的数字积分过程中将在每个采样时刻减去当前的前一脉冲的幅值更新值，以得到后一脉冲的真实幅值；
e.如果在后一个脉冲的数字积分时间内再次遇到新的MVT触发信号，则将后一个脉冲视作前一个脉冲按照上述处理方法迭代进行。
为实现上述目的，本发明提供如下技术方案：
一种脉冲堆积事件实时处理系统，其包括：
触发模块，用以当MVT触发信号到来时，触发模块产生积分触发信号，并将前一脉冲已采样下降沿波形更新到当前基线值；
数字积分模块，用以当接收到积分触发信号后将RTS采样值减去基线校正值作为真实的脉冲幅值进行积分；
拟合能量值查询模块，用以根据当前基线值和已积分时间来给出脉冲未积分部分能量的拟合值。
拟合基线值查询模块，用以根据当前基线值和已积分时间来获取前一个脉冲因堆积而未采集的波形幅值，并作为后一个脉冲的基线校正值。
上述的脉冲堆积事件实时处理系统，优选地，所述数字积分模块中如果在数字积分周期内未发生第二次MVT触发信号，认为无堆积事件发生，数字积分周期完成时输出的脉冲能量值即为当前MVT触发信号对应的脉冲能量值。
上述的脉冲堆积事件实时处理系统，优选地，所述数字积分模块中如果在数字积分周期内发生第二次MVT触发信号，数字积分模块用拟合信息查询模块得到的未积分部分的拟合能量值与第二次MVT触发信号前的已完成积分部分能量值相加作为当前MVT触发信号对应的脉冲能量值。
从上述技术方案可以看出，本发明基于MVT技术的脉冲堆积事件实时处理方法，首先利用多电压阈值采样方法完成了高计数率下堆积脉冲的识别和脉冲信号处理过程的触发，接下来利用脉冲先验信息和已采集脉冲信息以查找表的方式获得因脉冲堆积未正确采样部分的信息，完成了堆积脉冲信息的实时恢复。本发明的方法简洁高效，易于在探测器级别的实时数据采集系统上实现，在低采样率的情况下仍然能够较优的堆积脉冲信息恢复效果。
附图说明
图1为本发明脉冲堆积事件实时处理方法的步骤示意图；
图2为本发明脉冲堆积事件实时处理方法的流程图，其中，图2(a)为两个伽马光子先后打入探测器单元的示意图，图2(b)为闪烁脉冲电信号分别被MVT方法和RTS方法采样的示意图，图2(c)为利用MVT方法得到的电平翻转信号作为有脉冲出现的触发信号，利用RTS方法得到的电压采样点进行数字积分求取脉冲能量，以得到堆积事件中各个脉冲的信息的示意图；
图3为本发明脉冲堆积事件实时处理方法在探测器中的实现框架的示意图；
图4为本发明脉冲堆积事件实时处理方法中闪烁脉冲下降沿衰减时间分布的示意图。
具体实施方式
本发明公开了一种基于多电压阈值采样(Multi-voltageThreshold，简称MVT)技术的脉冲堆积事件实时处理方法，该方法简洁高效，易于在探测器级别的实时数据采集系统上实现，在低采样率的情况下仍然能够较优的堆积脉冲信息恢复效果。
闪烁脉冲通常被认为是具有一个快速上升沿和持续较久的下降沿构成，其快速上升沿包含了脉冲的时间信息，而整个脉冲波形的幅值变化包含了脉冲的能量信息。由于绝大多数堆积事件中各个脉冲上升沿仍然完好存在，数据获取系统只需完成脉冲上升沿的采样即可判断脉冲是否出现以及相应的时间信息；另一方面闪烁脉冲波形符合一定的变化规律，因此对于特定的探测器和前端电子学而言，其脉冲数学模型相对恒定，可通过已经得到的脉冲采样点信息推算出剩余未采样部分脉冲的信息。因此可利用MVT数字化方法对快速的上升沿进行采样以识别脉冲有无和到来时刻，再对脉冲波形下降沿进行部分采样以恢复因脉冲堆积未正确采样的脉冲波形，从而实现了堆积脉冲的识别和信息恢复。
如图1所示，图1为本发明脉冲堆积事件实时处理方法的步骤示意图，本发明脉冲堆积事件实时处理方法，其包括步骤
S1：获取非堆积脉冲数字波形数据库，生成脉冲下降沿的数字波形查询表，并存储在数据采集系统中；
S2：在进行有堆积脉冲的数据采集时，利用多电压阈值采样方法完成堆积脉冲的识别和脉冲信号处理过程的触发；
S3：在脉冲信号处理过程中，利用已采集的未堆积部分脉冲信息查找数字波形表以获得因脉冲堆积未正确采样部分的信息，完成了堆积脉冲的实时处理。
本发明采取MVT方法进行堆积脉冲的识别。MVT方法低成本、高性能的实现了闪烁脉冲上升沿的采样，可用于脉冲堆积中脉冲的实时识别。此外该方法为同一伽马光子在位置敏感型探测器中产生的多个闪烁脉冲角信号的处理提供了同步触发信号，提高了位置信息的求取精度。
本发明基于查表的堆积脉冲信息实时处理。根据脉冲先验信息预先得到未积分部分的能量值和幅值变化的可能情况并写入查找表，在实时数据采集过程中根据脉冲已采样部分的信息进行查表，得到堆积脉冲中前一脉冲未准确采样部分的能量信息，以及前一脉冲在发生脉冲堆积后对后一脉冲采样波形的能量贡献。
所述步骤S1包括步骤：
P1：在无堆积脉冲的情况下，用多电压阈值采样方法作为采集正常脉冲的触发信号；
P2：触发后用RTS方法进行脉冲波形的采样，得到闪烁脉冲数字波形；
P3：提取闪烁脉冲数字波形库中各个闪烁脉冲的下降沿的幅值变化情况，并生成对应的拟合基线值查询表和拟合能量值查询表；
P4：将上述表存储在数据采集系统中。
所述步骤S2包括步骤：
P1：数据采集系统对高计数下的闪烁脉冲电信号进行采集，由于闪烁脉冲具有一定持续时间，且出现时间不固定，因此相距过近的两个或多个闪烁脉冲将发生脉冲堆积；
P2：闪烁脉冲电信号分别被MVT方法采样，MVT方法利用多个电压阈值与闪烁脉冲波形快速上升沿进行电压比较，无论该脉冲是否为堆积脉冲都能在其上升沿得到电平翻转信号；
P3：将MVT方法得到的电平翻转信号作为脉冲的识别信号，触发RTS方法对脉冲进行采样和数字积分求取脉冲能量信息。
所述步骤S3具体包括步骤：
a.在数据采集系统进行高计数下的数据采样过程中，当有堆积事件发生时，前一个MVT触发信号作为前一个脉冲的数字积分开始信号；
b.当第二个MVT触发信号发生时，利用堆积脉冲中前一脉冲未堆积部分脉冲下降沿作为输入参数查找拟合能量值查询表以得到前一脉冲堆积部分能量值，并与已采集的未堆积部分能量值加和作为前一脉冲的能量值；
c.同时利用前一脉冲未堆积部分脉冲下降沿作为输入参数查询拟合基线值查询表以得到前一脉冲在堆积发生后的脉冲幅值，并在每个采样时刻将该脉冲幅值按照脉冲下降沿衰减规律进行更新；
d.在后一脉冲的数字积分过程中将在每个采样时刻减去当前的前一脉冲的幅值更新值，以得到后一脉冲的真实幅值；
e.如果在后一个脉冲的数字积分时间内再次遇到新的MVT触发信号，则将后一个脉冲视作前一个脉冲按照上述处理方法迭代进行。
本发明具体描述如图2所示，图2为多电压阈值采样技术的脉冲堆积事件实时处理方法流程图。堆积事件产生及处理流程如下所述：
图2(a)：两个伽马光子200以较短的时间间隔先后打入探测器单元100。
图2(b)：探测器将伽马光子转换成闪烁脉冲电信号，由于闪烁脉冲具有一定持续时间，因此相距过近的两个闪烁脉冲将发生脉冲堆积。闪烁脉冲电信号分别被MVT方法和RTS方法采样，其中MVT方法利用电压阈值与脉冲波形快速上升沿进行电压比较，得到电平翻转信号；RTS方法以固定的时间间隔对脉冲进行连续的电压采样。
图2(c)：利用MVT方法得到的电平翻转信号作为有脉冲出现的触发信号，利用RTS方法得到的电压采样点进行数字积分求取脉冲能量，以得到堆积事件中各个脉冲的信息，具体过程如下：
a.数据采集系统获取低计数下的非堆积脉冲数字波形数据库，生成脉冲下降沿的拟合基线值查询表和拟合能量值查询表，并存储在数据采集系统中；
b.在数据采集系统进行高计数下的数据采样过程中，当有堆积事件发生时，前一个MVT触发信号作为前一个脉冲的数字积分开始信号；
c.当第二个MVT触发信号发生时，利用堆积脉冲中前一脉冲未堆积部分脉冲下降沿作为输入参数查找拟合能量值查询表以得到前一脉冲堆积部分能量值，并与已采集的未堆积部分能量值加合作为前一脉冲的能量值；
d.同时利用前一脉冲未堆积部分脉冲下降沿作为输入参数查询拟合基线值查询表以得到前一脉冲在堆积发生后的脉冲幅值，并在每个采样时刻将该脉冲幅值按照脉冲下降沿衰减规律进行更新；
e.在后一脉冲的数字积分过程中将在每个采样时刻减去对应的前一脉冲的幅值，以得到后一脉冲的真实幅值。
f.如果在后一个脉冲的数字积分时间内再次遇到新的MVT触发信号，则将后一个脉冲视作前一个脉冲按照上述处理方法迭代进行。
本发明方法在探测器中的实现范例如图3所示，即本发明脉冲堆积实时处理系统，其包括：
触发模块：当MVT触发信号到来时，触发模块产生积分触发信号，并将前一脉冲已采样信息更新。
数字积分模块：数字积分模块接受到积分触发信号后将RTS采样值减去基线校正值作为真实的脉冲幅值进行积分。如果在数字积分周期内未发生第二次MVT触发信号，认为无堆积事件发生，数字积分周期完成时输出的脉冲能量值即为当前MVT触发信号对应的脉冲能量值。如果在数字积分周期内发生第二次MVT触发信号，数字积分模块用拟合信息查询模块得到的未积分部分的拟合能量值与第二次MVT触发信号前的已完成积分部分能量值相加作为当前MVT触发信号对应的脉冲能量值。
拟合能量值查询模块，根据当前基线值和已积分时间来给出脉冲未积分部分能量的拟合值。
拟合基线值查询模块，根据当前基线值和已积分时间来获取前一个脉冲因堆积而未采集的波形幅值，并作为后一个脉冲的基线校正值。
由于闪烁脉冲波形符合一定的变化规律，因此对于特定的探测器和前端电子学而言，其脉冲数学模型相对恒定。如图4所示，对闪烁脉冲下降沿的数字波形采样点进行指数衰减拟合后，其衰减常数的分布集中在一个特定值附件。因此可预先对闪烁脉冲下降沿所有可能的已采样点数和幅值，根据下降沿指数衰减方程计算出对应的剩余未采样点数的幅值以及能量总和，分别写入拟合基线值查询模块和拟合能量值查询模块中，在实际数据采集过程中将已采样点数及幅值作为地址输入到上述查找表中即可实时获得未采样脉冲的幅值变化和能量信息。
与现有技术相比，本发明的优点包括：
1、MVT数字化方法对探测器输出的闪烁脉冲上升沿进行快速准确的采样，可实时分辨堆积脉冲中各个脉冲的发生时刻。
2、MVT数字化方法应用于位置敏感型探测器的脉冲堆积处理时，可通过对伽马光子产生的总闪烁脉冲信号进行采样，实现由总信号分别产生的多个位置信号的同步采样和处理，便于准确求取各个位置信号之间的比例以获取位置信息。
3、由于根据特定的探测器及前端电子学下的闪烁脉冲模型预先确定了闪烁脉冲未积分部分的能量值和幅值变化的可能情况，并以查找表的形式储存在探测器硬件中，在实际数据采样过程中可根据已采样脉冲信息查表迅速得到对应结果，不需要进行复杂的脉冲拟合方程计算，从而实现了实时、高通量的堆积脉冲信息处理。
4、实现了对两个及以上连续到来的堆积脉冲的迭代处理流程，且由于无论有无脉冲堆积，都会在每个脉冲到来时将基线校正值更新，因此能够对两个脉冲相距时间间隔大于积分周期但前一脉冲的下降沿幅值尚能对后续脉冲信息求取造成影响的情况进行处理。
本发明的方法与系统中，PET探测器所使用的闪烁晶体可以是闪烁晶体硅酸钇镥(简称LYSO)、晶体锗酸铋(简称BGO)、闪烁晶体硅酸钆(简称GSO)等任意可将伽马光子转换成可见光光子的器件。PET探测器所使用的光电探测器件可以是PMT、硅光电倍增器(SiliconPhotomultiplier，简称SiPM)等任意可将可见光信号转换成闪烁脉冲电信号的器件。脉冲能量信息的获取方法可以是ADC等RTS类型的采样方法，也可以是MVT，超越阈值时间采样方法(TimeOverThreshold，简称TOT)等电压域采样方法。根据已采样部分脉冲信息得到未采样部分的能量值和幅值变化的方法可采用脉冲下降沿指数拟合、脉冲分段积分值比例等方法。
Claims (7)
1.一种脉冲堆积事件实时处理方法，其特征在于：包括步骤：
S1：获取非堆积脉冲数字波形数据库，生成脉冲下降沿的数字波形查询表，并存储在数据采集系统中；
S2：在进行有堆积脉冲的数据采集时，利用多电压阈值采样方法完成堆积脉冲的识别和脉冲信号处理过程的触发；
S3：在脉冲信号处理过程中，利用已采集的未堆积部分脉冲信息查找数字波形表以获得因脉冲堆积未正确采样部分的信息，完成了堆积脉冲的实时处理。
2.根据权利要求1所述的脉冲堆积事件实时处理方法，其特征在于：所述步骤S1包括步骤：
P1：在无堆积脉冲的情况下，用多电压阈值采样方法作为采集正常脉冲的触发信号；
P2：触发后用RTS方法进行脉冲波形的采样，得到闪烁脉冲数字波形；
P3：提取闪烁脉冲数字波形库中各个闪烁脉冲的下降沿的幅值变化情况，并生成对应的拟合基线值查询表和拟合能量值查询表；
P4：将上述表存储在数据采集系统中。
3.根据权利要求1所述的脉冲堆积事件实时处理方法，其特征在于：所述步骤S2包括步骤：
P1：数据采集系统对高计数下的闪烁脉冲电信号进行采集，由于闪烁脉冲具有一定持续时间，且出现时间不固定，因此相距过近的两个或多个闪烁脉冲将发生脉冲堆积；
P2：闪烁脉冲电信号分别被MVT方法采样，MVT方法利用多个电压阈值与闪烁脉冲波形快速上升沿进行电压比较，无论该脉冲是否为堆积脉冲都能在其上升沿得到电平翻转信号；
P3：将MVT方法得到的电平翻转信号作为脉冲的识别信号，触发RTS方法对脉冲进行采样和数字积分求取脉冲能量信息。
4.根据权利要求1所述的脉冲堆积事件实时处理方法，其特征在于：所述步骤S3具体包括步骤：
a.在数据采集系统进行高计数下的数据采样过程中，当有堆积事件发生时，前一个MVT触发信号作为前一个脉冲的数字积分开始信号；
b.当第二个MVT触发信号发生时，利用堆积脉冲中前一脉冲未堆积部分脉冲下降沿作为输入参数查找拟合能量值查询表以得到前一脉冲堆积部分能量值，并与已采集的未堆积部分能量值加和作为前一脉冲的能量值；
c.同时利用前一脉冲未堆积部分脉冲下降沿作为输入参数查询拟合基线值查询表以得到前一脉冲在堆积发生后的脉冲幅值，并在每个采样时刻将该脉冲幅值按照脉冲下降沿衰减规律进行更新；
d.在后一脉冲的数字积分过程中将在每个采样时刻减去当前的前一脉冲的幅值更新值，以得到后一脉冲的真实幅值；
e.如果在后一个脉冲的数字积分时间内再次遇到新的MVT触发信号，则将后一个脉冲视作前一个脉冲按照上述处理方法迭代进行。
5.一种脉冲堆积事件实时处理系统，其特征在于：包括：
触发模块，用以当MVT触发信号到来时，触发模块产生积分触发信号，并将前一脉冲已采样下降沿波形更新到当前基线值；
数字积分模块，用以当接收到积分触发信号后将RTS采样值减去基线校正值作为真实的脉冲幅值进行积分；
拟合能量值查询模块，用以根据当前基线值和已积分时间来给出脉冲未积分部分能量的拟合值。
拟合基线值查询模块，用以根据当前基线值和已积分时间来获取前一个脉冲因堆积而未采集的波形幅值，并作为后一个脉冲的基线校正值。
6.根据权利要求5所述的脉冲堆积事件实时处理系统，其特征在于：所述数字积分模块中如果在数字积分周期内未发生第二次MVT触发信号，认为无堆积事件发生，数字积分周期完成时输出的脉冲能量值即为当前MVT触发信号对应的脉冲能量值。
7.根据权利要求5所述的脉冲堆积事件实时处理系统，其特征在于：所述数字积分模块中如果在数字积分周期内发生第二次MVT触发信号，数字积分模块用拟合信息查询模块得到的未积分部分的拟合能量值与第二次MVT触发信号前的已完成积分部分能量值相加作为当前MVT触发信号对应的脉冲能量值。
Priority Applications (6)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510751817.7A CN105212954B (zh) | 2015-11-05 | 2015-11-05 | 一种脉冲堆积事件实时处理方法与系统 |
| PCT/CN2016/104499 WO2017076312A1 (zh) | 2015-11-05 | 2016-11-04 | 一种脉冲堆积事件实时处理方法与系统 |
| JP2018522504A JP6657398B2 (ja) | 2015-11-05 | 2016-11-04 | パルスパイルアップイベントのリアルタイム処理の方法及びシステム |
| FIEP16861580.5T FI3361288T3 (fi) | 2015-11-05 | 2016-11-04 | Menetelmä ja järjestelmä pulssin kasaantumistapahtuman tosiaikaista käsittelyä varten |
| US15/773,997 US10228470B2 (en) | 2015-11-05 | 2016-11-04 | Method and system for real-time processing of pulse pile-up event |
| EP16861580.5A EP3361288B1 (en) | 2015-11-05 | 2016-11-04 | Method and system for real-time processing of pulse pile-up event |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510751817.7A CN105212954B (zh) | 2015-11-05 | 2015-11-05 | 一种脉冲堆积事件实时处理方法与系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105212954A true CN105212954A (zh) | 2016-01-06 |
| CN105212954B CN105212954B (zh) | 2018-03-16 |
Family
ID=54982499
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510751817.7A Active CN105212954B (zh) | 2015-11-05 | 2015-11-05 | 一种脉冲堆积事件实时处理方法与系统 |
Country Status (6)
| Country | Link |
|---|---|
| US (1) | US10228470B2 (zh) |
| EP (1) | EP3361288B1 (zh) |
| JP (1) | JP6657398B2 (zh) |
| CN (1) | CN105212954B (zh) |
| FI (1) | FI3361288T3 (zh) |
| WO (1) | WO2017076312A1 (zh) |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106125127A (zh) * | 2016-06-18 | 2016-11-16 | 武汉京邦科技有限公司 | 一种硅光电倍增器的时间标记方法及其微元阵列编码系统 |
| WO2017076312A1 (zh) * | 2015-11-05 | 2017-05-11 | 苏州瑞派宁科技有限公司 | 一种脉冲堆积事件实时处理方法与系统 |
| CN106842277A (zh) * | 2017-01-22 | 2017-06-13 | 沈阳东软医疗系统有限公司 | 一种堆积事件处理方法及装置 |
| CN106821409A (zh) * | 2017-01-22 | 2017-06-13 | 沈阳东软医疗系统有限公司 | 堆积事件处理方法及装置 |
| CN107462917A (zh) * | 2017-09-19 | 2017-12-12 | 南京瑞派宁信息科技有限公司 | 一种电离辐射探测装置及方法 |
| CN107976706A (zh) * | 2016-10-25 | 2018-05-01 | 上海东软医疗科技有限公司 | 一种pet系统的计数丢失校正方法和装置 |
| CN109669206A (zh) * | 2019-03-03 | 2019-04-23 | 南昌华亮光电有限责任公司 | 环绕式液体闪烁体智能能谱仪系统及其信号处理方法 |
| WO2020042664A1 (zh) * | 2018-08-27 | 2020-03-05 | 苏州瑞派宁科技有限公司 | 脉冲信号的采样方法、装置和计算机程序介质 |
| CN110954934A (zh) * | 2019-10-24 | 2020-04-03 | 中国船舶重工集团公司第七一九研究所 | 一种自适应核脉冲堆积信号峰值提取方法 |
| CN111158039A (zh) * | 2020-01-02 | 2020-05-15 | 苏州瑞派宁科技有限公司 | 信号采样、重建方法及装置 |
| CN112587161A (zh) * | 2020-12-09 | 2021-04-02 | 明峰医疗系统股份有限公司 | Pet成像设备的堆积信号恢复方法、pet成像设备及计算机可读存储介质 |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111538067B (zh) * | 2020-05-06 | 2022-09-06 | 东华理工大学 | 一种数字化核脉冲直线成形方法 |
| CN112882081B (zh) * | 2021-01-13 | 2022-12-23 | 明峰医疗系统股份有限公司 | Pet探测器全通道增益校准方法 |
| CN115545064B (zh) * | 2021-12-09 | 2023-05-23 | 西安中核核仪器股份有限公司 | 一种脉冲信号堆积识别方法 |
| CN114791620B (zh) * | 2022-03-28 | 2024-08-16 | 西北核技术研究所 | 一种针对离散数字核脉冲信号的梯形成形方法 |
| CN115951392A (zh) * | 2022-12-17 | 2023-04-11 | 苏州瑞派宁科技有限公司 | 能量校正方法、装置、电子设备、探测设备及存储介质 |
| CN117235436B (zh) * | 2023-11-16 | 2024-03-22 | 苏州瑞派宁科技有限公司 | 脉冲数字化方法、校正方法、装置、电子设备及存储介质 |
| CN117972536B (zh) * | 2024-04-01 | 2024-07-09 | 成都大学 | 一种脉冲分类方法及系统 |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6903344B2 (en) * | 2003-03-25 | 2005-06-07 | Cti Pet Systems, Inc. | Baseline correction in PET utilizing continuous sampling ADCs to compensate for DC and count rate errors |
| US7552017B1 (en) * | 2005-10-06 | 2009-06-23 | Lawrence Livermore National Security, Llc | Tailpulse signal generator |
| CN101561507A (zh) * | 2009-05-21 | 2009-10-21 | 清华大学 | 射线探测器的射线能量检测方法 |
| CN102073059A (zh) * | 2010-12-31 | 2011-05-25 | 华中科技大学 | 一种数字化pileup波形处理方法及系统 |
| CN102262238A (zh) * | 2011-04-19 | 2011-11-30 | 苏州瑞派宁科技有限公司 | 一种提取闪烁脉冲信息的方法及装置 |
| CN102843139A (zh) * | 2012-09-20 | 2012-12-26 | 苏州瑞派宁科技有限公司 | 一种闪烁脉冲数字化的方法及装置 |
| CN103454671A (zh) * | 2013-08-21 | 2013-12-18 | 中国人民解放军第二炮兵工程大学 | 一种基于高速数字采样的核辐射脉冲堆积判断与校正方法 |
| CN103969675A (zh) * | 2013-02-05 | 2014-08-06 | 苏州瑞派宁科技有限公司 | 数字化闪烁脉冲的基线校正方法及系统 |
| CN104101894A (zh) * | 2013-04-07 | 2014-10-15 | 苏州瑞派宁科技有限公司 | 一种堆积事件侦测与重建的方法及系统 |
| WO2017076312A1 (zh) * | 2015-11-05 | 2017-05-11 | 苏州瑞派宁科技有限公司 | 一种脉冲堆积事件实时处理方法与系统 |
Family Cites Families (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5210423A (en) * | 1989-09-22 | 1993-05-11 | Siemens Gammasonics, Inc. | Method and apparatus for unpiling pulses generated by piled-up scintillation events |
| EP0424651A1 (en) * | 1989-09-22 | 1991-05-02 | Siemens Aktiengesellschaft | Method and apparatus for unpiling pulses generated by piled-up scintillation events |
| JP3152476B2 (ja) * | 1992-01-28 | 2001-04-03 | 浜松ホトニクス株式会社 | 放射線検出装置 |
| US5270547A (en) * | 1992-05-07 | 1993-12-14 | Independent Scintillation Imaging Systems (Isis) Inc. | Scintillation camera valid event discrimination |
| US6252232B1 (en) * | 1998-06-29 | 2001-06-26 | General Electric Company | Digital integrator |
| US6525323B1 (en) * | 2000-04-18 | 2003-02-25 | Koninklijke Philips Electronics, N.V. | Method and apparatus for improved estimation of characteristics of pulses detected by a nuclear camera |
| JP4160275B2 (ja) | 2001-05-28 | 2008-10-01 | 浜松ホトニクス株式会社 | エネルギー測定方法及び測定装置 |
| CN101124725B (zh) * | 2004-09-16 | 2012-06-20 | 南方创新国际私人有限公司 | 用于从检测器输出数据中分解单个信号的方法和设备 |
| US7211800B2 (en) * | 2005-09-29 | 2007-05-01 | Siemens Medical Solutions Usa, Inc. | Afterglow DC-offset calibration in a nuclear imaging system |
| CN104586414B (zh) * | 2014-12-29 | 2017-03-29 | 沈阳东软医疗系统有限公司 | 一种堆积脉冲信号的处理方法和装置 |
-
2015
- 2015-11-05 CN CN201510751817.7A patent/CN105212954B/zh active Active
-
2016
- 2016-11-04 JP JP2018522504A patent/JP6657398B2/ja active Active
- 2016-11-04 EP EP16861580.5A patent/EP3361288B1/en active Active
- 2016-11-04 US US15/773,997 patent/US10228470B2/en active Active
- 2016-11-04 WO PCT/CN2016/104499 patent/WO2017076312A1/zh active Application Filing
- 2016-11-04 FI FIEP16861580.5T patent/FI3361288T3/fi active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6903344B2 (en) * | 2003-03-25 | 2005-06-07 | Cti Pet Systems, Inc. | Baseline correction in PET utilizing continuous sampling ADCs to compensate for DC and count rate errors |
| US7552017B1 (en) * | 2005-10-06 | 2009-06-23 | Lawrence Livermore National Security, Llc | Tailpulse signal generator |
| CN101561507A (zh) * | 2009-05-21 | 2009-10-21 | 清华大学 | 射线探测器的射线能量检测方法 |
| CN102073059A (zh) * | 2010-12-31 | 2011-05-25 | 华中科技大学 | 一种数字化pileup波形处理方法及系统 |
| CN102262238A (zh) * | 2011-04-19 | 2011-11-30 | 苏州瑞派宁科技有限公司 | 一种提取闪烁脉冲信息的方法及装置 |
| CN102843139A (zh) * | 2012-09-20 | 2012-12-26 | 苏州瑞派宁科技有限公司 | 一种闪烁脉冲数字化的方法及装置 |
| CN103969675A (zh) * | 2013-02-05 | 2014-08-06 | 苏州瑞派宁科技有限公司 | 数字化闪烁脉冲的基线校正方法及系统 |
| CN104101894A (zh) * | 2013-04-07 | 2014-10-15 | 苏州瑞派宁科技有限公司 | 一种堆积事件侦测与重建的方法及系统 |
| CN103454671A (zh) * | 2013-08-21 | 2013-12-18 | 中国人民解放军第二炮兵工程大学 | 一种基于高速数字采样的核辐射脉冲堆积判断与校正方法 |
| WO2017076312A1 (zh) * | 2015-11-05 | 2017-05-11 | 苏州瑞派宁科技有限公司 | 一种脉冲堆积事件实时处理方法与系统 |
Cited By (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2017076312A1 (zh) * | 2015-11-05 | 2017-05-11 | 苏州瑞派宁科技有限公司 | 一种脉冲堆积事件实时处理方法与系统 |
| US10228470B2 (en) | 2015-11-05 | 2019-03-12 | Raycan Technology Co., Ltd. (Su Zhou) | Method and system for real-time processing of pulse pile-up event |
| CN106125127A (zh) * | 2016-06-18 | 2016-11-16 | 武汉京邦科技有限公司 | 一种硅光电倍增器的时间标记方法及其微元阵列编码系统 |
| CN107976706B (zh) * | 2016-10-25 | 2019-10-29 | 上海东软医疗科技有限公司 | 一种pet系统的计数丢失校正方法和装置 |
| CN107976706A (zh) * | 2016-10-25 | 2018-05-01 | 上海东软医疗科技有限公司 | 一种pet系统的计数丢失校正方法和装置 |
| US10360699B2 (en) | 2016-10-25 | 2019-07-23 | Shenyang Neusoft Medical Systems Co., Ltd. | Correcting count loss |
| CN106842277A (zh) * | 2017-01-22 | 2017-06-13 | 沈阳东软医疗系统有限公司 | 一种堆积事件处理方法及装置 |
| CN106821409A (zh) * | 2017-01-22 | 2017-06-13 | 沈阳东软医疗系统有限公司 | 堆积事件处理方法及装置 |
| CN106821409B (zh) * | 2017-01-22 | 2020-10-13 | 东软医疗系统股份有限公司 | 堆积事件处理方法及装置 |
| CN107462917A (zh) * | 2017-09-19 | 2017-12-12 | 南京瑞派宁信息科技有限公司 | 一种电离辐射探测装置及方法 |
| WO2020042664A1 (zh) * | 2018-08-27 | 2020-03-05 | 苏州瑞派宁科技有限公司 | 脉冲信号的采样方法、装置和计算机程序介质 |
| US11944470B2 (en) | 2018-08-27 | 2024-04-02 | Raycan Technology Co., Ltd. (Suzhou) | Method and device for sampling a pulse signal, and computer program medium |
| CN109669206A (zh) * | 2019-03-03 | 2019-04-23 | 南昌华亮光电有限责任公司 | 环绕式液体闪烁体智能能谱仪系统及其信号处理方法 |
| CN110954934A (zh) * | 2019-10-24 | 2020-04-03 | 中国船舶重工集团公司第七一九研究所 | 一种自适应核脉冲堆积信号峰值提取方法 |
| CN111158039A (zh) * | 2020-01-02 | 2020-05-15 | 苏州瑞派宁科技有限公司 | 信号采样、重建方法及装置 |
| WO2021135338A1 (zh) * | 2020-01-02 | 2021-07-08 | 苏州瑞派宁科技有限公司 | 信号采样、重建方法及装置 |
| CN112587161A (zh) * | 2020-12-09 | 2021-04-02 | 明峰医疗系统股份有限公司 | Pet成像设备的堆积信号恢复方法、pet成像设备及计算机可读存储介质 |
Also Published As
| Publication number | Publication date |
|---|---|
| US20180321391A1 (en) | 2018-11-08 |
| EP3361288A1 (en) | 2018-08-15 |
| EP3361288A4 (en) | 2019-06-26 |
| EP3361288B1 (en) | 2024-04-03 |
| WO2017076312A1 (zh) | 2017-05-11 |
| CN105212954B (zh) | 2018-03-16 |
| JP2018538525A (ja) | 2018-12-27 |
| FI3361288T3 (fi) | 2024-04-17 |
| US10228470B2 (en) | 2019-03-12 |
| JP6657398B2 (ja) | 2020-03-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105212954A (zh) | 一种脉冲堆积事件实时处理方法与系统 | |
| CN102073059B (zh) | 一种数字化pileup波形处理方法及系统 | |
| US9772408B2 (en) | Method and device for extracting scintillation pulse information | |
| CN103969675A (zh) | 数字化闪烁脉冲的基线校正方法及系统 | |
| CN103858022B (zh) | 数据获取装置、对湮没伽马射线检测器及对湮没伽马射线检测方法 | |
| US6936822B2 (en) | Method and apparatus to prevent signal pile-up | |
| CN102141632B (zh) | 提高γ射线检测中的定时分辨率的装置和相关方法 | |
| US10175369B2 (en) | Signal data processing for radiation detection | |
| CN109171787B (zh) | 脉冲信号的采样方法、装置和计算机程序介质 | |
| CN102307046B (zh) | 一种时间分辨光子计数成像系统及方法 | |
| Streun et al. | Coincidence detection by digital processing of free-running sampled pulses | |
| CN102323959B (zh) | 用于时间分辨光子计数成像的采集卡 | |
| Wang et al. | Advantages of digitally sampling scintillation pulses in pileup processing in PET | |
| CN102313574B (zh) | 基于延迟线的信号测量通道合并方法及系统 | |
| CN104586414A (zh) | 一种堆积脉冲信号的处理方法和装置 | |
| CN1444733A (zh) | 具有高计数率的伽马照相机系统 | |
| CN106842277A (zh) | 一种堆积事件处理方法及装置 | |
| Annand | A fast module for pulse shape analysis | |
| CN110954934A (zh) | 一种自适应核脉冲堆积信号峰值提取方法 | |
| CN107767427B (zh) | 一种信号波形恢复方法及装置 | |
| Anderson et al. | A digital data acquisition system for the detectors at Gammasphere | |
| Deng et al. | Empirical bayesian energy estimation for multi-voltage threshold digitizer in PET | |
| CN111505696A (zh) | 一种双时步pileup波形处理系统及方法 | |
| CN109491960A (zh) | 一种减小图像畸变的位置读出电路 | |
| US9778383B2 (en) | Transmission of PET-signals by means of time division multiplexing |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |