CN102338880A - 核脉冲幅度数字化方法及系统 - Google Patents

核脉冲幅度数字化方法及系统 Download PDF

Info

Publication number
CN102338880A
CN102338880A CN2011101474232A CN201110147423A CN102338880A CN 102338880 A CN102338880 A CN 102338880A CN 2011101474232 A CN2011101474232 A CN 2011101474232A CN 201110147423 A CN201110147423 A CN 201110147423A CN 102338880 A CN102338880 A CN 102338880A
Authority
CN
China
Prior art keywords
signal
tot
threshold
circuit
dynamic 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.)
Granted
Application number
CN2011101474232A
Other languages
English (en)
Other versions
CN102338880B (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.)
University of Science and Technology of China USTC
Original Assignee
University of Science and Technology of China USTC
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 University of Science and Technology of China USTC filed Critical University of Science and Technology of China USTC
Priority to CN2011101474232A priority Critical patent/CN102338880B/zh
Publication of CN102338880A publication Critical patent/CN102338880A/zh
Priority to PCT/CN2012/075885 priority patent/WO2012163239A1/zh
Priority to JP2013517000A priority patent/JP5698352B2/ja
Priority to US13/807,765 priority patent/US9182500B2/en
Priority to EP12793792.8A priority patent/EP2592437B1/en
Application granted granted Critical
Publication of CN102338880B publication Critical patent/CN102338880B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01TMEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/16Measuring radiation intensity
    • G01T1/17Circuit arrangements not adapted to a particular type of detector

Abstract

本发明实施例公开了核脉冲幅度数字化方法及系统,该方法包括:得到本次核事例的到达时刻Td;获得待测的成形信号;从To时刻起,比较待测的成形信号与动态阈值信号在同时刻的取值,当某一时刻后者由小于前者变化为不小于前者时,将该时刻作为本次核事例的过阈点时刻Tot,Td与To相差ΔT;将Td与Tot之间的时间长度转化为数字量。所述动态阈值信号满足:
Figure DDA0000065929360000011
其中,Tot1为某一成形信号与动态阈值信号比较得到的过阈点时刻,Vot1为该成形信号的幅度,Tot2为另一成形信号与动态阈值信号比较得到的过阈点时刻,Vot2为所述另一成形信号的幅度。可见,在一次核事例探测中,待测的成形信号的幅度Vot与(Tot-Td+ΔT)或(Tot-Td-ΔT)成线性比例关系。将Td和Tot之间的时间长度数字化,即可实现对幅度的数字化。

Description

