CN109684908A - 一种用于血栓弹力图仪的信号滤波方法 - Google Patents

一种用于血栓弹力图仪的信号滤波方法 Download PDF

Info

Publication number
CN109684908A
CN109684908A CN201811109037.2A CN201811109037A CN109684908A CN 109684908 A CN109684908 A CN 109684908A CN 201811109037 A CN201811109037 A CN 201811109037A CN 109684908 A CN109684908 A CN 109684908A
Authority
CN
China
Prior art keywords
signal
digital signal
value
filtering method
elastic force
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
CN201811109037.2A
Other languages
English (en)
Other versions
CN109684908B (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.)
Shenzhen Ward Life Technology Co Ltd
Original Assignee
Shenzhen Ward Life 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 Shenzhen Ward Life Technology Co Ltd filed Critical Shenzhen Ward Life Technology Co Ltd
Priority to CN201811109037.2A priority Critical patent/CN109684908B/zh
Publication of CN109684908A publication Critical patent/CN109684908A/zh
Application granted granted Critical
Publication of CN109684908B publication Critical patent/CN109684908B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/02—Preprocessing
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H40/00—ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Business, Economics & Management (AREA)
  • Biomedical Technology (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Signal Processing (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Artificial Intelligence (AREA)
  • Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Epidemiology (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Primary Health Care (AREA)
  • Public Health (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
  • Investigating Or Analysing Biological Materials (AREA)

Abstract

本发明提出一种用于血栓弹力图仪的信号滤波方法，通过本发明的用于血栓弹力图仪的信号滤波方法，可以消除采样噪声信号的干扰、缓变信号中的尖脉冲干扰、大幅度短时间非周期信号干扰，且可以屏蔽信号采集探针高频振动带来的信号震荡，消除长周期大脉冲干扰，达到数据筛选的目的，且在整个凝血纤溶周期，能够很好消除噪声，使得输出的检测数据更加准确。

Description

一种用于血栓弹力图仪的信号滤波方法
技术领域
本发明涉及信号处理领域，尤其一种用于血栓弹力图仪的信号滤波方法。
背景技术
血栓弹力图是通过动力学传感器检测技术，采集血液在凝固到纤溶工程中的动力学变化信号，通过数据处理计算出参数结果，用于病人的凝血特性评估，应用于LVD血液检测领域的凝血分析，凝血特性评估常用于临床评价，比如术后出血、心外手术中或术后的血栓、器官移植、外伤和心脏手术过程中或术后发生的血栓症。
目前美国血液公司TEG弹力图仪能够监测病人的凝血状况，主要基于以下两个事实：(1)血液凝固过程的最终结果是形成血凝块；(2)血凝块的物理性质(速率、硬度、稳定性)将决定病人是否具有正常的凝血功能，是否会出血或形成血栓。Teg数据分析软件对采样的动力学信号进行滤波识别，数据筛选，输出参数结果。
然而，上述方法却存在以下问题：(1)Teg数据分析软件对通过动力学信号进行了单一简单滤波，仪器操作工程中往往会引入机械振动，机械振动产生的噪声信号有变化的频率，滤波器根本无法消除噪声信号；(2)数据筛选算法只能在一定程度上消除小振动信号的干扰，且在凝血初期，对凝血快速变化的信号的筛选效果不是很明显。
发明内容
为了解决上述问题，本发明提出一种用于血栓弹力图仪的信号滤波方法，能够消除噪声信号和小振动信号的干扰，使得弹力图输出的数据更加准确。
本发明通过以下技术方案实现的：
本发明提出一种用于血栓弹力图仪的信号滤波方法，包括以下步骤：
(1)利用动力学信号采样获得动力学数字信号S1；
(2)将动力学数字信号S1进行均值滤波，获取数字信号S2；
(3)将数字信号S2进行限幅滤波，输出数字信号S3；
(4)对数字信号S3进行中值滤波，获得数字信号S4；
(5)对数字信号S4进行低通滤波处理，获得数字信号S5；
(6)对数字信号S5进行识别处理，计算出多个稳定的极大值点和多个稳定的极小值点；
(7)把所有极大值点组合，生成弹力图上包络线，把所有极小值点组合，生成弹力图下包络线，得到弹力图G1；
(8)通过高阶多项式拟合技术，原始数据拟合窗口的大小依据弹力图G1的幅度变化率进行自适应调整，消除凝血快速期出现的尖峰干扰，得到平滑后的弹力图G2；
(9)对平滑后的弹力图G2进行计算，可以输出相应的参数结果。
优选地，其中，步骤(2)～步骤(5)可以任意调整步骤的执行顺序。
优选地，在步骤(1)中，感应探针浸入血液中，血液凝固和纤溶的强度传递给感应探针转变成模拟电信号，经信号放大器放大后，被处理器AD采样成动力学数字信号S1。
优选地，在步骤(2)中，对动力学数字信号S1进行N点均值滤波器消除采样噪声，同时保留了动力学数字信号S1中动力学信号的细节,滤波后得到数字信号S2。
优选地，在步骤(3)中，首先确定N次采样之间最大偏差限幅水平值(THR)，每次检测到新值时进行判断；如果本次与N次采样均值之差小于THR，则本次值有效，如果本次值与前面N次采样值均值大于等于THR，则本次值无效；根据采样值一定窗口大小的前后值设置限幅频率，通过限幅频率来调整限幅水平THR 的大小，从而消除数字信号S2中的小幅尖脉冲干扰。
优选地，在步骤(4)中，设置一个自适应大小窗口N，N的确定可以根据数字信号S3本来的变化快慢来确定大小，通过当前信号的差分量来确定N值，当前信号在窗口的中间，获取N个信号值进行排序，取N个排序后的信号的中间值作为当前信号值，即数字信号S4。
优选地，在步骤(5)中，低通滤波的截止频率设置为10HZ，滤除数字信号 S4中大于10HZ的高频噪声信号。
优选地，在步骤(6)中，对数字信号S5进行正态分布脉冲识别，通过爬坡算法识别到有效脉冲的初始极大值和初始极小值点，把极值点周围一定窗口大小的信号输入高斯正态分布拟合函数，计算出多个稳定的极大值点和多个稳定的极小值点。
本发明的有益效果：
通过本发明的用于血栓弹力图仪的信号滤波方法，可以消除采样噪声信号的干扰、缓变信号中的尖脉冲干扰、大幅度短时间非周期信号干扰，且可以屏蔽信号采集探针高频振动带来的信号震荡，消除长周期大脉冲干扰，达到数据筛选的目的，且在整个凝血纤溶周期，能够很好消除噪声，使得输出的检测数据更加准确。
附图说明
图1为本发明的用于血栓弹力图仪的信号滤波方法的步骤流程示意图；
图2为本发明的用于血栓弹力图仪的信号滤波方法的动力学数字信号的波形图；
图3为本发明的用于血栓弹力图仪的信号滤波方法的滤波后的动力学数字信号的波形图；
图4为本发明的用于血栓弹力图仪的信号滤波方法的C弹力包络线；
图5为本发明的用于血栓弹力图仪的信号滤波方法的数据筛选后的弹力包络线。
具体实施方式
为了更加清楚、完整的说明本发明的技术方案，下面结合附图对本发明作进一步说明。
请参考图1，本发明一种用于血栓弹力图仪的信号滤波方法，包括以下步骤：
(1)利用动力学信号采样获得动力学数字信号S1；
(2)将动力学数字信号S1进行均值滤波，获取数字信号S2；
(3)将数字信号S2进行限幅滤波，输出数字信号S3；
(4)对数字信号S3进行中值滤波，获得数字信号S4；
(5)对数字信号S4进行低通滤波处理，获得数字信号S5；
(6)对数字信号S5进行识别处理，计算出多个稳定的极大值点和多个稳定的极小值点；
(7)把所有极大值点组合，生成弹力图上包络线，把所有极小值点组合，生成弹力图下包络线，得到弹力图G1；
(8)通过高阶多项式拟合技术，原始数据拟合窗口的大小依据弹力图G1的幅度变化率进行自适应调整，消除凝血快速期出现的尖峰干扰，得到平滑后的弹力图G2；
(9)对平滑后的弹力图G2进行计算，可以输出相应的参数结果。
优选地，其中，步骤(2)～步骤(5)可以任意调整步骤的执行顺序。
在本发明中，步骤(2)～步骤(5)可以根据实际需求任意调整步骤的执行顺序，如按照步骤(5)-步骤(2)-步骤(3)-步骤(4)的顺序执行；且也可以选择性删除步骤(2)～步骤(5)中任意一个步骤或者多个步骤，这些改变或变形均在本发明的保护范围内。
优选地，在步骤(1)中，感应探针浸入血液中，血液凝固和纤溶的强度传递给感应探针转变成模拟电信号，经信号放大器放大后，被处理器AD采样成动力学数字信号S1，如图2所示。
优选地，在步骤(2)中，对动力学数字信号S1进行N点均值滤波器消除采样噪声，同时保留了动力学数字信号S1中动力学信号的细节,滤波后得到数字信号S2，如图3所示。
在本实施方式中，其中N的取值根据AD采样频率来选定，均值滤波器采用特有的自适应加权算法，运用模糊数学的思想，对滤波器权值优化，不仅抑制噪声而且很好的保留了动力学信号的细节，观察数字信号S2，相比于原来的数字信号S1，噪声点被消除了。
优选地，在步骤(3)中，首先确定N次采样之间最大偏差限幅水平值(THR)，每次检测到新值时进行判断；如果本次与N次采样均值之差小于THR，则本次值有效，如果本次值与前面N次采样值均值大于等于THR，则本次值无效；根据采样值一定窗口大小的前后值设置限幅频率，通过限幅频率来调整限幅水平THR 的大小，从而消除数字信号S2中的小幅尖脉冲干扰。
在本实施方式中，与传统滤波相比，本发明采用可变的限幅滤波，通过对限幅水平、限幅频率等参数进行优化设计，信号滤波后具有更高的信噪比。
优选地，在步骤(4)中，设置一个自适应大小窗口N，N的确定可以根据数字信号S3本来的变化快慢来确定大小，通过当前信号的差分量来确定N值，当前信号在窗口的中间，获取N个信号值进行排序，取N个排序后的信号的中间值作为当前信号值，即数字信号S4。
在本实施方式中，与传统中值滤波相比，本文采用自适应中值滤波，通过计算噪声信号和有效信号的变化率来计算窗口尺寸N，在信号变化快速的时间段保持信号平滑。
优选地，在步骤(5)中，低通滤波的截止频率设置为10HZ，滤除数字信号 S4中大于10HZ的高频噪声信号。
在本实施方式中，由于数字信号S4的频率在0.1HZ左右，选择低通滤波器，截止频率设置10HZ，滤除大于10HZ的高频噪声信号，数字低通滤波器类型可以选用巴特沃斯低通滤波器，利用双线性变化法，按指标进行频率变化，计算出低通滤波器系数，本发明的滤波器通带衰减特性平坦，衰减幅宽窄，对于0.1HZ 的信号来说，能很精准的滤波，滤波后得到数字信号S5，达到脉冲识别的准入条件。
优选地，在步骤(6)中，对数字信号S5进行正态分布脉冲识别，通过爬坡算法识别到有效脉冲的初始极大值和初始极小值点，把极值点周围一定窗口大小的信号输入高斯正态分布拟合函数，计算出多个稳定的极大值点和多个稳定的极小值点。
在本实施方式，如图4所示，把所有极大值点组合，生成弹力图上包络线，把所有极小值点组合，生成弹力图下包络线，上下两条包络线组成血栓弹力图 G1曲线。
在本实施方式，由于仪器设备在长期振动干扰下，带来长周期的振动噪声，弹力图G1会出现比较大的尖峰，通过高阶多项式拟合技术，原始数据拟合窗口的大小依据弹力图G1的幅度变化率进行自适应调整，在凝血快速期，能很好的保留信号的细节，也能消除快速期出现的尖峰干扰，多项式拟合后，得到平滑后的弹力图G2，如图5所示。
在本实施方式，比较弹力图G1和G2，弹力图G1中长周期的大脉冲干扰被消除了，形成了平滑后的弹力图G2。
通过本发明的用于血栓弹力图仪的信号滤波方法，可以消除采样噪声信号的干扰、缓变信号中的尖脉冲干扰、大幅度短时间非周期信号干扰，且可以屏蔽信号采集探针高频振动带来的信号震荡，消除长周期大脉冲干扰，达到数据筛选的目的，且在整个凝血纤溶周期，能够很好消除噪声，使得输出的检测数据更加准确。
当然，本发明还可有其它多种实施方式，基于本实施方式，本领域的普通技术人员在没有做出任何创造性劳动的前提下所获得其他实施方式，都属于本发明所保护的范围。

Claims (8)

1.一种用于血栓弹力图仪的信号滤波方法，其特征在于，包括以下步骤：
(1)利用动力学信号采样获得动力学数字信号S1；
(2)将动力学数字信号S1进行均值滤波，获取数字信号S2；
(3)将数字信号S2进行限幅滤波，输出数字信号S3；
(4)对数字信号S3进行中值滤波，获得数字信号S4；
(5)对数字信号S4进行低通滤波处理，获得数字信号S5；
(6)对数字信号S5进行识别处理，计算出多个稳定的极大值点和多个稳定的极小值点；
(7)把所有极大值点组合，生成弹力图上包络线，把所有极小值点组合，
生成弹力图下包络线，得到弹力图G1；
(8)通过高阶多项式拟合技术，原始数据拟合窗口的大小依据弹力图G1的幅度变化率进行自适应调整，消除凝血快速期出现的尖峰干扰，得到平滑后的弹力图G2；
(9)对平滑后的弹力图G2进行计算，可以输出相应的参数结果。
2.根据权利要求1所述的用于血栓弹力图仪的信号滤波方法，其特征在于，其中，步骤(2)～步骤(5)可以任意调整步骤的执行顺序。
3.根据权利要求1所述的用于血栓弹力图仪的信号滤波方法，其特征在于，在步骤(1)中，感应探针浸入血液中，血液凝固和纤溶的强度传递给感应探针转变成模拟电信号，经信号放大器放大后，被处理器AD采样成动力学数字信号S1。
4.根据权利要求1所述的用于血栓弹力图仪的信号滤波方法，其特征在于，在步骤(2)中，对动力学数字信号S1进行N点均值滤波器消除采样噪声，同时保留了动力学数字信号S1中动力学信号的细节,滤波后得到数字信号S2。
5.根据权利要求1所述的用于血栓弹力图仪的信号滤波方法，其特征在于，在步骤(3)中，首先确定N次采样之间最大偏差限幅水平值(THR)，每次检测到新值时进行判断；如果本次与N次采样均值之差小于THR，则本次值有效，如果本次值与前面N次采样值均值大于等于THR，则本次值无效；根据采样值一定窗口大小的前后值设置限幅频率，通过限幅频率来调整限幅水平THR的大小，从而消除数字信号S2中的小幅尖脉冲干扰。
6.根据权利要求1所述的用于血栓弹力图仪的信号滤波方法，其特征在于，在步骤(4)中，设置一个自适应大小窗口N，N的确定可以根据数字信号S3本来的变化快慢来确定大小，通过当前信号的差分量来确定N值，当前信号在窗口的中间，获取N个信号值进行排序，取N个排序后的信号的中间值作为当前信号值，即数字信号S4。
7.根据权利要求1所述的用于血栓弹力图仪的信号滤波方法，其特征在于，在步骤(5)中，低通滤波的截止频率设置为10HZ，滤除数字信号S4中大于10HZ的高频噪声信号。
8.根据权利要求1所述的用于血栓弹力图仪的信号滤波方法，其特征在于，在步骤(6)中，对数字信号S5进行正态分布脉冲识别，通过爬坡算法识别到有效脉冲的初始极大值和初始极小值点，把极值点周围一定窗口大小的信号输入高斯正态分布拟合函数，计算出多个稳定的极大值点和多个稳定的极小值点。
CN201811109037.2A 2018-09-21 2018-09-21 一种用于血栓弹力图仪的信号滤波方法 Active CN109684908B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811109037.2A CN109684908B (zh) 2018-09-21 2018-09-21 一种用于血栓弹力图仪的信号滤波方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811109037.2A CN109684908B (zh) 2018-09-21 2018-09-21 一种用于血栓弹力图仪的信号滤波方法

Publications (2)

Publication Number Publication Date
CN109684908A true CN109684908A (zh) 2019-04-26
CN109684908B CN109684908B (zh) 2023-05-02

Family

ID=66184521

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811109037.2A Active CN109684908B (zh) 2018-09-21 2018-09-21 一种用于血栓弹力图仪的信号滤波方法

Country Status (1)

Country Link
CN (1) CN109684908B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110221544A (zh) * 2019-06-14 2019-09-10 四川烟叶复烤有限责任公司 一种烟叶复烤入口水分数据滤波方法
CN110991275A (zh) * 2019-11-19 2020-04-10 宁波展海电子科技有限公司 红外线信号快速防抖滤波方法
CN111693685A (zh) * 2020-05-18 2020-09-22 中国科学院苏州生物医学工程技术研究所 血栓弹力信号提取方法、存储介质及计算机设备
CN113962252A (zh) * 2021-09-16 2022-01-21 深圳市国赛生物技术有限公司 凝血时间计算方法、装置、系统与可读存储介质
CN117233364A (zh) * 2023-11-16 2023-12-15 苏州思迈德生物科技有限公司 一种血栓弹力图的参数提取方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103398922A (zh) * 2013-07-09 2013-11-20 广东石油化工学院 一种血栓弹力测量装置及其测量方法
US20150024473A1 (en) * 2012-08-22 2015-01-22 Biomedica Usa, Llc Device and Method for Performing Blood Thromboelastographic Assays by Magnetic Sensing
CN107231142A (zh) * 2017-06-21 2017-10-03 北京乐普医疗科技有限责任公司 一种血栓弹力图仪自适应信号处理算法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150024473A1 (en) * 2012-08-22 2015-01-22 Biomedica Usa, Llc Device and Method for Performing Blood Thromboelastographic Assays by Magnetic Sensing
CN103398922A (zh) * 2013-07-09 2013-11-20 广东石油化工学院 一种血栓弹力测量装置及其测量方法
CN107231142A (zh) * 2017-06-21 2017-10-03 北京乐普医疗科技有限责任公司 一种血栓弹力图仪自适应信号处理算法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
SHANSHAN LI等: "Photoelectronic device for measuring of small rotary angles in thromboelastography applications", 《2014 7TH INTERNATIONAL CONFERENCE ON BIOMEDICAL ENGINEERING AND INFORMATICS》 *
罗强等: "血栓弹力图仪恒温控制系统研究", 《自动化应用》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110221544A (zh) * 2019-06-14 2019-09-10 四川烟叶复烤有限责任公司 一种烟叶复烤入口水分数据滤波方法
CN110991275A (zh) * 2019-11-19 2020-04-10 宁波展海电子科技有限公司 红外线信号快速防抖滤波方法
CN110991275B (zh) * 2019-11-19 2024-02-06 宁波展海电子科技有限公司 红外线信号快速防抖滤波方法
CN111693685A (zh) * 2020-05-18 2020-09-22 中国科学院苏州生物医学工程技术研究所 血栓弹力信号提取方法、存储介质及计算机设备
CN111693685B (zh) * 2020-05-18 2022-06-17 中国科学院苏州生物医学工程技术研究所 血栓弹力信号提取方法、存储介质及计算机设备
CN113962252A (zh) * 2021-09-16 2022-01-21 深圳市国赛生物技术有限公司 凝血时间计算方法、装置、系统与可读存储介质
CN113962252B (zh) * 2021-09-16 2023-07-18 深圳市国赛生物技术有限公司 凝血时间计算方法、装置、系统与可读存储介质
CN117233364A (zh) * 2023-11-16 2023-12-15 苏州思迈德生物科技有限公司 一种血栓弹力图的参数提取方法
CN117233364B (zh) * 2023-11-16 2024-02-02 苏州思迈德生物科技有限公司 一种血栓弹力图的参数提取方法

Also Published As

Publication number Publication date
CN109684908B (zh) 2023-05-02

Similar Documents

Publication Publication Date Title
CN109684908A (zh) 一种用于血栓弹力图仪的信号滤波方法
JP4881574B2 (ja) 経穴位置評価装置
US8979760B2 (en) Method and apparatus for automatic optimization of doppler imaging parameters
CN108764073B (zh) 一种结合频谱能量形态拟合的加速度滤噪和积分方法
JP5409907B2 (ja) 機械発振器の共振パラメータの特定
EP2085028A1 (en) Processing low signal-to-noise ratio signals
CN106691376B (zh) 心电信号自适应滤波方法及装置
Maniruzzaman et al. Least-Mean-Square algorithm based adaptive filters for removing power line interference from ECG signal
WO2018082190A1 (zh) 一种ecg信号处理方法及装置
CN106706320B (zh) 一种基于前馈控制随机共振的轴承早期故障诊断方法
CN108427031A (zh) 基于多项式拟合及非干扰区域划分的间谐波检测方法
CN109009086A (zh) 一种心电信号r波检测方法及系统
KR101276973B1 (ko) 맥박수 측정 방법 및 장치
CN111736222A (zh) 单炮数据信噪比确定方法及装置
CN106559055B (zh) 一种用于调制域分析仪中兼容连续波和脉冲调制信号的自动增益控制电路及方法
US8483984B2 (en) Method and apparatus for testing the operating conditions of an electric network
CN106214143B (zh) 瞬时运动干扰识别方法及心率置信度计算方法
CN111839847A (zh) 肌电假肢控制方法、装置、电子设备和计算机可读介质
CN109682958B (zh) 一种用于血栓弹力图仪的加速度传感器信号补偿方法
CN107270954B (zh) 一种振弦式传感器性能分析装置及方法
JP6793299B2 (ja) 信号検出装置及び信号検出方法
CN107582068A (zh) 一种基于周期自相干原理的畸变产物耳声发射探测方案
CN111693685B (zh) 血栓弹力信号提取方法、存储介质及计算机设备
CN109068992A (zh) 睡眠信号调节装置和方法
JPS60210245A (ja) 超音波ドプラ診断装置

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