CN1264474C - 剔除反常数据方法和应用该方法的血液成分光谱分析系统 - Google Patents

剔除反常数据方法和应用该方法的血液成分光谱分析系统 Download PDF

Info

Publication number
CN1264474C
CN1264474C CNB2004100283223A CN200410028322A CN1264474C CN 1264474 C CN1264474 C CN 1264474C CN B2004100283223 A CNB2004100283223 A CN B2004100283223A CN 200410028322 A CN200410028322 A CN 200410028322A CN 1264474 C CN1264474 C CN 1264474C
Authority
CN
China
Prior art keywords
wavelength
light
ratio
individual
parameter
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.)
Expired - Fee Related
Application number
CNB2004100283223A
Other languages
English (en)
Other versions
CN1519557A (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.)
Samsung Electronics Co Ltd
Original Assignee
Samsung Electronics 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 Samsung Electronics Co Ltd filed Critical Samsung Electronics Co Ltd
Publication of CN1519557A publication Critical patent/CN1519557A/zh
Application granted granted Critical
Publication of CN1264474C publication Critical patent/CN1264474C/zh
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/145—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
    • A61B5/1455—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/145—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
    • A61B5/14532—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue for measuring glucose, e.g. by tissue impedance measurement
    • 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
    • Y10—TECHNICAL SUBJECTS COVERED BY FORMER USPC
    • Y10T—TECHNICAL SUBJECTS COVERED BY FORMER US CLASSIFICATION
    • Y10T436/00—Chemistry: analytical and immunological testing
    • Y10T436/10—Composition for standardization, calibration, simulation, stabilization, preparation or preservation; processes of use in preparation for chemical testing
    • Y10T436/102499—Blood gas standard or control
    • 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
    • Y10—TECHNICAL SUBJECTS COVERED BY FORMER USPC
    • Y10T—TECHNICAL SUBJECTS COVERED BY FORMER US CLASSIFICATION
    • Y10T436/00—Chemistry: analytical and immunological testing
    • Y10T436/10—Composition for standardization, calibration, simulation, stabilization, preparation or preservation; processes of use in preparation for chemical testing
    • Y10T436/105831—Protein or peptide standard or control [e.g., hemoglobin, etc.]
    • 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
    • Y10—TECHNICAL SUBJECTS COVERED BY FORMER USPC
    • Y10T—TECHNICAL SUBJECTS COVERED BY FORMER US CLASSIFICATION
    • Y10T436/00—Chemistry: analytical and immunological testing
    • Y10T436/10—Composition for standardization, calibration, simulation, stabilization, preparation or preservation; processes of use in preparation for chemical testing
    • Y10T436/106664—Blood serum or blood plasma standard or control

Abstract

本发明提供了一种剔除反常数据的方法和一种利用光谱学技术作血液成分分析的系统。在利用光谱学技术作血液成分分析的系统中剔除反常数据的方法包括：(a)在预定的单位时间周期内采集与第一和第二波长相对应的PPG信号；(b)计算“n”个参数，“n”为所采集的PPG信号中包含的脉冲数据的数目，其中n为正整数；(c)计算“n”个参数的平均值；和(d)将对平均值的偏差大于预定标准偏差的参数数目与“n”个参数的比率与预定的剔除参考值相比较，以便确定“n”个脉冲数据是否有效。

Description

