CN113892931A - 一种基于深度学习的fmcw雷达提取分析腹内压力方法 - Google Patents

一种基于深度学习的fmcw雷达提取分析腹内压力方法 Download PDF

Info

Publication number
CN113892931A
CN113892931A CN202111199566.8A CN202111199566A CN113892931A CN 113892931 A CN113892931 A CN 113892931A CN 202111199566 A CN202111199566 A CN 202111199566A CN 113892931 A CN113892931 A CN 113892931A
Authority
CN
China
Prior art keywords
signal
formula
abdominal
intra
abdominal pressure
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
CN202111199566.8A
Other languages
English (en)
Other versions
CN113892931B (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.)
Chongqing University
Original Assignee
Chongqing University
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 Chongqing University filed Critical Chongqing University
Priority to CN202111199566.8A priority Critical patent/CN113892931B/zh
Publication of CN113892931A publication Critical patent/CN113892931A/zh
Application granted granted Critical
Publication of CN113892931B publication Critical patent/CN113892931B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/05Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/03Detecting, measuring or recording fluid pressure within the body other than blood pressure, e.g. cerebral pressure; Measuring pressure in body tissues or organs
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/103Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/11Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
    • A61B5/113Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb occurring during breathing
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7203Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/725Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7253Details of waveform analysis characterised by using transforms
    • A61B5/7257Details of waveform analysis characterised by using transforms using Fourier transforms
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7264Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Pathology (AREA)
  • Biophysics (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Veterinary Medicine (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Public Health (AREA)
  • Artificial Intelligence (AREA)
  • Physiology (AREA)
  • Signal Processing (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Psychiatry (AREA)
  • Mathematical Physics (AREA)
  • Dentistry (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Hematology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Evolutionary Computation (AREA)
  • Fuzzy Systems (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)

Abstract

本发明涉及生命体特征采集计算领域,尤其涉及一种基于深度学习的FMCW雷达提取分析腹内压力方法,通过FMCW雷达来采集测试对象含腹部呼吸信号的生命体特征信号,将该信号进行特征提取,而后输入至预先搭建好的神经网络中进行计算,通过引入注意力机制,对多通道的腹部呼吸信号进行重要性计算,从不同重要性的通道自适应的提取出不同量的特征,并与通过离散小波变换所提取出的隐藏特征相结合,提高该神经网络算法的能力,获取更多的腹部呼吸信号中的时域信息,将时域信息与测试对象身体腹内压力相关联,从而推算出测试对象身体的腹内压力;相对与传统侵入式或接触式的测量腹内压力的方法,该方法是一种全新的无接触式的测量腹内压力的方法。

Description

一种基于深度学习的FMCW雷达提取分析腹内压力方法
技术领域
发明涉及生命体特征采集计算领域,尤其涉及一种基于深度学习的FMCW雷达提取分析腹内压力方法。
背景技术
现有的腹内压力在临床上有两种测量方式:
1.直接测压,置管于腹腔内,然后连接压力传感器或者腹腔镜手术中通过自动气腹机对压力进行连续监测,直接腹腔内穿刺测压是将针头刺入或导管置入腹腔,接生理盐水并将通过三通与换能器或输液器连接,测量时以腋中线为零点测压。或是通过微导管测压,采用尖端带有微电极的导管置入腹腔内,连接监测仪进行测压,但是这些方法由于是直接测压,采用了侵入式的测压方法,会有一定的创伤,导管易堵塞和有腹腔感染的风险;
2.间接测压,通过测量下腔静脉压力,膀胱压力及胃内压力间接反映腹腔内压力,可以经过股静脉置管测定下腔静脉压力,其与腹内压力变化有较好的相关性,放置股静脉时导管尖端达腹腔下腔静脉位置,通过三通连接股静脉置管和压力换能器进行测压。经膀胱测压在危重病人置有尿管时,不需要重新置尿管,通过连接测压管或者传感器,通过水柱或者传感器的监护仪数据测量膀胱内压,通过已经研究发现的相关性可以推出腹内压。同样也是侵入式,当不是危重患者时,无法插入尿管至膀胱进行测量,并且测量膀胱内压,对病人测量时的舒适度会有很大的影响。
发明内容
针对现有技术中所存在的不足,发明提供了一种基于深度学习的FMCW雷达提取分析腹内压力方法,通过该方法可以避免现有临床上的接触式测压方法,提升患者在采集腹内压力时的舒适度以及避免采集时意外症状的发生。
根据发明的实施例,一种基于深度学习的FMCW雷达提取分析腹内压力方法,包括下述步骤:
S1:将FMCW毫米波雷达系统,放置在观测对象旁边,朝向观测对象腹部位置发射毫米级别的线性调频连续波信号,采集得到从观测对象腹部反射回来的雷达波;将反射回来的雷达波与原始发射信号相混频得到原始数据信号;
S2:所述原始数据信号通过快速傅里叶变换提取出包含有观测对象腹部呼吸信号的生命体征信号,再通过巴特沃斯滤波器滤去噪声;
S3:将包含有观测对象腹部呼吸信号的生命体征信号输入搭建好的神经网络中,对包含有观测对象腹部呼吸信号的生命体征信号进行特征提取,并通过回归分析与腹内压力进行关联,进而分析出腹内压力。
进一步的,在S2中,将取得的原始数据信号定义为X1(t),对X1(t)进行相位特征
Figure BDA0003304414880000026
的提取,将X1(t)的相位
Figure BDA0003304414880000027
定义为:X1(t)的虚部I与实部R之比;拟出公式(1):
Figure BDA0003304414880000021
原始数据信号X1(t)可以被分解为原始信号实部R(t),原始信号虚部I(t),原始数据信号 X1(t)还包含有实部的直流信号dcr,虚部的直流信号dci;将原始信号实部R(t),原始信号虚部I(t)、实部的直流信号dcr、虚部的直流信号dci代入公式(1),可以得到相位特征
Figure BDA0003304414880000028
的计算公式(2):
Figure BDA0003304414880000022
进一步的,在S2中,将得到的
Figure BDA0003304414880000029
重新定义为:包含有腹部呼吸的时域信号X2(t);而后通过快速傅里叶变换将X2(t)变换为腹部呼吸的频域信号X3(k),对X2(t)的频域进行分析;计算公式为:
Figure BDA0003304414880000023
式(3)中:k=0,1,......N-1;
Figure BDA0003304414880000024
得到包含有腹部呼吸的频域信号X3(k)后,由于实际腹部呼吸信号频率范围为0.2Hz~0.4Hz,在X3(k)中观察到符合实际范围的腹部呼吸频域信号X4(k),将X3(k)与设计的巴特沃斯滤波器频域相乘,对X4(k)进行滤噪;
巴特沃斯滤波器公式为:
|Ha(jω)|2=1/1+(ω/ωc)2N (4)
式(4)中:N-滤波器阶数,ωc-滤波器的截止频率,Ha(jω)-设计的滤波器幅频响应;
滤噪腹部呼吸频域信号X4(k)计算公式为:
X4(k)=X3(k)*Ha(jω) (5)
式(5)中:Ha(jω)-设计的巴特沃斯滤波器幅频响应;
再通过快速傅里叶逆变换将滤噪后的腹部呼吸频域信号X4(k)变换为腹部呼吸时域信号X4(t);
计算公式为:
Figure BDA0003304414880000025
式(6)中:k=0,1,......N-1;
Figure BDA0003304414880000031
进一步的,在S3中:将获取的腹部时域呼吸信号X4(t)输入到预先搭建的深度学习神经网络,进行进一步的特征提取以及信号映射到腹内压力的分析,其具体步骤为:
S3.1:将腹部时域呼吸信号X4(t)通过卷积神经网络进行一次特征提取,得到提取的特征定义为信号X5(t);X4(t)为输入的四通道二维矩阵,设进行卷积提取特征的卷积核矩阵为 k,阶矩阵为n×n;
提取特征X5(t)的计算公式为:
Figure BDA0003304414880000032
式(7)中:
Figure BDA0003304414880000033
ki,j-卷积核;
S3.2:再通过线性整流函数加入非线性因素,增强模型表达力,
线性整流函数公式为:
X6(t)=max(0,X5(t)) (8)
S3.3:对S3.1中的卷积神经网络以及S3.2的线性整流函数进行五层堆叠计算(即将式 (7)与式(8)堆叠计算5次),目的是提取四个通道的数据特征,之后引入非线性因素,得到特征信号X7(t);此时,设X7(t)中大小为H*W的矩阵为C个;加入压缩-激活神经网络,引入注意力机制,对每个通道进行特征权重的计算,其具体步骤为:
S3.3.1:压缩-激活神经网络首先需要进行对C个特征图进行压缩操作,获取上述提取特征图的全局信息;
计算公式为:
Figure BDA0003304414880000034
式(9)中:W-特征矩阵的宽,H-特征矩阵的高;
S3.3.2:之后进行激活操作,计算每个特征图的重要性,通过全连接层和非线性层计算得到C个特征图的重要性,计算公式为:
S=Fex(Z,W)=sigmoid(W2×relu(W1×Z)) (10)
式(10)中:W1,W2-进行全连接操作的隐藏层矩阵;Z-式(9)中提取的全局信息;
S3.3.3:得到重要性权重S,将权重乘回特征图可得到学习到重要性的特征信息X8(t);
X8(t)=S×X7(t) (11)
S3.4:通过批归一化层将X8(t)进行归一化处理,得到X9(t),使X9(t)符合标准正态分布,均值为0,标准差为1,加快收敛速度,防止梯度消失;
批归一化层的计算公式为:
Figure BDA0003304414880000041
式(12)中:γ,β-两个可学习重构参数;μ-均值;(σ)2+ε-方差;
S3.5:通过长短记忆神经网络对X9(t)进行时间序列数据关联;
计算公式为:
X10(t)=sigmoid(W×X9(t)+V×h(t-1)+U×c(t-1) (13)
式(13)中:h(t)-计算的传递状态;c(t)-选择记忆的状态;W,V,U-计算系数;
S3.6:最后使用多层感知机对进行X10(t)回归计算,得到腹内压力的数据y,
计算公式为:
y=W3×{W2×(W1×X10(t)+B1)+B2}+B3 (14)
式(14)中:W1,W2,W3,-计算系数;B1,B2,B3-计算偏置。
相比于现有技术,发明具有如下有益效果:
与传统的临床采用的接触式不同,毫米波雷达可不接触测试对象身体,隔一定距离探测,避免了穿刺测量导致的感染风险,并且提升了采集时的舒适性;本发明的关键点在于:使用雷达进行无接触式的生命体特征信号收集,对采集的生命体特征信号使用基于深度学习的神经网络进行计算,通过引入注意力机制,对多通道的腹部呼吸信号进行重要性计算,从不同重要性的通道自适应的提取出不同量的特征,提高该神经网络算法的能力,获取更多的腹部呼吸信号中的时域信息,将时域信息与测试对象身体腹内压力相关联,从而推算出测试对象身体的腹内压力;相对与传统侵入式或接触式的测量腹内压力的方法,该方法是一种全新的无接触式的测量腹内压力的方法。
附图说明
图1为发明实施例的逻辑路线示意图;
图2为发明实施例的实验数据图。
具体实施方式
下面结合附图及实施例对发明中的技术方案进一步说明。
一种基于深度学习的FMCW雷达提取分析腹内压力方法,包括下述步骤:
S1:将FMCW毫米波雷达系统(具有发射天线(TX),接收天线(RX)),放置在观测对象旁边,朝向观测对象腹部位置发射毫米级别的线性调频连续波信号,采集得到从观测对象腹部反射回来的雷达波;将反射回来的雷达波与原始发射信号相混频得到原始数据信号;
S2:所述原始数据信号通过快速傅里叶变换提取出包含有观测对象腹部呼吸信号的生命体征信号,再通过巴特沃斯滤波器滤去噪声;
S3:将包含有观测对象腹部呼吸信号的生命体征信号输入搭建好的神经网络中,对包含有观测对象腹部呼吸信号的生命体征信号进行特征提取,并通过回归分析与腹内压力进行关联,进而分析出腹内压力。
在S2中,将取得的原始数据信号定义为X1(t),对X1(t)进行相位特征
Figure BDA0003304414880000055
的提取,将X1(t) 的相位
Figure BDA0003304414880000056
定义为:X1(t)的虚部I与实部R之比;拟出公式(1):
Figure BDA0003304414880000051
原始数据信号X1(t)可以被分解为原始信号实部R(t),原始信号虚部I(t),原始数据信号 X1(t)还包含有实部的直流信号dcr,虚部的直流信号dci;将原始信号实部R(t),原始信号虚部I(t)、实部的直流信号dcr、虚部的直流信号dci代入公式(1),可以得到相位特征
Figure BDA0003304414880000057
的计算公式(2):
Figure BDA0003304414880000052
在S2中,将得到的
Figure BDA0003304414880000058
重新定义为:包含腹部呼吸的时域信号X2(t),而后通过快速傅里叶变换将X2(t)变换为包含腹部呼吸的频域信号X3(k),对X2(t)的频域进行分析,计算公式为:
Figure BDA0003304414880000053
式(3)中:k=0,1,......N-1;
Figure BDA0003304414880000054
得到包含腹部呼吸的频域信号X3(k)后,(X3(k)在率噪前还含有许多噪声)由于实际腹部呼吸信号频率范围为0.2Hz~0.4Hz,在X3(k)中观察到符合实际范围的腹部呼吸频域信号X4(k),此时要对X4(k)进行滤噪,具体做法为:将X3(k)与设计的巴特沃斯滤波器频域相乘,X4(k)进行滤噪;
巴特沃斯滤波器公式为:
|Ha(jω)|2=1/1+(ω/ωc)2N (4)
式(4)中:N-滤波器阶数,ωc-滤波器的截止频率,Ha(jω)-设计的滤波器幅频响应;
滤噪腹部呼吸频域信号X4(k)计算公式为:
X4(k)=X3(k)*Ha(jω) (5)
式(5)中:Ha(jω)-设计的巴特沃斯滤波器幅频响应;
滤噪后的X4(k)即为纯净的腹部呼吸频域信号,再通过快速傅里叶逆变换将滤噪后的腹部呼吸频域信号X4(k)变换为腹部呼吸时域信号X4(t)(X4(t)为四通道信号);
计算公式为:
Figure BDA0003304414880000061
式(6)中:k=0,1,......N-1;
Figure BDA0003304414880000062
在S3中:将获取的腹部时域呼吸信号X4(t)输入到预先搭建的深度学习神经网络(具体为卷积神经网络(CNN)),进行进一步的特征提取以及信号映射到腹内压力的分析,其具体步骤为:
S3.1:将腹部时域呼吸信号X4(t)通过卷积神经网络进行一次特征提取,得到提取的特征定义为信号X5(t);X4(t)为四通道二维矩阵,设进行卷积提取特征的卷积核矩阵为k,阶矩阵为n×n;
提取特征X5(t)的计算公式为:
Figure BDA0003304414880000063
式(7)中:
Figure BDA0003304414880000064
ki,j-卷积核;
S3.2:再通过线性整流函数加入非线性因素,增强模型的表达力,如果没有该函数,每层的神经网络计算都相当于矩阵相乘,都为原数据的线性组合,无法表征出比较复杂的情况,而复杂的模型都为非线性的,加入了非线性因素使的神经网络能更好的逼近实际的模型;
线性整流函数公式为:
X6(t)=max(0,X5(t)) (8)
S3.3:对S3.1中的卷积神经网络以及S3.2的线性整流函数进行五层堆叠计算(即将式 (7)与式(8)堆叠计算5次),目的是提取四个通道的数据特征,之后引入非线性因素,得到特征信号X7(t);此时,设X7(t)中大小为H*W的矩阵为C个;加入压缩-激活神经网络(squeeze and excitation networks)引入注意力机制,进行对每个通道的特征权重计算操作,依照计算的重要程度去提升有用的特征并抑制对当前任务用处不大的特征;压缩-激活神经网络首先需要进行对C个特征图进行压缩操作,获取上述提取特征图的全局信息,其具体步骤为:
S3.3.1:压缩-激活神经网络首先需要进行对C个特征图进行压缩操作,获取上述提取特征图的全局信息;
计算公式为:
Figure BDA0003304414880000071
式(9)中:W-特征矩阵的宽,H-特征矩阵的高;
S3.3.2:之后进行激活操作,计算每个特征图的重要性,通过全连接层和非线性层计算得到C个特征图的重要性,计算公式为:
S=Fex(Z,W)=sigmoid(W2×relu(W1×Z)) (10)
式(10)中:W1,W2-进行全连接操作的隐藏层矩阵;Z-式(9)中提取的全局信息;
S3.3.3:得到重要性权重S,将权重乘回特征图可得到学习到重要性的特征信息X8(t);
X8(t)=S×X7(t) (11)
S3.4:通过批归一化层将X8(t)进行归一化处理,得到X9(t),使X9(t)符合标准正态分布,均值为0,标准差为1,加快收敛速度,防止梯度消失,当数据输出较大时,学习梯度很小,浅层网络基本不学习,只有较后层的网络表征整个网络,失去了深度的意义,而加入了批归一化层后,网络输出不会很大,即梯度不会很小;
批归一化层的计算公式为:
Figure BDA0003304414880000072
式(12)中:γ,β-两个可学习重构参数;μ-均值;(σ)2+ε-方差;
S3.5:通过长短记忆神经网络对X9(t)进行时间序列数据关联;
计算公式为:
X10(t)=sigmoid(W×X9(t)+V×h(t-1)+U×c(t-1) (13)
式(13)中:h(t)-计算的传递状态;c(t)-选择记忆的状态;W,V,U-计算系数;
S3.6:最后使用多层感知机对进行X10(t)回归计算,得到腹内压力的数据y,
计算公式为:
y=W3×{W2×(W1×X10(t)+B1)+B2}+B3 (14)
式(14)中:W1,W2,W3-计算系数;B1,B2,B3-计算偏置;
接下来对照实验组进一步说明本发明的实际效果:
参数设定:
雷达频率:60-64ghz;
每次测量时间:1min;
雷达测量位置:测试对象的腹部;
雷达放置距离:垂直距离测量对象位置50cm;
雷达天线数量:4根接收(RX)天线,3根发射(TX)天线;
雷达扫描范围:120°方位角视野,30°仰角视野;
对照组测量方式:直接测压(将传感器插入测量对象膀胱)
折线图中,横坐标为时间h,纵坐标为腹压cmH2O,▲:雷达测量后的计算值;●:雷达测量后的计算值的平均值;▼:临床测量值;
正三角形的标点为对一次测量的数据所进行分段处理得到的预测值,由于腹部压力会随着呼吸而变换,临床(选用侵入式)测量的时间为一次呼吸的末尾,本发明测量的一次数据为1分钟,由于不知道1分钟内的何时为呼吸末尾,因此将1分钟数据分为多段进行数据预测,并将多次预测的值求计算后的平均值作为预测值,即圆形标点折线,临床测量值折线为倒三角形标点折线;
由图(2)可以看出,随着时间的推移,1h后,测试对象腹部会出现明显感染导致腹压升高,预测值折线与临床测量值折线升高趋势一致,两者最高相差大约1cm,逐渐趋近。
最后说明的是,以上实施例仅用以说明发明的技术方案而非限制,尽管参照较佳实施例对发明进行了详细说明,本领域的普通技术人员应当理解,可以对发明的技术方案进行修改或者等同替换,而不脱离发明技术方案的宗旨和范围,其均应涵盖在发明的权利要求范围当中。

Claims (4)

1.一种基于深度学习的FMCW雷达提取分析腹内压力方法,其特征在于,包括下述步骤:
S1:将FMCW毫米波雷达系统,放置在观测对象旁边,朝向观测对象腹部位置发射毫米级别的线性调频连续波信号,采集得到从观测对象腹部反射回来的雷达波;将反射回来的雷达波与原始发射信号相混频得到原始数据信号;
S2:所述原始数据信号通过快速傅里叶变换提取出包含有观测对象腹部呼吸信号的生命体征信号,再通过巴特沃斯滤波器滤去噪声;
S3:将包含有观测对象腹部呼吸信号的生命体征信号输入搭建好的神经网络中,对包含有观测对象腹部呼吸信号的生命体征信号进行特征提取,并通过回归分析与腹内压力进行关联,进而分析出腹内压力。
2.如权利要求1所述的一种基于深度学习的FMCW雷达提取分析腹内压力方法,其特征在于,在S2中,将取得的原始数据信号定义为X1(t),对X1(t)进行相位特征
Figure FDA0003304414870000011
的提取,将X1(t)的相位
Figure FDA0003304414870000012
定义为:X1(t)的虚部I与实部R之比;拟出公式(1):
Figure FDA0003304414870000013
原始数据信号X1(t)可以被分解为原始信号实部R(t),原始信号虚部I(t),原始数据信号X1(t)还包含有实部的直流信号dcr,虚部的直流信号dci;将原始信号实部R(t),原始信号虚部I(t)、实部的直流信号dcr、虚部的直流信号dci代入公式(1),可以得到相位特征
Figure FDA0003304414870000014
的计算公式(2):
Figure FDA0003304414870000015
3.如权利要求2所述的一种基于深度学习的FMCW雷达提取分析腹内压力方法,其特征在于,在S2中,将得到的
Figure FDA0003304414870000016
重新定义为:包含有腹部呼吸的时域信号X2(t);而后通过快速傅里叶变换将X2(t)变换为腹部呼吸的频域信号X3(k),对X2(t)的频域进行分析;计算公式为:
Figure FDA0003304414870000017
式(3)中:k=0,1,......N-1;
Figure FDA0003304414870000018
得到包含有腹部呼吸的频域信号X3(k)后,由于实际腹部呼吸信号频率范围为0.2Hz~0.4Hz,在X3(k)中观察到符合实际范围的腹部呼吸频域信号X4(k),将X3(k)与设计的巴特沃斯滤波器频域相乘,对X4(k)进行滤噪;
巴特沃斯滤波器公式为:
|Ha(jω)|2=1/1+(ω/ωc)2N (4)
式(4)中:N-滤波器阶数,ωc-滤波器的截止频率,Ha(jω)-设计的滤波器幅频响应;
滤噪腹部呼吸频域信号X4(k)计算公式为:
X4(k)=X3(k)*Ha(jω) (5)
式(5)中:Ha(jω)-设计的巴特沃斯滤波器幅频响应;
再通过快速傅里叶逆变换将滤噪后的腹部呼吸频域信号X4(k)变换为腹部呼吸时域信号X4(t);
计算公式为:
Figure FDA0003304414870000021
式(6)中:k=0,1,......N-1;
Figure FDA0003304414870000022
4.如权利要求3所述的一种基于深度学习的FMCW雷达提取分析腹内压力方法,其特征在于,在S3中:将获取的腹部时域呼吸信号X4(t)输入到预先搭建的深度学习神经网络,进行进一步的特征提取以及信号映射到腹内压力的分析,其具体步骤为:
S3.1:将腹部时域呼吸信号X4(t)通过卷积神经网络进行一次特征提取,得到提取的特征定义为信号X5(t);X4(t)为输入的四通道二维矩阵,设进行卷积提取特征的卷积核矩阵为k,阶矩阵为n×n;
提取特征X5(t)的计算公式为:
Figure FDA0003304414870000023
式(7)中:
Figure FDA0003304414870000024
ki,j-卷积核;
S3.2:再通过线性整流函数加入非线性因素,增强模型表达力,
线性整流函数公式为:
X6(t)=max(0,X5(t)) (8)
S3.3:对S3.1中的卷积神经网络以及S3.2的线性整流函数进行五层堆叠计算,目的是提取四个通道的数据特征,之后引入非线性因素,得到特征信号X7(t);此时,设X7(t)中大小为H*W的矩阵为C个;加入压缩-激活神经网络,引入注意力机制,对每个通道进行特征权重的计算,其具体步骤为:
S3.3.1:压缩-激活神经网络首先需要进行对C个特征图进行压缩操作,获取上述提取特征图的全局信息;
计算公式为:
Figure FDA0003304414870000025
式(9)中:W-特征矩阵的宽,H-特征矩阵的高;
S3.3.2:之后进行激活操作,计算每个特征图的重要性,通过全连接层和非线性层计算得到C个特征图的重要性,计算公式为:
S=Fex(Z,W)=sigmoid(W2×relu(W1×Z)) (10)
式(10)中:W1,W2-进行全连接操作的隐藏层矩阵;Z-式(9)中提取的全局信息;
S3.3.3:得到重要性权重S,将权重乘回特征图可得到学习到重要性的特征信息X8(t);
X8(t)=S×X7(t) (11)
S3.4:通过批归一化层将X8(t)进行归一化处理,得到X9(t),使X9(t)符合标准正态分布,均值为0,标准差为1,加快收敛速度,防止梯度消失;
批归一化层的计算公式为:
Figure FDA0003304414870000031
式(12)中:γ,β-两个可学习重构参数;μ-均值;(σ)2+ε-方差;
S3.5:通过长短记忆神经网络对X9(t)进行时间序列数据关联;
计算公式为:
X10(t)=sigmoid(W×X9(t)+V×h(t-1)+U×c(t-1) (13)
式(13)中:h(t)-计算的传递状态;c(t)-选择记忆的状态;W,V,U-计算系数;
S3.6:最后使用多层感知机对进行X10(t)回归计算,得到腹内压力的数据y,
计算公式为:
y=W3×{W2×(W1×x10(t)+B1)+B2}+B3 (14)
式(14)中:W1,W2,W3-计算系数;B1,B2,B3-计算偏置。
CN202111199566.8A 2021-10-14 2021-10-14 一种基于深度学习的fmcw雷达提取分析腹内压力方法 Active CN113892931B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111199566.8A CN113892931B (zh) 2021-10-14 2021-10-14 一种基于深度学习的fmcw雷达提取分析腹内压力方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111199566.8A CN113892931B (zh) 2021-10-14 2021-10-14 一种基于深度学习的fmcw雷达提取分析腹内压力方法

Publications (2)

Publication Number Publication Date
CN113892931A true CN113892931A (zh) 2022-01-07
CN113892931B CN113892931B (zh) 2023-08-22

Family

ID=79192067

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111199566.8A Active CN113892931B (zh) 2021-10-14 2021-10-14 一种基于深度学习的fmcw雷达提取分析腹内压力方法

Country Status (1)

Country Link
CN (1) CN113892931B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115381467A (zh) * 2022-10-31 2022-11-25 浙江浙大西投脑机智能科技有限公司 一种基于注意力机制的时频信息动态融合解码方法及装置

Citations (26)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050163349A1 (en) * 2003-06-06 2005-07-28 Daniela Brunner System and method for assessing motor and locomotor deficits and recovery therefrom
US20080001735A1 (en) * 2006-06-30 2008-01-03 Bao Tran Mesh network personal emergency response appliance
US20130157881A1 (en) * 2010-02-05 2013-06-20 Astute Medical, Inc. Methods and compositions for diagnosis and prognosis of renal injury and renal failure
JP2014207935A (ja) * 2013-04-16 2014-11-06 富士通株式会社 生体情報取得装置、方法及びプログラム
US20150164375A1 (en) * 2012-05-30 2015-06-18 Resmed Sensor Technologies Limited Method and apparatus for monitoring cardio-pulmonary health
US20150237834A1 (en) * 2014-02-24 2015-08-27 Protequus LLC Mobile animal surveillance and distress monitoring
US9610443B1 (en) * 2013-02-27 2017-04-04 Steven G Dean Methods to trigger high amplitude oscillations or resonance in the cardiovascular system of a patient using electrical stimulation
CN110113988A (zh) * 2016-11-17 2019-08-09 美第奇科技有限责任公司 用于中心静脉压测量状态的方法和设备
CN110638482A (zh) * 2019-10-18 2020-01-03 北京大学第三医院(北京大学第三临床医学院) 一种肠鸣音和腹腔压力的实时监护系统及方法
US20200022637A1 (en) * 2017-03-24 2020-01-23 The Regents Of The University Of California System, device and method for bladder volume sensing
CN111603138A (zh) * 2020-05-19 2020-09-01 杭州电子科技大学 一种基于毫米波雷达的睡眠呼吸暂停监测系统
CN111965636A (zh) * 2020-07-20 2020-11-20 重庆大学 一种基于毫米波雷达和视觉融合的夜间目标检测方法
US20200367810A1 (en) * 2017-12-22 2020-11-26 Resmed Sensor Technologies Limited Apparatus, system, and method for health and medical sensing
JP2021023615A (ja) * 2019-08-06 2021-02-22 学校法人 関西大学 脈圧推定装置、脈圧推定システム、脈圧推定方法、及び制御プログラム
US20210089131A1 (en) * 2019-09-23 2021-03-25 Apple Inc. Finger-Mounted Input Devices
US20210093203A1 (en) * 2019-09-30 2021-04-01 DawnLight Technologies Systems and methods of determining heart-rate and respiratory rate from a radar signal using machine learning methods
US20210106275A1 (en) * 2019-10-11 2021-04-15 Boe Technology Group Co., Ltd. Method, device and storage medium for sleep monitoring
CN112754431A (zh) * 2020-12-31 2021-05-07 杭州电子科技大学 基于毫米波雷达和轻量级神经网络的呼吸心跳监测系统
CN112998701A (zh) * 2021-03-27 2021-06-22 复旦大学 基于毫米波雷达的生命体征检测与身份识别系统和方法
CN113143270A (zh) * 2020-12-02 2021-07-23 长春理工大学 一种基于生物雷达和语音信息的双模态融合情绪识别方法
US20210353195A1 (en) * 2018-10-30 2021-11-18 Renovia Inc. Devices, systems, and methods for monitoring bladder function
CN114246563A (zh) * 2021-12-17 2022-03-29 重庆大学 基于毫米波雷达的心肺功能智能监测设备
US20220133222A1 (en) * 2019-02-19 2022-05-05 Koninklijke Philips N.V. A sleep monitoring system and method
US20220404483A1 (en) * 2020-01-07 2022-12-22 The Regents Of The University Of California Integrated wideband stepped-chirp radar sensor
CN115758122A (zh) * 2022-11-22 2023-03-07 江西朴拙医疗设备有限公司 基于多尺度卷积神经网络的睡眠呼吸事件定位方法及装置
CN115736888A (zh) * 2022-11-24 2023-03-07 南京邮电大学 基于fmcw雷达的呼吸信号提取方法

Patent Citations (26)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050163349A1 (en) * 2003-06-06 2005-07-28 Daniela Brunner System and method for assessing motor and locomotor deficits and recovery therefrom
US20080001735A1 (en) * 2006-06-30 2008-01-03 Bao Tran Mesh network personal emergency response appliance
US20130157881A1 (en) * 2010-02-05 2013-06-20 Astute Medical, Inc. Methods and compositions for diagnosis and prognosis of renal injury and renal failure
US20150164375A1 (en) * 2012-05-30 2015-06-18 Resmed Sensor Technologies Limited Method and apparatus for monitoring cardio-pulmonary health
US9610443B1 (en) * 2013-02-27 2017-04-04 Steven G Dean Methods to trigger high amplitude oscillations or resonance in the cardiovascular system of a patient using electrical stimulation
JP2014207935A (ja) * 2013-04-16 2014-11-06 富士通株式会社 生体情報取得装置、方法及びプログラム
US20150237834A1 (en) * 2014-02-24 2015-08-27 Protequus LLC Mobile animal surveillance and distress monitoring
CN110113988A (zh) * 2016-11-17 2019-08-09 美第奇科技有限责任公司 用于中心静脉压测量状态的方法和设备
US20200022637A1 (en) * 2017-03-24 2020-01-23 The Regents Of The University Of California System, device and method for bladder volume sensing
US20200367810A1 (en) * 2017-12-22 2020-11-26 Resmed Sensor Technologies Limited Apparatus, system, and method for health and medical sensing
US20210353195A1 (en) * 2018-10-30 2021-11-18 Renovia Inc. Devices, systems, and methods for monitoring bladder function
US20220133222A1 (en) * 2019-02-19 2022-05-05 Koninklijke Philips N.V. A sleep monitoring system and method
JP2021023615A (ja) * 2019-08-06 2021-02-22 学校法人 関西大学 脈圧推定装置、脈圧推定システム、脈圧推定方法、及び制御プログラム
US20210089131A1 (en) * 2019-09-23 2021-03-25 Apple Inc. Finger-Mounted Input Devices
US20210093203A1 (en) * 2019-09-30 2021-04-01 DawnLight Technologies Systems and methods of determining heart-rate and respiratory rate from a radar signal using machine learning methods
US20210106275A1 (en) * 2019-10-11 2021-04-15 Boe Technology Group Co., Ltd. Method, device and storage medium for sleep monitoring
CN110638482A (zh) * 2019-10-18 2020-01-03 北京大学第三医院(北京大学第三临床医学院) 一种肠鸣音和腹腔压力的实时监护系统及方法
US20220404483A1 (en) * 2020-01-07 2022-12-22 The Regents Of The University Of California Integrated wideband stepped-chirp radar sensor
CN111603138A (zh) * 2020-05-19 2020-09-01 杭州电子科技大学 一种基于毫米波雷达的睡眠呼吸暂停监测系统
CN111965636A (zh) * 2020-07-20 2020-11-20 重庆大学 一种基于毫米波雷达和视觉融合的夜间目标检测方法
CN113143270A (zh) * 2020-12-02 2021-07-23 长春理工大学 一种基于生物雷达和语音信息的双模态融合情绪识别方法
CN112754431A (zh) * 2020-12-31 2021-05-07 杭州电子科技大学 基于毫米波雷达和轻量级神经网络的呼吸心跳监测系统
CN112998701A (zh) * 2021-03-27 2021-06-22 复旦大学 基于毫米波雷达的生命体征检测与身份识别系统和方法
CN114246563A (zh) * 2021-12-17 2022-03-29 重庆大学 基于毫米波雷达的心肺功能智能监测设备
CN115758122A (zh) * 2022-11-22 2023-03-07 江西朴拙医疗设备有限公司 基于多尺度卷积神经网络的睡眠呼吸事件定位方法及装置
CN115736888A (zh) * 2022-11-24 2023-03-07 南京邮电大学 基于fmcw雷达的呼吸信号提取方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
SHARANYA MANGA 等: "Estimation of physiologic pressures:invasive and non-invasive techniques, AI models, and future perspectives", 《SENSORS》, vol. 23, no. 2, pages 5744 *
TANG H*: "Deep domain adaptation ofr predicting intra-abdominal pressure with multichannel attention fusion radar", 《ADCANCED INTELLIGENT SYSTEMS》, vol. 4, no. 5, pages 2100209 *
刘声 等: "双基地MIMO雷达的快速测向算法", 吉林大学学报(工学版), vol. 46, no. 6, pages 1675 - 1680 *
康瑞;杨千朋;: "大黄红藤消炎汤配合超声引导下经皮穿刺置管引流术治疗急性胰腺炎合并胰腺周围组织坏死感染的研究", 现代中西医结合杂志, no. 34, pages 65 - 69 *
韩伟奇 等: "丙泊酚对腹腔镜手术气腹应激反应的影响", 《中国社区医师》, no. 18, pages 46 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115381467A (zh) * 2022-10-31 2022-11-25 浙江浙大西投脑机智能科技有限公司 一种基于注意力机制的时频信息动态融合解码方法及装置

Also Published As

Publication number Publication date
CN113892931B (zh) 2023-08-22

Similar Documents

Publication Publication Date Title
US10825569B2 (en) Universal non-invasive blood glucose estimation method based on time series analysis
CN103385702B (zh) 一种无创血压连续检测装置及方法
CN104173043B (zh) 一种适合于移动平台的心电数据分析方法
EP0327645B1 (en) Apparatus for measuring intracranial pressure
CN111568399A (zh) 一种基于雷达的呼吸和心跳信号检测方法及系统
CN110464320A (zh) 多目标人体心率和呼吸频率测量系统及方法
CN104000573B (zh) 基于体表两点脉搏波的中心动脉脉搏波监测系统及方法
CN105930665B (zh) 一种基于超声衰减系数的肝脏脂肪定量测量方法
CN102499669B (zh) 心脏参数的测量装置
CN110742593B (zh) 一种基于线谱跟踪的生命信号特征提取方法
CN112806977B (zh) 基于多尺度融合网络的生理参数测量方法
CN108175384A (zh) 基于子宫肌电信号识别宫缩的方法及装置
CN104644151B (zh) 一种基于光电容积脉搏信号的压力脉搏波波形传播预测方法
CN113892931A (zh) 一种基于深度学习的fmcw雷达提取分析腹内压力方法
CN111839488B (zh) 基于脉搏波的无创连续血压测量装置和方法
CN108670238A (zh) 基于脑氧信号与血压信号的脑血流自主调节能力评价方法
CN111887828B (zh) 围术期患者非接触式生理信息监测装置、计算机设备和存储介质
CN114652351A (zh) 基于超声多普勒的连续血压测量方法、装置和电子设备
CN115736888A (zh) 基于fmcw雷达的呼吸信号提取方法
CN102579008A (zh) 用于监测脑积水和脑水肿的装置及其监测方法
CN117958761A (zh) 一种基于毫米波雷达的睡眠呼吸暂停检测方法及系统
CN114145725B (zh) 一种基于无创连续血压测量的ppg采样率估算方法
CN116725507A (zh) 基于拟合匹配滤波与Double-CZT的雷达高精度心率检测方法
CN108836316B (zh) 一种基于bp神经网络的心电信号r波提取方法
Shi et al. Systolic blood pressure measurement algorithm with mmWave radar sensor

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