CN109171787B - 脉冲信号的采样方法、装置和计算机程序介质 - Google Patents

脉冲信号的采样方法、装置和计算机程序介质 Download PDF

Info

Publication number
CN109171787B
CN109171787B CN201810982519.2A CN201810982519A CN109171787B CN 109171787 B CN109171787 B CN 109171787B CN 201810982519 A CN201810982519 A CN 201810982519A CN 109171787 B CN109171787 B CN 109171787B
Authority
CN
China
Prior art keywords
sampling
time
preset
pulse signal
threshold
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
CN201810982519.2A
Other languages
English (en)
Other versions
CN109171787A (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.)
Raycan Technology Co Ltd
Original Assignee
Raycan Technology 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 Raycan Technology Co Ltd filed Critical Raycan Technology Co Ltd
Priority to CN201810982519.2A priority Critical patent/CN109171787B/zh
Publication of CN109171787A publication Critical patent/CN109171787A/zh
Priority to PCT/CN2019/086239 priority patent/WO2020042664A1/zh
Priority to US17/256,900 priority patent/US11944470B2/en
Priority to JP2020568696A priority patent/JP7265791B2/ja
Priority to EP19855793.6A priority patent/EP3796048B1/en
Priority to HUE19855793A priority patent/HUE063109T2/hu
Application granted granted Critical
Publication of CN109171787B publication Critical patent/CN109171787B/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
    • 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)
    • 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
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R19/00—Arrangements for measuring currents or voltages or for indicating presence or sign thereof
    • G01R19/04—Measuring peak values or amplitude or envelope of ac or of pulses
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R29/00—Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
    • G01R29/02—Measuring characteristics of individual pulses, e.g. deviation from pulse flatness, rise time or duration
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R29/00—Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
    • G01R29/02—Measuring characteristics of individual pulses, e.g. deviation from pulse flatness, rise time or duration
    • G01R29/027—Indicating that a pulse characteristic is either above or below a predetermined value or within or beyond a predetermined range of values
    • 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
    • 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
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10072—Tomographic images
    • G06T2207/10104—Positron emission tomography [PET]

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Medical Informatics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Molecular Biology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Public Health (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • Optics & Photonics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Pathology (AREA)
  • Veterinary Medicine (AREA)
  • Biophysics (AREA)
  • Theoretical Computer Science (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Power Engineering (AREA)
  • Measurement Of Radiation (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
  • Nuclear Medicine (AREA)

Abstract

本申请公开了脉冲信号的采样方法、装置和计算机程序介质。该采样方法可以包括：步骤S1，根据多个预设采样阈值在脉冲信号的上升沿部分上采集多个第一采样点，每个第一采样点均由一个预设采样阈值以及与该预设采样阈值对应的第一时间来表征；步骤S2，根据多个预设采样阈值在脉冲信号的下降沿部分上采集多个第二采样点，每个第二采样点均由一个预设采样阈值以及与一个该预设采样阈值对应的第二时间来表征；以及步骤S3，在与第一时间或第二时间相隔预设时间间隔的第三时间处，在脉冲信号上采集第三采样点，第三采样点由第三时间以及与第三时间对应的响应幅度来表征。通过本申请提供的技术方案，可以提高脉冲信号的采样效率并且降低能耗。

Description

脉冲信号的采样方法、装置和计算机程序介质
技术领域
本申请涉及信号处理技术领域，特别涉及一种脉冲信号的采样方法、装置和计算机程序介质。
背景技术
正电子发射断层成像(Positron Emission Tomography，简称PET)是一种利用放射性元素进行临床显像的技术，其过程为：将发射正电子的放射性核素标记到能够参与活体组织血流或代谢过程的化合物上，再将标有放射性核素的化合物注射到受检者体内。放射性核素在体内发射出的正电子移动大约1mm后与受检者体内的负电子结合发生电子对的湮灭，产生γ光子，释放出的γ光子可以被闪烁晶体接受转换为可见光，再通过光电倍增管元件转换为脉冲信号以进行重建，从而帮助确定放射性核素的富集部位，并且帮助定位代谢旺盛区域并进行活度评估。
在PET或其它相关领域中，在对脉冲信号进行重建之前，需要对脉冲信号进行采样处理。现有技术中主要利用时间间隔采样方法来进行采样处理，具体地，按照相同的间隔时间，记录脉冲信号在每个时间点对应的幅值，每个时间点和与其对应的幅值便构成了一个采样点。
发明内容
在实现本申请过程中，发明人发现现有技术中至少存在如下问题：
现有技术中的时间间隔采样方法采样点多、记录信息多，所以其采样时间长，这导致其采样效率较低并且能耗较大。
为解决上述技术问题，本申请实施例提供的一种脉冲信号的采样方法、装置和计算机程序介质是这样实现的：
一种脉冲信号的采样方法，包括以下步骤：
步骤S1，根据多个预设采样阈值在脉冲信号的上升沿部分上采集多个第一采样点，每个所述第一采样点均由一个所述预设采样阈值以及与一个所述预设采样阈值对应的第一时间来表征；
步骤S2，根据多个所述预设采样阈值在所述脉冲信号的下降沿部分上采集多个第二采样点，每个所述第二采样点均由一个所述预设采样阈值以及与一个所述预设采样阈值对应的第二时间来表征；
步骤S3，在与所述第一时间或所述第二时间相隔预设时间间隔的第三时间处，在所述脉冲信号上采集第三采样点，所述第三采样点由所述第三时间以及与所述第三时间对应的响应幅度来表征。
优选地，所述步骤S1包括：
记录在所述脉冲信号的上升沿部分上信号幅度达到多个所述预设采样阈值时的多个所述第一时间，以确定出多个所述第一采样点。
优选地，所述步骤S1包括：
当一个所述预设采样阈值对应有多个所述第一时间时，计算多个所述第一时间中的最大时间和最小时间的平均值，或者计算多个所述第一时间的平均值，并将计算得到的平均值作为与一个所述预设采样阈值对应的所述第一时间。
优选地，所述预设采样阈值包括电压阈值、电流阈值或磁场强度阈值。
优选地，所述预设采样阈值是根据所述脉冲信号的触发特性或预先统计的幅度经验值来设定的。
优选地，所述步骤S3包括：
在与所述第一时间相隔预设时间间隔的第三时间处，在所述脉冲信号中多个所述第一采样点后面的上升沿部分、波峰或下降沿部分采集所述第三采样点，或者，
在与所述第二时间相隔预设时间间隔的第三时间处，在所述脉冲信号中多个所述第二采样点后面的下降沿部分或波谷采集所述第三采样点。
优选地，所述步骤S3包括：
在所述第三时间处，确定所述脉冲信号上与所述第三时间对应的响应幅度以确定出所述第三采样点。
优选地，所述预设时间间隔根据以下方式来确定：
分别计算预先采集的多个脉冲信号中的脉冲幅值点的横坐标与所选择的标定点的横坐标之间的差值，所述标定点为多个所述第一采样点中的一个或者多个所述第二采样点中的一个；
对计算得到的多个所述差值进行高斯拟合，得到高斯曲线和拟合参数；
比较所得到的高斯曲线的离散程度，将符合最佳正态分布的高斯曲线所对应的拟合参数作为所述预设时间间隔。
一种脉冲信号的采样装置，包括：
第一采样单元，其被配置为根据多个预设采样阈值在脉冲信号的上升沿部分上采集多个第一采样点，每个所述第一采样点均由一个所述预设采样阈值以及与一个所述预设采样阈值对应的第一时间来表征；
第二采样单元，其被配置为根据多个所述预设采样阈值在所述脉冲信号的下降沿部分上采集多个第二采样点，每个所述第二采样点均由一个所述预设采样阈值以及与一个所述预设采样阈值对应的第二时间来表征；以及
第三采样单元，其被配置为在与所述第一时间或所述第二时间相隔预设时间间隔的第三时间处，在所述脉冲信号上采集第三采样点，所述第三采样点由所述第三时间以及与所述第三时间对应的响应幅度来表征。
一种计算机程序介质，包括：
存储器，其上存储有指令，
处理器，其与所述存储器联接，并且被配置为根据所述存储器存储的指令来执行以下操作：
根据多个预设采样阈值在脉冲信号的上升沿部分上采集多个第一采样点，每个所述第一采样点均由一个所述预设采样阈值以及与一个所述预设采样阈值对应的第一时间来表征；
根据多个所述预设采样阈值在所述脉冲信号的下降沿部分上采集多个第二采样点，每个所述第二采样点均由一个所述预设采样阈值以及与一个所述预设采样阈值对应的第二时间来表征；以及
在与所述第一时间或所述第二时间相隔预设时间间隔的第三时间处，在所述脉冲信号上采集第三采样点，所述第三采样点由所述第三时间以及与所述第三时间对应的响应幅度来表征。
由以上本申请实施例提供的技术方案可见，本申请实施例通过在脉冲信号的上升沿和下降沿部分进行多阈值采样，并且在多阈值采样的基础上对脉冲信号进行时间间隔采样，这可以实现提高脉冲信号的采样效率并且降低能耗的目的。另外，通过利用本申请实施例所提供的技术方案，可以很好地捕捉到脉冲信号的幅值变化，从而可以提高后续还原脉冲信号的波形的准确性。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1是根据本申请实施例的脉冲信号的采样方法的流程图；
图2是图1中的采样方法的具体操作流程图；
图3是利用图2中的采样方法进行采样的结果示意图；
图4是现有技术中的利用时间间隔方法来对脉冲信号进行采样的结果示意图；
图5是现有技术中的对脉冲信号进行多阈值采样后经过拟合处理所得到的拟合曲线与实际测量曲线的对比示意图；
图6是利用本申请实施例所提供的采样方法对脉冲信号进行采样后经过拟合处理所得到的拟合曲线与实际测量曲线的对比示意图；
图7是根据本申请实施例的脉冲信号的采样装置的结构示意图；
图8是根据本申请实施例的计算机程序介质的结构示意图。
具体实施方式
下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是用于解释说明本申请的一部分实施例，而不是全部的实施例，并不希望限制本申请的范围或权利要求书。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其它实施例，都应当属于本申请保护的范围。
需要说明的是，当元件被称为“设置在”另一个元件上，它可以直接设置在另一个元件上或者也可以存在居中的元件。当元件被称为“连接/联接”至另一个元件，它可以是直接连接/联接至另一个元件或者可能同时存在居中元件。本文所使用的术语“连接/联接”可以包括电气和/或机械物理连接/联接。本文所使用的术语“包括/包含”指特征、步骤或元件的存在，但并不排除一个或更多个其它特征、步骤或元件的存在或添加。本文所使用的术语“和/或”包括一个或多个相关所列项目的任意的和所有的组合。
除非另有定义，本文所使用的所有的技术和科学术语与属于本申请的技术领域的技术人员通常理解的含义相同。本文中所使用的术语只是为了描述具体实施例的目的，而并不是旨在限制本申请。
另外，在本申请的描述中，术语“第一”、“第二”、“第三”等仅用于描述目的和区别类似的对象，两者之间并不存在先后顺序，也不能理解为指示或暗示相对重要性。此外，在本申请的描述中，除非另有说明，“多个”的含义是两个或两个以上。
下面结合附图对本申请实施例所提供的脉冲信号的采样方法、装置和计算机程序介质进行详细说明。
如图1所示，本申请实施例提供了一种脉冲信号的采样方法，其包括以下步骤：
S1：根据多个预设采样阈值在脉冲信号的上升沿部分上采集多个第一采样点。
在主动获取或接收到来自闪烁晶体探测器或气体电离探测器等辐射探测器或其它装置的脉冲信号之后，可以对脉冲信号进行多阈值采样，例如，可以根据多个预设采样阈值在脉冲信号的上升沿部分上依次采集多个第一采样点。具体地，首先，可以记录在脉冲信号的上升沿部分上信号幅度达到多个预设采样阈值时的多个第一时间，从而可以确定出脉冲信号的上升沿部分上的多个第一采样点，进而可以记录所述多个第一采样点的相关信息，例如，幅度、时间等。每个第一采样点均可以由一个预设采样阈值以及与该预设采样阈值对应的第一时间来表征，例如，第一采样点的坐标可以用(Ti,Mi)来表示，其中，T表示时间，M表示幅度，i为正整数。
在脉冲信号的上升沿部分上，如果一个预设采样阈值对应有多个第一时间，则可以计算所述多个第一时间中的最大时间和最小时间的平均值，或者计算所述多个第一时间的平均值，并将计算得到的平均值作为与该预设采样阈值对应的第一时间。当然，也可以按照其它方式来确定第一时间，在此不进行限制。脉冲信号可以是电信号(例如，电脉冲信号)、光信号、声音信号等。预设采样阈值可以是电性阈值，例如电压阈值或电流阈值，也还可以是其它阈值，例如磁场强度等磁性阈值。预设采样阈值可以是根据脉冲信号的特性(例如，触发电平、脉冲幅值)来设定的，也可以是根据预先统计的脉冲信号的幅度经验值来设定的。对于预设采样阈值的个数，可以根据实际情况来选择，在此不进行限制。
S2：根据多个预设采样阈值在脉冲信号的下降沿部分上采集多个第二采样点。
在脉冲信号的上升沿部分上确定出与多个预设采样阈值对应的多个第一时间后，可以按照相同的方式对脉冲信号的下降沿部分进行多阈值采样。具体地，可以记录在脉冲信号的下降沿部分上信号幅度达到多个预设采样阈值时的多个第二时间，从而可以确定出脉冲信号的下降沿部分上的多个第二采样点，进而可以记录所述多个第二采样点的相关信息，例如，幅度、时间等。每个第二采样点均可以由一个预设采样阈值以及与该预设采样阈值对应的第二时间来表征。
同样，如果一个预设采样阈值对应有多个第二时间，则可以计算所述多个第二时间中的最大时间和最小时间的平均值，或者计算所述多个第二时间的平均值，将计算得到的平均值作为与该预设采样阈值对应的第二时间。当然，也可以按照其它方式来确定第二时间，在此不进行限制。
关于该步骤的详细描述，可以参照上述步骤S1，在此不再赘叙。
S3：在与第一时间或第二时间相隔预设时间间隔的第三时间处，在脉冲信号上采集第三采样点。
在脉冲信号的上升沿部分上确定出与多个预设采样阈值对应的多个第一时间或者在其下降沿部分上确定出与多个预设采样阈值对应的多个第二时间后，可以对脉冲信号再进行时间采样，以提高采样的精准度，并且有助于确定出脉冲信号上的最大幅值点或最低幅值点。
在具体实施方式中，可以以任意一个或多个第一采样点(例如，所记录的最后一个第一采样点，也可以是其它采样点)作为标定点，在与其对应的第一时间相隔预设时间间隔(△)的第三时间处，确定脉冲信号中与第三时间对应的响应幅度，从而确定出脉冲信号中的多个第一采样点后面的上升沿部分、波峰或下降沿部分上的第三采样点，进而可以记录第三采样点的相关信息，例如，幅度、时间等。第三时间可以是第一时间与预设时间间隔(△)之和。
在另一具体实施方式中，也可以以任意一个或多个第二采样点(例如，所记录的第三个第二采样点，也可以是其它采样点)作为标定点，在与其对应的第二时间相隔预设时间间隔(△)的第三时间处，确定脉冲信号中与第三时间对应的响应幅度，从而确定出脉冲信号中的多个第二采样点后面的波谷或下降沿部分上的第三采样点。第三时间可以是第二时间与预设时间间隔(△)之和。
第三采样点优选的为脉冲幅值点(例如，波峰)，但也可以是其它数据点，其可以是一个，也可以是多个。每个第三采样点均可以由一个第三时间以及与该第三时间对应的响应幅度来表征。
预设时间间隔可以是任意设置的，也可以是预先通过实验来确定的，以提高采样效果。例如，以波峰为第三采样点时，可以通过以下方式来确定预设时间间隔，但不限于以下这种方式：
(1)预先采集多个(例如，10万个)脉冲，记录每个脉冲信号的上升沿部分的第一采样点和脉冲幅值点(即，波峰)；
(2)分别计算所述多个脉冲信号中的脉冲幅值点的横坐标与所选择的标定点的横坐标之间的差值，该标定点为多个第一采样点中的一个；
(3)对计算得到的多个差值进行高斯拟合(具体过程可以参照现有技术)，得到高斯曲线和拟合参数；
(4)比较所得到的高斯曲线的离散程度，将符合最佳正态分布的高斯曲线所对应的拟合参数作为所述预设时间间隔。
同样，也可以以波谷为第三采样点，按照上述方式来确定预设时间间隔，此时标定点为多个第二采样点中的一个。
需要说明的是，上述步骤S1至步骤S3并不一定需要严格按照上述顺序来执行。步骤S3也可以与步骤S1交叉执行，例如，在步骤S1中采集到第一个第一采样点后，便可以执行步骤S3。另外，对于第三时间与第一时间相隔预设时间间隔的情况，步骤S3可以在步骤S2之前执行，也可以与步骤S2交叉执行。
所采集的多个第一采样点、多个第二采样点和第三采样点构成了脉冲信号的所有采样点。通过这些记录采样点的相关信息，便可以得知脉冲信号的特征信息，例如，幅度大小、上升斜率、上升时间、下降时间等，从而可以为后续进行分类拟合提供依据。
下面以预设采样阈值为电压阈值并且采集9个采样点为例来说明上述步骤的执行过程。
如图2所示，按照电压阈值来进行采样：可以利用预先设置的4个电压阈值(V1，V2，V3，V4)，在脉冲信号的上升沿部分确定出对应的第一时间(T1，T2，T3，T4)，从而得到4个第一采样点(T1，V1)、(T2，V2)、(T3，V3)和(T4，V4)，在脉冲信号的下降沿部分确定出对应的第二时间(T5，T6，T7，T8)，从而得到4个第二采样点(T5，V4)、(T6，V3)、(T7，V2)和(T8，V1)；按照时间间隔来进行采样：可以在第4个第一采样点后通过预设时间间隔来采样得到第三采样点(T9，V9)，如图3所示。上述9个采样点构成了脉冲信号的所有采样点。
通过以上描述可以看出，本申请实施例通过在脉冲信号的上升沿和下降沿部分进行多阈值采样，并且在多阈值采样的基础上对脉冲信号进行时间间隔采样，这可以提高采样效率并且降低硬件的能耗，从而降低成本，而且还可以为后续对采样点进行分类提供依据。另外，在多阈值采样的基础上对脉冲信号进行时间间隔采样，可以让原脉冲信号保留的信息更多，从而降低了后续重建脉冲信号的难度，提高了后续重建脉冲信号的精度。另外，通过利用本申请实施例所提供的技术方案，可以很好地捕捉到脉冲信号的幅值变化，从而可以提高后续重建脉冲信号的波形的准确性。
下面以具体实例来说明本申请实施例的有益效果。
图4为现有技术中按照时间间隔来对脉冲信号进行采样的结果示意图，图5为现有技术中对脉冲信号进行多阈值采样后经过拟合处理所得到的拟合曲线与实际测量曲线的对比示意图；图6为利用本申请实施例所提供的采样方法对脉冲信号进行采样后经过拟合处理所得到的拟合曲线与实际测量曲线的对比示意图。其中，实线表示实际测量曲线，虚线表示拟合曲线，并且图5和图6中的拟合方法相同。
通过对比图4与图6可以看出，对于在同一时间段内对同一脉冲信号进行采样，本申请实施例所提供的采样方法所需采集的采样点的点数远远少于现有技术中的时间间隔采样方法所需采集的采样点的点数。本申请实施例所提供的采样方法只需要采集数个(例如，9个)采样点，而现有技术中的时间间隔采样方法起码需要采集几百个采样点，甚至成千上万个采样点。由此可见，利用本申请实施例所提供的采样方法可以减少采样点数，从而可以提高采样效率并且降低成本。
通过对比图5与图6可以看出，相对于现有技术中的多阈值采样方法，利用本申请实施例所提供的采样方法能使后续拟合得到的拟合曲线与实际测量曲线更加接近。由此可见，利用本申请实施例所提供的采样方法可以提高后续重建脉冲信号的波形的准确性，并且提高后续重建信号的精度。
本申请实施例还提供了一种脉冲信号的采样装置，如图7所示。该采样装置可以包括：
第一采样单元101，其可以被配置为根据多个预设采样阈值在脉冲信号的上升沿部分上采集多个第一采样点，每个第一采样点均由一个预设采样阈值以及与该预设采样阈值对应的第一时间来表征；
第二采样单元102，其可以被配置为根据多个预设采样阈值在脉冲信号的下降沿部分上采集多个第二采样点，每个第二采样点均由一个预设采样阈值以及与该预设采样阈值对应的第二时间来表征；
第三采样单元103，其可以被配置为在与第一时间或第二时间相隔预设时间间隔的第三时间处，在脉冲信号上采集第三采样点，第三采样点由第三时间以及与该第三时间对应的响应幅度来表征。
在主动获取或接收到来自闪烁晶体探测器或气体电离探测器等辐射探测器的脉冲信号之后，第一采样单元101可以记录在脉冲信号的上升沿部分上信号幅度达到多个预设采样阈值时的多个第一时间，从而可以确定出脉冲信号的上升沿部分上的多个第一采样点，并记录所述多个第一采样点的相关信息。在第一采样单元101完成第一采样点的采集之后，第二采样单元102可以记录在脉冲信号的下降沿部分上信号幅度达到多个预设采样阈值时的多个第二时间，从而可以确定出脉冲信号的下降沿部分上的多个第二采样点，并记录所述多个第二采样点的相关信息。第三采样单元103也可以在与第一时间或第二时间相隔预设时间间隔的第三时间处，确定脉冲信号中与第三时间对应的响应幅度，从而可以确定出脉冲信号上的第三采样点，并记录第三采样点的相关信息。
关于该实施例的详细描述，可以参照图1中的采样方法来进行描述，在此不再赘叙。
通过利用以上装置，可以实现在脉冲信号的上升沿和下降沿部分进行多阈值采样，并且在多阈值采样的基础上对脉冲信号进行时间间隔采样，这可以提高采样效率并且降低硬件的能耗，从而降低成本，而且还可以为后续对采样点进行分类提供依据。另外，在多阈值采样的基础上对脉冲信号进行时间间隔采样，可以让原脉冲信号保留的信息更多，从而降低了后续重建脉冲信号的难度，也提高了后续重建脉冲信号的精度。另外，通过利用本申请实施例所提供的技术方案，可以很好地捕捉到脉冲信号的幅值变化，从而可以提高后续重建脉冲信号的波形的准确性。
本申请实施例还提供了一种用于对脉冲信号进行采样的计算机程序介质，如图8所示。该计算机程序介质可以包括：
存储器201，其上存储有指令，
处理器202，其与存储器201联接，并且可以被配置为根据存储器201存储的指令来执行以下操作：
根据多个预设采样阈值在脉冲信号的上升沿部分上采集多个第一采样点，其中，每个第一采样点均由其中一个预设采样阈值以及与该预设采样阈值对应的第一时间来表征；
根据多个预设采样阈值在脉冲信号的下降沿部分上采集多个第二采样点，其中，每个第二采样点均由其中一个预设采样阈值以及与该预设采样阈值对应的第二时间来表征；以及
在与第一时间或第二时间相隔预设时间间隔的第三时间处，在脉冲信号上采集第三采样点，其中，第三采样点由第三时间以及与该第三时间对应的响应幅度来表征。
关于该实施例的详细描述，可以参照图1中的采样方法来进行描述，在此不再赘叙。
需要说明的是，本申请实施例所提供的上述采样方法、装置及计算机程序介质并不限于应用于PET、多电压阈值采样(Multi-Voltage Threshold，简称MVT)等领域，其也可以应用于需要对脉冲信号进行采样的任何领域。
上述实施例阐明的装置、单元等，具体可以由计算机芯片和/或实体实现，或者由具有某种功能的产品来实现。为了描述的方便，描述以上装置时以功能分为各种单元分别描述。当然，在实施本申请时可以把各单元的功能在同一个或多个计算机芯片中实现。
虽然本申请提供了如上述实施例或流程图所述的方法操作步骤，但基于常规或者无需创造性的劳动在所述方法中可以包括更多或者更少的操作步骤。在逻辑性上不存在必要因果关系的步骤中，这些步骤的执行顺序不限于本申请实施例提供的执行顺序。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其它实施例的不同之处。
上述实施例是为便于该技术领域的普通技术人员能够理解和使用本申请而描述的。熟悉本领域技术的人员显然可以容易地对这些实施例做出各种修改，并把在此说明的一般原理应用到其它实施例中而不必经过创造性的劳动。因此，本申请不限于上述实施例，本领域技术人员根据本申请的揭示，不脱离本申请范畴所做出的改进和修改都应该在本申请的保护范围之内。

Claims (10)

1.一种脉冲信号的采样方法，其特征在于，所述采样方法包括以下步骤：
步骤S1，根据多个预设采样阈值在脉冲信号的上升沿部分上采集多个第一采样点，每个所述第一采样点均由一个所述预设采样阈值以及与一个所述预设采样阈值对应的第一时间来表征；
步骤S2，根据多个所述预设采样阈值在所述脉冲信号的下降沿部分上采集多个第二采样点，每个所述第二采样点均由一个所述预设采样阈值以及与一个所述预设采样阈值对应的第二时间来表征；
步骤S3，在与所述第一时间或所述第二时间相隔预设时间间隔的第三时间处，在所述脉冲信号上采集第三采样点，所述第三采样点由所述第三时间以及与所述第三时间对应的响应幅度来表征。
2.根据权利要求1所述的采样方法，其特征在于，所述步骤S1包括：
记录在所述脉冲信号的上升沿部分上信号幅度达到多个所述预设采样阈值时的多个所述第一时间，以确定出多个所述第一采样点。
3.根据权利要求1或2所述的采样方法，其特征在于，所述步骤S1包括：
当一个所述预设采样阈值对应有多个所述第一时间时，计算多个所述第一时间中的最大时间和最小时间的平均值，或者计算多个所述第一时间的平均值，并将计算得到的平均值作为与一个所述预设采样阈值对应的所述第一时间。
4.根据权利要求1所述的采样方法，其特征在于，所述预设采样阈值包括电压阈值、电流阈值或磁场强度阈值。
5.根据权利要求1所述的采样方法，其特征在于，所述预设采样阈值是根据所述脉冲信号的触发特性或预先统计的脉冲信号的幅度经验值来设定的。
6.根据权利要求1所述的采样方法，其特征在于，所述步骤S3包括：
在与所述第一时间相隔预设时间间隔的第三时间处，在所述脉冲信号中的多个所述第一采样点后面的上升沿部分、波峰或下降沿部分采集所述第三采样点，或者，
在与所述第二时间相隔预设时间间隔的第三时间处，在所述脉冲信号中的多个所述第二采样点后面的下降沿部分或波谷采集所述第三采样点。
7.根据权利要求1或6所述的采样方法，其特征在于，所述步骤S3包括：
在所述第三时间处，确定所述脉冲信号上与所述第三时间对应的响应幅度以确定出所述第三采样点。
8.根据权利要求1所述的采样方法，其特征在于，所述预设时间间隔根据以下方式来确定：
分别计算预先采集的多个脉冲信号中的脉冲幅值点的横坐标与所选择的标定点的横坐标之间的差值，所述标定点为多个所述第一采样点中的一个或者或多个所述第二采样点中的一个；
对计算得到的多个所述差值进行高斯拟合，得到高斯曲线和拟合参数；
比较所得到的高斯曲线的离散程度，将符合最佳正态分布的高斯曲线所对应的拟合参数作为所述预设时间间隔。
9.一种脉冲信号的采样装置，其特征在于，所述采样装置包括：
第一采样单元，其被配置为根据多个预设采样阈值在脉冲信号的上升沿部分上采集多个第一采样点，每个所述第一采样点均由一个所述预设采样阈值以及与一个所述预设采样阈值对应的第一时间来表征；
第二采样单元，其被配置为根据多个所述预设采样阈值在所述脉冲信号的下降沿部分上采集多个第二采样点，每个所述第二采样点均由一个所述预设采样阈值以及与一个所述预设采样阈值对应的第二时间来表征；以及
第三采样单元，其被配置为在与所述第一时间或所述第二时间相隔预设时间间隔的第三时间处，在所述脉冲信号上采集第三采样点，所述第三采样点由所述第三时间以及与所述第三时间对应的响应幅度来表征。
10.一种计算机程序介质，其特征在于，包括：
存储器，其上存储有指令，
处理器，其与所述存储器联接，并且被配置为根据所述存储器存储的指令来执行以下操作：
根据多个预设采样阈值在脉冲信号的上升沿部分上采集多个第一采样点，每个所述第一采样点均由一个所述预设采样阈值以及与一个所述预设采样阈值对应的第一时间来表征；
根据多个所述预设采样阈值在所述脉冲信号的下降沿部分上采集多个第二采样点，每个所述第二采样点均由一个所述预设采样阈值以及与一个所述预设采样阈值对应的第二时间来表征；以及
在与所述第一时间或所述第二时间相隔预设时间间隔的第三时间处，在所述脉冲信号上采集第三采样点，所述第三采样点由所述第三时间以及与所述第三时间对应的响应幅度来表征。
CN201810982519.2A 2018-08-27 2018-08-27 脉冲信号的采样方法、装置和计算机程序介质 Active CN109171787B (zh)

Priority Applications (6)

Application Number Priority Date Filing Date Title
CN201810982519.2A CN109171787B (zh) 2018-08-27 2018-08-27 脉冲信号的采样方法、装置和计算机程序介质
PCT/CN2019/086239 WO2020042664A1 (zh) 2018-08-27 2019-05-09 脉冲信号的采样方法、装置和计算机程序介质
US17/256,900 US11944470B2 (en) 2018-08-27 2019-05-09 Method and device for sampling a pulse signal, and computer program medium
JP2020568696A JP7265791B2 (ja) 2018-08-27 2019-05-09 パルス信号をサンプリングするための方法及びデバイス並びにコンピュータプログラム媒体
EP19855793.6A EP3796048B1 (en) 2018-08-27 2019-05-09 Method and device for sampling pulse signal, and computer program medium
HUE19855793A HUE063109T2 (hu) 2018-08-27 2019-05-09 Impulzusjel-mintavételi módszer és eszköz, valamint számítógépes programhordozó

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810982519.2A CN109171787B (zh) 2018-08-27 2018-08-27 脉冲信号的采样方法、装置和计算机程序介质

Publications (2)

Publication Number Publication Date
CN109171787A CN109171787A (zh) 2019-01-11
CN109171787B true CN109171787B (zh) 2021-02-26

Family

ID=64916375

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810982519.2A Active CN109171787B (zh) 2018-08-27 2018-08-27 脉冲信号的采样方法、装置和计算机程序介质

Country Status (6)

Country Link
US (1) US11944470B2 (zh)
EP (1) EP3796048B1 (zh)
JP (1) JP7265791B2 (zh)
CN (1) CN109171787B (zh)
HU (1) HUE063109T2 (zh)
WO (1) WO2020042664A1 (zh)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109171787B (zh) 2018-08-27 2021-02-26 苏州瑞派宁科技有限公司 脉冲信号的采样方法、装置和计算机程序介质
CN111158039B (zh) * 2020-01-02 2022-01-04 苏州瑞派宁科技有限公司 信号采样、重建方法及装置
CN113406877B (zh) * 2021-05-20 2022-09-23 西北核技术研究所 基于特征点的脉冲信号高时间精度测量方法和系统
CN113375765B (zh) * 2021-07-06 2024-06-11 金卡水务科技有限公司 计量仪表的检测方法、装置和存储介质
CN113906299A (zh) * 2021-09-03 2022-01-07 长江存储科技有限责任公司 用于自动波形分析的方法和系统
CN115113512B (zh) * 2022-05-26 2024-03-29 合肥综合性国家科学中心人工智能研究院(安徽省人工智能实验室) 脉冲的采样方法、采样系统、装置及计算机可读存储介质
CN115032877B (zh) * 2022-05-26 2024-03-29 合肥综合性国家科学中心人工智能研究院(安徽省人工智能实验室) 脉冲的采样方法、采样系统、装置及计算机可读存储介质
WO2024032432A1 (zh) * 2022-08-08 2024-02-15 深圳市恒运昌真空技术有限公司 一种射频电源信号采集方法及装置
CN115015797B (zh) * 2022-08-08 2022-10-14 深圳市恒运昌真空技术有限公司 一种射频电源信号采集方法及装置

Family Cites Families (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US1022847A (en) 1911-11-21 1912-04-09 Gustave J Jaccard Game.
DE3207556C2 (de) * 1982-03-03 1983-12-22 Vierling, Oskar, Prof. Dr.Phil.Habil., 8553 Ebermannstadt Anordnung zum Messen der charakteristischen Zeiten von Impulsen und Impulsserien
US4495586A (en) * 1982-07-29 1985-01-22 Tektronix, Inc. Waveform acquisition apparatus and method
DE4100250A1 (de) * 1991-01-07 1991-05-23 Klaus Brinker Verfahren zur messung des tastverhaeltnisses von taktsignalen
US5274569A (en) * 1991-10-15 1993-12-28 International Business Machines Corporation Dual sense non-differencing digital peak detector
US5446650A (en) * 1993-10-12 1995-08-29 Tektronix, Inc. Logic signal extraction
DE19948892C2 (de) * 1999-10-11 2002-07-18 Asm Automation Sensorik Messte Impulsdetektor und Verfahren zur Detektion von sinusförmigen Impulsen
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
WO2008054883A2 (en) * 2006-06-07 2008-05-08 Xrf Corporation Devices and methods for detecting and analyzing radiation
CN102688155A (zh) 2011-03-23 2012-09-26 日本乐敦制药株式会社 外用组合物
CN102262238B (zh) * 2011-04-19 2014-07-23 苏州瑞派宁科技有限公司 一种提取闪烁脉冲信息的方法及装置
CN102789952A (zh) * 2011-05-19 2012-11-21 安徽中科大建成海晟科技有限责任公司 一种基于脉冲前沿时间测量和幅度修正算法的飞行时间质谱仪电子学读出方法
WO2013085923A1 (en) * 2011-12-06 2013-06-13 The Board Of Trustees Of The Leland Stanford Junior University Time-based digitizer for pet photodetector
CN202801645U (zh) * 2012-09-20 2013-03-20 苏州瑞派宁科技有限公司 一种闪烁脉冲数字化的装置
CN102843139B (zh) * 2012-09-20 2015-10-28 苏州瑞派宁科技有限公司 一种闪烁脉冲数字化的方法及装置
CN103961126B (zh) * 2013-02-05 2016-07-06 苏州瑞派宁科技有限公司 一种多阈值采样数字化器件的阈值校正方法
CN104639123B (zh) 2013-11-14 2017-08-25 苏州瑞派宁科技有限公司 闪烁脉冲越过阈值的时间点获取方法及装置
CN104656119B (zh) * 2013-11-19 2018-06-05 苏州瑞派宁科技有限公司 一种闪烁脉冲信息复原的方法及系统
US11131756B2 (en) * 2015-09-29 2021-09-28 Qualcomm Incorporated LIDAR system with reflected signal strength measurement
CN105212954B (zh) 2015-11-05 2018-03-16 苏州瑞派宁科技有限公司 一种脉冲堆积事件实时处理方法与系统
CN107024711B (zh) * 2017-04-17 2019-02-26 苏州瑞派宁科技有限公司 一种闪烁脉冲数字化信号的拟合方法
KR101930402B1 (ko) * 2017-04-17 2018-12-18 서울대학교산학협력단 톱니모양 문턱전압을 이용한 시간 기반의 신호 획득 장치 및 방법
CN109171787B (zh) 2018-08-27 2021-02-26 苏州瑞派宁科技有限公司 脉冲信号的采样方法、装置和计算机程序介质

Also Published As

Publication number Publication date
CN109171787A (zh) 2019-01-11
US20210275114A1 (en) 2021-09-09
EP3796048A4 (en) 2021-09-08
EP3796048A1 (en) 2021-03-24
JP2021535357A (ja) 2021-12-16
EP3796048B1 (en) 2023-08-30
HUE063109T2 (hu) 2024-01-28
US11944470B2 (en) 2024-04-02
WO2020042664A1 (zh) 2020-03-05
EP3796048C0 (en) 2023-08-30
JP7265791B2 (ja) 2023-04-27

Similar Documents

Publication Publication Date Title
CN109171787B (zh) 脉冲信号的采样方法、装置和计算机程序介质
CN102073059B (zh) 一种数字化pileup波形处理方法及系统
CN109350098B (zh) 信号的拟合方式的确定方法、重建方法和装置
US9279888B2 (en) System for online/offline singles-pairing with keeping/rejecting multi-coincidences for positron emission tomography
US20130240721A1 (en) Pet calibrations with varying coincidence windows
US20150069255A1 (en) Noise discrimination in photon counting system
EP3427667B1 (en) Time calibration in pet device
US10890671B2 (en) Time-based signal acquisition apparatus and method using sawtooth-shaped threshold voltage
CN110415311B (zh) Pet图像重建方法、系统、可读存储介质和设备
WO2021135337A1 (zh) 信号采样电路、探测装置及成像系统
CN110432922B (zh) 一种提高pet系统时间校准精度的方法
CN106842277B (zh) 一种堆积事件处理方法及装置
CN107767427B (zh) 一种信号波形恢复方法及装置
US4812656A (en) Processing of radioisotope-distribution imaging signals in a scintillation camera apparatus
CN110988961B (zh) 信号处理方法、装置和探测系统
CN111158039B (zh) 信号采样、重建方法及装置
US10663607B2 (en) Apparatuses for processing signals for a plurality of energy regions, and systems and methods for detecting radiation of a plurality of energy regions
CN114660650A (zh) 一种成形时间自适应的数字脉冲滤波方法、装置及电子设备
US20170074996A1 (en) Transmission Of PET-Signals By Means Of Time Division Multiplexing
CN115113512B (zh) 脉冲的采样方法、采样系统、装置及计算机可读存储介质
US20240188918A1 (en) Timer-based amplitude correction method for photon counting computed tomography
US11733405B1 (en) System for timing pickoff of qualified signals
CN117043641A (zh) 脉冲整形器电路
Popescu Image quality evaluation using automatic image scanning and a novel nonparametric free-response data analysis method. Application to PET energy based scatter correction evaluation

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
GR01 Patent grant
GR01 Patent grant