剔除反常数据方法和应用该方法的 血液成分光谱分析系统
技术领域
本发明涉及一种血液成分分析系统，更具体地说，涉及一种剔除由于人体内部和外部生命要素而产生的反常数据的方法，以及应用了此方法的血液成分光谱学分析系统。
背景技术
在现代医学工程领域中，用非损伤性的方法测量生理信号是一个重要的问题，人们已经研制出一些测量各种生理变量的方法。特别是已经发表了多种无需采集血样就能分析血液成分的方法，还进行了有关无需采血而作血液成分分析的研究。例如，在一种用光谱学分析血液成分的方法中，无需采集血样，病人就不必忍受疼痛，也不会受到感染，而分析结果也能够实时得到。由于这些优点以及相关技术的发展，用于测量氧饱和度的传统仪器的常规技术已应用于分析诸如血红蛋白和葡萄糖等血液成分的医疗器械当中。这些医疗器械利用具有特定波长的光对血液成分发生反应，通过测量由被人体吸收的不同波长的光量的差值的方法来分析血液成分。
当使用光谱学技术(spectroscopy)时会产生一种光体积描记(photo-plethysmographic)(PPG)信号，其包含能脉动的成分和非脉动的成分。将结合附图1对所述PPG信号作进一步的描述。在附图1中，参考标号lo表示放射到人体的光的量，参考标号lt表示透过人体的光的量，参考标号la表示被人体吸收的光的量，参考标号To表示一个心跳周期，参考标号lp表示脉动成分的最大值点，参考标号lv表示脉动成分的最小值点，参考标号P1表示由脉动成分即交流(AC)成分造成的光强度的变化，而参考标号P2表示由非脉动成分即直流(DC)成分造成的光强度的变化。吸收发射到人体的光的成分大部分属于非脉动的成分，即P2成分，例如骨骼以及不会随时间而改变的重要器官组织(vital tissue)，而主要的脉动成分，即P1成分由于心脏的跳动会随时间而改变。
对于用比率计来分析血液成分的方法来说，必须精确地测量随着时间而改变的脉动成分吸收的光量。然而，由于内部生命要素，如呼吸、血压、脉搏率、体温、血管状态或者植物性神经系统，和外部生命要素，如自发或非自发的运动，PPG信号的振幅和基线会经常变化。在定量条件下，与非脉动成分吸收的光量相比，脉动成分的振幅变化范围在2至5％，PPG信号的基线变化范围在3至5％。由于这些内部生命要素造成的脉动成分的振幅和PPG信号的基线的变化导致了反常数据。结果，在利用光谱学技术所作的血液成分分析中就产生了误差。
美国专利号5025791、5517988、5645060以及6393311中披露了PPG信号中产生反常数据的传统解决方法。在美国专利号5025791中，在放置探针的身体部位连接运动传感器和加速计，超出预定极限的运动会被检测到并形成与之对应的信号，将其传输至控制系统的微处理器上，于是在测出有运动的时间内测得的数据会在测量氧饱和度的分析中被剔除。在所述技术中，由于运动传感器和加速计是作为光学单元的附件使用的，该系统的结构复杂，因而提高了测量的成本。
美国专利5517988中，提出了表明那种数据是由非生理因素引起的实例，这是基于有关PPG信号的先前预定的标准，例如AC成分振幅的预定变化或者DC成分的预定斜率，这样就可以避免由于反常数据而造成的错误分析。然而，即使当PPG信号的AC和DC成分是稳定的，由于内部生命要素的影响，在很短时间内变量仍会变化很大，所述变量例如用于回归方程中，并在单位测量时间内从脉动成分中测得的比率的比值。
在美国专利号5645060中，包含运动和噪声项的Beer-Lambert’s方程中比率的比值被用在了比率计中以除去噪声的影响。在假设每个波长信号中运动/噪声值是成比例的，并且脉动成分不会受运动影响的情况下，这项技术是有效的。然而，脉动成分根据呼吸变化很大并影响了用于计算血液成分浓度的PPG信号，因此计算结果并不精确。
美国专利号6393311中，在单个的数据段上使用了傅立叶变换来找出具有代表性的频率，设计了一个带通滤波器来提取只与代表性频率相对应的信号，并且提取出从那些运动引起的噪声中已经去除的信号。在这项技术中，傅立叶变换中的计算量和用于窄带过滤的滤波器的抽头的数目都增加了，因此需要处理大量的数据。结果，很难实时地实现连续的信号处理。
正如上面所描述的，传统技术的方法仅局限于排除由于外部生命要素(例如在测量氧饱和度系统中由运动引起的噪声和探针的不稳定接触)产生的反常数据。再者，在传统技术中，诸如氧饱和信号等生理信号的测量是在预定时间过后直到病人生理上稳定下来才测得的，因此需要很长测量时间。
发明内容
本发明提供了一种在血液成分光谱分析系统中，剔除由于内部和外部生命要素(例如运动引起的噪声或探针的不稳定接触)所产生的反常数据的方法。
本发明还提供了应用上述方法剔除反常数据以便增加分析精确度并减少测量时间的一种血液成分光谱分析系统。
根据本发明的一个方面，提供了一种利用光谱学技术在估测血液组份浓度的血液成分光谱分析系统中剔除反常数据的方法，所述方法通过对通过将对血液成分作出反应的第一和第二波长的光发射到试样上和对穿过试样的所述第一和第二波长的光进行光电转换所获得的光体积描记(PPG)信号进行分析，来估测血液成分的浓度，信号，所述方法包括：(a)在预定的单位时间周期内采集与第一和第二波长相对应的PPG信号；(b)计算“n”个参数，“n”为所采集PPG信号中包含的脉冲数据的数目，为正整数；(c)计算“n”个参数的平均值；(d)将对其平均值的偏差大于预定标准偏差的参数数目与“n”个参数的比率与预定的剔除参考值相比较，以便确定“n”个脉冲数据是否有效，其中，所述参数为第一波长光的每个脉冲的最大值点和最小值点的比率对数与第二波长光的每个脉冲的最大值点和最小值点的比率对数的比值。
最好地，所述方法还包括：(e)更新全部或“n”个脉冲数据中的一些，并且当作为步骤(d)中的比较结果的比率等于或大于预定的剔除参考值时重复步骤(b)至(d)。
根据本发明的另一个方面，提供了一个利用光谱技术的血液成分分析系统，包括：光发射装置，其包含至少两个能发射对特定的血液成分作出反应的波长的光的光源，并且能根据预定的定时向试样发射特定波长的光；光接收装置，用于接收从试样发出的光并将所接收的光转换成电信号；以及信号处理器，用于从由光接收装置所接收的电信号中提取与对特定的血液成分作出反应的特定的波长的光的体积描记(PPG)信号，对于在预定单元时间周期内采集的与所述波长相对应的PPG信号中包括的“n”个脉冲，来计算“n”个参数的平均值，以及将对其平均值的偏差大于预定标准偏差的参数的数目与“n”个参数的比率与预定的剔除参考值相比较，以便确定“n”个脉冲数据是否有效，这里n为正整数，并且其中，所述参数为第一波长光的每个脉冲的最大值点和最小值点的比率对数与第二波长光的每个脉冲的最大值点和最小值点的比率对数的比值。
所述信号处理器最好包括用于存储“n”个脉冲数据的先入先出(FIFO)数据缓冲器，和当所述“n”个脉冲数据被确定为无效时，将所述“n”个脉冲数据划分成预定的单元组，用与第一预定单元组相同数目的新数据来取代由初始输入数据组成的所述第一预定单元组，和确定新的“n”个脉冲数据是否有效。
附图说明
通过结合其附图详细描述优选实施例，本发明上述的和其它的特征和优点将会变得更清楚，其中：
图1说明了一个典型的光体积描记(photo-plethysmographic)(PPG)信号的曲线图；
图2是根据本发明的一个实施例的血液成分光谱分析系统的方框图；
图3是根据本发明的一个实施例中的剔除反常数据的方法的流程图；以及
图4展示了第一波长和第二波长的光强的曲线图的例子，其中详细注明了图3中所示的方法的每一阶段中使用的参数。
具体实施方式
以下将根据本发明的优选实施例，结合附图详细描述一种剔除反常数据的方法，以及应用了该方法的一种血液成分光谱学分析系统。
图2是根据本发明一个实施例中的血液成分光谱分析系统方框图。该血液成分光谱分析系统包括光发射装置21、受实验者22、光接收装置23、放大和过滤单元24、信号处理器25、存储单元26、以及显示单元27。
参见图2光发射装置21至少包含两个光源，例如发光二极管，其中每个都能发射一种能与特定血液成起作用的波长的光。根据信号处理器25控制的时间，每个光源都向受实验者22发射特定波长的光。
光接收装置23接收从放置在光发射装置21和光接收装置23之间的受实验者22传播、散射或者反射来的光，将接收光转化为电信号，并将该电信号提供给放大和过滤单元24。放大和过滤单元24将该电信号放大至预定的标准，然后过滤该放大的电信号使除去该电信号中的噪声成分。
信号处理器25从放大和过滤单元24提供的电信号中提取与特定血液成分作用的光体积描记(PPG)信号，将该PPG信号转换成数字数据，进行数据分析并对该数字数据进行处理以便除去其中的反常数据，然后将信号处理得到的信号提供给存储单元26和显示单元27。信号处理器25包含可机读的记录介质，其记录了实现根据本发明中的剔除反常数据的方法的程序。该信号处理器25优选含有一个具有先入先出(FIFO)结构的数据缓冲装置以储存“n”个数据，其中“n”为正整数。这里，当“n”个数据作为信号分析和处理的结果被确定为应被剔除时，“n”个数据就被划分到预定量的组中，例如，有五个组时，只有最先输入的n/5个数据能得到更新并且应用到剔除反常数据的信号分析中。
存储单元26储存那些如信号处理器25处理结果那样被确定为有效的数字数据，而显示单元27显示出信号处理器25处理的结果以便用户得知。
图3是根据本发明一个实施例中的剔除反常数据的方法流程图。该方法包括采集数据(31)，计算波长之间参数的平均值(32至35)，以及利用该参数的平均值、标准偏差和预定的剔除参考值来确定该参数是否有效(36至38)。图3所示的方法将结合图4来描述。
参见图3，步骤31中，在预定的单位时间周期内，例如T1到T2，由放大和过滤单元24提供的PPG信号的“n”个脉冲数据根据它们的波长进行采集。“n”个脉冲数据通常可包含一个心跳周期，并且可以由调制信号产生，该调制信号由于外部压力引起的血流变化而具有较大的振幅。
在步骤32中，预定的单位时间周期内每一脉冲的光强的最大值点Ip和最小值点Iv都是从根据它们波长采集的脉冲数据组中获得的。关于第一波长和第二波长的光强曲线G1和G2，当从T1到T2的时间内出现“n”个脉冲时，在关于第一波长的曲线G1上可得到“n”个脉冲中的最大值点Ip11至Ip1n和最小值点Iv11至Iv1n，在关于第二波长的曲线G2上可获得“n”个脉冲中的最大值点Ip21至Ip2n和从最小值点Iv21至Iv2n。
在步骤33中，计算出每种波长光的每一脉冲的最大值点Ip和最小值点Iv的比率对数ln(Ip/Iv)。也就是说，在第一波长光的曲线G1上，计算出第一到第n个脉冲的比率对数为ln(Ip11/Iv11)到ln(Ip1n/I1n)。在第二波长光的曲线G2上，计算出第一到第n个脉冲的比率对数为ln(Ip21/Iv21)到ln(Ip2n/Iv2n)。
在步骤34中，每种波长的光得到的比率对数ln(Ip/Iv)被用来计算参数，例如波长之间比率的比值(ROR)，即比率R12为第一波长光的每一个脉冲的比率对数ln(Ip1k/Iv1k)(其中k为1到n)与第二波长光的每一个脉冲的比率对数ln(Ip2k/Iv2k)(其中k为1到n)的比值。此第一和第二波长光间的ROR值R12可用公式(1)表示。
R 12 = ln ( I p 1 k / I v 1 k ) ln ( I p 2 k / I v 2 k ) · · · ( 1 )
在步骤35中，第一和第二波长间光的ROR值R12的平均值M12可根据公式(2)计算。
M 12 = Σ k = 1 n ln ( I p 1 k / I v 1 k ) ln ( I P 2 k / I v 2 k ) n · · · ( 2 )
在步骤36中，在由步骤34计算出的“n”个参数中，计算与步骤35中计算出的平均值M12偏差大于标准偏差的参数数目，即参数“m”。这里，标准偏差可以通过实验或模拟而设置一个最佳值。例如，标准偏差可依据用于获得不同波长之间ROR值的两种波长类型而设置不同的值。当两种波长对应于等消光点时，标准偏差设置为±3.5％至±4.5％，优选设置为±4％。当两种波长在红光或红外范围时，标准偏差应设置为±1.5％至±2.5％，优选设置为±2％。
在步骤37中，从步骤34得出的所有参数数目“n”与从步骤36得出的平均偏差大于标准偏差的参数数目“m”的比值，即m/n，与预设的剔除参考值相比较。这里，为保证分析的精确度，剔除参考值可以通过实验或模拟而设置一个最佳值。步骤37的比较结果中，如果比值m/n大于或等于预设的剔除参考值，当前采集的“n”个数据将被丢弃并不用于分析，然后该方法退回到步骤31。这里步骤31中所有的“n”个数据将更新为新的数据。或者，“n”个数据可以被划入预定的单元组中，相同数量的新数据将代替由前次输入数据组成的初始预设单元组，并且成为初始预设单元组。
如果m/n小于预设的剔除参考值，步骤34中得到的第一和第二波长光之间“n”个ROR值R12就被确定为有效数据并被储存到存储单元26之中。储存在存储单元26中有效数据的平均值，即第一和第二波长光之间的RORs值R12的平均值，可被运用到例如通过多元线性回归分析确定的回归方程中以便计算特定血液成分的浓度。
例如，当计算血红蛋白量(Hb)的回归方程用公式(3)表示时，在步骤31至38中确定为有效的第一和第二波长光之间的RORs平均值举例如下表1：
Hb＝50.3-21.5×R13-1.30×R15-9.00×R45                 …(3)
                            表1
  脉冲索引   R13   R15   R45
  1   1.11935   1.56065   1.37540
  2   1.13040   1.51785   1.28079
  3   1.11839   1.55493   1.32513
  4   1.15221   1.51489   1.29946
  5   1.16072   1.57932   1.33420
  6   1.15716   1.49149   1.25697
  7   1.19721   1.52200   1.26801
  8   1.09604   1.45250   1.30688
  9   1.08103   1.56253   1.32151
  10   1.11863   1.55737   1.32719
  11   1.11456   1.50028   1.31471
  12   1.17718   1.59165   1.36132
  13   1.13893   1.59417   1.40585
  14   1.11350   1.53898   1.36942
  15   1.11424   1.52678   1.30636
  16   1.18262   1.60972   1.34542
  17   1.18803   1.66496   1.45322
  18   1.17468   1.59366   1.36323
  19   1.17791   1.63599   1.38658
  20   1.12564   1.63940   1.38449
  21   1.15660   1.59757   1.38367
  22   1.12855   1.58233   1.35407
  23   1.11076   1.62814   1.42659
  24   1.19795   1.59025   1.35692
  25   1.16404   1.59834   1.37347
  26   1.12490   1.55290   1.28225
  27   1.12226   1.50934   1.34515
  28   1.16658   1.55562   1.32509
  29   1.13476   1.58237   1.35572
  30   1.10259   1.62487   1.39802
  平均值   1.14158   1.56769   1.34624
