CN111870252A - 基于激光传感的浅表血管震颤测量方法与装置 - Google Patents
基于激光传感的浅表血管震颤测量方法与装置 Download PDFInfo
- Publication number
- CN111870252A CN111870252A CN202010707385.0A CN202010707385A CN111870252A CN 111870252 A CN111870252 A CN 111870252A CN 202010707385 A CN202010707385 A CN 202010707385A CN 111870252 A CN111870252 A CN 111870252A
- Authority
- CN
- China
- Prior art keywords
- superficial
- tremor
- laser
- laser sensor
- detected
- 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.)
- Pending
Links
- 206010044565 Tremor Diseases 0.000 title claims abstract description 117
- 230000002792 vascular Effects 0.000 title claims abstract description 37
- 238000000691 measurement method Methods 0.000 title claims abstract description 9
- 238000001914 filtration Methods 0.000 claims abstract description 53
- 210000004204 blood vessel Anatomy 0.000 claims abstract description 44
- 238000000034 method Methods 0.000 claims description 26
- 238000005259 measurement Methods 0.000 claims description 24
- 230000003044 adaptive effect Effects 0.000 claims description 15
- 239000010410 layer Substances 0.000 claims description 11
- 239000012790 adhesive layer Substances 0.000 claims description 8
- 238000004458 analytical method Methods 0.000 claims description 6
- 230000009467 reduction Effects 0.000 claims description 5
- 238000013016 damping Methods 0.000 claims description 3
- 238000002955 isolation Methods 0.000 claims description 3
- 230000002093 peripheral effect Effects 0.000 claims description 2
- 230000001678 irradiating effect Effects 0.000 claims 1
- 206010003226 Arteriovenous fistula Diseases 0.000 abstract description 11
- 230000007310 pathophysiology Effects 0.000 abstract description 5
- 206010016717 Fistula Diseases 0.000 description 16
- 230000003890 fistula Effects 0.000 description 16
- 238000010586 diagram Methods 0.000 description 10
- 230000008859 change Effects 0.000 description 8
- 230000001133 acceleration Effects 0.000 description 7
- 230000000694 effects Effects 0.000 description 7
- 239000008280 blood Substances 0.000 description 6
- 210000004369 blood Anatomy 0.000 description 6
- 210000003462 vein Anatomy 0.000 description 6
- 201000010099 disease Diseases 0.000 description 5
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 5
- 230000017531 blood circulation Effects 0.000 description 4
- 238000003745 diagnosis Methods 0.000 description 4
- 238000006073 displacement reaction Methods 0.000 description 4
- 238000013186 photoplethysmography Methods 0.000 description 4
- 230000004044 response Effects 0.000 description 4
- 210000001367 artery Anatomy 0.000 description 3
- 238000001514 detection method Methods 0.000 description 3
- 239000000284 extract Substances 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000033001 locomotion Effects 0.000 description 2
- 239000004814 polyurethane Substances 0.000 description 2
- 229920002635 polyurethane Polymers 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000001028 reflection method Methods 0.000 description 2
- 108010076504 Protein Sorting Signals Proteins 0.000 description 1
- 206010070834 Sensitisation Diseases 0.000 description 1
- 208000007536 Thrombosis Diseases 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 1
- 229910052782 aluminium Inorganic materials 0.000 description 1
- 238000002555 auscultation Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000000747 cardiac effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000003203 everyday effect Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000009532 heart rate measurement Methods 0.000 description 1
- 239000000017 hydrogel Substances 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000002559 palpation Methods 0.000 description 1
- 231100000915 pathological change Toxicity 0.000 description 1
- 230000036285 pathological change Effects 0.000 description 1
- 230000001575 pathological effect Effects 0.000 description 1
- 230000004802 pathophysiological signal Effects 0.000 description 1
- 230000002980 postoperative effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000008313 sensitization Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 208000019553 vascular disease Diseases 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
- A61B5/11—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
- A61B5/1101—Detecting tremor
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0059—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
- A61B5/0082—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence adapted for particular medical purposes
- A61B5/0084—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence adapted for particular medical purposes for introduction into the body, e.g. by catheters
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
- A61B5/11—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
- A61B5/1126—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb using a particular sensing technique
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Public Health (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Physiology (AREA)
- Dentistry (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
本发明提出一种基于激光传感的浅表血管震颤测量方法与装置,通过激光传感器非接触采集浅表血管震颤信号,再通过中值滤波和/或自适应滤波的方法还原获得纯净的震颤信号。本发明测量的浅表血管震颤可以是脉搏引起的震颤和动静脉瘘震颤等。该装置有单个或者双个激光传感器可供选择进行工作,浅表血管震颤信号可通过单激光配合中值滤波器进行提取,也可以采用双激光配合自适应滤波滤除人体自身震颤信号,从而获得纯净的病生理引起的浅表血管震颤信号。
Description
技术领域
本发明涉及信号检测方法和设备领域,尤其涉及一种基于基于激光传感的浅表血管震颤测量方法与装置,可测量浅表血管,尤其是动静脉内瘘震颤时传导到皮肤表面的震动信号。
背景技术
术语解释:
浅表血管震颤:浅表血管震颤为触诊时指尖感到的一种细微震动感,比如脉搏引起的震颤;也可能为疾病的特征性体征,比如动静脉瘘引起的浅表震颤等。震颤的发生是血液流经狭窄的口径或向异常的方向流动时,形成漩涡造成血管壁震动传至皮肤表面所致。
动静脉瘘:通过手术在皮下将距离相近的某一动脉和浅表静脉血管联通,人工形成体内动静脉间直通管道。
临床医生通过触摸患者,感受浅表血管震颤进行疾病判断,然而一般情况下只能主观感受到其震颤,而无法定量的解释震颤中包含的病生理信号。因此,能准确测量震颤信号的装置成为迫切需求。
如,当患者进行动静脉造瘘术后,会形成一个动静脉内瘘。动静脉内瘘形成后,压力较高的动脉血快速经过压力较低侧的静脉会形成湍流而产生震颤和杂音。手术成功的标志是术后静脉侧能触及震颤,听到血管杂音。术后医嘱普遍包含要求病人每天自我检查震颤及杂音,以便能及时发现血栓等病变及时就医处理。
目前临床过程中,医生用听诊器听诊患者的血管杂音是比较常规的疾病诊断手段,但只能通过主观感受进行判断,并且出院后患者通常无法自主进行听诊。还有就是由于个体性差异,对相同的信号,每个人主观听到的杂音都略有差异,无法客观准确反映动静脉内瘘的状态。
在现有技术中主要提供了以下几种方法:
(1)脉搏测量
光电容积脉搏波描记法(PPG)是借光电手段在活体组织中检测血液容积变化的一种红外无损检测技术。当一定波长的光束照射到指端皮肤表面时,光束将通过透射或反射方式传送到光电接收器。利用光电传感器检测经过人体血液和组织吸收后的反射光强度的不同,描记出血管容积在心动周期内的变化。
(2)动静脉瘘震颤的测量
利用听诊器收集内瘘的声学特性。将听诊器的一端放置于内瘘区域,另一端利用声学传感器将内瘘的杂音转化成电信号,再将信号进行放大再进行收集,最后将震颤信号以声学的方式传递给医生。
(3)利用加速度传感器收集内瘘的震颤特性。
将加速度传感器固定到内瘘的静脉侧,使得内瘘的震颤能带动加速度传感器运动,再将加速度传感器传递出来的信号经过放大后进行收集。但该方法依然需要紧贴皮肤进行测量,可能会对内瘘造成压迫。
在上述测量脉搏引起的震颤方案中,光电容积脉搏波描记法本质上时记录的血液容积变化,不是震颤,该方法不能测量血流冲击血管引起的震颤。此外由于测量部位不可避免的移动,周围环境中自然光、日光灯等光源的干扰,会导致结果进一步地不准确。
听诊器测量方式虽然能采集到内瘘的声学特性,但是受环境影响比较明显,容易将外界环境的嘈杂信号掺杂在正常信号中,影响对内瘘状态的判断,听诊方式不能有效提取信号中次声波成分。
采用加速度传感器,对于震颤信号中的匀速部分或者加速度小的部分无法收集到,信号不能正确反映出测量点包含的血管病生理信息。此外,以动静脉瘘为例,临床上建议除了穿刺和必要检查外尽可能不要在瘘上施加外力。
发明内容
针对现有技术存在的缺陷和不足,本发明提出了一种基于激光传感的浅表血管震颤测量方法与装置,通过激光传感器非接触采集浅表血管震颤信号,再通过中值滤波或自适应滤波的方法还原获得纯净的震颤信号。本发明测量的浅表血管震颤可以是脉搏引起的震颤和动静脉瘘震颤等。该装置有单个或者双个激光传感器可供选择进行工作,浅表血管震颤信号可通过单激光配合中值滤波器进行提取,也可以采用双激光配合自适应滤波滤除人体自身震颤信号,从而获得纯净的病生理引起的浅表血管震颤信号。提出的装置不仅能够在输出信号的基础上以可视化的方式进行展示,还能有效地测量震颤信号,为医生进行病例状态评估提供准确的原始数据来源。
本发明提出的基于激光传感的浅表血管震颤测量方法与装置,能够检测浅表血管震颤,包括脉搏引起的震颤和动静脉瘘震颤等。被试者将手臂放在本装置的手臂放置区,移动可伸缩支架调整激光测量点位置。当实用双激光时,主激光测出的信号将包含手臂自身震颤信号(人体本身含有低频震颤,在本发明中属于噪声)。辅助激光用于人体自身低频震颤噪声信号测量。调整主激光传感器,使其射出的激光对准测量区域,调整辅助激光传感器使其射出的激光避开测量区域,落在其他区域的皮肤上。血液流经血管产生震颤,传递到皮肤表面,导致激光发射口与皮肤表面之间的距离随震颤信号而发生改变,产生的位移变化将被传感器记录下来。最后将采集到的两路信号传递到主机进行分析,通过自适应滤波的方法去除手臂自身发生的位移带来的干扰,得到干净的浅表血管震颤信号。或者,主机上采用主传感器信号结合中值滤波滤除噪声,获得浅表血管震颤信号。
与现有装置相比,本发明的方案是采用非接触式设备,且具备大响应范围(0-2000Hz),能够有效覆盖所测量的病理和生理信号。发明采用激光测距原理测量血液流经血管产生震颤导致的皮肤表面的距离变化,浅表血管震颤信号可通过单激光配合中值滤波器进行提取,也可以采用双激光配合自适应滤波滤除人体自身震颤信号,从而获得纯净的病生理引起的浅表血管震颤信号。本发明属于无创应用,同时由于采用激光测距原理间接测量震颤,测距精度可以达到0.1um,所以测量结果准确。发明提供的方法不需要与皮肤直接接触,不会给检查部位施加外力,安全性高。同时与光电容积脉搏波描记法相比,激光光源受外界环境影响小,不受环境中光源的干扰,使用场景比较广。与听诊器方案相比,使用激光三角反射法采集震颤信号,采集过程干扰小,精度高,能够有效的收集0-2000Hz范围内,特别是含次声波频段内的信号。与加速度传感器方案相比,由于激光的单色性好、方向性强,所以测距精度高,对微小信号变化敏感性高,且不需要与皮肤直接接触,安全性也更高。
其具体采用以下技术方案:
一种基于激光传感的浅表血管震颤测量方法,其特征在于:通过激光传感器非接触采集浅表血管震颤信号,再通过中值滤波和/或自适应滤波还原获得纯净的震颤信号。
优选地,采用单激光传感器照射皮肤上待测浅表血管所在区域,将采集到的反射光对应的数字信号进行中值滤波,还原获得纯净的震颤信号。
优选地,所述中值滤波包括以下步骤:
步骤A1:定义一个长度为奇数的K长窗口,K=2N+1,N为正整数;
步骤A2:设在某一个时刻,窗口内的信号样本为x(i-N),…,x(i),…,x(i+N),其中x(i)为位于窗口中心的信号样本值;
步骤A3:对K个信号样本值按从小到大的顺序排列后,取其中值为在i处的样值,定义为中值滤波的输出值。
优选地,采用双激光传感器,其中,主激光传感器照射皮肤上待测浅表血管所在区域,辅助激光传感器照射皮肤上非待测浅表血管所在区域,将采集到的双激光传感器的反射光对应的数字信号进行自适应滤波,还原获得纯净的震颤信号。
优选地,采用最小均方差准则LMS作为自适应滤波的最优准则。
一种基于激光传感的浅表血管震颤测量装置,其特征在于,包括:主激光传感器、辅助激光传感器、U型手臂放置槽和主机;所述U型手臂放置槽用于放置待检测人员的手臂;所述主激光传感器对准皮肤上待测浅表血管所在区域;所述辅助激光传感器对准皮肤上非待测浅表血管所在区域;所述主激光传感器和辅助激光传感器分别连接主机。
优选地,所述主激光传感器和辅助激光传感器分别通过带阻尼的多自由度激光传感器支撑单元安装在与U型手臂放置槽集成的底座上。
优选地,还包括可供选择的反光片,用于敷设于皮肤上待测浅表血管所在区域及其周边区域,用于增强激光测量信号的效果。
优选地,所述反光片是由反光层、粘贴层和隔离纸组成的方形贴片,尺寸至少为2mm x 2mm。
以及根据以上优选装置的使用方法,其特征在于,包括以下步骤:
步骤B1:将待检测人员的手臂放在底座上的U型手臂放置槽,手臂内侧朝上;
步骤B2:在待测浅表血管所在区域贴上反光片:撕开反光片上的隔离纸,反光层向上,粘贴层贴在待测浅表血管所在区域,与使用者的皮肤表面接触;
步骤B3:调节两侧的激光传感器支撑单元的位置和角度,使激光传感器位于待测浅表血管所在区域上方;
步骤B4:调整主激光传感器的支撑单元,使其射出的激光对准待测浅表血管所在区域的反光片上,采集浅表血管震颤信号;
步骤B5:若要进行自适应滤波,则再调整辅助激光传感器的支撑单元,使其射出的激光避开待测浅表血管所在区域,落在待测浅表血管所在区域附近皮肤的反光片上,采集人体自身的震颤信号;若进行中值滤波,则跳过此步骤;
步骤B6:将主激光传感器和辅助激光传感器采集到的信号传递到主机进行分析;
步骤B7:单激光配合中值滤波器提取震颤信号,双激光配合自适应滤波滤除人体自身震颤信号,从而获得纯净的浅表血管震颤信号。
本发明及其优选方案相比于现有技术有几大显著优点:
1、将震颤转化为位移进行测量。提出了一种基于激光传感的浅表血管震颤测量方法与装置,测量血液流经血管产生的浅表震颤。选择的激光光源测量方法,受外界环境影响小,抗干扰性好,传感响应速度快,最高精度可以达到1um,能够在自然光下做到高频震颤测量(可实现0-2000Hz振动分析)。
2、本发明可通过单激光配合中值滤波器提取信号,在保护信号轮廓的同时对信号进行平滑处理,剔除脉冲型噪声又可以使震颤信号的参数不受影响。
3、本发明可采用自适应滤波法进行信号的增强和提取,将采集到的浅表血管震颤信号和人体自身的震颤信号,利用自适应抵消模型提取干净、高质量的浅表血管震颤信号。
4、本发明属于无创技术,不需要与皮肤直接接触,安全性高。
上述提及的指标参数受传感器性能决定,选用更好的传感器可以达到更高的性能。
需要特别强调的是,本发明方案并不属于疾病的诊断和治疗方案,因采集震颤信息仅为对人体信息的一种采集,可以作为辅助医生进行判断的参考,但并不能直接形成对疾病的诊断判断及治疗。
附图说明
下面结合附图和具体实施方式对本发明进一步详细的说明:
图1是本发明实施例装置整体结构示意图;
图2是本发明实施例采用的体表震颤的激光三角反射法原理示意图;
图3是本发明实施例采用的窗长K=1.5s的中值滤波结果示意图;
图4是本发明实施例中值滤波后提取的效果示意图;
图5是本发明实施例自适应滤波原理示意图;
图6是本发明实施例采用自适应滤波测量得到的震颤信号示意图(瘘口信号);
图7是本发明实施例采用自适应滤波测量得到的震颤信号示意图(无震颤皮肤信号)
图8是本发明实施例自适应滤波提取结果示意图。
图中:1-激光传感器支撑单元;2-主激光传感器;3-辅助激光传感器;4-底座;5-U型手臂放置槽;6-反光片;7-主机;8-反光层;9-粘贴层;10-隔离纸。
具体实施方式
为让本专利的特征和优点能更明显易懂,下文特举实施例,并配合附图,作详细说明如下:
1.装置总述
本实施例装置整体结构如图1所示,被试者将手臂放在底座4上的U型手臂放置槽5,调节激光传感器支撑单元1,该支撑单元是带阻尼的多自由度激光传感器支撑单元,能够固定并调整激光与测量点之间的位置和距离。调整主激光传感2器,使其射出的激光对准测量区域,调整辅助激光传感器3使其射出的激光避开测量区域,落在其他区域的皮肤上。血液流经血管产生震颤,传递到皮肤表面,导致激光发射口与皮肤表面之间的距离随震颤信号而发生改变,产生的位移变化被主激光传感器2和辅助激光传感器3记录。最后将采集到的两路信号传递到主机7进行分析,浅表血管震颤信号可通过单激光配合中值滤波器进行提取,也可以采用双激光配合自适应滤波滤除人体自身震颤信号,从而获得纯净的病生理引起的浅表血管震颤信号。
为了增强激光的反射效果,根据使用装置时的具体情况,还可以在测量区域贴上反光片6。反光片6是由反光层8、粘贴层9和隔离纸10组成的方形贴片,考虑到反光片6需覆盖激光到达测量点的所有光束,因此建议反光片6的尺寸至少为2mm x 2mm以上。反光层8采用铝模结构,具有良好的光反射性能。粘贴层为聚氨酯水凝9,便于将反光片6贴在测量区域。聚氨酯水凝胶有生物相容、无致敏、不粘黏人体毛发、使用周期长等优点。使用时,撕开反光片6上的隔离纸10,将粘贴层9贴在测量区域,与使用者的皮肤表面接触,反光层8则向上,确保激光光束落在反光层8上。需要注意,在本装置中反光片并不是必须的,但是它的使用能够增强测量的效果。
2.测量原理
激光距离传感器发射出的激光经患者皮肤表面的反射后又被测距仪接收,传感器同时记录下了激光往返的时间。假设光在空气中传播的速度为c,在传感器和被测量点之间往返一次所需的时间为t,则传感器和被测量点之间的距离L可用以下公式表示:
本实施例的激光传感器采用的是激光三角反射法原理,原理图如图2所示。半导体激光器通过镜头将可见的红色激光射向被测物体表面,经被测物体反射的激光通过接收器镜头,被内部的线性CCD阵列接收,根据不同的距离,线性CCD阵列可以在不同的角度下“看见”这个光点。当物体发生位移,由位置A移动到位置B时,激光被反射到线性CCD阵列上的不同位置。根据这个角度及已知的激光和相机之间的距离,数字信号处理器通过三角函数计算阵列上的光点位置,就能计算出传感器和被测物体之间的距离。同时,光束在接收元件的位置通过模拟电路和数字电路的处理,并通过微处理器分析,计算出相应的输出值,并在设定的模拟量窗口内,按比例输出标准数据信号。本实施例使用的高精度激光三角测量传感器,传感响应速度快,精度可以达到0.1um,能提供高精度的测量检测。
3.中值滤波法
生物信号一般都要用数字滤波器进行处理以消除由于生理或物理原因产生的噪声对信号质量的影响。线性滤波在不同程度上都能用于信号的处理,但线性滤波无法完全抑制脉冲型噪声,而本实施例提供的非线性的中值滤波方法可以在保护信号轮廓的同时对信号进行平滑处理。特别是既可以剔除脉冲型噪声又可以使震颤信号的参数不受影响。
中值滤波是基于排序统计理论的一种能有效抑制噪声的非线性信号处理技术。方法是通过采样窗口取出奇数个数据进行排序,用排序后的中值取代要处理的数据。中值滤波的算法比较简单,也易于用硬件实现。对一个数字信号序列进行滤波处理时,本实施例提供的具体步骤如下:
(1)定义一个长度为奇数的K长窗口,K=2N+1,N为正整数。
(2)设在某一个时刻,窗口内的信号样本为x(i-N),…,x(i),…,x(i+N),其中x(i)为位于窗口中心的信号样本值。
对这K个信号样本值按从小到大的顺序排列后,取其中值为在i处的样值,即定义为中值滤波的输出值。
其实施的效果如图3、图4所示。
4.自适应滤波法
自适应滤波器具有滤波效果好,可操作性强等优点,属于数字滤波器的一种,但它不同于一般的数字滤波器,它不采用固定的滤波系数,能够跟踪和适应系统或环境的变化,在工程应用中不需要预先知道干扰噪声的统计特性,通过逐次迭代能够自适应地将其自身的工作状态调整到最佳。自适应滤波器能够随着外界信号特性动态地改变参数,保持最佳滤波状态,参数的调整是由自适应算法来决定的。本实施例选用最小均方差准则(LMS)作为自适应滤波的最优准则,该算法具有计算量小、稳定性好且易于实现等优点。
最小均方差准则(LMS)是一种迭代算法,具体流程为:设定初始的滤波器权重系数ωi(n),读取自适应滤波器的输入u(n)和期待响应信号d(n)。计算自适应滤波器的输出y(n),公式如下:
滤波器输出y(n)相对于d(n)的误差为e(n),有:
e(n)=d(n)-y(n)
得到计算因子μe(n),更新滤波器权重系数,进行迭代:
ωi(n+1)=ωi(n)+2μe(n)u(n-i)
LMS算法收敛的充分条件是:
其中:μ表示步长因子;M表示FIR滤波器阶数,有0<i<M-1;λmax表示自相关矩阵R的最大特征值。
在本实施例中,两个激光传感器分别采集到浅表血管震颤信号s和人体自身的震颤信号NNL。浅表血管震颤信号s包含有效的浅表血管震颤信号X和人体自身的震颤噪声信号NL。由于两个激光传感器不是测量同一个位置,有NNL≠NL,但是他们共同包含有人体自身震颤的信号成分。为获取有效的浅表血管震颤信号X,本发明采用最小均方误差准则(LMS)对收集的信号进行滤波处理,滤除人体自身的震颤信号NNL。将带有噪声信号NL的浅表血管震颤信号s作为主通道的信号源,人体自身的震颤信号NNL作为参考通道的信号源。自适应滤波原理框图如图5所示。
其中,主通道包含有效的浅表血管震颤信号X和人体自身的震颤噪声信号NL:
d(n)=S=X+NL
e(n)=d(n)-y(n)
=X+NL-NNL
计算均方差:
ξ=E[e2(n)]
=E{[X+NL-NNL]2}
=E[X2]+E{[NL-NNL]2}+2E{X[NL-NNL]}
已知NNL≠NL,要想输出稳定的浅表血管震颤信号,方差ξ需要达到最小值,令:
E{[v0(n)-y(n)]2}→0
可通过调节自适应滤波器的权重系数ω(n),使得:
ω(n)→ω*(n)
此时方差ξ最小,采集到的人体自身的震颤信号NNL将十分接近人体自身的震颤噪声信号NL,取系统的误差信号作为输出信号,输出信号为仅含有效的浅表血管震颤信号X。
其最终的信号提取结果如图6-图8所示。
5.测量流程
根据本实施例以上提供的装置和算法,即可以实现对浅表血管震颤信号的采集,以对内瘘部位的血管震颤为例,包括以下步骤:
第一步:将待检测人员的手臂放在底座上的U型手臂放置槽,手臂内侧朝上使得动静脉内瘘朝上。
第二步:在测量区域贴上反光片。撕开反光片上的隔离纸,反光层向上,粘贴层贴在测量区域,与使用者的皮肤表面接触。
第三步:调节两侧的激光传感器支撑单元的位置和角度,使激光传感器位于内瘘上方。
第四步:调整主激光传感器的支撑单元,使其射出的激光对准内瘘测量区域的反光片上,采集浅表血管震颤信号。
第五步:若要进行自适应滤波,则再调整辅助激光传感器的支撑单元,使其射出的激光避开内瘘测量区域,落在内瘘附近皮肤的反光片上,采集人体自身的震颤信号。若进行中值滤波,则跳过此步骤。
第六步:将主激光传感器和辅助激光传感器采集到的信号传递到主机进行分析。
第七步:单激光配合中值滤波器提取震颤信号,双激光配合自适应滤波滤除人体自身震颤信号,从而获得纯净的病生理引起的浅表血管震颤信号。
本专利不局限于上述最佳实施方式,任何人在本专利的启示下都可以得出其它各种形式的基于激光传感的浅表血管震颤测量方法与装置,凡依本发明申请专利范围所做的均等变化与修饰,皆应属本专利的涵盖范围。
Claims (10)
1.一种基于激光传感的浅表血管震颤测量方法,其特征在于:通过激光传感器非接触采集浅表血管震颤信号,再通过中值滤波和/或自适应滤波还原获得纯净的震颤信号。
2.根据权利要求1所述的基于激光传感的浅表血管震颤测量方法,其特征在于:采用单激光传感器照射皮肤上待测浅表血管所在区域,将采集到的反射光对应的数字信号进行中值滤波,还原获得纯净的震颤信号。
3.根据权利要求2所述的基于激光传感的浅表血管震颤测量方法,其特征在于:所述中值滤波包括以下步骤:
步骤A1:定义一个长度为奇数的K长窗口,K=2N+1,N为正整数;
步骤A2:设在某一个时刻,窗口内的信号样本为x(i-N),…,x(i),…,x(i+N),其中x(i)为位于窗口中心的信号样本值;
步骤A3:对K个信号样本值按从小到大的顺序排列后,取其中值为在i处的样值,定义为中值滤波的输出值。
4.根据权利要求1所述的基于激光传感的浅表血管震颤测量方法,其特征在于:采用双激光传感器,其中,主激光传感器照射皮肤上待测浅表血管所在区域,辅助激光传感器照射皮肤上非待测浅表血管所在区域,将采集到的双激光传感器的反射光对应的数字信号进行自适应滤波,还原获得纯净的震颤信号。
5.根据权利要求4所述的基于激光传感的浅表血管震颤测量方法,其特征在于:采用最小均方差准则LMS作为自适应滤波的最优准则。
6.一种基于激光传感的浅表血管震颤测量装置,其特征在于,包括:主激光传感器、辅助激光传感器、U型手臂放置槽和主机;所述U型手臂放置槽用于放置待检测人员的手臂;所述主激光传感器对准皮肤上待测浅表血管所在区域;所述辅助激光传感器对准皮肤上非待测浅表血管所在区域;所述主激光传感器和辅助激光传感器分别连接主机。
7.根据权利要求6所述的基于激光传感的浅表血管震颤测量装置,其特征在于:所述主激光传感器和辅助激光传感器分别通过带阻尼的多自由度激光传感器支撑单元安装在与U型手臂放置槽集成的底座上。
8.根据权利要求7所述的基于激光传感的浅表血管震颤测量装置,其特征在于:还包括反光片,用于敷设于皮肤上待测浅表血管所在区域及其周边区域。
9.根据权利要求8所述的基于激光传感的浅表血管震颤测量装置,其特征在于:所述反光片是由反光层、粘贴层和隔离纸组成的方形贴片,尺寸至少为2mm x 2mm。
10.根据权利要求9所述的基于激光传感的浅表血管震颤测量装置的使用方法,其特征在于,包括以下步骤:
步骤B1:将待检测人员的手臂放在底座上的U型手臂放置槽,手臂内侧朝上;
步骤B2:在待测浅表血管所在区域贴上反光片:撕开反光片上的隔离纸,反光层向上,粘贴层贴在待测浅表血管所在区域,与使用者的皮肤表面接触;
步骤B3:调节两侧的激光传感器支撑单元的位置和角度,使激光传感器位于待测浅表血管所在区域上方;
步骤B4:调整主激光传感器的支撑单元,使其射出的激光对准待测浅表血管所在区域的反光片上,采集浅表血管震颤信号;
步骤B5:若要进行自适应滤波,则再调整辅助激光传感器的支撑单元,使其射出的激光避开待测浅表血管所在区域,落在待测浅表血管所在区域附近皮肤的反光片上,采集人体自身的震颤信号;若进行中值滤波,则跳过此步骤;
步骤B6:将主激光传感器和辅助激光传感器采集到的信号传递到主机进行分析;
步骤B7:单激光配合中值滤波器提取震颤信号,双激光配合自适应滤波滤除人体自身震颤信号,从而获得纯净的浅表血管震颤信号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010707385.0A CN111870252A (zh) | 2020-07-21 | 2020-07-21 | 基于激光传感的浅表血管震颤测量方法与装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010707385.0A CN111870252A (zh) | 2020-07-21 | 2020-07-21 | 基于激光传感的浅表血管震颤测量方法与装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111870252A true CN111870252A (zh) | 2020-11-03 |
Family
ID=73156252
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010707385.0A Pending CN111870252A (zh) | 2020-07-21 | 2020-07-21 | 基于激光传感的浅表血管震颤测量方法与装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111870252A (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120283802A1 (en) * | 2007-01-17 | 2012-11-08 | Chih-Ching Lin | Far-infrared radiation in use of improving patency of arteriovenous fistula, decreasing failure of arteriovenous fistula maturation, and preventing and/or ameliorating peripheral artery diseases |
CN103479332A (zh) * | 2013-09-03 | 2014-01-01 | 河海大学 | 一种动静脉内瘘脉动测试系统 |
CN104116510A (zh) * | 2014-08-11 | 2014-10-29 | 西华大学 | 一种用于帕金森病人震颤的传感装置及检测方法 |
CN107949349A (zh) * | 2015-03-10 | 2018-04-20 | 斯滕特泰克有限公司 | 手术系统、及其用于经皮动静脉瘘(avf)创建的装置和方法 |
CN110179494A (zh) * | 2019-05-31 | 2019-08-30 | 福建工程学院 | 基于加速度传感器的动静脉内瘘震颤测量装置及测量方法 |
CN210963605U (zh) * | 2019-06-26 | 2020-07-10 | 上海军慧医疗管理有限公司 | 一种用于改善动静脉内瘘损伤的便携式可调节低强度激光治疗仪 |
-
2020
- 2020-07-21 CN CN202010707385.0A patent/CN111870252A/zh active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120283802A1 (en) * | 2007-01-17 | 2012-11-08 | Chih-Ching Lin | Far-infrared radiation in use of improving patency of arteriovenous fistula, decreasing failure of arteriovenous fistula maturation, and preventing and/or ameliorating peripheral artery diseases |
CN103479332A (zh) * | 2013-09-03 | 2014-01-01 | 河海大学 | 一种动静脉内瘘脉动测试系统 |
CN104116510A (zh) * | 2014-08-11 | 2014-10-29 | 西华大学 | 一种用于帕金森病人震颤的传感装置及检测方法 |
CN107949349A (zh) * | 2015-03-10 | 2018-04-20 | 斯滕特泰克有限公司 | 手术系统、及其用于经皮动静脉瘘(avf)创建的装置和方法 |
CN110179494A (zh) * | 2019-05-31 | 2019-08-30 | 福建工程学院 | 基于加速度传感器的动静脉内瘘震颤测量装置及测量方法 |
CN210963605U (zh) * | 2019-06-26 | 2020-07-10 | 上海军慧医疗管理有限公司 | 一种用于改善动静脉内瘘损伤的便携式可调节低强度激光治疗仪 |
Non-Patent Citations (3)
Title |
---|
THAD B.WELCH 等: "《实时数字信号处理(从MATLAB到C)-基于TMS320C6x DSP(第3版)》", 30 June 2020 * |
张建忠 等: "《传感器与检测技术》", 31 August 2014 * |
曾海峰 等: "《机电一体化系统设计》", 31 May 2018 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8419643B2 (en) | Ultrasonic method and apparatus for assessment of bone | |
US5727561A (en) | Method and apparatus for non-invasive detection and analysis of turbulent flow in a patient's blood vessels | |
US20220378299A1 (en) | Noninvasive method for measuring sound frequencies created by vortices in a carotid artery, visualization of stenosis, and ablation means | |
CN110141196B (zh) | 基于双三角血流模型的外周动脉血管弹性评价系统 | |
KR20100065084A (ko) | 움직임 잡음에 강인한 맥파 측정 장치 및 그 방법 | |
CN109640828B (zh) | 超声波血流量监测 | |
JP2009527270A (ja) | 内皮機能の非侵襲的評価のための方法および装置 | |
Saito et al. | Noninvasive assessment of arterial stiffness by pulse wave analysis | |
US20180289297A1 (en) | System and method for monitoring glucose level | |
JP5534674B2 (ja) | 非正弦波周期性の流動作用を測定又は検出する超音波法及び装置 | |
JP2009525117A5 (zh) | ||
CN113598728A (zh) | 生理信号的降噪方法、监测方法、监测装置及可穿戴设备 | |
Joseph et al. | ARTSENS® Pen: A portable, image-free device for automated evaluation of vascular stiffness | |
CN111870252A (zh) | 基于激光传感的浅表血管震颤测量方法与装置 | |
Arathy et al. | An accelerometer probe for local pulse wave velocity measurement | |
CN105725983B (zh) | 一种周围动脉硬化的早期筛查方法及系统 | |
CN112367912A (zh) | 改进的个人健康数据收集 | |
US20150305622A1 (en) | Optical Low Coherence Microphone | |
CN111870251B (zh) | 电容传感原理的动静脉内瘘震颤测量方法及装置 | |
Pratiwi et al. | A review of equipment and signal processing of the automated auscultation for blood pressure measurement | |
KR100908998B1 (ko) | 원격 맥진 장치 및 방법 | |
KR100401213B1 (ko) | 심음과 맥파를 이용한 맥파전달속도 측정시스템 | |
JP2000175911A (ja) | 骨強度測定方法及び装置 | |
JP2002524181A (ja) | 超音波伝搬プローブの較正のための装置及び方法 | |
RU2302199C1 (ru) | Способ исследования биомеханических свойств суставов |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20201103 |
|
RJ01 | Rejection of invention patent application after publication |