CN103961126A - 一种多阈值采样数字化器件的阈值校正方法 - Google Patents
一种多阈值采样数字化器件的阈值校正方法 Download PDFInfo
- Publication number
- CN103961126A CN103961126A CN201310045681.9A CN201310045681A CN103961126A CN 103961126 A CN103961126 A CN 103961126A CN 201310045681 A CN201310045681 A CN 201310045681A CN 103961126 A CN103961126 A CN 103961126A
- Authority
- CN
- China
- Prior art keywords
- threshold value
- comparator
- reference voltage
- dot
- many
- 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
- 238000005070 sampling Methods 0.000 title claims abstract description 101
- 238000000034 method Methods 0.000 title claims abstract description 24
- 230000000630 rising effect Effects 0.000 claims abstract description 14
- 238000005259 measurement Methods 0.000 claims 1
- 238000002600 positron emission tomography Methods 0.000 description 4
- 230000000052 comparative effect Effects 0.000 description 3
- 230000002123 temporal effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000009191 jumping Effects 0.000 description 1
- 238000002372 labelling Methods 0.000 description 1
- 230000004060 metabolic process Effects 0.000 description 1
- 238000009206 nuclear medicine Methods 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 230000001991 pathophysiological effect Effects 0.000 description 1
- 239000000700 radioactive tracer Substances 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03M—CODING; DECODING; CODE CONVERSION IN GENERAL
- H03M1/00—Analogue/digital conversion; Digital/analogue conversion
- H03M1/06—Continuously compensating for, or preventing, undesired influence of physical parameters
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03M—CODING; DECODING; CODE CONVERSION IN GENERAL
- H03M1/00—Analogue/digital conversion; Digital/analogue conversion
- H03M1/10—Calibration or testing
- H03M1/1009—Calibration
- H03M1/1033—Calibration over the full range of the converter, e.g. for correcting differential non-linearity
- H03M1/1057—Calibration over the full range of the converter, e.g. for correcting differential non-linearity by trimming, i.e. by individually adjusting at least part of the quantisation value generators or stages to their nominal values
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03M—CODING; DECODING; CODE CONVERSION IN GENERAL
- H03M1/00—Analogue/digital conversion; Digital/analogue conversion
- H03M1/10—Calibration or testing
- H03M1/1009—Calibration
- H03M1/1014—Calibration at one point of the transfer characteristic, i.e. by adjusting a single reference value, e.g. bias or gain error
- H03M1/1019—Calibration at one point of the transfer characteristic, i.e. by adjusting a single reference value, e.g. bias or gain error by storing a corrected or correction value in a digital look-up table
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03M—CODING; DECODING; CODE CONVERSION IN GENERAL
- H03M1/00—Analogue/digital conversion; Digital/analogue conversion
- H03M1/12—Analogue/digital converters
- H03M1/34—Analogue value compared with reference values
- H03M1/36—Analogue value compared with reference values simultaneously only, i.e. parallel type
- H03M1/361—Analogue value compared with reference values simultaneously only, i.e. parallel type having a separate comparator and reference value for each quantisation level, i.e. full flash converter type
Abstract
一种多阈值采样数字化器件的阈值校正方法,其步骤为:生成一三角波,测量该三角波波形的上升沿部分的斜率k1、下降沿部分的斜率k2以及其峰值幅度Vpeak,该脉冲高于实际工作阈值Vt部分的宽度DOT表示为DOT(Vt)=(Vpeak-Vt)/k1–(Vpeak–Vt)/k2;设置n组阈值对;根据公式利用测得的脉冲宽度DOT计算出实际工作状态下的阈值;根据实际工作状态阈值与实际设定参考电压之间的对应关系建立阈值校正函数,进而根据该函数对设定的阈值进行校正。该方法可以有效的解决比较器在实际工作状态下因实际设定的设定参考电压与实际工作状态阈值不一致带来的采样精度下降的问题,由于本方法的运用,可以避免MVT器件中使用高精度的比较器,从而有效的降低了MVT器件的成本。
Description
技术领域
本发明涉及正电子发射断层成像设备领域,尤其涉及一种正电子发射断层成像设备中多阈值采样数字化器件的阈值校正方法。
背景技术
正电子发射断层成像(Positron Emission Tomography,全文均简称PET)通过捕捉人体内因正子湮灭而发出的gamma光子获取以正电子核素为标记的示踪剂在人体内的分布情况,进而获取脏器功能,代谢等病理生理特征。获取gamma光子所携带的能量、位置以及时间信息的准确性直接影响到系统成像的性能。闪烁脉冲采集与处理单元是PET系统中的关键核心部件,主要功能为处理前端探测器形成的闪烁脉冲,获取gamma光子所携带的能量、位置以及时间信息。为了保证PET系统的性能,希望所使用的闪烁脉冲采集与处理单元,具有精度高,性能稳定,便于实时校正、集成度高等特点。
多阈值采样(Multi-Voltage Threshold,以下可简称MVT)方法为PET等核医学设备中闪烁脉冲的数字化提供了一套低廉而高效的解决方案。该方法根据闪烁脉冲的特征合理的设置多个电压阈值,通过对闪烁脉冲越过电压阈值的时间的数字化最终实现闪烁脉冲的数字化采样。根据该方法,发展出的MVT数字化器件已经开始应用于PET系统中,这些MVT数字化器件多由比较器和时间数字转换器构成。其中,比较器用于实现闪烁脉冲与电压阈值间的比较,在闪烁脉冲跃过电压阈值的时刻输出跳变信号,时间数字转换器用于对比较器输出的跳变信号的时间进行数字化采样。
现有的这些MVT数字化器件为避免因比较器精度所带来的比较电压与实际电压阈值之间的偏差而引起的采样误差,多采用高精密的比较器。这种高精度比较器的应用虽然在一定程度上解决了上述问题,但也使得整个数字化器件的成本急剧上升。同时这种高精度比较器对工作环境要求苛刻,器件的功耗极高,使得整个MVT数字化器件的成本过于昂贵,通道密度较低。
因此,针对现有技术中存在的问题,有必要提供一种新的多阈值采样数字化器件的阈值校正方法,以降低MVT数字化器件在实现过程中对比较器精度的要求。
发明内容
有鉴于此,本发明的目的在于提供一种多阈值采样数字化器件的阈值校正方法,利用该方法可以获取多阈值采样数字化器件中比较器设置的比较电压与实际工作时比较电压之间的关系,利用该关系可以有效的提高多阈值采样数字化器件的采样精度,降低其对比较器精度的要求。
为实现上述目的,本发明提供如下技术方案:
一种多阈值采样数字化器件的阈值校正方法,其步骤如下:
(1)生成一用于多阈值采样数字化器件的阈值校正的三角波,测量该三角波波形的k1、k2以及Vpeak,该三角波波形的DOT表示为:
其中,k1为三角波波形的上升沿部分的斜率,k2为三角波波形的下降沿部分的斜率,Vpeak为三角波波形的峰值幅度,Vt为实际工作阈值,DOT(Vt)为脉冲高于实际工作阈值部分的宽度;
(2)配置多阈值采样数字化器件中各个比较器的参考电压,并利用该数字化器件对步骤(1)中的脉冲进行数字化采样获取脉冲的采样点;
(3)根据步骤(2)中获得的采样点,计算出每一个比较器在设定的参考电压下的DOT,并根据步骤(1)中所给出的公式计算出每一个比较器的实际工作阈值,将计算出的每个实际工作阈值与设定的参考电压相对应;
(4)重新设定数字化器件中各个比较器的参考电压并重复步骤(2)和(3)m次,针对多阈值采样数字化器件中的每一个比较器获得n组实际工作阈值与设定参考电压的对应关系,其中n为大于等于1的正整数,m为大于等于0的正整数;
(5)将步骤(4)所获得的实际工作阈值及其对应的设定参考电压按其相关的比较器分为k组,k为正整数,k表示多阈值采样数字化器件中比较器的个数;
(6)采用线性拟合的方式获得每组数据实际工作阈值与设定参考电压之间的函数关系,具体如下:
y(a,i)=k(a)×x(a,i)+offset(a),
其中,y(a,i)表示多阈值采样数字化器件中第a个比较器对应的第i对的实际工作阈值,x(a,i)表示多阈值采样数字化器件中第a个比较器对应的第i对的设定的参考电压,k(a)与offset(a)表示与第a个比较器对应的数据组内的实际工作阈值与设定参考电压经过线性拟和得到的系数;
(7)根据步骤(6)中获得的阈值校正函数获得每个比较器在设定的参考电压下的实际阈值电压,达成阀值的校正。
优选的,在上述多阈值采样数字化器件的阈值校正方法中,所述步骤(3)中所述DOT的具体计算方式为:从所获取的脉冲采样点中将经由同一个比较器获得的两个采样点提取出来,计算两个采样点之间的时间间隔,该时间间隔的时间数值即为该比较器在该次采样时与其设置参考电压所对应的实际工作阈值电压下的DOT。
优选的,在上述多阈值采样数字化器件的阈值校正方法中,所述步骤(3)中所述DOT为在相同条件下多次测量后的平均值。
优选的,在上述多阈值采样数字化器件的阈值校正方法中,所述步骤(6)中的阈值校正函数的另一种表达式为:
y(a,i)=x(a,i)+offset(a),
当步骤(4)中的n等于1时,offset(a)为与第a个比较器对应的工作阈值与设定的参考电压之间的差值;当步骤(4)中的n大于1时,offset(a)为与第a个比较器对应的多组工作阈值与设定的参考电压之间差值的平均值。
优选的,在上述多阈值采样数字化器件的阈值校正方法中,步骤(4)中n的值大于2时,步骤(6)中的校正函数对步骤(4)中获得与该比较器对应的多组工作阈值与设定的参考电压进行样条拟合或样条插值获得实际工作阈值与实际设定参考电压之间的关系。
优选的,在上述多阈值采样数字化器件的阈值校正方法中,所述步骤(4)中n的值大于2时,步骤(6)中的校正函数可以对步骤(4)中获得与该比较器对应的多组工作阈值与设定的参考电压进行线性插值的方式获得实际工作阈值与实际设定参考电压之间的关系。
优选的,在上述多阈值采样数字化器件的阈值校正方法中,所述步骤(7)中利用阈值校正函数或查找表根据每个比较器欲设定的阈值电压设置对应的参考电压。
优选的,在上述多阈值采样数字化器件的阈值校正方法中,所述步骤(1)中的三角波可用正弦波替换,该正弦波的公式如下:
测量该正弦波的峰值为2k,周期T和均值电压b,则正弦波高于电压Vt的DOT可以表示为:
优选的,在上述多阈值采样数字化器件的阈值校正方法中,所述步骤(1)中的三角波形可替换为任意波形,该波形在t时刻其电压值可以由确定的数学表达式描述:
其中,t0为该波形中峰值电压Vpeak所对应的时间;则对于电压Vt所对应的DOT可以表示为:
DOT(Vt)=g-1(Vt)-f-1(Vt)。
优选的,在上述多阈值采样数字化器件的阈值校正方法中,所述步骤(1)中的三角波形可替换为任意波形,针对该任意波形其Vt和DOT之间的关系是通过测量该脉冲多个电压下的DOT后建立查找表的方式获得。
从上述技术方案可以看出,本发明的多阈值采样数字化器件的阈值校正方法可以有效的解决比较器在实际工作状态下因实际设定的设定参考电压与实际工作状态阈值不一致带来的采样精度下降的问题,由于本方法的运用,可以避免MVT器件中使用高精度的比较器,从而有效的降低了MVT器件的成本;同时由于比较器的精度的降低,整个器件的功耗将会得到有效的改善,高密度的集成也变为可能。
与现有技术相比,本发明的有益效果是:
(1)通过获取比较器设置阈值与实际工作阈值之间的关系进行阈值校正,避免高精度MVT器件中使用高精度的比较器。
(2)通过利用三角波等特殊波形的脉冲宽度DOT与实际阈值之间的关系完成校正函数的获取。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的有关本发明的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明多阈值采样数字化器件的阈值校正方法的流程图;
图2为MVT采样方法的基本原理图。
具体实施方式
本发明公开了一种能有效的提高多阈值采样数字化器件的采样精度,降低其对比较器精度的要求的多阈值采样数字化器件的阈值校正方法。
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行详细地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图2所示,图2为MVT采样方法的基本原理图。图中A表示闪烁脉冲。图中设置了4个阈值(也可以为其他大于1的值,如无说明下文中所述情况均针对该值为4的情况)分别为V1、V2、V3、V4。脉冲越过这些阈值的时间将会被识别并转换为数字信号。比较器将会被用来识别脉冲越过这些阈值的时刻,时间数字转换器将会被用来将相应得时刻数字化。采用4个阈值的MVT将会产生8个采样点其中4个在脉冲的上升沿,4个在脉冲的下降沿。通常情况下这样的MVT器件需要4个比较器实现,每个比较器对应一个阈值用来识别脉冲超过或低于该阈值的时刻,因此一个比较器将会对应一个阈值同时会涉及到2个采样点,一个在脉冲的上升沿一个在脉冲的下降沿。图中B表示了与阈值V3所对应的比较器识别出的脉冲低于该阈值的时刻。
根据MVT采样原理,对于一个脉冲,由一个比较器决定的两个采样点的时间间隔就表示了这个脉冲在超过该比较器实际工作阈值下的脉冲宽度(Duration over threshold,以下简称DOT)。实际的脉冲都是连续而光滑的,因此上文所说的DOT与比较器实际工作阈值将会是一一对应的关系。当我们确切的知道输入脉冲的形状特点的时候就可以通过MVT器件中每个比较器测量出的DOT计算出该比较器的实际工作阈值,从而建立该比较器的设置的参考电压与实际工作的阈值电压之间的关系,进而达到校正的目的。
如图1所示,根据上述原理,本发明公开的多阈值采样数字化器件的阈值校正方法,其步骤如下:
一种多阈值采样数字化器件的阈值校正方法,其步骤如下:
(1)生成一用于多阈值采样数字化器件的阈值校正的三角波,测量该三角波波形的k1、k2以及Vpeak,该三角波波形的DOT表示为:
其中,k1为三角波波形的上升沿部分的斜率,k2为三角波波形的下降沿部分的斜率,Vpeak为三角波波形的峰值幅度,Vt为实际工作阈值,DOT(Vt)为脉冲高于实际工作阈值部分的宽度;
(2)配置多阈值采样数字化器件中各个比较器的参考电压,并利用该数字化器件对步骤(1)中的脉冲进行数字化采样获取脉冲的采样点;
(3)根据步骤(2)中获得的采样点,计算出每一个比较器在设定的参考电压下的DOT,并根据步骤(1)中所给出的公式计算出每一个比较器的实际工作阈值,将计算出的每个实际工作阈值与设定的参考电压相对应;
(4)重新设定数字化器件中各个比较器的参考电压并重复步骤(2)和(3)m次,针对多阈值采样数字化器件中的每一个比较器获得n组实际工作阈值与设定参考电压的对应关系,其中n为大于等于1的正整数,m为大于等于0的正整数;
(5)将步骤(4)所获得的实际工作阈值及其对应的设定参考电压按其相关的比较器分为k组,k为正整数,k表示多阈值采样数字化器件中比较器的个数;
(6)采用线性拟合的方式获得每组数据实际工作阈值与设定参考电压之间的函数关系,具体如下:
y(a,i)=k(a)×x(a,i)+offset(a),
其中,y(a,i)表示多阈值采样数字化器件中第a个比较器对应的第i对的实际工作阈值,x(a,i)表示多阈值采样数字化器件中第a个比较器对应的第i对的设定的参考电压,k(a)与offset(a)表示与第a个比较器对应的数据组内的实际工作阈值与设定参考电压经过线性拟和得到的系数;
(7)根据步骤(6)中获得的阈值校正函数获得每个比较器在设定的参考电压下的实际阈值电压,达成阀值的校正。
所述步骤(3)中所述DOT的具体计算方式为:从所获取的脉冲采样点中将经由同一个比较器获得的两个采样点提取出来,计算两个采样点之间的时间间隔,该时间间隔的时间数值即为该比较器在该次采样时与其设置参考电压所对应的实际工作阈值电压下的DOT。
所述步骤(3)中所述DOT为在相同条件下多次测量后的平均值。
所述步骤(6)中的阈值校正函数可以表示为:
y(a,i)=x(a,i)+offset(a),
当步骤(4)中的n等于1时,offset(a)为与第a个比较器对应的工作阈值与设定的参考电压之间的差值;当步骤(4)中的n大于1时,offset(a)为与第a个比较器对应的多组工作阈值与设定的参考电压之间差值的平均值。
所述步骤(4)中n的值大于2时,步骤(6)中的校正函数对步骤(4)中获得与该比较器对应的多组工作阈值与设定的参考电压进行样条拟合或样条插值获得实际工作阈值与实际设定参考电压之间的关系。插值和拟合是两个不同的概念。拟合是得到一组数据的函数关系y=f(x),插值是根据该组数据利用一定的模型和假设获得该数据中不存在的数据对的数值。
所述步骤(4)中n的值大于2时,步骤(6)中的校正函数可以对步骤(4)中获得与该比较器对应的多组工作阈值与设定的参考电压进行线性插值的方式获得实际工作阈值与实际设定参考电压之间的关系。
所述步骤(7)中利用阈值校正函数或查找表根据每个比较器欲设定的阈值电压设置对应的参考电压。
所述步骤(1)中的三角波可用正弦波替换,该正弦波的公式如下:
测量该正弦波的峰值为2k,周期T和均值电压b,则正弦波高于电压Vt的DOT可以表示为:
所述步骤(1)中的三角波形可替换为任意波形,该波形在t时刻其电压值可以由确定的数学表达式描述:
其中,t0为该波形中峰值电压Vpeak所对应的时间;则对于电压Vt所对应的DOT可以表示为:
DOT(Vt)=g-1(Vt)-f-1(Vt)。
所述步骤(1)中的三角波形可替换为任意波形,针对该任意波形其Vt和DOT之间的关系是通过测量该脉冲多个电压下的DOT后建立查找表的方式获得。
下面将结合具体的实施例对本发明的方法进行阐述。
实施例一:三角波
(1)利用任意波形发生器产生一三角波,该三角波的上升沿部分的斜率k1=10mV/ns,下降沿部分的斜率k2=-10mV/ns,峰值幅度为500mV。则该波形的DOT可以表示为:
(2)利用多阈值采样数字化器件对该三角波重复采样100次,其中,比较器1、2、3和4的参考电压均设置为40mV。则对于一次采样我们将会获得8个采样时间点,每个比较器对应2个采样点,一个在脉冲的上升沿部分,一个在脉冲的下降沿部分。为方便表述对于第n次采样时比较器i在上升沿获得的采样时间点记录为Tr(n,i),在下降沿获得的采样时间点记录为Tf(n,i),
(3)对于比较器i在参考电压设置为40mV时,其对应的DOT可以表示为:
根据实际测量结果如下表:
比较器1 | 比较器2 | 比较器3 | 比较器4 | |
参考电压 | 40mV | 40mV | 40mV | 40mV |
DOT | 95.2ns | 93.2ns | 96.6ns | 88.6ns |
实际工作阈值 | 24mV | 34mV | 17mV | 57mV |
根据步骤(1)中所给公式可以计算出4个比较器各自的实际工作阈值,如上表。
(4)将比较器1、2、3和4的参考电压均设置为200mV并重复步骤(2)和(3)得到各个比较器在该参考电压下测量出的DOT以及其对应的实际工作阈值,结果如下表
比较器1 | 比较器2 | 比较器3 | 比较器4 | |
参考电压 | 200mV | 200mV | 200mV | 200mV |
DOT | 63.2ns | 61.2ns | 64.6ns | 56.6ns |
实际工作阈值 | 184mV | 194mV | 177mV | 217mV |
(5)对于各个比较器其参考电压与实际工作阈值电压如下表表示
(6)根据步骤(5)中得到的数据采用线性拟合的方式得到各个参考电压与工作阈值之间的函数关系如下:
y(a,i)=k(a)×x(a,i)+offset(a)
其中k(1),k(2),k(3),k(4)均为1,offset(1)、offset(2)、offset(3)、offset(4)分别为16mV,6mV,23mV,-17mV。
(7)当多阈值采样数字化器件中比较器1、2、3和4的参考电压分别设置为40mV,80mV,120mV,160mV的时候根据步骤(6)中的公式可以知道该多阈值采样数字化器件实际工作的阈值电压为24mV、74mV、97mV、177mV。
实施例二:锯齿波(锯齿波是三角波的一种特殊情况)
(1)利用任意波形发生器产生一锯齿波,该锯齿波的上升沿部分的斜率k1=1mV/ns,下降沿部分的斜率k2=-500mV/ns,峰值幅度为200mV。则该波形的DOT可以表示为:
(2)利用多阈值采样数字化器件对该三角波重复采样100次,其中比较器1、2、3和4的参考电压均设置为40mV。则对于一次采样我们将会获得8个采样时间点,每个比较器对应2个采样点一个在脉冲的上升沿部分,一个在脉冲的下降沿部分。为方便表述对于第n次采样时比较器i在上升沿获得的采样时间点记录为Tr(n,i),在下降沿获得的采样时间点记录为Tf(n,i),
(3)对于比较器i在参考电压设置为40mV时,其对应的DOT可以表示为:
根据实际测量结果如下表:
比较器1 | 比较器2 | 比较器3 | 比较器4 | |
参考电压 | 40mV | 40mV | 40mV | 40mV |
DOT | 176.2ns | 166.1ns | 183.4ns | 143.6ns |
实际工作阈值 | 24.2mV | 34.3mV | 17.0mV | 56.8mV |
根据步骤(1)中所给公式可以计算出4个比较器各自的实际工作阈值,如上表。
(4)将比较器1、2、3和4的参考电压均设置为200mV并重复步骤(2)
和(3)得到各个比较器在该参考电压下测量出的DOT以及其对应的实际工作阈值,结果如下表
比较器1 | 比较器2 | 比较器3 | 比较器4 | |
参考电压 | 140mV | 140mV | 140mV | 140mV |
DOT | 76.8ns | 67.2ns | 83.3ns | 43.1ns |
实际工作阈值 | 123.6mV | 133.2mV | 117.1mV | 157.3mV |
(5)对于各个比较器其参考电压与实际工作阈值电压如下表表示
(6)根据步骤(5)中得到的数据采用线性拟合的方式得到各个参考电压与工作阈值之间的函数关系如下:
y(a,i)=k(a)×x(a,i)+offset(a)
其中k(1),k(2),k(3),k(4)均为1,offset(1)、offset(2)、offset(3)、offset(4)分别为16mV,6mV,23mV,-17mV。
(7)当多阈值采样数字化器件中比较器1、2、3和4的参考电压分别设置为40mV,80mV,120mV,160mV的时候根据步骤(6)中的公式可以知道该多阈值采样数字化器件实际工作的阈值电压为24mV、74mV、97mV、177mV。
实施例三:正弦波
(1)利用任意波形发生器产生正弦波,该正弦波的周期T=200ns,均值电b=200mV,峰峰值为400mV即k=200mV。则该波形的DOT可以表示为
(2)利用多阈值采样数字化器件对该三角波重复采样100次,其中比较器1、2、3和4的参考电压均设置为40mV。则对于一次采样我们将会获得8个采样时间点,每个比较器对应2个采样点一个在脉冲的上升沿部分,一个在脉冲的下降沿部分。为方便表述对于第n次采样时比较器i在上升沿获得的采样时间点记录为Tr(n,i),在下降沿获得的采样时间点记录为Tf(n,i),
(3)对于比较器i在参考电压设置为40mV时,其对应的DOT可以表示为:
根据实际测量结果如下表:
比较器1 | 比较器2 | 比较器3 | 比较器4 | |
参考电压 | 40mV | 40mV | 40mV | 40mV |
DOT | 168.6ns | 162.4ns | 173.6ns | 150.7ns |
实际工作阈值 | 23.9mV | 33.8mV | 16.9mV | 57.0mV |
根据步骤(1)中所给公式可以计算出4个比较器各自的实际工作阈值,如上表。
(4)将比较器1、2、3和4的参考电压均设置为200mV并重复步骤(2)和(3)得到各个比较器在该参考电压下测量出的DOT以及其对应的实际工作阈值,结果如下表
比较器1 | 比较器2 | 比较器3 | 比较器4 | |
参考电压 | 140mV | 140mV | 140mV | 140mV |
DOT | 124.8mV | 121.3mV | 127.2mV | 113.8mV |
实际工作阈值 | 124.0mV | 134.2mV | 117.0mV | 157.1mV |
(5)对于各个比较器其参考电压与实际工作阈值电压如下表表示
(6)根据步骤(5)中得到的数据采用线性拟合的方式得到各个参考电压与工作阈值之间的函数关系如下:
y(a,i)=k(a)×x(a,i)+offset(a)
其中k(1),k(2),k(3),k(4)均为1,offset(1)、offset(2)、offset(3)、offset(4)分别为16mV,6mV,23mV,-17mV。
(7)当多阈值采样数字化器件中比较器1、2、3和4的参考电压分别设置为40mV,80mV,120mV,160mV的时候根据步骤(6)中的公式可以知道该多阈值采样数字化器件实际工作的阈值电压为24mV、74mV、97mV、177mV。
本发明实施例的多阈值采样数字化器件的阈值校正方法可以有效的解决比较器在实际工作状态下因实际设定的设定参考电压与实际工作状态阈值不一致带来的采样精度下降的问题,由于本方法的运用,可以避免MVT器件中使用高精度的比较器,从而有效的降低了MVT器件的成本。同时由于比较器的精度的降低整个器件的功耗将会得到有效的改善,高密度的集成也变为可能。
本发明与现有技术的区别在于:
(1)通过获取比较器设置阈值与实际工作阈值之间的关系进行阈值校正,避免高精度MVT器件中使用高精度的比较器。
(2)通过利用三角波等特殊波形的DOT与实际阈值之间的关系完成校正函数的获取。
对于本领域技术人员而言,显然本发明不限于上述示范性实施例的细节,而且在不背离本发明的精神或基本特征的情况下,能够以其他的具体形式实现本发明。因此,无论从哪一点来看,均应将实施例看作是示范性的,而且是非限制性的,本发明的范围由所附权利要求而不是上述说明限定,因此旨在将落在权利要求的等同要件的含义和范围内的所有变化囊括在本发明内。不应将权利要求中的任何附图标记视为限制所涉及的权利要求。
此外,应当理解,虽然本说明书按照实施方式加以描述,但并非每个实施方式仅包含一个独立的技术方案,说明书的这种叙述方式仅仅是为清楚起见,本领域技术人员应当将说明书作为一个整体,各实施例中的技术方案也可以经适当组合,形成本领域技术人员可以理解的其他实施方式。
Claims (10)
1.一种多阈值采样数字化器件的阈值校正方法,其特征在于:步骤如下:
(1)生成一用于多阈值采样数字化器件的阈值校正的三角波,测量该三角波波形的k1、k2以及Vpeak,该三角波波形的DOT表示为:
其中,k1为三角波波形的上升沿部分的斜率,k2为三角波波形的下降沿部分的斜率,Vpeak为三角波波形的峰值幅度,Vt为实际工作阈值,DOT(Vt)为脉冲高于实际工作阈值部分的宽度;
(2)配置多阈值采样数字化器件中各个比较器的参考电压,并利用该数字化器件对步骤(1)中的脉冲进行数字化采样获取脉冲的采样点;
(3)根据步骤(2)中获得的采样点,计算出每一个比较器在设定的参考电压下的DOT,并根据步骤(1)中所给出的公式计算出每一个比较器的实际工作阈值,将计算出的每个实际工作阈值与设定的参考电压相对应;
(4)重新设定数字化器件中各个比较器的参考电压并重复步骤(2)和(3)m次,针对多阈值采样数字化器件中的每一个比较器获得n组实际工作阈值与设定参考电压的对应关系,其中n为大于等于1的正整数,m为大于等于0的正整数;
(5)将步骤(4)所获得的实际工作阈值及其对应的设定参考电压按其相关的比较器分为k组,k为正整数,k表示多阈值采样数字化器件中比较器的个数;
(6)采用线性拟合的方式获得每组数据实际工作阈值与设定参考电压之间的函数关系,具体如下:
y(a,i)=k(a)×x(a,i)+offset(a),
其中,y(a,i)表示多阈值采样数字化器件中第a个比较器对应的第i对的实际工作阈值,x(a,i)表示多阈值采样数字化器件中第a个比较器对应的第i对的设定的参考电压,k(a)与offset(a)表示与第a个比较器对应的数据组内的实际工作阈值与设定参考电压经过线性拟和得到的系数;
(7)根据步骤(6)中获得的阈值校正函数获得每个比较器在设定的参考电压下的实际阈值电压,达成阀值的校正。
2.根据权利要求1所述的多阈值采样数字化器件的阈值校正方法,其特征在于:所述步骤(3)中所述DOT的具体计算方式为:从所获取的脉冲采样点中将经由同一个比较器获得的两个采样点提取出来,计算两个采样点之间的时间间隔,该时间间隔的时间数值即为该比较器在该次采样时与其设置参考电压所对应的实际工作阈值电压下的DOT。
3.根据权利要求1所述的多阈值采样数字化器件的阈值校正方法,其特征在于:所述步骤(3)中所述DOT为在相同条件下多次测量后的平均值。
4.根据权利要求1所述的多阈值采样数字化器件的阈值校正方法,其特征在于:所述步骤(6)中的阈值校正函数的另一种表达式为:
y(a,i)=x(a,i)+offset(a),
当步骤(4)中的n等于1时,offset(a)为与第a个比较器对应的工作阈值与设定的参考电压之间的差值;当步骤(4)中的n大于1时,offset(a)为与第a个比较器对应的多组工作阈值与设定的参考电压之间差值的平均值。
5.根据权利要求1所述的多阈值采样数字化器件的阈值校正方法,其特征在于:所述步骤(4)中n的值大于2时,步骤(6)中的校正函数对步骤(4)中获得与该比较器对应的多组工作阈值与设定的参考电压进行样条拟合或样条插值获得实际工作阈值与实际设定参考电压之间的关系。
6.根据权利要求1所述的多阈值采样数字化器件的阈值校正方法,其特征在于:所述步骤(4)中n的值大于2时,步骤(6)中的校正函数可以对步骤(4)中获得与该比较器对应的多组工作阈值与设定的参考电压进行线性插值的方式获得实际工作阈值与实际设定参考电压之间的关系。
7.根据权利要求1所述的多阈值采样数字化器件的阈值校正方法,其特征在于:所述步骤(7)中利用阈值校正函数或查找表根据每个比较器欲设定的阈值电压设置对应的参考电压。
8.根据权利要求1所述的多阈值采样数字化器件的阈值校正方法,其特征在于:所述步骤(1)中的三角波可用正弦波替换,该正弦波的公式如下:
测量该正弦波的峰值为2k,周期T和均值电压b,则正弦波高于电压Vt的DOT可以表示为:
9.根据权利要求1所述的多阈值采样数字化器件的阈值校正方法,其特征在于:所述步骤(1)中的三角波形可替换为任意波形,该波形在t时刻其电压值可以由确定的数学表达式描述:
其中,t0为该波形中峰值电压Vpeak所对应的时间;则对于电压Vt所对应的DOT可以表示为:
DOT(Vt)=g-1(Vt)-f-1(Vt)
10.根据权利要求1所述的多阈值采样数字化器件的阈值校正方法,其特征在于:所述步骤(1)中的三角波形可替换为任意波形,针对该任意波形其DOT和Vt之间的关系是通过实际测量该脉冲多个电压下的DOT后建立查找表的方式获得。
Priority Applications (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310045681.9A CN103961126B (zh) | 2013-02-05 | 2013-02-05 | 一种多阈值采样数字化器件的阈值校正方法 |
EP13874702.7A EP2954838A4 (en) | 2013-02-05 | 2013-02-28 | THRESHOLD CORRECTION FOR MULTI-VOLTAGE-SWITCHING DIGITIZATION DEVICE |
PCT/CN2013/072006 WO2014121530A1 (zh) | 2013-02-05 | 2013-02-28 | 一种多阈值采样数字化器件的阈值校正方法 |
US14/765,551 US9344102B2 (en) | 2013-02-05 | 2013-02-28 | Threshold correction method for multi-voltage threshold sampling digitization device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310045681.9A CN103961126B (zh) | 2013-02-05 | 2013-02-05 | 一种多阈值采样数字化器件的阈值校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103961126A true CN103961126A (zh) | 2014-08-06 |
CN103961126B CN103961126B (zh) | 2016-07-06 |
Family
ID=51231476
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310045681.9A Active CN103961126B (zh) | 2013-02-05 | 2013-02-05 | 一种多阈值采样数字化器件的阈值校正方法 |
Country Status (4)
Country | Link |
---|---|
US (1) | US9344102B2 (zh) |
EP (1) | EP2954838A4 (zh) |
CN (1) | CN103961126B (zh) |
WO (1) | WO2014121530A1 (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107422364A (zh) * | 2017-06-26 | 2017-12-01 | 沈阳东软医疗系统有限公司 | 一种检测器的时间定时校正方法及装置 |
CN109171787A (zh) * | 2018-08-27 | 2019-01-11 | 苏州瑞派宁科技有限公司 | 脉冲信号的采样方法、装置和计算机程序介质 |
CN109444559A (zh) * | 2018-10-26 | 2019-03-08 | 苏州瑞迈斯医疗科技有限公司 | 脉冲信号的采样方法、重建方法和装置 |
CN110811659A (zh) * | 2019-10-17 | 2020-02-21 | 湖北锐世数字医学影像科技有限公司 | 时间数据处理方法、装置及系统 |
CN110954938A (zh) * | 2019-12-11 | 2020-04-03 | 湖北锐世数字医学影像科技有限公司 | 比较器的阈值校正方法、装置和采样电路 |
CN110988961A (zh) * | 2019-12-11 | 2020-04-10 | 湖北锐世数字医学影像科技有限公司 | 信号处理方法、装置和探测系统 |
CN111158039A (zh) * | 2020-01-02 | 2020-05-15 | 苏州瑞派宁科技有限公司 | 信号采样、重建方法及装置 |
US10755452B2 (en) | 2017-07-10 | 2020-08-25 | Shanghai Neusoft Medical Technology Co., Ltd. | Calibrating time in PET device |
CN111669182A (zh) * | 2020-06-24 | 2020-09-15 | 苏州瑞迈斯医疗科技有限公司 | 信号采样装置、系统及方法 |
CN111929569A (zh) * | 2020-09-18 | 2020-11-13 | 深圳英集芯科技有限公司 | Ic芯片的校准方法、系统及装置 |
CN112290918A (zh) * | 2020-10-30 | 2021-01-29 | 湖北锐世数字医学影像科技有限公司 | 用于对比较器进行阈值校正的方法、装置及系统 |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9843337B1 (en) | 2017-03-16 | 2017-12-12 | Analog Devices Global | Background flash offset calibration in continuous-time delta-sigma ADCS |
KR101930402B1 (ko) * | 2017-04-17 | 2018-12-18 | 서울대학교산학협력단 | 톱니모양 문턱전압을 이용한 시간 기반의 신호 획득 장치 및 방법 |
CN110336313B (zh) * | 2019-07-11 | 2022-07-12 | 南方电网科学研究院有限责任公司 | 一种触发角的校正方法、装置、设备及介质 |
CN113533950B (zh) * | 2021-06-11 | 2023-12-26 | 国网浙江省电力有限公司嘉兴供电公司 | 一种基于BiLSTM的高压断路器操作电压采样方法 |
CN115113512B (zh) * | 2022-05-26 | 2024-03-29 | 合肥综合性国家科学中心人工智能研究院(安徽省人工智能实验室) | 脉冲的采样方法、采样系统、装置及计算机可读存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040049120A1 (en) * | 2002-09-11 | 2004-03-11 | Jian Cao | Methods and apparatus for cardiac R-wave sensing in a subcutaneous ECG waveform |
US20060055985A1 (en) * | 2004-09-13 | 2006-03-16 | Canon Kabushiki Kaisha | Color-noise reduction circuit and image-sensing apparatus therewith |
US20090103361A1 (en) * | 2007-05-04 | 2009-04-23 | Lee Wang | Level verification and adjustment for multi-level cell (mlc) non-volatile memory (nvm) |
CN102843139A (zh) * | 2012-09-20 | 2012-12-26 | 苏州瑞派宁科技有限公司 | 一种闪烁脉冲数字化的方法及装置 |
Family Cites Families (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6822485B2 (en) * | 2002-08-14 | 2004-11-23 | Guide Technology | Method for calibrating threshold levels on comparators with dithered DC signals |
US7046014B2 (en) * | 2004-02-27 | 2006-05-16 | Broadcom Corporation | ATE measurement technique for comparator threshold voltage |
JP2006042050A (ja) * | 2004-07-28 | 2006-02-09 | Nec Corp | 送信装置およびピーク低減方法 |
JP2006121387A (ja) * | 2004-10-21 | 2006-05-11 | Nec Corp | 識別再生方法および識別再生装置 |
JP2007085932A (ja) * | 2005-09-22 | 2007-04-05 | Verigy (Singapore) Pte Ltd | 入力閾値レベル測定方法および入力閾値レベル測定装置 |
CN101454977B (zh) * | 2006-03-31 | 2011-06-15 | 安立股份有限公司 | 波形整形装置以及差错测定装置 |
TWI331455B (en) * | 2007-01-30 | 2010-10-01 | Univ Nat Chiao Tung | A precisely self-calibrating, low power and high speed analog to digital converter |
US7750732B1 (en) * | 2007-12-04 | 2010-07-06 | Fairchild Semiconductor Corporation | Adaptive rail amplifier (ARA) technology |
JP4439572B2 (ja) * | 2008-07-11 | 2010-03-24 | 任天堂株式会社 | ディジタルデータ補正プログラムおよびディジタルデータ補正装置 |
JP5427658B2 (ja) * | 2010-03-16 | 2014-02-26 | パナソニック株式会社 | コンパレータのオフセット補正装置 |
FR2966663B1 (fr) * | 2010-10-26 | 2013-07-05 | St Microelectronics Grenoble 2 | Dispositif et procede de traitement d'un signal analogique |
JP2012191412A (ja) * | 2011-03-10 | 2012-10-04 | Advantest Corp | A/dコンバータの試験装置および試験方法 |
US8649471B1 (en) * | 2011-12-08 | 2014-02-11 | Exelis Inc. | Application of crest factor reduction to a signal |
JP5807549B2 (ja) * | 2012-01-10 | 2015-11-10 | 富士通株式会社 | 比較回路およびa/d変換回路 |
US8779952B1 (en) * | 2012-04-12 | 2014-07-15 | Pmc-Sierra Us, Inc. | Background calibration of threshold errors in analog to digital converters |
JP5942798B2 (ja) * | 2012-11-12 | 2016-06-29 | 富士通株式会社 | 比較回路およびa/d変換回路 |
JP5965825B2 (ja) * | 2012-11-30 | 2016-08-10 | 株式会社ソシオネクスト | コンパレータ及びその補正方法 |
-
2013
- 2013-02-05 CN CN201310045681.9A patent/CN103961126B/zh active Active
- 2013-02-28 EP EP13874702.7A patent/EP2954838A4/en active Pending
- 2013-02-28 US US14/765,551 patent/US9344102B2/en active Active
- 2013-02-28 WO PCT/CN2013/072006 patent/WO2014121530A1/zh active Application Filing
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040049120A1 (en) * | 2002-09-11 | 2004-03-11 | Jian Cao | Methods and apparatus for cardiac R-wave sensing in a subcutaneous ECG waveform |
US20060055985A1 (en) * | 2004-09-13 | 2006-03-16 | Canon Kabushiki Kaisha | Color-noise reduction circuit and image-sensing apparatus therewith |
US20090103361A1 (en) * | 2007-05-04 | 2009-04-23 | Lee Wang | Level verification and adjustment for multi-level cell (mlc) non-volatile memory (nvm) |
CN102843139A (zh) * | 2012-09-20 | 2012-12-26 | 苏州瑞派宁科技有限公司 | 一种闪烁脉冲数字化的方法及装置 |
Non-Patent Citations (1)
Title |
---|
QINGGUO XIE 等: "An investigation of baseline calibration method for digitally sampling scintillation pulses in PET", 《NUCLEAR SCIENCE SYMPOSIUM CONFERENCE RECORD,2010IEEE》 * |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10162068B1 (en) | 2017-06-26 | 2018-12-25 | Shenyang Neusoft Medical Systems Co., Ltd. | Calibrating time in PET device |
CN107422364A (zh) * | 2017-06-26 | 2017-12-01 | 沈阳东软医疗系统有限公司 | 一种检测器的时间定时校正方法及装置 |
US10755452B2 (en) | 2017-07-10 | 2020-08-25 | Shanghai Neusoft Medical Technology Co., Ltd. | Calibrating time in PET device |
CN109171787A (zh) * | 2018-08-27 | 2019-01-11 | 苏州瑞派宁科技有限公司 | 脉冲信号的采样方法、装置和计算机程序介质 |
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 |
CN109444559B (zh) * | 2018-10-26 | 2021-02-05 | 苏州瑞迈斯医疗科技有限公司 | 脉冲信号的采样方法、重建方法和装置 |
CN109444559A (zh) * | 2018-10-26 | 2019-03-08 | 苏州瑞迈斯医疗科技有限公司 | 脉冲信号的采样方法、重建方法和装置 |
CN110811659A (zh) * | 2019-10-17 | 2020-02-21 | 湖北锐世数字医学影像科技有限公司 | 时间数据处理方法、装置及系统 |
CN110811659B (zh) * | 2019-10-17 | 2022-11-18 | 湖北锐世数字医学影像科技有限公司 | 时间数据处理方法、装置及系统 |
CN110954938B (zh) * | 2019-12-11 | 2021-08-31 | 湖北锐世数字医学影像科技有限公司 | 比较器的阈值校正方法、装置和采样电路 |
CN110988961B (zh) * | 2019-12-11 | 2021-09-07 | 湖北锐世数字医学影像科技有限公司 | 信号处理方法、装置和探测系统 |
CN110988961A (zh) * | 2019-12-11 | 2020-04-10 | 湖北锐世数字医学影像科技有限公司 | 信号处理方法、装置和探测系统 |
CN110954938A (zh) * | 2019-12-11 | 2020-04-03 | 湖北锐世数字医学影像科技有限公司 | 比较器的阈值校正方法、装置和采样电路 |
WO2021135338A1 (zh) * | 2020-01-02 | 2021-07-08 | 苏州瑞派宁科技有限公司 | 信号采样、重建方法及装置 |
CN111158039A (zh) * | 2020-01-02 | 2020-05-15 | 苏州瑞派宁科技有限公司 | 信号采样、重建方法及装置 |
CN111669182A (zh) * | 2020-06-24 | 2020-09-15 | 苏州瑞迈斯医疗科技有限公司 | 信号采样装置、系统及方法 |
CN111669182B (zh) * | 2020-06-24 | 2023-11-03 | 苏州瑞迈斯科技有限公司 | 信号采样装置、系统及方法 |
CN111929569A (zh) * | 2020-09-18 | 2020-11-13 | 深圳英集芯科技有限公司 | Ic芯片的校准方法、系统及装置 |
CN112290918A (zh) * | 2020-10-30 | 2021-01-29 | 湖北锐世数字医学影像科技有限公司 | 用于对比较器进行阈值校正的方法、装置及系统 |
WO2022089080A1 (zh) * | 2020-10-30 | 2022-05-05 | 湖北锐世数字医学影像科技有限公司 | 用于对比较器进行阈值校正的方法、装置及系统 |
CN112290918B (zh) * | 2020-10-30 | 2024-03-26 | 湖北锐世数字医学影像科技有限公司 | 用于对比较器进行阈值校正的方法、装置及系统 |
Also Published As
Publication number | Publication date |
---|---|
US9344102B2 (en) | 2016-05-17 |
EP2954838A4 (en) | 2016-11-02 |
US20150372689A1 (en) | 2015-12-24 |
EP2954838A1 (en) | 2015-12-16 |
CN103961126B (zh) | 2016-07-06 |
WO2014121530A1 (zh) | 2014-08-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103961126A (zh) | 一种多阈值采样数字化器件的阈值校正方法 | |
CN105824817B (zh) | 一种闪烁脉冲的数字化方法 | |
CN102736511B (zh) | 时间测量系统及时间测量方法 | |
CN109752584B (zh) | 一种周期信号有效值测量方法 | |
CN109581461B (zh) | 核脉冲能量测量方法及系统 | |
CN104808056A (zh) | 一种基于比较器转换的频率特性测试方法与装置 | |
Larios et al. | A computational investigation of the finite-time blow-up of the 3D incompressible Euler equations based on the Voigt regularization | |
CN102843139A (zh) | 一种闪烁脉冲数字化的方法及装置 | |
CN106771456A (zh) | 基于功率谱熵的数字示波器触发方法及系统 | |
CN104991118A (zh) | 一种高分辨异频信号频率测量系统和测量方法 | |
EP2044692B1 (en) | Undersampling of a repetitive signal for measuring transistion times to reconstruct an analog waveform | |
CN103135650B (zh) | 电流/频率变换电路线性度及对称性数字补偿方法 | |
CN103870238A (zh) | 基于fpga的伪随机序列发生器及其生成方法 | |
CN110988961B (zh) | 信号处理方法、装置和探测系统 | |
JP7418871B2 (ja) | 信号サンプリング方法及び再構成方法、並びに装置 | |
CN111562537A (zh) | 电能表动态误差的同步测量方法 | |
RU58825U1 (ru) | Аналого-цифровой преобразователь | |
CN106483546A (zh) | 信号处理装置及放射线测定装置 | |
CN110954938B (zh) | 比较器的阈值校正方法、装置和采样电路 | |
CN205388624U (zh) | 基于集成鉴相鉴频器的相位差测量电路 | |
CN104849551B (zh) | 一种谐相角分析方法 | |
CN209841927U (zh) | 带触发的表源一体同步动态畸变功率源 | |
CN102739200B (zh) | 信号发生器 | |
JP2010160093A (ja) | 波形測定装置 | |
Liu et al. | Enhanced MVT digitizer in PET |
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 | ||
EE01 | Entry into force of recordation of patent licensing contract | ||
EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20140806 Assignee: Suzhou Ruimaisi Medical Technology Co Ltd Assignor: Huazhong University of Science & Technology Contract record no.: 2017320010042 Denomination of invention: Threshold value rectifying method of multi-threshold value sampling digitalization device Granted publication date: 20160706 License type: Common License Record date: 20171018 |