这里，预定单位时间周期内包含的脉冲数“n”为30，R13表示569nm波长和805nm波长光之间的ROR值，R15表示569nm波长和970nm波长光之间的ROR值，R45表示940nm波长和970nm波长光之间的ROR值。血红蛋白浓度可通过将两种波长光之间RORs的平均值应用于R13，R15或R45中估测出来，根据本发明此RORs值被确定为有效。
本发明可以通过记录在可机读的记录介质上并能被计算机识别的编码来实现。例如，实现剔除反常数据的方法可以通过以下方式实现：在可机读的记录介质上记录第一程序，用于计算出与第一和第二波长光相对应的PPG信号中包含的有关“n”个脉冲的“n”个参数，其中PPG信号是在预定的单位时间周期内采集到的；以及记录第二程序，用来计算“n”个参数的平均值，并将对平均值的偏差大于所预定的标准偏差的参数数目与“n”个参数的比率与预定的剔除参考值相比较，以便确定“n”个脉冲数据是否有效。同时，可机读的记录介质可以是能够在其上记录计算机可读数据的任何类型的介质，例如，ROM、RAM、CD-ROM、磁带、软盘或者光学数据存储装置。本发明还可以通过载波(例如通过互联网传播)实现。或者，可机读的记录介质分配到通过网络连接的计算机系统中，这样本发明就能作为存储在记录介质上并能被计算机读取并执行的编码来实现。本发明领域的程序员很容易推导出可实现本发明的功能程序、编码和编码片段。
<实验>
根据本发明的一个实施例，为了获得剔除反常数据方法的定量结果，要连续10分钟从30个血红蛋白浓度在10.5-16.5g/dl的仰卧受试者身上采集数据。30个被试者以2比1的比例分组，通过对作为标准模型的20个受试者试样作多元线性回归分析获得回归方程。并将该回归方程应用于从余下10个受试者采集的数据以计算出结果作为预测模型。作为参考的所有受试者的血红蛋白值是用HemoCue AB(瑞典)来测量的，它是一种用小试管作非损伤性测量的装置。
为了确定参数，即不同波长光之间的RORs值是否是反常数据，将有关569nm波长光的对数比与有关805nm波长光的对数比之比值R13的标准偏差，和有关940nm波长光的对数比与有关970nm波长光的对数比之比值R45的标准偏差分别增加2％和1％。剔除参考值设置为20％，用于测定的单位时间周期设置为1分钟。由于该标准偏差小，所以相关系数就大，但是满足标准偏差的试样数目N减少了，能够用来分析的数据就不多了。相关系数与满足标准偏差的试样数目N成反比。因此，R13和R45的标准偏差最好分别设为4％和2％。表2给出了有关R13和R45的标准模型和预测模型中的相关系数和被分析试样数目N。
                                  表2
  标准模型             校准模型             预测模型
  R13(％)   相关系   标准偏差   试样的   相关系   标准偏差   试样的
  R45(％)   数(R)   (g/dl)  数目(N)   数(R)   (g/dl)  数目(N)
  2   1   0.71   1.2   15   0.65   1.4   6
  4   2   0.71   1.3   20   0.61   1.7   10
  6   3   0.59   1.5   20   0.39   2.05   10
  8   4   0.53   1.82   20   0.3   2.34   10
  没有应用本发明   0.51   1.9   20   0.3   2.4   10