核脉冲幅度数字化方法及系统
技术领域
本发明涉及核信号处理技术领域,更具体地说,涉及核脉冲幅度数字化方法及系统。
背景技术
核事例能量探测器的基本原理都是将入射射线转化为可测量的电信号进行测量的,这些电信号是一系列分离的脉冲信号,其中脉冲信号的幅度值代表着被测射线的能量。现在,需要有一种技术方案来对上述脉冲信号的幅度进行数字化,以便于后续分析和处理。
发明内容
有鉴于此,本发明目的在于提供一种核脉冲幅度数字化方法及系统。
为实现上述目的,本发明提供如下技术方案:
一种核脉冲幅度数字化方法,基于核事例能量探测器、恒比定时电路和成形电路,所述方法包括:
利用恒比定时电路对所述核事例能量探测器输出的待测脉冲信号进行恒比定时处理,输出恒比定时信号,所述恒比定时信号用于标志本次核事例的到达时刻Td;
利用成形电路对所述待测脉冲信号进行成形处理,获得待测的成形信号;
从To时刻起,比较待测的成形信号与动态阈值信号在相同时刻的取值,当某一时刻后者由小于前者变化为不小于前者时,将该时刻作为本次核事例的过阈点时刻Tot,所述Td与To相差预设时长ΔT;
将所述Td与所述Tot之间的时间长度转化为数字量;
所述动态阈值信号满足下述公式:
Tot 2 - To Vot 2 = Tot 1 - To Vot 1 , 其中,所述Tot1为某一成形信号与所述动态阈值信号比较得到的过阈点时刻,Vot1为该成形信号的幅度,Tot2为另一成形信号与所述动态阈值信号比较得到的过阈点时刻,所述Vot2为所述另一成形信号的幅度。
一种核脉冲幅度数字化系统,其特征在于,包括核事例能量探测器、成形电路、恒比定时单元、比较器电路和时间数字转换TDC单元,所述核事例能量探测器的输出端分别与成形电路和恒比定时电路的输入端相连接,所述恒比定时单元包括恒比定时电路,所述时间数字转换TDC单元包括TDC电路;
其中:
所述成形电路用于,对所述核事例能量探测器输出的待测脉冲信号进行成形处理,输出待测的成形信号;
所述恒比定时电路用于,对所述核事例能量探测器输出的待测脉冲信号进行恒比定时处理,输出恒比定时信号,所述恒比定时信号用于标志本次核事例的到达时刻Td;
所述比较器电路用于,从To时刻起,比较待测的成形信号与动态阈值信号在相同时刻的取值,当某一时刻后者由小于前者变化为不小于前者时,输出脉冲信号,所述比较器电路输出的脉冲信号用于标志本次核事例的过阈点时刻Tot,所述Td与To相差预设时长;
所述TDC电路,记录所述Td,记录所述Tot,以及将Td与Tot之间的时间长度转化为数字量;
所述动态阈值信号满足下述公式:
Tot 2 - To Vot 2 = Tot 1 - To Vot 1 , 其中,所述Tot1为某一成形信号与所述动态阈值信号比较得到的过阈点时刻,Vot1为该成形信号的幅度,Tot2为另一成形信号与所述动态阈值信号比较得到的过阈点时刻,所述Vot2为所述另一成形信号的幅度。
可见,在本发明提供的技术方案,由于动态阈值信号满足 Tot 2 - To Vot 2 = Tot 1 - To Vot 1 , 则可知,对于待测脉冲信号代表的核事例探测而言,假设待测的成形信号的幅度为Vot,则同样有 Tot - To Vot = Tot 2 - To Vot 2 = Tot 1 - To Vot 1 , 又知输出时刻Td与To相差预设时长ΔT(可以是To=Td-ΔT,或To=ΔT+Td),而则 Tot - To Vot = Tot 2 - To Vot 2 = Tot 1 - To Vot 1 可进一步变化为 Tot - Td - ΔT Vot = Tot 2 - To Vot 2 = Tot 1 - To Vot 1 或者 Tot - Td + ΔT Vot = Tot 2 - To Vot 2 = Tot 1 - To Vot 1 , 也即,在一次核事例探测中,根据本发明所提供的技术方案,待测的成形信号的幅度Vot与(Tot-Td+ΔT)或(Tot-Td-ΔT)成线性比例关系。由于,Td和Tot均可测,而ΔT为一预设的已知值,将Td和Tot之间的时间长度数字化,即可实现对信号幅度的数字化。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的511KeVγ射线的电流波形;
图2为本发明实施例提供的不同能量核事例的电流波形经积分后的电压波形;
图3为本发明实施例提供的恒比定时电路结构示意图;
图4为本发明实施例提供的核脉冲幅度数字化方法流程图;
图5为本发明实施例提供的核脉冲幅度数字化系统结构示意图;
图6为本发明实施例提供的有源RC积分电路的结构示意图;
图7为本发明实施例提供的有源RC积分电路的输出信号的波形;
图8为本发明实施例提供的核脉冲幅度数字化系统另一结构示意图;
图9a为本发明实施例提供的数字阈值电压值的产生图;
图9b为本发明实施例提供的数字阈值电压值的另一产生图;
图9c为本发明实施例提供的根据图9b所示的D1-D3三条动态阈值曲线而测量到的过阈时间与核事例能量探测器输出信号的幅度值之间的对应关系;
图10为本发明实施例提供的核脉冲幅度数字化系统又一结构示意图;
图11为本发明实施例提供的基于FPGA芯片的动态阈值电压产生单元的结构示意图;
图12为本发明实施例提供的恒比定时甑别电路的结构图;
图13为本发明实施例提供的恒比定时甑别电路的又一结构图;
图14为本发明实施例提供的计时示意图。
具体实施方式
为了引用和清楚起见,下文中使用的技术名词、简写或缩写总结如下:
恒比定时,是取信号上升沿上高度相当于脉冲幅度固定比例的点作为定时瞬间的定时方法。
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
为便于理解,先对核事例能量探测器、成形电路、恒比定时电路进行简单介绍。
本发明中的核事例能量探测器主要用于探测γ射线的能量,核事例能量探测器一般由闪烁晶体,比如LYSO(硅酸钇镥)闪烁晶体、碘化钠单晶,和光电转换控制器(比如光电倍增管-PMT)组合而成,由光电转换器(本实施例用PMT)输出脉冲信号。图1示出了PMT输出的、511KeVγ射线的电流波形。该电流波形前沿的上升时间小于1ns,后沿的下降时间由LYSO的衰减时间决定,在本发明某一实施例中可将其设定为42ns,而电流波形的面积代表着该次核事例被探测到的能量值。
成形电路主要对上述核事例能量探测器输出的电流波形(也即上述脉冲信号)进行积分处理,得到电压波形(或将其称为成形信号,其幅度与入射的γ射线的能量成正比),当然在积分处理后,也可对电压波形再进行低通滤波或放大。参见图2(图2中,横轴表示时间,单位为纳秒,纵轴表示电压,单位为伏),可知,不同能量核事例的电流波形经积分后的电压波形的差异仅是电压波形的峰值(即幅度值),它们的达峰时间Tf、上升时间和下降时间均相同,其中,上述成形信号的上升时间主要取决于核事例能量探测器的发光衰减特性,下降时间取决于成形电路的时间常数。也即,在选定核事例能量探测器和成形电路后,成形信号的上升时间、达峰时间以及下降时间即为固定值。本发明后续记载中,采用恒比定时方法确定信号到达的时刻T1,由图2可知,对于不同幅度的信号,取信号上升沿上高度相当于脉冲幅度固定比例的点所对应的时刻作为定时时刻T1(即后续的Td),T1和信号的达峰时刻Tf之间的时间差是一个固定的值Δt1。
恒比定时电路的一种结构如图3所示:核事例能量探测器输出的脉冲信号分成两路,一路经延时电路31延迟,一路经衰减电路32衰减,经比较电路B1进行比较,比较电路B1输出的恒比定时信号的上升沿所发生时刻在上述电流波形上的取值与电流波形的最大幅度值为固定比例,该比例的大小由延迟电路31的延迟量、衰减电路32的衰减系数和电流波形的上升时间共同决定,可通过改变延迟量来对上述比例进行调置,达到调节恒比定时信号的上升沿的发生时刻的目的。
基于上述介绍,本发明实施例提供了核脉冲幅度数字化方法,以实现核事例能量探测器所输出脉冲信号幅度的数字化。参见图4,该方法至少包括如下步骤:
S1、利用恒比定时电路对上述核事例能量探测器输出的脉冲信号进行恒比定时处理,输出恒比定时信号。
S2、利用成形电路对上述核事例能量探测器输出的待测脉冲信号进行成形处理,获得待测的成形信号。
S3、从To时刻起,比较待测的成形信号与动态阈值信号在相同时刻的取值,当某一时刻后者由小于前者变化为不小于前者时,将该时刻作为过阈点时刻Tot,上述Td与To相差预设时长ΔT;
S4、将上述Td与Tot之间的时间长度转化为数字量。
上述动态阈值信号是根据探测器特性和所用成形电路参数预先计算好并存储起来,在被测信号到来时的T0时刻由动态阈值产生电路产生。
对于任意两个幅度为Vot1和Vot2的成形信号,上述动态阈值信号要满足下述公式:
Tot 2 - To Vot 2 = Tot 1 - To Vot 1 (公式一)
其中,Tot1为某一成形信号与上述动态阈值信号比较得到的过阈点时刻(也即两波形相交叉的时刻点),Tot2为另一成形信号与该动态阈值信号比较得到的过阈点时刻(也即两波形相交叉的时刻点)。
在本发明中,将上述Td与Tot之间的时间长度称为过阈时间。
本发明同时还公开了一种能实施上述方法的核脉冲幅度数字化系统。图5示出了该系统的一种结构,包括:
核事例能量探测器1、成形电路2、恒比定时单元3、比较器电路4、时间数字转换(TDC)单元5,核事例能量探测器1的输出端分别与成形电路2和恒比定时电路3的输入端相连接,恒比定时单元包括上述恒比定时电路,而时间数字转换TDC单元5包括TDC电路。
其中:
成形电路2用于,对核事例能量探测器1输出的待测脉冲信号进行成形处理,输出待测的成形信号;
恒比定时电路用于,对核事例能量探测器1输出的待测脉冲信号进行恒比定时处理,输出恒比定时信号;
比较器电路4用于,从To时刻起,比较待测的成形信号与上述动态阈值信号在相同时刻的取值,当某一时刻后者由小于前者变化为不小于前者时,输出脉冲信号;
TDC电路,记录接收到所述恒比定时信号的到达时刻Td,记录接收到所述比较器电路输出的脉冲信号的时刻作为过阈点时刻Tot,以及将Td与Tot之间的时间长度转化为数字量。
由于动态阈值信号满足 Tot 2 - To Vot 2 = Tot 1 - To Vot 1 , 则可知,对于待测脉冲信号代表的核事例探测而言,假设待测的成形信号的幅度为Vot,则同样有 Tot - To Vot = Tot 2 - To Vot 2 = Tot 1 - To Vot 1 , 又知输出时刻Td与To相差预设时长ΔT(可以是To=Td-ΔT,或To=ΔT+Td),而则 Tot - To Vot = Tot 2 - To Vot 2 = Tot 1 - To Vot 1 可进一步变化为 Tot - Td - ΔT Vot = Tot 2 - To Vot 2 = Tot 1 - To Vot 1 或者 Tot - Td + ΔT Vot = Tot 2 - To Vot 2 = Tot 1 - To Vot 1 , 也即,在一次核事例探测中,根据本发明所提供的技术方案,待测的成形信号的幅度Vot与(Tot-Td+ΔT)或(Tot-Td-ΔT)成线性比例关系。由于,Td和Tot均可测,而ΔT为一预设的已知值,将Td和Tot之间的时间长度数字化,即可实现对信号幅度的数字化。
在本发明其他实施例中,上述成形处理具体可为有源积分成形处理。上述成形电路具体可为有源积分电路(例如有源RC积分电路)。图6则示出了有源RC积分电路的一种结构。该电路将核事例能量探测器1输出的电流信号积分,输出的电压波形可参见图7。
在本发明其他实施例中,上述有源积分电路中的电阻R=5KΩ,电容C=18pf,输出电压波形的达峰时间为60ns,峰值后恢复到基线的时间大于300ns。不同能量核事例积分后电压波形的差异仅是电压波形的峰值,即电压波形幅度,它们的达峰时间和下降时间均相同。
在本发明其他实施例中,参见图8,在比较成形信号与动态阈值信号在相同时刻的取值之前,还包括如下步骤:
S5、生成动态阈值信号。
动态阀值信号只要在To时刻以前生成,并在To时刻送入上述比较器电路作为输入之一即可。
在本发明其他实施例中,上述动态阈值信号可根据预先存储的一系列数字阈值电压值生成。由于根据上述公式一可以预先计算出任一幅度的成形信号所对应的阈值的取值,这些取值可以组成阈值序列(即上述一系列数字阈值电压值)加以存储。
下面将对如何产生上述一系列数字阈值电压值进行介绍。
方式一:
请参见图9a,上述一系列数字阈值电压值的产生流程如下:
一,确定时间间隔Δt和上述一系列数字阈值电压值的个数N,N为不小于1的正整数,在图9a中,N等于5;
二,选取N个成形信号S1-SN,其中,SN的幅度为A,S1的幅度为A/N,所述S2至SN-1的幅度值以A/N递增,例如在图9a中,S5-S1的幅度分别为2.5V,2.0V,1.5V,1.0V和0.5V;
前已述及,在选定核事例能量探测器和成形电路后,成形信号的上升时间、达峰时间以及下降时间即为固定值。因此,可以很容易地使用信号发生器或软件仿真生成N个成形信号。
三,设定成形信号的某一时刻为时间原点Ts1,TS1也即上述提及的To。在本实施例中选择成形信号的达峰时刻为时间原点;
四,将成形信号SM在Ts1+MΔt时刻对应的电压取值作为数字阈值电压值进行存储,所述M为正整数,并且1≤M≤N,即可得到5个数字阈值电压值,上述5个数字阈值电压值可构成离散的动态阈值曲线。
本领域技术人员完全可根据精度及其他需要对时间间隔Δt和N的取值进行灵活设计以获取不同的动态阀值曲线(如图9a中的D1-D3曲线),但需要保证Ts1+NΔt不超出上述TDC电路所能测量的最大时间范围,此不再赘述。
方式二:
请参见图9b,上述一系列数字阈值电压值的产生流程如下:
一,确定最大时间测量值Tmax,图9b示出三条动态阈值曲线D1-D3对应的最大时间测量值分别为156ns、312.5ns和625ns;
二,设定成形信号的某一时刻为时间原点Ts1,在本实施例中选择成形信号的达峰时刻为时间原点;
三,确定时间间隔Δt;
四,根据公式(Tmax-Ts1)/Δt可确定数字阈值电压值的个数N,N为不小于1的正整数。在本实施例中,N等于5,当然,本领域技术人员完全可根据精度及其他需要对时间间隔Δt和N的取值进行灵活设计,在此不再赘述;
五,选取N个成形信号S1-SN,所述SN的幅度为A,所述S1的幅度为A/N,所述S2至SN-1的幅度值以A/N递增;
六,将成形信号SM在Ts1+MΔt时刻(Ts1+NΔt=Tmax)对应的电压取值作为数字阈值电压值进行存储,所述M为正整数,并且1≤M≤N,即可得到N个数字阈值电压值,上述N个数字阈值电压值可构成离散的动态阈值曲线。
为了提高测量精度,Δt可选取得小一点,在本发明其他实施例中,可采用12bitDAC,其变换时钟为160MHz,时钟周期是6.25ns。仍以上述156ns、312.5ns和625ns为例,其分别被6.25ns等间隔划分为25、50和100个数据计算点(为方便观看,图9b仅示出了5个计算机点,与之对应,也需要选取25、50和100个幅度不同的成形信号。
为证明本发明所提供技术方案中,过阈时间与成形电路输出信号的幅度值呈良好的线性关系。可使用任意波形信号发生器产生和上述核事例能量探测器输出波形一样的信号波形,并控制该波形的幅度以得到不同的测量值。图9c示出了根据图9b所示的D1-D3三条动态阈值曲线而测量到的过阈时间与核事例能量探测器输出信号的幅度值之间的对应关系。可观察到,使用本发明提供的技术方案,过阈时间与成形电路输出信号的幅度值确实呈良好的线性关系。
上述方式一和方式二均可通过PC机完成。
在本发明其他实施例中,参见图10,上述系统还可包括生成上述动态阈值信号的动态阈值产生单元6。
作为本发明的一个实施例,图11是一种基于FPGA(现场可编程门阵列)芯片的动态阈值电压产生单元的结构示意图。在图11中,由PC机计算好的一系列数字阈值电压值存储在FPAG的静态存储器(RAM)中,写入读出逻辑块按照预设的时间间隔(6.25ns)将上述数字阈值电压值逐个读出,经过数字模拟变换电路(DAC)和低通滤波电路后,一个连续时间和连续幅度的动态阈值电压信号就产生了。为了提高阈值电压信号的精确度,可以选择较高分辨率的DAC芯片、较高的变换时钟速率。
在本发明的一个实施例中,参见图12,可使用预甑别电路配合恒比定时电路共同完成一个有效事例的恒比定时输出。预甑别电路是一种简单的前沿定时电路,它将输入信号和固定的前沿定时阈值(即一直流电)进行比较,只有当输入信号具有超过该前沿定时阈值的取值时,才会有预甑别信号输出。引入预甑别电路的目的是避免噪声引起恒比定时信号的输出,预甑别电路的前沿定时阈值由一片数模转化(DAC)芯片的输出来给定。需要注意的是,预甑别信号的输出先于恒比定时信号的输出。在图12中,恒比定时电路的输出只有在预甑别电路有输出时才能通过与门输出。
在本发明的一个实施例中,恒比定时电路和预甑别电路的输出均与FPGA芯片的输入端相连(见图13),在FPGA内部实现必要的逻辑判选(该逻辑判选与上述与门的功能相类似),只有在收到预甑别信号时,恒比定时电路的输出才被认为是一次核事例信号的到来,才记录恒比定时信号的到达时刻。此时,上述TDC单元5还可包括指示单元,指示单元可以在接收到上述预甑别信号时,指示上述TDC电路工作。
目前常用的时间测量技术的基本原理是用计数器计算从开始时刻到结束时刻之间的时钟个数(在本发明中,开始时刻即为恒比定时信号的到达时刻Td,结束时刻即为过阈点时刻Tot),提高时钟频率是提高测量精度的一种方法。为进一步提高测量精度,在使用FPGA实现TDC电路的功能时,可采用四相位时钟将时间测量误差减小为一个时钟周期的四分之一。图14中的恒比定时信号到来时(图14中用定时信号代表恒比定时信号和下述比较器电路输出的脉冲信号),9-bit计数器的计算值连同四相位的状态一并被锁存记录(即记录上述到达时刻),表明起始时刻的计数值和相位,等到比较器电路输出的脉冲信号到来时,同样记录下当时的计数器值和四相位时钟的状态(即记录上述过阈点时刻Tot),由两个时刻的计数值和四相位时钟的状态,计算机可以根据将过阈时间减去ΔT,而最终实现幅度的数字化。以上述160MHz时钟周期为例,采用四相位时钟得到的时间测量精度是160MHz时钟周期的四分之一,即1.5625ns。
本发明在具体实施过程中,影响最终测量精度的主要因素有恒比定时精度、TDC测量精度、动态阈值产生电路的精度等,还取决于TDC测量范围和被测信号幅度范围的对应关系,将核事例能量探测器的特性与这些因素作综合考虑,才能够实现高精度的核脉冲幅度数字化测量。
综上可知,本发明在不明显增加核事例能量探测器的前端读出电路的基础上,利用时间数字转换技术,将核脉冲信号的电压幅度值,转换为过阈时间量,从而实现能量信号的数字化。
在本发明某些实施例中,还利用了FPGA芯片技术或半导体集成电路技术,利用上述技术最突出优点是可以实现大量通道数的并行读出,可以实现多通道核事例能量信号测量的高集成度。因此,其尤其适用于核成像探测器(核成像探测器要求同时读出几十乃至几百或几万个通道所探测到的核事例能量信息),可以满足核成像探测器的实验室研究和实际应用的需要。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的装置而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下,在其它实施例中实现。因此,本发明将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。

Claims (10)

1.一种核脉冲幅度数字化方法,其特征在于,基于核事例能量探测器、恒比定时电路和成形电路,所述方法包括:
利用恒比定时电路对所述核事例能量探测器输出的待测脉冲信号进行恒比定时处理,输出恒比定时信号,所述恒比定时信号用于标志本次核事例的到达时刻Td;
利用成形电路对所述待测脉冲信号进行成形处理,获得待测的成形信号;
从To时刻起,比较待测的成形信号与动态阈值信号在相同时刻的取值,当某一时刻后者由小于前者变化为不小于前者时,将该时刻作为本次核事例的过阈点时刻Tot,所述Td与To相差预设时长ΔT;
将所述Td与所述Tot之间的时间长度转化为数字量;
所述动态阈值信号满足下述公式:
Tot 2 - To Vot 2 = Tot 1 - To Vot 1 , 其中,所述Tot1为某一成形信号与所述动态阈值信号比较得到的过阈点时刻,Vot1为该成形信号的幅度,Tot2为另一成形信号与所述动态阈值信号比较得到的过阈点时刻,所述Vot2为所述另一成形信号的幅度。
2.如权利要求1所述的方法,其特征在于,在比较待测的成形信号与动态阈值信号在相同时刻的取值之前,还包括:生成动态阈值信号。
3.如权利要求2所述的方法,其特征在于,所述生成所述动态阈值信号的具体实施步骤包括:
根据预先存储的一系列数字阈值电压值生成动态阈值信号。
4.如权利要求3所述的方法,其特征在于,
在所述根据预先存储一系列数字阈值电压值生成动态阈值信号的步骤之前,还包括:
确定时间间隔Δt和所述一系列数字阈值电压值的个数N,N为不小于1的正整数;
选取N个成形信号S1-SN,所述SN的幅度为A,所述S1的幅度为A/N,所述S2至SN-1的幅度值以A/N递增;
设定成形信号的某一时刻为时间原点Ts1;
将成形信号SM在Ts1+MΔt时刻对应的电压取值作为数字阈值电压值进行存储,所述M为正整数,并且1≤M≤N。
5.如权利要求4所述的方法,其特征在于,所述根据预先存储一系列数字阈值电压值生成动态阈值信号的具体实施方式包括:
按所述预设时间间隔依次读出预先存储的一系列数字阈值电压值,对其进行数模转换,转换为模拟信号;
对所述模拟信号进行低通滤波和放大,得到所述动态阈值信号。
6.如权利要求1-5任一项所述的方法,其特征在于,将所述Td与Tot之间的时间长度转化为数字量的具体实施方式包括:
当所述核事例能量探测器输出的待测脉冲信号具有不低于固定阈值的取值时,将所述Td与Tot之间的时间长度转化为数字量。
7.如权利要求1-5任一项所述的方法,其特征在于,所述成形处理具体为有源积分成形处理。
8.一种核脉冲幅度数字化系统,其特征在于,包括核事例能量探测器、成形电路、恒比定时单元、比较器电路和时间数字转换TDC单元,所述核事例能量探测器的输出端分别与成形电路和恒比定时电路的输入端相连接,所述恒比定时单元包括恒比定时电路,所述时间数字转换TDC单元包括TDC电路;
其中:
所述成形电路用于,对所述核事例能量探测器输出的待测脉冲信号进行成形处理,输出待测的成形信号;
所述恒比定时电路用于,对所述核事例能量探测器输出的待测脉冲信号进行恒比定时处理,输出恒比定时信号,所述恒比定时信号用于标志本次核事例的到达时刻Td;
所述比较器电路用于,从To时刻起,比较待测的成形信号与动态阈值信号在相同时刻的取值,当某一时刻后者由小于前者变化为不小于前者时,输出脉冲信号,所述比较器电路输出的脉冲信号用于标志本次核事例的过阈点时刻Tot,所述Td与To相差预设时长;
所述TDC电路,记录所述Td,记录所述Tot,以及将Td与Tot之间的时间长度转化为数字量;
所述动态阈值信号满足下述公式:
Tot 2 - To Vot 2 = Tot 1 - To Vot 1 , 其中,所述Tot1为某一成形信号与所述动态阈值信号比较得到的过阈点时刻,Vot1为该成形信号的幅度,Tot2为另一成形信号与所述动态阈值信号比较得到的过阈点时刻,所述Vot2为所述另一成形信号的幅度。
9.如权利要求8所述的系统,其特征在于,还包括生成动态阈值信号的动态阈值产生单元。
10.如权利要求8所述的系统,其特征在于,所述恒比定时单元还包括预甑别电路,所述TDC单元还包括指示单元,其中:
预甑别电路用于当所述核事例能量探测器输出的待测脉冲信号具有高于固定的前沿定时阈值的取值时,输出预甑别信号;
所述指示单元用于在接收到所述预甑别信号时,指示所述TDC电路工作。
CN2011101474232A 2011-06-02 2011-06-02 核脉冲幅度数字化方法及系统 Expired - Fee Related CN102338880B (zh)

Priority Applications (5)

Application Number Priority Date Filing Date Title
CN2011101474232A CN102338880B (zh) 2011-06-02 2011-06-02 核脉冲幅度数字化方法及系统
PCT/CN2012/075885 WO2012163239A1 (zh) 2011-06-02 2012-05-22 核辐射脉冲幅度的数字化方法及系统
JP2013517000A JP5698352B2 (ja) 2011-06-02 2012-05-22 核放射パルス幅のディジタル化方法及びシステム
US13/807,765 US9182500B2 (en) 2011-06-02 2012-05-22 Method and system for amplitude digitization of nuclear radiation pulses
EP12793792.8A EP2592437B1 (en) 2011-06-02 2012-05-22 Method and system for digitalizing nuclear radiation pulse amplitude

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2011101474232A CN102338880B (zh) 2011-06-02 2011-06-02 核脉冲幅度数字化方法及系统

Publications (2)

Publication Number Publication Date
CN102338880A true CN102338880A (zh) 2012-02-01
CN102338880B CN102338880B (zh) 2012-12-26

Family

ID=45514707

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2011101474232A Expired - Fee Related CN102338880B (zh) 2011-06-02 2011-06-02 核脉冲幅度数字化方法及系统

Country Status (5)

Country Link
US (1) US9182500B2 (zh)
EP (1) EP2592437B1 (zh)
JP (1) JP5698352B2 (zh)
CN (1) CN102338880B (zh)
WO (1) WO2012163239A1 (zh)

Cited By (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012163239A1 (zh) * 2011-06-02 2012-12-06 中国科学技术大学 核辐射脉冲幅度的数字化方法及系统
CN102830420A (zh) * 2012-08-31 2012-12-19 中国原子能科学研究院 核辐射探测器工作电压标定方法及系统
CN102841368A (zh) * 2012-08-31 2012-12-26 中国原子能科学研究院 气体核辐射探测器收集的电荷数与外加电压的关系曲线测量方法及系统
CN102841366A (zh) * 2012-08-31 2012-12-26 中国原子能科学研究院 脉冲幅度甄别器甄别阈测定方法及系统
CN103675891A (zh) * 2013-12-11 2014-03-26 成都理工大学 基于双线性变换法的数字核脉冲高斯成形方法
CN104155680A (zh) * 2014-08-28 2014-11-19 东华理工大学 基于ads1605的多道能谱分析仪数据采集装置
CN104242877A (zh) * 2014-07-16 2014-12-24 成都理工大学 基于模拟cr-rc电路的核脉冲信号数字高斯成形方法
CN104777508A (zh) * 2015-03-30 2015-07-15 东南大学 一种基于模型库的数字脉冲重叠峰分离算法
CN105141294A (zh) * 2015-08-17 2015-12-09 西北核技术研究所 双预置甄别门控型恒比定时甄别器及数字恒比定时甄别方法
CN102916683B (zh) * 2012-10-18 2016-09-14 成都理工大学 一种可调参数核脉冲模拟方法
CN106160706A (zh) * 2015-04-24 2016-11-23 南京理工大学 一种双通道时刻鉴别电路
CN110133708A (zh) * 2019-04-24 2019-08-16 中国科学技术大学 一种核脉冲信号的测量装置及测量方法
CN110531404A (zh) * 2019-06-03 2019-12-03 中国科学技术大学 核脉冲电荷时间转换方法与系统
CN110568468A (zh) * 2019-08-27 2019-12-13 福州智元仪器设备有限公司 一种辐射脉冲计数突变算法
CN110954938A (zh) * 2019-12-11 2020-04-03 湖北锐世数字医学影像科技有限公司 比较器的阈值校正方法、装置和采样电路
CN111308223A (zh) * 2019-12-12 2020-06-19 江苏科技大学 基于传输响应的脉冲信号幅度测量方法
CN112255662A (zh) * 2020-10-14 2021-01-22 南昌华亮光电有限责任公司 基于智能时间采样的伽马相机及探测方法
CN112327347A (zh) * 2020-10-29 2021-02-05 中广核久源(成都)科技有限公司 曲率可调的数字化核脉冲成形系统
CN113126140A (zh) * 2019-12-30 2021-07-16 中核控制系统工程有限公司 一种高速数字n/γ波形实时甄别的系统及方法

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8981756B2 (en) * 2012-04-11 2015-03-17 Rockwell Automation Technologies, Inc. Adaptive threshold voltage for frequency input modules
JP6214031B2 (ja) * 2013-07-19 2017-10-18 国立研究開発法人理化学研究所 放射線検出器のための信号データ処理方法、信号データ処理装置、および放射線検出システム
FR3058230B1 (fr) * 2016-10-27 2019-03-15 Detection Technology Sas Dispositif de spectrometrie
JP7057630B2 (ja) 2017-06-23 2022-04-20 浜松ホトニクス株式会社 放射線検出装置
CN109959962B (zh) * 2017-12-14 2022-07-26 中国核动力研究设计院 基于脉冲型中子探测器信号特性的核信号发生器
RU2705933C1 (ru) * 2018-12-26 2019-11-12 Федеральное государственное автономное образовательное учреждение высшего образования "Национальный исследовательский Томский политехнический университет" Двухканальный сцинтилляционный счетчик ионизирующего излучения
CN110347965B (zh) * 2019-07-17 2023-02-03 四川新先达测控技术有限公司 核脉冲信号处理方法及装置
RU198513U1 (ru) * 2019-12-30 2020-07-14 Федеральное государственное автономное образовательное учреждение высшего образования "Новосибирский национальный исследовательский государственный университет" (Новосибирский государственный университет, НГУ) Стенд для проведения испытаний сцинтилляционных датчиков
CN111487663B (zh) * 2020-03-16 2022-06-03 中国石油天然气集团有限公司 一种核测井仪器能量电流谱获取方法
WO2022004639A1 (ja) * 2020-07-03 2022-01-06 国立大学法人東京大学 検出回路
CN113359181B (zh) * 2021-07-01 2022-11-04 成都理工大学 一种新型平头尖顶脉冲成形系统及方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6085382A (ja) * 1983-10-17 1985-05-14 Hitachi Ltd 広レンジ信号処理方法
CN1732847A (zh) * 2004-03-24 2006-02-15 株式会社东芝 磁共振成像装置和磁共振成像数据处理方法

Family Cites Families (35)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4217497A (en) * 1978-03-14 1980-08-12 The United States Of America As Represented By The Department Of Health, Education And Welfare Portable instrument for measuring neutron energy spectra and neutron dose in a mixed n-γ field
DE4223773C2 (de) * 1992-07-18 1997-07-10 Berthold Lab Prof Dr Verfahren zur Unterscheidung und gleichzeitigen oder getrennten Messung von Einzel- und Mehrelektronenereignissen in einem optoelektronischen Detektor
JPH10197643A (ja) * 1997-01-06 1998-07-31 Toshiba Corp 放射線計測装置及び放射線計測システム
JP3103047B2 (ja) * 1997-09-01 2000-10-23 シュランベルジェ、ホールディング、リミテッド 核分光システムにおけるパルス波形調整及び弁別方法
US6252232B1 (en) * 1998-06-29 2001-06-26 General Electric Company Digital integrator
US6160259A (en) * 1998-10-23 2000-12-12 Adac Laboratories Channel-specific control of pulse integration in a gamma camera system in response to pulse pile-up
US6232604B1 (en) * 1998-11-06 2001-05-15 General Electric Company Analog time adjustment for coincidence detection electronics
JP4064009B2 (ja) * 1999-07-30 2008-03-19 株式会社東芝 線種弁別型放射線検出装置
US6459086B1 (en) * 2000-11-24 2002-10-01 Koninklijke Philips Electronics, N.V. Digital peak detector for radiation detection systems
JP4160275B2 (ja) * 2001-05-28 2008-10-01 浜松ホトニクス株式会社 エネルギー測定方法及び測定装置
US6781426B2 (en) * 2001-07-26 2004-08-24 The Regents Of The University Of California Nanosecond monolithic CMOS readout cell
US7818047B2 (en) * 2001-11-09 2010-10-19 Nova R&D, Inc. X-ray and gamma ray detector readout system
US6664543B2 (en) * 2002-01-28 2003-12-16 Cti Pet Systems, Inc. Continuous sampling and digital integration for PET scintillation
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
US6822506B2 (en) * 2003-04-17 2004-11-23 Concorde Microsystems, Inc. Amplitude and rise-time sensitive timing-shaping filters with built-in pulse-tail cancellation for high count-rate operation
US7126126B2 (en) * 2003-10-16 2006-10-24 Brookhaven Science Associates, Llc Compact conscious animal positron emission tomography scanner
US7045802B2 (en) * 2003-11-26 2006-05-16 General Electric Company Method and apparatus for coincidence imaging digital triggering
DE102005015018B4 (de) * 2004-04-02 2009-01-02 National Health Research Institutes, Zhunan Vorrichtung und Verfahren zum Digitalisieren von Gammastrahlen-Energie und zum Charakterisieren von Spitzen-Zeit und Zerfallszeit-Konstanten ohne die Verwendung von ADCs
JP4611106B2 (ja) * 2005-03-11 2011-01-12 住友重機械工業株式会社 放射線検出回路及び放射線検査装置
JP4545628B2 (ja) * 2005-04-15 2010-09-15 住友重機械工業株式会社 放射線検出回路、放射線検出器および放射線検査装置
US8598536B2 (en) * 2006-10-04 2013-12-03 CERN—European Organization for Nuclear Research Apparatus and method for medical imaging
US7800439B2 (en) * 2006-10-27 2010-09-21 Ut-Battelle, Llc High speed preamplifier circuit, detection electronics, and radiation detection systems therefrom
CN101034120A (zh) * 2007-03-09 2007-09-12 中国科学院上海光学精密机械研究所 脉冲形状的测量装置及测量方法
US7605375B2 (en) * 2007-04-26 2009-10-20 Oy Ajat Ltd. Multi-functional radiation/photon identifying and processing application specific integrated circuit and device
CA2732756C (en) * 2007-08-03 2017-06-13 Pulsetor, Llc Pileup rejection in an energy-dispersive radiation spectrometry system
US7683334B2 (en) * 2007-08-07 2010-03-23 The State Of Oregon Acting By And Through The State Board Of Higher Education On Behalf Of Oregon State University Simultaneous beta and gamma spectroscopy
US8340377B2 (en) * 2007-09-17 2012-12-25 Siemens Medical Solutions Usa, Inc. Method for energy calculation and pileup determination for continuously sampled nuclear pulse processing
JP2009079969A (ja) * 2007-09-26 2009-04-16 Toshiba Corp 放射線スペクトル計測システム
EA016862B1 (ru) * 2007-11-02 2012-08-30 Юниверсити Оф Вашингтон Способ и устройство для оценки начального момента времени электронного импульса, предназначенного для использования в позитронно-эмиссионной томографии
US8258480B2 (en) * 2008-03-03 2012-09-04 The Board Of Trustees Of The Leland Stanford Junior University High energy photon detection using pulse width modulation
US8384037B2 (en) * 2008-09-11 2013-02-26 Siemens Medical Solutions Usa, Inc. Use of crystal location in nuclear imaging apparatus to minimize timing degradation in a photodetector array
US8492705B2 (en) * 2009-10-23 2013-07-23 Saint-Gobain Ceramics & Plastics, Inc. Radiation detector and method of using a radiation detector
CN101852861A (zh) * 2010-04-09 2010-10-06 罗康林 基于峰值侦测的多道脉冲幅度分析仪
CN102338880B (zh) * 2011-06-02 2012-12-26 中国科学技术大学 核脉冲幅度数字化方法及系统
US9223031B2 (en) * 2012-01-04 2015-12-29 Kabushiki Kaisha Toshiba Method for timing-pick-off of undersampled pulses from radiation detectors using a photosensor-based radiation detection system applying a prototype waveform and a weighting function

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6085382A (ja) * 1983-10-17 1985-05-14 Hitachi Ltd 広レンジ信号処理方法
CN1732847A (zh) * 2004-03-24 2006-02-15 株式会社东芝 磁共振成像装置和磁共振成像数据处理方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
薛?H栋等,1: "恒比定时单道分析器", 《核技术》 *

Cited By (29)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012163239A1 (zh) * 2011-06-02 2012-12-06 中国科学技术大学 核辐射脉冲幅度的数字化方法及系统
CN102841366B (zh) * 2012-08-31 2014-11-19 中国原子能科学研究院 脉冲幅度甄别器甄别阈测定方法及系统
CN102830420A (zh) * 2012-08-31 2012-12-19 中国原子能科学研究院 核辐射探测器工作电压标定方法及系统
CN102841368A (zh) * 2012-08-31 2012-12-26 中国原子能科学研究院 气体核辐射探测器收集的电荷数与外加电压的关系曲线测量方法及系统
CN102841366A (zh) * 2012-08-31 2012-12-26 中国原子能科学研究院 脉冲幅度甄别器甄别阈测定方法及系统
CN102830420B (zh) * 2012-08-31 2014-11-26 中国原子能科学研究院 核辐射探测器工作电压标定方法及系统
CN102841368B (zh) * 2012-08-31 2014-11-19 中国原子能科学研究院 气体核辐射探测器收集的电荷数与外加电压的关系曲线测量方法及系统
CN102916683B (zh) * 2012-10-18 2016-09-14 成都理工大学 一种可调参数核脉冲模拟方法
CN103675891B (zh) * 2013-12-11 2016-03-02 成都理工大学 基于双线性变换法的数字核脉冲高斯成形方法
CN103675891A (zh) * 2013-12-11 2014-03-26 成都理工大学 基于双线性变换法的数字核脉冲高斯成形方法
CN104242877A (zh) * 2014-07-16 2014-12-24 成都理工大学 基于模拟cr-rc电路的核脉冲信号数字高斯成形方法
CN104242877B (zh) * 2014-07-16 2017-04-05 成都理工大学 基于模拟cr‑rc电路的核脉冲信号数字高斯成形方法
CN104155680A (zh) * 2014-08-28 2014-11-19 东华理工大学 基于ads1605的多道能谱分析仪数据采集装置
CN104155680B (zh) * 2014-08-28 2016-05-04 东华理工大学 基于ads1605的多道能谱分析仪数据采集装置
CN104777508A (zh) * 2015-03-30 2015-07-15 东南大学 一种基于模型库的数字脉冲重叠峰分离算法
CN106160706A (zh) * 2015-04-24 2016-11-23 南京理工大学 一种双通道时刻鉴别电路
CN105141294A (zh) * 2015-08-17 2015-12-09 西北核技术研究所 双预置甄别门控型恒比定时甄别器及数字恒比定时甄别方法
CN105141294B (zh) * 2015-08-17 2018-03-30 西北核技术研究所 双预置甄别门控型恒比定时甄别器及数字恒比定时甄别方法
CN110133708A (zh) * 2019-04-24 2019-08-16 中国科学技术大学 一种核脉冲信号的测量装置及测量方法
CN110531404A (zh) * 2019-06-03 2019-12-03 中国科学技术大学 核脉冲电荷时间转换方法与系统
CN110568468A (zh) * 2019-08-27 2019-12-13 福州智元仪器设备有限公司 一种辐射脉冲计数突变算法
CN110568468B (zh) * 2019-08-27 2023-01-10 福州智元仪器设备有限公司 一种辐射脉冲计数突变算法
CN110954938A (zh) * 2019-12-11 2020-04-03 湖北锐世数字医学影像科技有限公司 比较器的阈值校正方法、装置和采样电路
CN110954938B (zh) * 2019-12-11 2021-08-31 湖北锐世数字医学影像科技有限公司 比较器的阈值校正方法、装置和采样电路
CN111308223A (zh) * 2019-12-12 2020-06-19 江苏科技大学 基于传输响应的脉冲信号幅度测量方法
CN111308223B (zh) * 2019-12-12 2022-05-13 江苏科技大学 基于传输响应的脉冲信号幅度测量方法
CN113126140A (zh) * 2019-12-30 2021-07-16 中核控制系统工程有限公司 一种高速数字n/γ波形实时甄别的系统及方法
CN112255662A (zh) * 2020-10-14 2021-01-22 南昌华亮光电有限责任公司 基于智能时间采样的伽马相机及探测方法
CN112327347A (zh) * 2020-10-29 2021-02-05 中广核久源(成都)科技有限公司 曲率可调的数字化核脉冲成形系统

Also Published As

Publication number Publication date
US9182500B2 (en) 2015-11-10
CN102338880B (zh) 2012-12-26
EP2592437B1 (en) 2019-07-10
US20130161528A1 (en) 2013-06-27
WO2012163239A1 (zh) 2012-12-06
EP2592437A1 (en) 2013-05-15
JP5698352B2 (ja) 2015-04-08
EP2592437A4 (en) 2015-02-11
JP2013534629A (ja) 2013-09-05

Similar Documents

Publication Publication Date Title
CN102338880B (zh) 核脉冲幅度数字化方法及系统
CN104502684B (zh) 一种全数字化峰值到达时刻鉴别方法
CN101561507B (zh) 射线探测器的射线能量检测方法
US9344102B2 (en) Threshold correction method for multi-voltage threshold sampling digitization device
US9910167B2 (en) Method for digitalizing scintillation pulse
CN104360376A (zh) 具有放射源核素识别功能的伽马相机及核素识别方法
CN105022082B (zh) 光子测量前端电路
CN103913749A (zh) 一种基于激光脉冲飞行时间测量的测距方法
CN102323959B (zh) 用于时间分辨光子计数成像的采集卡
CN106227026B (zh) 一种双延迟内插法的时间间隔计数器
Ruegg et al. Drift‐Free Mössbauer Spectrometer
CN105301627B (zh) 一种能谱分析方法、能谱分析系统及伽马射线探测系统
CN109581461A (zh) 核脉冲能量测量方法及系统
CN104251986A (zh) 一种全波形数字化检测装置
CN104121956B (zh) 一种时差式超声波流量计时差测量方法
CN102307046A (zh) 一种时间分辨光子计数成像系统及方法
CN103199870A (zh) 一种触发点快速定位装置
CN110703583A (zh) 基于soc的多通道高精度大量程时间数字转换器
US5371362A (en) Nuclear detection process with base potential correction and correspnding apparatus (particularly a gamma-camera)
CN106527099B (zh) 一种时间数字转换器及其时间测量电路与方法
CN111413725B (zh) 一种利用虚拟仪器技术实现γ-γ数字符合测量的系统及方法
CN202334490U (zh) 用于时间分辨光子计数成像的采集卡
Kong et al. An FPGA-based fast linear discharge readout scheme enabling simultaneous time and energy measurements for TOF-PET detectors
CN105720950A (zh) 一种核脉冲信号的产生方法及装置
CN104977556B (zh) 基于死区时间测量的平均波形捕获率测试方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20121226

Termination date: 20190602

CF01 Termination of patent right due to non-payment of annual fee