CN108618778A - 一种麻深指数计算方法及麻醉深度监护仪 - Google Patents
一种麻深指数计算方法及麻醉深度监护仪 Download PDFInfo
- Publication number
- CN108618778A CN108618778A CN201810427545.9A CN201810427545A CN108618778A CN 108618778 A CN108618778 A CN 108618778A CN 201810427545 A CN201810427545 A CN 201810427545A CN 108618778 A CN108618778 A CN 108618778A
- Authority
- CN
- China
- Prior art keywords
- module
- electric signal
- fabo
- signal
- energy 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.)
- Granted
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/369—Electroencephalography [EEG]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/48—Other medical applications
- A61B5/4821—Determining level or depth of anaesthesia
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7225—Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7253—Details of waveform analysis characterised by using transforms
- A61B5/7257—Details of waveform analysis characterised by using transforms using Fourier transforms
Abstract
本发明公开了一种麻醉深度监护仪,包括信号输入模块、放大模块、模数转换模块,信号输入模块用于接收模拟脑皮层电信号,放大模块和信号输入模块相连用于将脑皮层电信号放大,模数转换模块和放大模块相连用于将模拟脑皮层电信号转换为数字脑皮层电信号,还包括计算模块,计算模块和模数转换模块相连用于将获得的数字脑皮层电信号进行计算,计算模块包括傅里叶变换系统、加权计算模块或者自适应神经模糊系统,傅里叶变换系统用于分析脑电信号频谱的数字化能量参数加权计算模块用于计算脑皮层活动电信号的综合能量参数D和BIS值。本发明能够获得较为准确的麻醉深度的监护值,使得BIS的值更加准确。
Description
技术领域
本发明涉及麻醉领域,具体涉及一种麻深指数计算方法及麻醉深度监护仪。
背景技术
现代镇痛实践中由于强效吸入麻醉药、阿片类药、肌松药和静脉麻醉药的使用,使得镇痛的定义不能简单、统一化。麻醉现今的另一种定义就是,麻醉是一种药物诱导的无意识状态,认为意识消失是域值性的(全或无现象),故麻醉不存在深度,而只有合适的麻醉蛇毒。合适的镇痛是一种或几种麻药的浓度达到足以满足手术并使病人舒适的效应。联合用药时由于药理学作用的多样性,用一种测度方法确定不同作用的强度几乎是不可能的。
现在医学认为麻醉也可能是由于中枢神经系统不同水平的功能5甚或是兴奋的结果。麻醉药对网状结构诱发的、自发的神经活动依据麻醉药和选择的神经元不同产生各种影响,有的增强,有的减弱,有的无影响。许多麻醉药是增强自发的网状结构神经活动,而不是简单的抑制它。脑电是来自脑的神经组织的电活动,一切活组织在兴奋过程中都有电位的变化。研究表明,当大脑处于觉醒和睡眠等不同状态时,参加电活动神经单元的数量、每个神经元放电的频率、动作电位的神经传导速度都会有所不同,脑电信号是由脑的神经系统产生的一种可测的电生理反应,由脑电的电位与时间关系组成的图形曲线就是脑电图EEG。脑电图代表大脑皮层的自发神经电活动。而全身麻醉能引起可逆的中枢神经系统的抑制和兴奋,从而达到意识消失和止痛的目的。脑电波的分类,主要根据其频率的不同来人为划分,通常有如下划分方式。
α波:是正常成人安静闭眼时主要的脑电活动,大脑各区均有,以顶枕部最为明显。一般认为α波是皮质处在安静状态时的主要脑电活动表现,为8~13Hz。
β波:以额区和中央区为最明显。当正常人从安静闭眼状态下被唤醒时,会出现波。一般波是新皮质处在特殊紧张活动状态时的主要脑电活动表现,14~30Hz。
BIS EEG分析是应用非线性相位锁定原理对原始EEG波形进行处理的一种方法,属于一种回归的处理方法。是在功率谱分析的基础上又加入了相关函数谱的分析,既测定EEG的线性成分,又分析EEG成分之间的非线性关系。BIS可测定麻醉的催眠部分,对几种临床目标和几种麻醉药有着很好的敏感度和特异度。但是不同数据具有很大的个体差异化,需要将个体差异化考虑进去以提高测试的准确性。因而采用目前的BIS计算方法和装置具有一定的局限性。
发明内容
本发明的目的在于提供一种麻深指数计算方法及麻醉深度监护仪,解决现有差异化不高引起的BIS计算具有局限性的问题。本发明采用将EEG进行能量变换获得数量化能量参数,再对数量化能量参数进行加权和法计算获得综合能量参数,再由综合能量参数和BSR进行加权算法计算或者输入自适应神经模糊逻辑运算系统进行计算获得了较为准确的麻醉深度的监护值,使得BIS的值更加准确。
本发明通过下述技术方案实现:
一种麻深指数计算方法,将得EEG信号进行放大、A/D转换及消除噪声处理,还包括对脑皮层活动性计算:
S1、对脑皮层电信号进行FFT变换获得脑皮层电信号的能量参数
S2、将脑皮层活动电信号的能量参数用加权算法计算得到脑电信号的综合能量参数D。
S3、将综合能量参数D、爆发抑制比BSR%进行加权算法或者输入自适应神经模糊系统计算获得BIS值。
FFT变换为短时FFT变换。
一种麻醉深度监护仪,包括信号输入模块、放大模块、模数转换模块,信号输入模块用于接收模拟脑皮层电信号,放大模块和信号输入模块相连用于将脑皮层电信号放大,模数转换模块和放大模块相连用于将模拟脑皮层电信号转换为数字脑皮层电信号,还包括计算模块,计算模块和模数转换模块相连用于将获得的数字脑皮层电信号进行计算,计算模块包括傅里叶变换系统、加权计算模块或者计算模块包括傅里叶变换系统、加权计算模块、自适应神经模糊系统,傅里叶变换系统用于分析脑电信号频谱的数字化能量参数 加权计算模块用于计算脑皮层活动电信号的综合能量参数D和BIS值;,自适应神经模糊系统用于将综合能量参数D、BSR%输入至自适应神经模糊系统中获得BIS值。
还包括脑电信号采集模块,脑电信号采集模块和信号输入模块相连,用于采集脑电信号并将脑电信号传递给脑电活动电信号计算模块。
还包括脑电信号偏重输入判定模块,脑电信号偏重输入判定模块用于将阿法波出现率输入并判端其大小,依照判定结果归属于对应类型脑电图,对综合能量参数进行加权算法计算,所述类型脑电图分为:阿法波优势型,阿法波出现率在75%以上;阿法波准优势型,阿法波出现率50%~75%;阿法波混合型,阿法波出现率在25%~50%;阿法波劣势型,阿法波出现率为0~25%。脑电信号偏重输入判定模块和计算模块的加权算法模块相连。
本发明充分利用个体差异,利用个体的阿法波型,根据事先测定的脑电图类型差异,将个体差异代入整个监护测试中,使得测试能加客观合理。在实施中,例如如果输入的类型是阿法波优势型,则能量参数中的D1权重势必增大,反之则反。
本发明与现有技术相比,具有如下的优点和有益效果:
本发明采用将EEG进行能量变换获得数量化能量参数,再对数量化能量参数进行加权和法计算获得综合能量参数,再由综合能量参数和BSR进行加权算法计算或者输入自适应神经模糊逻辑运算系统进行计算获得了较为准确的麻醉深度的监护值,使得BIS的值更加准确。
附图说明
此处所说明的附图用来提供对本发明实施例的进一步理解,构成本申请的一部分,并不构成对本发明实施例的限定。在附图中:
图1为本发明结构示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,下面结合实施例和附图,对本发明作进一步的详细说明,本发明的示意性实施方式及其说明仅用于解释本发明,并不作为对本发明的限定。
实施例1
如图1所示,一种麻醉深度监护仪,包括信号输入模块、放大模块、模数转换模块,信号输入模块用于接收模拟脑皮层电信号,放大模块和信号输入模块相连用于将脑皮层电信号放大,模数转换模块和放大模块相连用于将模拟脑皮层电信号转换为数字脑皮层电信号,还包括计算模块,计算模块和模数转换模块相连用于将获得的数字脑皮层电信号进行计算,计算模块包括傅里叶变换系统、加权计算模块,傅里叶变换系统用于分析脑电信号频谱的数字化能量参数加权计算模块用于计算脑皮层活动电信号的综合能量参数D和BIS值。
还包括脑电信号采集模块,脑电信号采集模块和信号输入模块相连,用于采集脑电信号并将脑电信号传递给脑电活动电信号计算模块。
还包括脑电信号偏重输入判定模块,脑电信号偏重输入判定模块用于将阿法波出现率输入并判端其大小,依照判定结果归属于对应类型脑电图,对综合能量参数进行加权算法计算,所述类型脑电图分为:阿法波优势型,阿法波出现率在75%以上;阿法波准优势型,阿法波出现率50%~75%;阿法波混合型,阿法波出现率在25%~50%;阿法波劣势型,阿法波出现率为0~25%。脑电信号偏重输入判定模块和计算模块的加权算法模块相连。
本发明中,脑电信号采集模块为脑电传感器。信号输入模块和模数转换模块均为现有技术。
实施例2
一种麻醉深度监护仪,包括信号输入模块、放大模块、模数转换模块,信号输入模块用于接收模拟脑皮层电信号,放大模块和信号输入模块相连用于将脑皮层电信号放大,模数转换模块和放大模块相连用于将模拟脑皮层电信号转换为数字脑皮层电信号,还包括计算模块,计算模块和模数转换模块相连用于将获得的数字脑皮层电信号进行计算,计算模块包括傅里叶变换系统、加权算法模块、自适应神经模糊系统,傅里叶变换系统用于分析脑电信号频谱的数字化能量参数加权计算模块用于计算脑皮层活动电信号的综合能量参数D,自适应神经模糊系统用于将综合能量参数D、BSR%输入至自适应神经模糊系统中获得BIS值。
还包括脑电信号采集模块,脑电信号采集模块和信号输入模块相连,用于采集脑电信号并将脑电信号传递给脑电活动电信号计算模块。
还包括脑电信号偏重输入判定模块,脑电信号偏重输入判定模块用于将阿法波出现率输入并判端其大小,依照判定结果归属于对应类型脑电图,对综合能量参数进行加权算法计算,所述类型脑电图分为:阿法波优势型,阿法波出现率在75%以上;阿法波准优势型,阿法波出现率50%~75%;阿法波混合型,阿法波出现率在25%~50%;阿法波劣势型,阿法波出现率为0~25%。脑电信号偏重输入判定模块和计算模块的加权算法模块相连。
本发明中,脑电信号采集模块为脑电传感器。信号输入模块和模数转换模块均为现有技术。
实施例3
一种麻深指数计算方法,将得EEG信号进行放大、A/D转换及消除噪声处理,还包括对脑皮层活动性计算:
S1、对脑皮层电信号进行FFT变换获得脑皮层电信号的能量参数
S2、将脑皮层活动电信号的能量参数用加权算法计算得到脑电信号的综合能量参数D。
S3、将综合能量参数D、爆发抑制比BSR%进行加权算法或者输入自适应神经模糊系统计算获得BIS值。
FFT变换为短时FFT变换。
例如若输入的类型是阿法波优势型,则能量参数中的D1权重势必增大;若输入的类型是阿法波混合型,则D1和D2权重相当;若输入的类型是阿法波劣势型,则D2权重增大。由此可针对个体化的不同进行有针对的参数修正。从而获得检测较为客观准确的结果。
以上所述的具体实施方式,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施方式而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (5)
1.一种麻深指数计算方法,将得EEG信号进行放大、A/D转换及消除噪声处理,其特征在于,还包括以下计算:
S1、对脑皮层电信号进行FFT变换获得脑皮层电信号的能量参数
S2、将脑皮层活动电信号的能量参数用加权算法计算得到脑电信号的综合能量参数D。
S3、将综合能量参数D、爆发抑制比BSR%进行加权算法或者输入自适应神经模糊系统计算获得BIS值。
2.根据权利要求1所述的麻深指数计算方法,其特征在于,FFT变换为短时FFT变换。
3.一种麻醉深度监护仪,包括信号输入模块、放大模块、模数转换模块,信号输入模块用于接收模拟脑皮层电信号,放大模块和信号输入模块相连用于将脑皮层电信号放大,模数转换模块和放大模块相连用于将模拟脑皮层电信号转换为数字脑皮层电信号,其特征在于,还包括计算模块,计算模块和模数转换模块相连用于将获得的数字脑皮层电信号进行计算,计算模块包括傅里叶变换系统、加权计算模块或者计算模块包括傅里叶变换系统、加权计算模块、自适应神经模糊系统,傅里叶变换系统用于分析脑电信号频谱的数字化能量参数加权计算模块用于计算脑皮层活动电信号的综合能量参数D和BIS值。
4.根据权利要求3所述的麻醉深度监护仪,其特征在于,还包括脑电信号采集模块,脑电信号采集模块和信号输入模块相连,用于采集脑电信号并将脑电信号传递给脑电活动电信号计算模块。
5.根据权利要求3所述的麻醉深度监护仪,其特征在于,还包括脑电信号偏重输入判定模块,脑电信号偏重输入判定模块用于将阿法波出现率输入并判端其大小,依照判定结果归属于对应类型脑电图,对综合能量参数进行加权算法计算,所述类型脑电图分为:阿法波优势型,阿法波出现率在75%以上;阿法波准优势型,阿法波出现率50%~75%;阿法波混合型,阿法波出现率在25%~50%;阿法波劣势型,阿法波出现率为0~25%。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810427545.9A CN108618778B (zh) | 2018-05-07 | 2018-05-07 | 一种麻深指数计算方法及麻醉深度监护仪 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810427545.9A CN108618778B (zh) | 2018-05-07 | 2018-05-07 | 一种麻深指数计算方法及麻醉深度监护仪 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108618778A true CN108618778A (zh) | 2018-10-09 |
CN108618778B CN108618778B (zh) | 2019-07-02 |
Family
ID=63695579
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810427545.9A Active CN108618778B (zh) | 2018-05-07 | 2018-05-07 | 一种麻深指数计算方法及麻醉深度监护仪 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108618778B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109893127A (zh) * | 2019-03-25 | 2019-06-18 | 惠州善雅医疗器械有限公司 | 麻醉深度监护仪 |
CN113100784A (zh) * | 2021-03-29 | 2021-07-13 | 合肥诺和电子科技有限公司 | 一种csi麻醉深度指数信号发生器 |
CN113116301A (zh) * | 2021-04-01 | 2021-07-16 | 合肥贺康医疗科技有限公司 | 一种麻醉深度多参数监护仪 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5458117A (en) * | 1991-10-25 | 1995-10-17 | Aspect Medical Systems, Inc. | Cerebral biopotential analysis system and method |
US20030167019A1 (en) * | 2000-10-16 | 2003-09-04 | Viertio-Oja Hanna E. | Method and apparatus for determining the cerebral state of a patient with fast response |
CN101301196A (zh) * | 2007-05-11 | 2008-11-12 | 上海世珈医疗设备有限责任公司 | 麻醉镇静深度的脑电图监测装置 |
US7783343B2 (en) * | 2005-07-07 | 2010-08-24 | The General Electric Company | Monitoring of the cerebral state of a subject |
CN103637798A (zh) * | 2013-12-17 | 2014-03-19 | 山东大学齐鲁医院 | 一种麻醉深度监测装置 |
CN104188627A (zh) * | 2014-08-27 | 2014-12-10 | 王远志 | 一种信息化麻醉深度监护仪 |
-
2018
- 2018-05-07 CN CN201810427545.9A patent/CN108618778B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5458117A (en) * | 1991-10-25 | 1995-10-17 | Aspect Medical Systems, Inc. | Cerebral biopotential analysis system and method |
US20030167019A1 (en) * | 2000-10-16 | 2003-09-04 | Viertio-Oja Hanna E. | Method and apparatus for determining the cerebral state of a patient with fast response |
US7783343B2 (en) * | 2005-07-07 | 2010-08-24 | The General Electric Company | Monitoring of the cerebral state of a subject |
CN101301196A (zh) * | 2007-05-11 | 2008-11-12 | 上海世珈医疗设备有限责任公司 | 麻醉镇静深度的脑电图监测装置 |
CN103637798A (zh) * | 2013-12-17 | 2014-03-19 | 山东大学齐鲁医院 | 一种麻醉深度监测装置 |
CN104188627A (zh) * | 2014-08-27 | 2014-12-10 | 王远志 | 一种信息化麻醉深度监护仪 |
Non-Patent Citations (2)
Title |
---|
RAMPIL,ET AL: "A Primer for EEG Signal Processing inAnesthesia", 《ANESTHESIOLOGY》 * |
朱政廉等: "麻醉深度神经检测模型研究现状综述", 《系统仿真技术》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109893127A (zh) * | 2019-03-25 | 2019-06-18 | 惠州善雅医疗器械有限公司 | 麻醉深度监护仪 |
CN113100784A (zh) * | 2021-03-29 | 2021-07-13 | 合肥诺和电子科技有限公司 | 一种csi麻醉深度指数信号发生器 |
CN113116301A (zh) * | 2021-04-01 | 2021-07-16 | 合肥贺康医疗科技有限公司 | 一种麻醉深度多参数监护仪 |
Also Published As
Publication number | Publication date |
---|---|
CN108618778B (zh) | 2019-07-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Denver et al. | Methodological issues in the quantification of respiratory sinus arrhythmia | |
Muthuswamy et al. | Higher-order spectral analysis of burst patterns in EEG | |
Wysocki et al. | Chaotic dynamics of resting ventilatory flow in humans assessed through noise titration | |
Colyer et al. | Textile electrodes embedded in clothing: A practical alternative to traditional surface electromyography when assessing muscle excitation during functional movements | |
CN100482155C (zh) | 基于脑机交互的注意力状态即时检测系统及检测方法 | |
CN101529429A (zh) | 预测心理和神经障碍治疗期间的不利事件的系统和方法 | |
US8280501B2 (en) | Systems and methods for quantitatively characterizing slow wave activities and states in sleep | |
CN101449974A (zh) | 一种自动实时估计麻醉深度的方法 | |
CN103989485A (zh) | 基于脑电波的人体疲劳度评价方法 | |
CN108618778B (zh) | 一种麻深指数计算方法及麻醉深度监护仪 | |
CN102178514A (zh) | 一种基于非线性与复杂性多重指标的昏迷程度评价方法 | |
Namazi et al. | Mathematical modeling of human brain neuronal activity in the absence of external stimuli | |
Ibáñez et al. | Only the fastest corticospinal fibers contribute to β corticomuscular coherence | |
Tiwari et al. | A comparative study of stress and anxiety estimation in ecological settings using a smart-shirt and a smart-bracelet | |
Zhang et al. | Real-time mental stress detection using multimodality expressions with a deep learning framework | |
KR101248118B1 (ko) | 뇌파의 뇌활성도 및 마취심도 분석장치 및 분석방법 | |
Zalewska et al. | Evaluation of MUAP shape irregularity-a new concept of quantification | |
Yang et al. | Identification of healthy and pathological heartbeat dynamics based on ECG-waveform using multifractal spectrum | |
Liu et al. | A novel delay estimation method for improving corticomuscular coherence in continuous synchronization events | |
Namazi et al. | The analysis of the influence of odorant’s complexity on fractal dynamics of human respiration | |
Kew et al. | Wearable patch-type ECG using ubiquitous wireless sensor network for healthcare monitoring application | |
Jelinek et al. | Reducing suicidal ideation by biofeedback-guided respiration–heart rate coherence | |
Murali et al. | Analysis of fractal and fast Fourier transform spectra of human electroencephalograms induced by odors | |
Halliday et al. | Neural spike train synchronization indices: definitions, interpretations, and applications | |
Liu | Model Construction of Using Physiological Signals to Detect Mental Health Status |
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 | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20221207 Address after: Room 103, Building 1, Military and Civilian Integration Science and Technology City, No. 769, Qingzhuhu Road, Qingzhuhu Street, Kaifu District, Changsha, Hunan 410000 Patentee after: Hunan Kemai Technology Co.,Ltd. Address before: 516000 4th floor, No.4 factory building, 2 Baini Road, he'nan'an, Huizhou City, Guangdong Province Patentee before: HUIZHOU SHANYA MEDICAL EQUIPMENT Co.,Ltd. |