参见表2，根据本发明，当将剔除反常数据的方法应用于30个被试者时，该方法用最佳标准偏差和预定的剔除参考值确定所采集的数据是否能用于非损伤性测量血液中血红蛋白浓度应用到，从标准模型中获得的相关系数为0.71，预测模型中获得的相关系数为0.61。与没有用本发明的方法所得到的分析结果相比，其中标准模型中获得的相关系数为0.51，预测模型中获得的相关系数为0.3，本发明改善了分析结果。因此证明了本发明有效地减小了由于内部生命要素如(呼吸、血压、脉搏速度、体温、血管状态或是植物性神经系统)和外部生命要素(如运动引起的噪声或与探针的不稳定接触)，导致的脉动成分生理信号变化的影响。
正如上面所描述的，根据本发明，利用光谱学技术的血液成分分析系统从被试者产生的PPG信号中计算回归方程的参数，在预定单位时间周期内分析参数的反常数据，并在分析中排除包含那些比率大于或等于预定的剔除参考值的反常数据的数据组，从而减小由于内部生命要素使脉动成分变化而造成的影响，以及外部生命要素造成的影响，并且无需单独增加硬件系统如运动传感器和加速计。结果能在无需增加血液成分分析系统生产成本的情况下提高分析的精确度。
另外，根据本发明，对于普通的生理信号测试，用于使病人达到稳定状态所需的准备时间消除或者减少了，这样用于分析的总时间就能减少。而且，与传统的由试验人员主观决定的分析相比，剔除或者不剔除反常数据可以自动和客观地决定，这样可以达到稳定的分析精确度。
尽管这里展示和记载了本发明的一些实施例，本领域技术人员应该知道，可以在不背离本发明原理和实质的情况下对这些元件做出改变，这些范围在所附的权利要求及其等价物中做出了限定。

Claims (14)

1.一种在利用光谱学技术的血液成分分析系统中剔除反常数据的方法，其通过对通过将对血液成分作出反应的第一和第二波长的光发射到试样上和对穿过试样的所述第一和第二波长的光进行光电转换所获得的光体积描记信号进行分析，来估测血液成分的浓度，所述方法包括：
(a)在预定的单位时间周期内采集与第一和第二波长相对应的PPG信号；
(b)计算“n”个参数，其中“n”为所采集的PPG信号中包含的脉冲数据的数目，这里n为正整数；
(c)计算“n”个参数的平均值；以及
(d)将对平均值的偏差大于所预定的标准偏差的参数的数目与“n”个参数的比率与预定的剔除参考值相比较，以便确定“n”个脉冲数据是否有效，
其中，所述参数为第一波长光的每个脉冲的最大值点和最小值点的比率对数与第二波长光的每个脉冲的最大值点和最小值点的比率对数的比值。
2.根据权利要求1所述的方法，还包括：更新全部“n”个脉冲数据，并且当作为步骤(d)中的比较结果的比率等于或大于预定的剔除参考值时，重复步骤(b)至(d)。
3.根据权利要求1所述的方法，还包括：更新“n”个脉冲数据中的一些，并且当作为步骤(d)中的比较结果的比率等于或大于预定的剔除参考值时，重复步骤(b)至(d)。
4.根据权利要求1所述的方法，其中所述标准偏差是根据第一和第二波长来确定。
5.根据权利要求4所述的方法，其中当所述第一和第二波长与等消光点相对应时，将所述标准偏差设置在±3.5％至±4.5％。
6.根据权利要求4所述的方法，其中当第一和第二波长在红光范围和红外光范围之一中时，将所述标准偏差设置在±1.5％至±2.5％。
7.根据权利要求1所述的方法，其中从由外部压力引起的血流中的变化所影响的调制信号来产生所述“n”个脉冲数据。
8.一种在利用光谱学技术估测血液成分浓度的血液成分光谱分析系统中使用的计算机可读记录介质，所述估测血液成分浓度是通过分析通过将对血液成分作出反应的第一和第二波长的光发射到试样上和对穿过试样的所述第一和第二波长的光进行光电转换所获得的光体积描记信号来进行的，所述计算机可读记录介质包括：
第一程序，用来对应于与第一和第二波长相对应的PPG信号中包含的“n”个脉冲来计算“n”个参数，其中PPG信号是在预定的单位时间周期内采集到的，n为正整数；以及
第二程序，用来计算“n”个参数的平均值，并将对其平均值的偏差大于预定的标准偏差的参数的数目与“n”个参数的比率与预定的剔除参考值相比较，以便确定“n”个脉冲数据是否有效，
其中，所述参数为第一波长光的每个脉冲的最大值点和最小值点的比率对数与第二波长光的每个脉冲的最大值点和最小值点的比率对数的比值。
9.一种利用光谱技术的血液成分分析系统，包括：
光发射装置，其包含至少两个能发射对特定的血液成分作出反应的波长的光的光源，并且能根据预定的定时向试样发射特定波长的光；
光接收装置，用于接收从试样发出的光并将所接收的光转换成电信号；以及
信号处理器，用于从由光接收装置所接收的电信号中提取与对特定的血液成分作出反应的特定的波长的光的体积描记信号，对于在预定单元时间周期内采集的与所述波长相对应的PPG信号中包括的“n”个脉冲，来计算“n”个参数的平均值，以及将对其平均值的偏差大于预定标准偏差的参数的数目与“n”个参数的比率与预定的剔除参考值相比较，以便确定“n”个脉冲数据是否有效，这里n为正整数，
并且其中，所述参数为第一波长光的每个脉冲的最大值点和最小值点的比率对数与第二波长光的每个脉冲的最大值点和最小值点的比率对数的比值。
10.根据权利要求9所述的血液成分分析系统，其中所述信号处理器包括用于存储“n”个脉冲数据的先入先出数据缓冲器，和当所述“n”个脉冲数据被确定为无效时，将所述“n”个脉冲数据划分成预定的单元组，用与第一预定单元组相同数目的新数据来取代由初始输入数据组成的所述第一预定单元组，和确定新的“n”个脉冲数据是否有效。
11.根据权利要求9所述的血液成分分析系统，其中所述标准偏差是根据第一和第二波长来确定。
12.根据权利要求11所述的血液成分分析系统，其中当所述第一和第二波长与等消光点相对应时，将所述标准偏差设置在±3.5％至±4.5％。
13.根据权利要求11所述的血液成分分析系统，其中当第一和第二波长是在红光范围或红外光范围中时，将所述标准偏差设置在±1.5％至±2.5％。
14.根据权利要求9所述的血液成分分析系统，其中从由外部压力引起的血流中的变化所影响的调制信号中来产生所述“n”个脉冲数据。
CNB2004100283223A 2003-01-07 2004-01-07 剔除反常数据方法和应用该方法的血液成分光谱分析系统 Expired - Fee Related CN1264474C (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
KR10-2003-0000779A KR100499139B1 (ko) 2003-01-07 2003-01-07 이상데이터 소거방법 및 이를 적용한 분광학을 이용한혈액성분분석시스템
KR779/03 2003-01-07
KR779/2003 2003-01-07

Publications (2)

Publication Number Publication Date
CN1519557A CN1519557A (zh) 2004-08-11
CN1264474C true CN1264474C (zh) 2006-07-19

Family

ID=36217582

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB2004100283223A Expired - Fee Related CN1264474C (zh) 2003-01-07 2004-01-07 剔除反常数据方法和应用该方法的血液成分光谱分析系统

Country Status (6)

Country Link
US (1) US8507281B2 (zh)
EP (1) EP1437087B1 (zh)
JP (1) JP4050706B2 (zh)
KR (1) KR100499139B1 (zh)
CN (1) CN1264474C (zh)
DE (1) DE602004000513T2 (zh)

Families Citing this family (81)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6675031B1 (en) 1999-04-14 2004-01-06 Mallinckrodt Inc. Method and circuit for indicating quality and accuracy of physiological measurements
US7190986B1 (en) 2002-10-18 2007-03-13 Nellcor Puritan Bennett Inc. Non-adhesive oximeter sensor for sensitive skin
US7006856B2 (en) 2003-01-10 2006-02-28 Nellcor Puritan Bennett Incorporated Signal quality metrics design for qualifying data for a physiological monitor
US7194293B2 (en) 2004-03-08 2007-03-20 Nellcor Puritan Bennett Incorporated Selection of ensemble averaging weights for a pulse oximeter based on signal quality metrics
US7590439B2 (en) 2005-08-08 2009-09-15 Nellcor Puritan Bennett Llc Bi-stable medical sensor and technique for using the same
US7657294B2 (en) 2005-08-08 2010-02-02 Nellcor Puritan Bennett Llc Compliant diaphragm medical sensor and technique for using the same
US7657295B2 (en) 2005-08-08 2010-02-02 Nellcor Puritan Bennett Llc Medical sensor and technique for using the same
US20070060808A1 (en) 2005-09-12 2007-03-15 Carine Hoarau Medical sensor for reducing motion artifacts and technique for using the same
US8092379B2 (en) 2005-09-29 2012-01-10 Nellcor Puritan Bennett Llc Method and system for determining when to reposition a physiological sensor
US7899510B2 (en) 2005-09-29 2011-03-01 Nellcor Puritan Bennett Llc Medical sensor and technique for using the same
US7904130B2 (en) 2005-09-29 2011-03-08 Nellcor Puritan Bennett Llc Medical sensor and technique for using the same
US7869850B2 (en) 2005-09-29 2011-01-11 Nellcor Puritan Bennett Llc Medical sensor for reducing motion artifacts and technique for using the same
US7483731B2 (en) 2005-09-30 2009-01-27 Nellcor Puritan Bennett Llc Medical sensor and technique for using the same
US8062221B2 (en) 2005-09-30 2011-11-22 Nellcor Puritan Bennett Llc Sensor for tissue gas detection and technique for using the same
US7555327B2 (en) 2005-09-30 2009-06-30 Nellcor Puritan Bennett Llc Folding medical sensor and technique for using the same
US8233954B2 (en) 2005-09-30 2012-07-31 Nellcor Puritan Bennett Llc Mucosal sensor for the assessment of tissue and blood constituents and technique for using the same
US7881762B2 (en) 2005-09-30 2011-02-01 Nellcor Puritan Bennett Llc Clip-style medical sensor and technique for using the same
US7486979B2 (en) 2005-09-30 2009-02-03 Nellcor Puritan Bennett Llc Optically aligned pulse oximetry sensor and technique for using the same
US8073518B2 (en) 2006-05-02 2011-12-06 Nellcor Puritan Bennett Llc Clip-style medical sensor and technique for using the same
US8145288B2 (en) 2006-08-22 2012-03-27 Nellcor Puritan Bennett Llc Medical sensor for reducing signal artifacts and technique for using the same
US8219170B2 (en) 2006-09-20 2012-07-10 Nellcor Puritan Bennett Llc System and method for practicing spectrophotometry using light emitting nanostructure devices
US8175671B2 (en) 2006-09-22 2012-05-08 Nellcor Puritan Bennett Llc Medical sensor for reducing signal artifacts and technique for using the same
US8396527B2 (en) 2006-09-22 2013-03-12 Covidien Lp Medical sensor for reducing signal artifacts and technique for using the same
US8195264B2 (en) 2006-09-22 2012-06-05 Nellcor Puritan Bennett Llc Medical sensor for reducing signal artifacts and technique for using the same
US7869849B2 (en) 2006-09-26 2011-01-11 Nellcor Puritan Bennett Llc Opaque, electrically nonconductive region on a medical sensor
US7574245B2 (en) 2006-09-27 2009-08-11 Nellcor Puritan Bennett Llc Flexible medical sensor enclosure
US7890153B2 (en) 2006-09-28 2011-02-15 Nellcor Puritan Bennett Llc System and method for mitigating interference in pulse oximetry
US7796403B2 (en) 2006-09-28 2010-09-14 Nellcor Puritan Bennett Llc Means for mechanical registration and mechanical-electrical coupling of a faraday shield to a photodetector and an electrical circuit
US7684842B2 (en) 2006-09-29 2010-03-23 Nellcor Puritan Bennett Llc System and method for preventing sensor misuse
US7476131B2 (en) 2006-09-29 2009-01-13 Nellcor Puritan Bennett Llc Device for reducing crosstalk
US8068891B2 (en) 2006-09-29 2011-11-29 Nellcor Puritan Bennett Llc Symmetric LED array for pulse oximetry
US8175667B2 (en) 2006-09-29 2012-05-08 Nellcor Puritan Bennett Llc Symmetric LED array for pulse oximetry
US7680522B2 (en) 2006-09-29 2010-03-16 Nellcor Puritan Bennett Llc Method and apparatus for detecting misapplied sensors
US8280469B2 (en) 2007-03-09 2012-10-02 Nellcor Puritan Bennett Llc Method for detection of aberrant tissue spectra
US7894869B2 (en) 2007-03-09 2011-02-22 Nellcor Puritan Bennett Llc Multiple configuration medical sensor and technique for using the same
US8265724B2 (en) 2007-03-09 2012-09-11 Nellcor Puritan Bennett Llc Cancellation of light shunting
US8352004B2 (en) 2007-12-21 2013-01-08 Covidien Lp Medical sensor and technique for using the same
US8346328B2 (en) 2007-12-21 2013-01-01 Covidien Lp Medical sensor and technique for using the same
US8366613B2 (en) 2007-12-26 2013-02-05 Covidien Lp LED drive circuit for pulse oximetry and method for using same
US8577434B2 (en) 2007-12-27 2013-11-05 Covidien Lp Coaxial LED light sources
US8442608B2 (en) 2007-12-28 2013-05-14 Covidien Lp System and method for estimating physiological parameters by deconvolving artifacts
US8452364B2 (en) 2007-12-28 2013-05-28 Covidien LLP System and method for attaching a sensor to a patient's skin
US8092993B2 (en) 2007-12-31 2012-01-10 Nellcor Puritan Bennett Llc Hydrogel thin film for use as a biosensor
US8897850B2 (en) 2007-12-31 2014-11-25 Covidien Lp Sensor with integrated living hinge and spring
US8070508B2 (en) 2007-12-31 2011-12-06 Nellcor Puritan Bennett Llc Method and apparatus for aligning and securing a cable strain relief
US8199007B2 (en) 2007-12-31 2012-06-12 Nellcor Puritan Bennett Llc Flex circuit snap track for a biometric sensor
WO2009120600A2 (en) 2008-03-25 2009-10-01 The Curators Of The University Of Missouri Method and system for non-invasive blood glucose detection utilizing spectral data of one or more components other than glucose
US8437822B2 (en) 2008-03-28 2013-05-07 Covidien Lp System and method for estimating blood analyte concentration
US8112375B2 (en) 2008-03-31 2012-02-07 Nellcor Puritan Bennett Llc Wavelength selection and outlier detection in reduced rank linear models
RU2566920C2 (ru) 2008-05-22 2015-10-27 Дзе Кьюрейторз Оф Дзе Юниверсити Оф Миссури Способ и система для неинвазивного оптического определения глюкозы крови, используя спектральный анализ данных
US8071935B2 (en) 2008-06-30 2011-12-06 Nellcor Puritan Bennett Llc Optical detector with an overmolded faraday shield
US7880884B2 (en) 2008-06-30 2011-02-01 Nellcor Puritan Bennett Llc System and method for coating and shielding electronic sensor components
US7887345B2 (en) 2008-06-30 2011-02-15 Nellcor Puritan Bennett Llc Single use connector for pulse oximetry sensors
US8364220B2 (en) 2008-09-25 2013-01-29 Covidien Lp Medical sensor and technique for using the same
US8417309B2 (en) 2008-09-30 2013-04-09 Covidien Lp Medical sensor
US8914088B2 (en) 2008-09-30 2014-12-16 Covidien Lp Medical sensor and technique for using the same
US8423112B2 (en) 2008-09-30 2013-04-16 Covidien Lp Medical sensor and technique for using the same
US8452366B2 (en) 2009-03-16 2013-05-28 Covidien Lp Medical monitoring device with flexible circuitry
US8221319B2 (en) 2009-03-25 2012-07-17 Nellcor Puritan Bennett Llc Medical device for assessing intravascular blood volume and technique for using the same
US8509869B2 (en) 2009-05-15 2013-08-13 Covidien Lp Method and apparatus for detecting and analyzing variations in a physiologic parameter
US8634891B2 (en) 2009-05-20 2014-01-21 Covidien Lp Method and system for self regulation of sensor component contact pressure
US8505821B2 (en) 2009-06-30 2013-08-13 Covidien Lp System and method for providing sensor quality assurance
US9010634B2 (en) 2009-06-30 2015-04-21 Covidien Lp System and method for linking patient data to a patient and providing sensor quality assurance
US8311601B2 (en) 2009-06-30 2012-11-13 Nellcor Puritan Bennett Llc Reflectance and/or transmissive pulse oximeter
US8391941B2 (en) 2009-07-17 2013-03-05 Covidien Lp System and method for memory switching for multiple configuration medical sensor
US8417310B2 (en) 2009-08-10 2013-04-09 Covidien Lp Digital switching in multi-site sensor
US8428675B2 (en) 2009-08-19 2013-04-23 Covidien Lp Nanofiber adhesives used in medical devices
CN102818769B (zh) * 2012-09-05 2015-04-22 江苏中烟工业有限责任公司 一种评价成品卷烟烟支内烟丝混合均匀性的方法
CN103018177A (zh) * 2012-12-06 2013-04-03 江苏易谱恒科技有限公司 基于随机采样一致集的谱图异常样本点检测方法
CN106999112A (zh) * 2014-10-10 2017-08-01 麦德托有限公司 用于无创医疗传感器的系统和方法
CN104535785A (zh) * 2014-12-24 2015-04-22 牛艳花 一种测量数据的降噪方法
JP6454211B2 (ja) * 2015-03-31 2019-01-16 シスメックス株式会社 検体分析装置、血液凝固分析装置、検体分析方法、及びコンピュータプログラム
KR20170008043A (ko) * 2015-07-13 2017-01-23 엘지전자 주식회사 이동 단말기의 심박/스트레스 측정회로 및 그 측정 방법
US10537253B2 (en) 2016-02-25 2020-01-21 Samsung Electronics Company, Ltd. Detecting live tissues using signal analysis
US10743821B2 (en) 2016-10-21 2020-08-18 Tata Consultancy Services Limited Anomaly detection by self-learning of sensor signals
CN110208248B (zh) * 2019-06-28 2021-11-19 南京林业大学 一种辨识拉曼光谱异常测量信号的方法
JP2021027933A (ja) * 2019-08-09 2021-02-25 Toto株式会社 生体情報取得装置
CN110547768B (zh) * 2019-08-30 2020-07-28 北京师范大学 一种近红外脑功能成像质量控制方法和控制系统
CN110974250B (zh) * 2019-12-27 2024-01-16 深圳市华讯方舟光电技术有限公司 基于太赫兹光谱的血糖检测方法、装置及计算机存储介质
CN111603151B (zh) * 2020-06-17 2023-05-16 深圳智领人工智能健康科技有限公司 一种基于时频联合分析的无创血液成分检测方法及系统
CN113080866A (zh) * 2021-04-09 2021-07-09 上海中医药大学附属曙光医院 一种针对实验动物的舌象采集装置和方法

Family Cites Families (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6334731A (ja) 1986-07-29 1988-02-15 Canon Inc 磁気記録媒体の製造方法
US4800495A (en) 1986-08-18 1989-01-24 Physio-Control Corporation Method and apparatus for processing signals used in oximetry
JPH0315502U (zh) 1989-06-28 1991-02-15
US5285782A (en) * 1992-01-17 1994-02-15 Physio-Control Corporation Method and apparatus for improving the accuracy of pulse transmittance oximeter
EP0658331B1 (en) 1993-12-11 1996-10-02 Hewlett-Packard GmbH A method for detecting an irregular state in a non-invasive pulse oximeter system
US5553615A (en) * 1994-01-31 1996-09-10 Minnesota Mining And Manufacturing Company Method and apparatus for noninvasive prediction of hematocrit
US5645060A (en) 1995-06-14 1997-07-08 Nellcor Puritan Bennett Incorporated Method and apparatus for removing artifact and noise from pulse oximetry
JP3420409B2 (ja) 1995-10-18 2003-06-23 セイコーエプソン株式会社 計測装置
US5842979A (en) * 1997-02-14 1998-12-01 Ohmeda Inc. Method and apparatus for improved photoplethysmographic monitoring of oxyhemoglobin, deoxyhemoglobin, carboxyhemoglobin and methemoglobin
US5983122A (en) * 1997-12-12 1999-11-09 Ohmeda Inc. Apparatus and method for improved photoplethysmographic monitoring of multiple hemoglobin species using emitters having optimized center wavelengths
US6393311B1 (en) * 1998-10-15 2002-05-21 Ntc Technology Inc. Method, apparatus and system for removing motion artifacts from measurements of bodily parameters
JP3835073B2 (ja) 1999-08-30 2006-10-18 株式会社デンソー 脈波計測装置
US6496723B1 (en) 1999-08-30 2002-12-17 Denso Corporation Method of obtaining information that corresponds to electrocardiogram of human body from pulse wave thereof
JP3415796B2 (ja) 2000-01-13 2003-06-09 日本コーリン株式会社 周期性生体情報測定装置

Also Published As

Publication number Publication date
CN1519557A (zh) 2004-08-11
DE602004000513T2 (de) 2007-03-15
US8507281B2 (en) 2013-08-13
JP2004209257A (ja) 2004-07-29
EP1437087B1 (en) 2006-03-22
DE602004000513D1 (de) 2006-05-11
EP1437087A1 (en) 2004-07-14
JP4050706B2 (ja) 2008-02-20
KR20040063383A (ko) 2004-07-14
US20040204865A1 (en) 2004-10-14
KR100499139B1 (ko) 2005-07-04

Similar Documents

Publication Publication Date Title
CN1264474C (zh) 剔除反常数据方法和应用该方法的血液成分光谱分析系统
US8073516B2 (en) Separating motion from cardiac signals using second order derivative of the photo-plethysmogram and fast fourier transforms
US9949675B2 (en) Noninvasive blood measurement platform
US11872060B2 (en) Methods and systems for calculating physiological parameters
EP0335357B1 (en) Improved method and apparatus for detecting optical pulses
US6094592A (en) Methods and apparatus for estimating a physiological parameter using transforms
US8442608B2 (en) System and method for estimating physiological parameters by deconvolving artifacts
US8478538B2 (en) Selection of signal regions for parameter extraction
US8594759B2 (en) Systems and methods for resolving the continuous wavelet transform of a signal
CN100515335C (zh) 能消除运动干扰的血氧测量方法及其装置
US20140378795A1 (en) Determination of a physiological parameter
US20130197328A1 (en) Signal processing apparatus and method
US20090326831A1 (en) Concatenated Scalograms
US20090326351A1 (en) Signal Processing Mirroring Technique
CN106137219B (zh) 双波长的绝对差值加和计算动脉血氧饱和度方法及其装置
US20120310051A1 (en) Systems And Methods For Signal Rephasing Using The Wavelet Transform
CN104224197A (zh) 双波长光电容积脉搏波特征量的动脉血氧饱和度计算方法
CA2728609A1 (en) Signal processing systems and methods for determining slope using an origin point
He et al. Spectral data quality assessment based on variability analysis: application to noninvasive hemoglobin measurement by dynamic spectrum
CN100589758C (zh) 交流分量的测量方法及测量装置
US20140187884A1 (en) Systems and methods for ensemble averaging in pulse oximetry
CN117017297A (zh) 驾驶员疲劳的预测和识别模型建立方法及其应用
US20120253147A1 (en) Calibration method and arrangement and sensor for non-invasively measuring blood characteristics of a subject
Habbu et al. Comparison of noninvasive blood glucose estimation using various regression models
Zhang et al. Non-invasive blood glucose detection using NIR based on GA and SVR

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
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20060719

Termination date: 20100208