CN113786177A - 生命体征信息提取方法、装置及电子设备 - Google Patents
生命体征信息提取方法、装置及电子设备 Download PDFInfo
- Publication number
- CN113786177A CN113786177A CN202110974749.6A CN202110974749A CN113786177A CN 113786177 A CN113786177 A CN 113786177A CN 202110974749 A CN202110974749 A CN 202110974749A CN 113786177 A CN113786177 A CN 113786177A
- Authority
- CN
- China
- Prior art keywords
- signal
- signals
- radar echo
- slice
- radar
- 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
- 238000000605 extraction Methods 0.000 title claims abstract description 43
- 238000000034 method Methods 0.000 claims abstract description 43
- 238000012545 processing Methods 0.000 claims abstract description 40
- 230000000241 respiratory effect Effects 0.000 claims abstract description 22
- 238000009825 accumulation Methods 0.000 claims abstract description 16
- 230000001427 coherent effect Effects 0.000 claims abstract description 16
- 238000005070 sampling Methods 0.000 claims description 21
- 238000001514 detection method Methods 0.000 claims description 17
- 238000012937 correction Methods 0.000 claims description 11
- 230000001186 cumulative effect Effects 0.000 claims description 10
- 238000004590 computer program Methods 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 7
- 238000001914 filtration Methods 0.000 claims description 6
- 230000005764 inhibitory process Effects 0.000 claims description 5
- 230000008030 elimination Effects 0.000 claims description 4
- 238000003379 elimination reaction Methods 0.000 claims description 4
- 238000003672 processing method Methods 0.000 claims description 4
- 230000006870 function Effects 0.000 description 7
- 230000029058 respiratory gaseous exchange Effects 0.000 description 7
- 239000011159 matrix material Substances 0.000 description 6
- 230000008859 change Effects 0.000 description 4
- 230000003068 static effect Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000000737 periodic effect Effects 0.000 description 3
- 239000000463 material Substances 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000035515 penetration Effects 0.000 description 2
- 210000000115 thoracic cavity Anatomy 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001629 suppression Effects 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/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/0205—Simultaneously evaluating both cardiovascular conditions and different types of body conditions, e.g. heart and respiratory condition
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/024—Detecting, measuring or recording pulse rate or heart rate
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/08—Detecting, measuring or recording devices for evaluating the respiratory organs
-
- 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/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/725—Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Public Health (AREA)
- Molecular Biology (AREA)
- Veterinary Medicine (AREA)
- General Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Animal Behavior & Ethology (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Physiology (AREA)
- Surgery (AREA)
- Signal Processing (AREA)
- Cardiology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Artificial Intelligence (AREA)
- Psychiatry (AREA)
- Pulmonology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Radar Systems Or Details Thereof (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Abstract
本申请提出一种生命体征信息提取方法、装置及电子设备，其中，方法包括：采用超宽带雷达探测得到雷达回波信号；确定所述雷达回波信号中慢时间方向的各切片信号对应的信号能量；基于预设宽度的能量窗以及所述各切片信号对应的信号能量，从所述雷达回波信号中提取人体体征信号；对所述人体体征信号进行相干累积处理，并从处理后的人体体征信号中提取呼吸信号和心跳信号。由此，可以实现根据雷达回波信号中慢时间方向的各切片信号对应的信号能量以及预设宽度的能量窗，从雷达回波信号中准确提取呼吸信号和心跳信号。
Description
技术领域
本申请涉及信号处理技术领域，尤其涉及一种生命体征信息提取方法、装置及电子设备。
背景技术
超宽带(Ultra-Wideband Radar，简称UWB)雷达具有距离分辨率高、穿透能力强、功耗低、抗干扰能力强等优点，因此可以非接触式、远距离检测人体对象的位置、呼吸和心跳等信息。
相关技术，由于复杂的噪声干扰或者其它外部原因导致的人体对象的身体小范围移动，导致很难准确提取人体对象的呼吸和心跳等生命体征信息。
发明内容
本申请旨在至少在一定程度上解决相关技术中的技术问题之一。
本申请提出一种生命体征信息提取方法及其装置，以实现根据雷达回波信号中慢时间方向的各切片信号对应的信号能量以及预设宽度的能量窗，从雷达回波信号中准确提取呼吸信号和心跳信号。
本申请第一方面实施例提出了一种生命体征信息提取方法，包括：采用超宽带雷达探测得到雷达回波信号；确定所述雷达回波信号中慢时间方向的各切片信号对应的信号能量；基于预设宽度的能量窗以及所述各切片信号对应的信号能量，从所述雷达回波信号中提取人体体征信号；对所述人体体征信号进行相干累积处理，并从处理后的人体体征信号中提取呼吸信号和心跳信号。
本申请实施例的生命体征信息提取方法，通过采用超宽带雷达探测得到雷达回波信号，确定雷达回波信号中慢时间方向的各切片信号对应的信号能量，基于预设宽度的能量窗以及各切片信号对应的信号能量，从雷达回波信号中提取人体体征信号，进而对人体体征信号进行相干累积处理，并从处理后的人体体征信号中提取呼吸信号和心跳信号。由此，可以实现根据雷达回波信号中慢时间方向的各切片信号对应的信号能量以及预设宽度的能量窗，从雷达回波信号中准确提取呼吸信号和心跳信号。
本申请第二方面实施例提出了一种生命体征信息提取装置，包括：探测模块，用于采用超宽带雷达探测得到雷达回波信号；确定模块，用于确定所述雷达回波信号中慢时间方向的各切片信号对应的信号能量；提取模块，用于基于预设宽度的能量窗以及所述各切片信号对应的信号能量，从所述雷达回波信号中提取人体体征信号；处理模块，用于对所述人体体征信号进行相干累积处理，并从处理后的人体体征信号中提取呼吸信号和心跳信号。
本申请实施例的生命体征信息提取装置，通过采用超宽带雷达探测得到雷达回波信号，确定雷达回波信号中慢时间方向的各切片信号对应的信号能量，基于预设宽度的能量窗以及各切片信号对应的信号能量，从雷达回波信号中提取人体体征信号，进而对人体体征信号进行相干累积处理，并从处理后的人体体征信号中提取呼吸信号和心跳信号。由此，可以实现根据雷达回波信号中慢时间方向的各切片信号对应的信号能量以及预设宽度的能量窗，从雷达回波信号中准确提取呼吸信号和心跳信号。
本申请第三方面实施例提出了一种电子设备，包括：存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时，实现如本申请第一方面实施例所述的方法。
本申请第四方面实施例提出了一种存储有计算机指令的非瞬时计算机可读存储介质，所述计算机指令用于使所述计算机执行本申请第一方面实施例所述的方法。
本申请第五方面实施例提出了一种计算机程序产品，包括计算机程序，所述计算机程序在被处理器执行时实现本申请第一方面实施例所述的方法。
本申请附加的方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本申请的实践了解到。
附图说明
本申请上述的和/或附加的方面和优点从下面结合附图对实施例的描述中将变得明显和容易理解，其中：
图1为本申请实施例所提供的一种生命体征信息提取方法的流程示意图；
图2为本申请实施例所提供的另一种生命体征信息提取方法的流程示意图；
图3为本申请实施例所提供的另一种生命体征信息提取方法的流程示意图；
图4为本申请实施例提供的一种生命体征信息提取装置的结构示意图；
图5是用来实现本公开实施例的生命体征信息提取方法的电子设备的框图。
具体实施方式
下面详细描述本申请的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，旨在用于解释本申请，而不能理解为对本申请的限制。
超宽带雷达具有距离分辨率高、穿透能力强、功耗低、抗干扰能力强等优点，因此可以非接触式、远距离检测人体的位置、呼吸和心跳等信息。
相关技术，可以利用雷达回波信号中慢时间方向的单一切片信号，获取人体对象的生命体征信息，但是，由于复杂的噪声干扰，慢时间方向的单一切片信号包含的生命体征信息较弱，很难被提取到，另外，由于其它外部原因导致人体对象的身体小范围移动时，会使慢时间方向的单一切片信号含有的生命体征信息缺失，从而会使得生命体征信息更难以提取。因此，相关技术中的生命体征信息提取方法无法准确提取人体对象的生命体征信息。
本申请实施例为了准确提取人体对象的呼吸和心跳等生命体征信息，提出一种生命体征信息提取方法、装置、电子设备、存储介质以及计算机程序产品。其中，生命体征信息提取方法，包括：采用超宽带雷达探测得到雷达回波信号，确定雷达回波信号中慢时间方向的各切片信号对应的信号能量，基于预设宽度的能量窗以及各切片信号对应的信号能量，从雷达回波信号中提取人体体征信号，进而对人体体征信号进行相干累积处理，并从处理后的人体体征信号中提取呼吸信号和心跳信号。由此，实现了根据雷达回波信号中慢时间方向的各切片信号对应的信号能量以及预设宽度的能量窗，从雷达回波信号中准确提取呼吸信号和心跳信号。
下面参考附图描述本申请实施例的生命体征信息提取方法及其装置。
图1为本申请实施例所提供的一种生命体征信息提取方法的流程示意图。
本申请实施例以该生命体征信息提取方法被配置于生命体征信息提取装置中来举例说明，该生命体征信息提取装置可以应用于任一具有计算能力的设备中，以使该设备可以执行生命体征信息提取功能。比如，电子设备可以为超宽带雷达生命探测仪。
如图1所示，该生命体征信息提取方法可以包括以下步骤：
步骤101，采用超宽带雷达探测得到雷达回波信号。
在本申请实施例中，可以通过超宽带雷达向检测对象发送检测信号，并接收检测对象返回的反馈信号，该反馈信号即为雷达回波信号。其中，雷达回波信号，包括无生命的检测对象返回的反馈信号，以及有生命的检测对象比如人体对象返回的反馈信号。
步骤102，确定雷达回波信号中慢时间方向的各切片信号对应的信号能量。
其中，雷达回波信号的快时间方向的快时间采样点可以为n，其中，n＝1,2,…,N。其中，N为快时间方向的最大采样点数，N的值可以根据需要设置，比如可以设置为512。
雷达回波信号的慢时间方向的慢时间采样点可以为m，其中，m＝1,2,…,M。其中，M为慢时间方向的最大采样点数，M的值可以根据需要设置。
雷达回波信号中慢时间方向的各切片信号，为雷达回波信号中快时间方向的各快时间采样点对应的信号。即，慢时间方向的每个切片信号，包括某快时间采样点下的各慢时间采样点对应的信号值。
在示例性实施例中，对于雷达回波信号中慢时间方向的每个切片信号，可以通过以下公式(1)，确定切片信号对应的信号能量。
其中，Ei为慢时间方向的第i个切片信号对应的信号能量，i＝1,2,...,N。fi(t)为雷达回波信号的慢时间方向的第i个切片信号，该切片信号为慢时间时刻t的函数。
步骤103，基于预设宽度的能量窗以及各切片信号对应的信号能量，从雷达回波信号中提取人体体征信号。
其中，能量窗的预设宽度，可以根据需要设置。比如，由于生命体征信息宽度一般为4-6个快时间采样点，则可以设置能量窗的预设宽度为5，从而利用宽度为5的能量窗进行人体对象的定位，从雷达回波信号中提取人体体征信号，以保证提取出的人体体征信号中包含生命体征信息。
在示例性实施例中，可以通过以下方法，基于预设宽度的能量窗以及各切片信号对应的信号能量，从雷达回波信号中提取人体体征信号：
将各切片信号分组置于能量窗内，并分别确定能量窗内的多个切片信号组的信号能量；其中，切片信号组包括雷达回波信号中慢时间方向的预设数量的相邻切片信号；
从多个切片信号组中，提取对应信号能量满足第一预设条件的目标切片信号组；
将目标切片信号组中的各切片信号作为人体体征信号。
其中，预设数量，与能量窗的宽度有关。比如，能量窗的宽度为5时，表示能量窗可以容纳5个快时间采样点对应的慢时间方向的切片信号，则切片信号组中包括的慢时间方向的切片信号的数量可以对应设置为5。
第一预设条件，可以根据需要设置。比如，可以设置第一预设条件为信号能量最大，即可以从多个切片信号组中，提取对应信号能量最大的目标切片信号组。
具体的，可以对雷达回波信号慢时间方向的各切片信号进行分组，每个切片信号组中包括雷达回波信号中慢时间方向的预设数量的相邻切片信号，并将分组后得到的多个切片信号组分别置于能量窗内，分别确定能量窗内的多个切片信号组的信号能量。比如，假设能量窗的宽度为5，切片信号组中包括的切片信号的数量为5，可以将快时间采样点n＝1，2，3，4，5分别对应的慢时间方向的切片信号置于能量窗内，确定这5个切片信号组成的切片信号组的信号能量，将快时间采样点n＝2，3，4，5，6分别对应的慢时间方向的切片信号置于能量窗内，确定这5个切片信号组成的切片信号组的信号能量，将快时间采样点n＝3，4，5，6，7分别对应的慢时间方向的切片信号置于能量窗内，确定这5个切片信号组成的切片信号组的信号能量，直至确定所有切片信号组的信号能量。
进而，可以对多个切片信号组的信号能量进行比较，以从多个切片信号组中，提取对应信号能量最大的目标切片信号组，将目标切片信号组中的各切片信号作为人体体征信号。
通过基于预设宽度的能量窗以及雷达回波信号慢时间方向的各切片信号对应的信号能量，可以自适应、准确的提取出静止人体对象以及微动人体对象的人体体征信号。
步骤104，对人体体征信号进行相干累积处理，并从处理后的人体体征信号中提取呼吸信号和心跳信号。
在示例性实施例中，可以对人体体征信号进行相干累积处理，使信号同相叠加以提高信噪比，完成相干累积处理后，可以对处理后的人体体征信号进行频谱分析，以提取呼吸信号和心跳信号。
其中，由于人体心跳引起的胸腔幅度变化约为0.3-0.7mm(毫米)，频率为0.9-1.6Hz(赫兹)，呼吸引起的胸腔幅度变化约为10-15mm，频率为0.15-0.5Hz，则可以从人体体征信号中提取幅度变化满足10-15mm、频率满足0.15-0.5Hz的信号作为呼吸信号，从人体体征信号中提取幅度变化满足0.3-0.7mm、频率满足0.9-1.6Hz的信号作为心跳信号。
本申请实施例的生命体征信息提取方法，通过采用超宽带雷达探测得到雷达回波信号，确定雷达回波信号中慢时间方向的各切片信号对应的信号能量，基于预设宽度的能量窗以及各切片信号对应的信号能量，从雷达回波信号中提取人体体征信号，进而对人体体征信号进行相干累积处理，并从处理后的人体体征信号中提取呼吸信号和心跳信号。由此，实现了根据雷达回波信号中慢时间方向的各切片信号对应的信号能量以及预设宽度的能量窗，从雷达回波信号中准确提取呼吸信号和心跳信号。
在示例性实施例中，在采用超宽带雷达探测得到雷达回波信号后，还可以先对雷达回波信号进行预处理，以去除雷达回波信号中的静态杂波、线性干扰以及随机噪声等，再利用处理后的雷达回波信号进行生命体征信息提取，以更准确的提取生命体征信息。下面结合图2，对本申请实施例所提供的生命体征信息提取方法进行进一步说明。
图2为本申请实施例所提供的另一种生命体征信息提取方法的流程示意图。
如图2所示，该生命体征信息提取方法可以包括以下步骤：
步骤201，采用超宽带雷达探测得到雷达回波信号。
在本申请实施例中，可以通过超宽带雷达向检测对象发送检测信号，并接收检测对象返回的反馈信号，该反馈信号即为雷达回波信号。
步骤202，对雷达回波信号进行数值偏差矫正处理，以矫正雷达回波信号中因雷达不稳定性导致的数值偏差。
在示例性实施例中，为解决雷达不稳定性导致的雷达回波信号中的数值偏差，可以对雷达回波信号进行数值偏差矫正处理。具体的，步骤202可以采用图3中的步骤301-305所示的方式实现：
步骤301，获取雷达回波信号的慢时间方向的各切片信号的累计概率分布。
步骤302，根据各切片信号的累计概率分布，将累计概率最大的切片信号确定为参考切片信号。
在示例性实施例中，可以获取雷达回波信号对应的信号矩阵R(m,n)，并获取信号矩阵R(m,n)的慢时间方向的各切片信号的累计概率分布，其中，对应的累计概率最大的慢时间方向的切片信号r(m,1)，为慢时间方向的最稳定的静止杂波切片信号，本申请实施例中，可以将对应的累计概率最大的慢时间方向的切片信号r(m,1)确定为参考切片信号，将其作为雷达稳定性参考。
步骤303，根据参考切片信号的累计概率分布，确定参考切片信号中概率满足第二预设条件的幅值。
步骤304，确定幅值对应的慢时间方向上的慢时间采样点。
其中，第二预设条件，可以根据需要设置，比如可以为概率小于0.5％。
在示例性实施例中，可以确定参考切片信号的累计概率分布，并认为在该参考切片信号内出现的概率满足第二预设条件的幅值为雷达不稳定性所致，进而确定雷达不稳定性所致的该幅值所在的慢时间方向上的慢时间采样点mx。
步骤305，对慢时间采样点对应的快时间方向的切片信号的幅值进行矫正，得到矫正后雷达回波信号。
在示例性实施例中，根据雷达不稳定性所致的幅值所在的慢时间方向上的慢时间采样点mx，可以定位mx位置的快时间方向的切片信号，并且，可以依据参考切片信号的变化趋势，确定稳定的快时间方向的切片信号mx-1，进而采用稳定的快时间方向的切片信号mx-1覆盖mx，直到完成所有偏差位置矫正，得到矫正后雷达回波信号对应的信号矩阵R1(m,n)。
步骤203，采用线性趋势抑制法，对雷达回波信号进行背景波消除，以消除雷达回波信号中的恒定成分。
其中，步骤203中的雷达回波信号，可以为步骤202进行数值偏差矫正处理后得到的矫正后雷达回波信号。
在示例性实施例中，雷达回波信号可能存在线性趋势问题，针对该问题，可以从慢时间方向对雷达回波信号进行线性趋势抑制，分别估计出慢时间方向的每个切片信号的线性趋势，并将线性趋势从原本慢时间方向的切片信号中减去，以对雷达回波信号进行背景波消除，消除雷达回波信号中的恒定成分，获得去除静态杂波和线性趋势项的雷达回波信号。
在示例性实施例中，可以采用下述公式(2)，确定消除恒定成分后的雷达回波信号：
其中，
是指消除恒定成分后的雷达回波信号对应的信号矩阵R2的转置，
是指未消除恒定成分之前的雷达回波信号对应的信号矩阵R1的转置，z＝[n”,1N]，n”＝[0,1,…,N-1]T，1N为大小为N×1维的向量。T表示矩阵的转置。
步骤204，对雷达回波信号进行信号增强处理，以增强雷达回波信号中的生命体征成分。
其中，步骤204中的雷达回波信号，可以为步骤203中消除雷达回波信号中的恒定成分后的雷达回波信号。
在实际应用中，因为噪声干扰以及多径效应等原因，会造成雷达回波信号中人体体征信号微弱，因此需对雷达回波信号进行增强处理，以增强雷达回波信号中的生命体征成分。
在示例性实施例中，可以采用巴特沃斯滤波方法，滤除雷达回波信号中的噪声成分，和/或，采用信号自相关处理方法，提高雷达回波信号的信噪比。
具体的，可以采用巴特沃斯滤波方法，设计带通滤波器，其中，带通频率可以为0.1Hz-5Hz，利用带通滤波器对雷达回波信号进行处理，以在保留雷达回波信号的原有频率特性基础上，滤除低频和高频干扰。
另外，由于经过上述处理的雷达回波信号主要包含生命体征信息以及零均值噪声，而生命体征信息具有一定的周期规律，零均值噪声不具备周期规律，可以通过自相关处理将零均值噪声减少到0，因此，本申请实施例中，可以采用信号自相关处理方法，提高雷达回波信号的信噪比。
具体的，可以采用以下公式(3)所示的方法，对雷达回波信号进行自相关处理：
R3(n')＝E[xm(n1)xm(n2)] (3)
其中，R3(n')＝E[xm(n1)xm(n2)]为进行自相关处理后的雷达回波信号，xm(n1)和xm(n2)分别为快时间采样点n1和n2对应的慢时间方向的切片信号，0≤n1，n2≤N-1，nˊ＝n1–n2。E为求平均。自相关处理，保留了慢时间方向的切片信号的周期分量和谐波分量。
通过上述过程，可以有效去除雷达回波信号中携带的复杂噪声干扰，提高雷达回波信号的信噪比。
需要说明的是，本申请实施例以对雷达回波信号进行数值偏差矫正处理、背景波消除以及信号增强处理为例进行说明，在实际应用中，可以仅进行上述三种处理中的任意一种或两种，本申请对此不作限制。
步骤205，确定雷达回波信号中慢时间方向的各切片信号对应的信号能量。
步骤206，基于预设宽度的能量窗以及各切片信号对应的信号能量，从雷达回波信号中提取人体体征信号。
在示例性实施例中，可以将各切片信号分组置于能量窗内，并分别确定能量窗内的多个切片信号组的信号能量，其中，切片信号组包括雷达回波信号中慢时间方向的预设数量的相邻切片信号，再从多个切片信号组中，提取对应信号能量满足第一预设条件的目标切片信号组，从而将目标切片信号组中的各切片信号作为人体体征信号。
步骤207，对人体体征信号进行相干累积处理，并从处理后的人体体征信号中提取呼吸信号和心跳信号。
其中，步骤205-207的具体实现过程及原理，可以参考上述实施例的描述，此处不再赘述。
本申请实施例的生命体征信息提取方法，通过采用超宽带雷达探测得到雷达回波信号，对雷达回波信号进行数值偏差矫正处理，以矫正雷达回波信号中因雷达不稳定性导致的数值偏差，采用线性趋势抑制法，对雷达回波信号进行背景波消除，以消除雷达回波信号中的恒定成分，对雷达回波信号进行信号增强处理，以增强雷达回波信号中的生命体征成分，确定雷达回波信号中慢时间方向的各切片信号对应的信号能量，基于预设宽度的能量窗以及各切片信号对应的信号能量，从雷达回波信号中提取人体体征信号，进而对人体体征信号进行相干累积处理，并从处理后的人体体征信号中提取呼吸信号和心跳信号。由此，有效的抑制了雷达回波信号中的复杂噪声干扰，提高了雷达回波信号的信噪比，且实现了根据雷达回波信号中慢时间方向的各切片信号对应的信号能量以及预设宽度的能量窗，自适应、准确地提取出静止人体对象和微动人体对象的呼吸信号和心跳信号。
与上述图1至图3实施例提供的生命体征信息提取方法相对应，本申请还提供一种生命体征信息提取装置，由于本申请实施例提供的生命体征信息提取装置与上述图1至图3实施例提供的生命体征信息提取方法相对应，因此在生命体征信息提取方法的实施方式也适用于本申请实施例提供的生命体征信息提取装置，在本申请实施例中不再详细描述。
图4为本申请实施例提供的一种生命体征信息提取装置的结构示意图。
如图4所示，该生命体征信息提取装置400可以包括：探测模块410、确定模块420、提取模块430以及处理模块440。
探测模块410，用于采用超宽带雷达探测得到雷达回波信号；
确定模块420，用于确定雷达回波信号中慢时间方向的各切片信号对应的信号能量；
提取模块430，用于基于预设宽度的能量窗以及各切片信号对应的信号能量，从雷达回波信号中提取人体体征信号；
处理模块440，用于对人体体征信号进行相干累积处理，并从处理后的人体体征信号中提取呼吸信号和心跳信号。
进一步地，在本申请实施例的一种可能的实现方式中，提取模块430，具体用于：
将各切片信号分组置于能量窗内，并分别确定能量窗内的多个切片信号组的信号能量；其中，切片信号组包括雷达回波信号中慢时间方向的预设数量的相邻切片信号；
从多个切片信号组中，提取对应信号能量满足第一预设条件的目标切片信号组；
将目标切片信号组中的各切片信号作为人体体征信号。
进一步地，在本申请实施例的另一种可能的实现方式中，该生命体征信息提取装置400还可以包括：
矫正模块，用于对雷达回波信号进行数值偏差矫正处理，以矫正雷达回波信号中因雷达不稳定性导致的数值偏差；
和/或，消除模块，用于采用线性趋势抑制法，对雷达回波信号进行背景波消除，以消除雷达回波信号中的恒定成分；
和/或，增强模块，用于对雷达回波信号进行信号增强处理，以增强雷达回波信号中的生命体征成分。
进一步地，在本申请实施例的另一种可能的实现方式中，矫正模块，具体用于：
获取雷达回波信号的慢时间方向的各切片信号的累计概率分布；
根据各切片信号的累计概率分布，将累计概率最大的切片信号确定为参考切片信号；
根据参考切片信号的累计概率分布，确定参考切片信号中概率满足第二预设条件的幅值；
确定幅值对应的慢时间方向上的慢时间采样点；
对慢时间采样点对应的快时间方向的切片信号的幅值进行矫正，得到矫正后雷达回波信号。
进一步地，在本申请实施例的另一种可能的实现方式中，增强模块，具体用于：
采用巴特沃斯滤波方法，滤除雷达回波信号中的噪声成分；
和/或，采用信号自相关处理方法，提高雷达回波信号的信噪比。
本申请实施例的生命体征信息提取装置，通过采用超宽带雷达探测得到雷达回波信号，确定雷达回波信号中慢时间方向的各切片信号对应的信号能量，基于预设宽度的能量窗以及各切片信号对应的信号能量，从雷达回波信号中提取人体体征信号，进而对人体体征信号进行相干累积处理，并从处理后的人体体征信号中提取呼吸信号和心跳信号。由此，实现了根据雷达回波信号中慢时间方向的各切片信号对应的信号能量以及预设宽度的能量窗，从雷达回波信号中准确提取呼吸信号和心跳信号。
为了实现上述实施例，如图5所示，本发明实施例提出了一种电子设备500，包括：存储器510、处理器520及存储在存储器510上并可在处理器520上运行的计算机程序，所述处理器520执行所述程序时，实现前述任一实施例提出的生命体征信息提取方法。
为了实现上述实施例，本申请还提出一种存储有计算机指令的非瞬时计算机可读存储介质，所述计算机指令用于使所述计算机执行如本申请前述任一实施例提出的生命体征信息提取方法。
为了实现上述实施例，本申请还提出一种计算机程序产品，包括计算机程序，所述计算机程序在被处理器执行时实现如本申请前述任一实施例提出的生命体征信息提取方法。
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本申请的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
此外，术语“第一”、“第二”仅用于描述目的，而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。在本申请的描述中，“多个”的含义是至少两个，例如两个，三个等，除非另有明确具体的限定。
流程图中或在此以其他方式描述的任何过程或方法描述可以被理解为，表示包括一个或更多个用于实现定制逻辑功能或过程的步骤的可执行指令的代码的模块、片段或部分，并且本申请的优选实施方式的范围包括另外的实现，其中可以不按所示出或讨论的顺序，包括根据所涉及的功能按基本同时的方式或按相反的顺序，来执行功能，这应被本申请的实施例所属技术领域的技术人员所理解。
在流程图中表示或在此以其他方式描述的逻辑和/或步骤，例如，可以被认为是用于实现逻辑功能的可执行指令的定序列表，可以具体实现在任何计算机可读介质中，以供指令执行系统、装置或设备(如基于计算机的系统、包括处理器的系统或其他可以从指令执行系统、装置或设备取指令并执行指令的系统)使用，或结合这些指令执行系统、装置或设备而使用。就本说明书而言，"计算机可读介质"可以是任何可以包含、存储、通信、传播或传输程序以供指令执行系统、装置或设备或结合这些指令执行系统、装置或设备而使用的装置。计算机可读介质的更具体的示例(非穷尽性列表)包括以下：具有一个或多个布线的电连接部(电子装置)，便携式计算机盘盒(磁装置)，随机存取存储器(RAM)，只读存储器(ROM)，可擦除可编辑只读存储器(EPROM或闪速存储器)，光纤装置，以及便携式光盘只读存储器(CDROM)。另外，计算机可读介质甚至可以是可在其上打印所述程序的纸或其他合适的介质，因为可以例如通过对纸或其他介质进行光学扫描，接着进行编辑、解译或必要时以其他合适方式进行处理来以电子方式获得所述程序，然后将其存储在计算机存储器中。
应当理解，本申请的各部分可以用硬件、软件、固件或它们的组合来实现。在上述实施方式中，多个步骤或方法可以用存储在存储器中且由合适的指令执行系统执行的软件或固件来实现。如，如果用硬件来实现和在另一实施方式中一样，可用本领域公知的下列技术中的任一项或他们的组合来实现：具有用于对数据信号实现逻辑功能的逻辑门电路的离散逻辑电路，具有合适的组合逻辑门电路的专用集成电路，可编程门阵列(PGA)，现场可编程门阵列(FPGA)等。
本技术领域的普通技术人员可以理解实现上述实施例方法携带的全部或部分步骤是可以通过程序来指令相关的硬件完成，所述的程序可以存储于一种计算机可读存储介质中，该程序在执行时，包括方法实施例的步骤之一或其组合。
此外，在本申请各个实施例中的各功能单元可以集成在一个处理模块中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个模块中。上述集成的模块既可以采用硬件的形式实现，也可以采用软件功能模块的形式实现。所述集成的模块如果以软件功能模块的形式实现并作为独立的产品销售或使用时，也可以存储在一个计算机可读取存储介质中。
上述提到的存储介质可以是只读存储器，磁盘或光盘等。尽管上面已经示出和描述了本申请的实施例，可以理解的是，上述实施例是示例性的，不能理解为对本申请的限制，本领域的普通技术人员在本申请的范围内可以对上述实施例进行变化、修改、替换和变型。
Claims (10)
1.一种生命体征信息提取方法，其特征在于，包括以下步骤：
采用超宽带雷达探测得到雷达回波信号；
确定所述雷达回波信号中慢时间方向的各切片信号对应的信号能量；
基于预设宽度的能量窗以及所述各切片信号对应的信号能量，从所述雷达回波信号中提取人体体征信号；
对所述人体体征信号进行相干累积处理，并从处理后的人体体征信号中提取呼吸信号和心跳信号。
2.根据权利要求1所述的方法，其特征在于，所述基于预设宽度的能量窗以及所述各切片信号对应的信号能量，从所述雷达回波信号中提取人体体征信号，包括：
将所述各切片信号分组置于所述能量窗内，并分别确定所述能量窗内的多个切片信号组的信号能量；其中，所述切片信号组包括所述雷达回波信号中慢时间方向的预设数量的相邻切片信号；
从多个所述切片信号组中，提取对应信号能量满足第一预设条件的目标切片信号组；
将所述目标切片信号组中的各切片信号作为所述人体体征信号。
3.根据权利要求1或2所述的方法，其特征在于，所述方法还包括：
对所述雷达回波信号进行数值偏差矫正处理，以矫正所述雷达回波信号中因雷达不稳定性导致的数值偏差；
和/或，采用线性趋势抑制法，对所述雷达回波信号进行背景波消除，以消除所述雷达回波信号中的恒定成分；
和/或，对所述雷达回波信号进行信号增强处理，以增强所述雷达回波信号中的生命体征成分。
4.根据权利要求3所述的方法，其特征在于，所述对所述雷达回波信号进行数值偏差矫正处理，以矫正所述雷达回波信号中因雷达不稳定性导致的数值偏差，包括：
获取所述雷达回波信号的慢时间方向的各所述切片信号的累计概率分布；
根据各所述切片信号的累计概率分布，将累计概率最大的切片信号确定为参考切片信号；
根据所述参考切片信号的累计概率分布，确定所述参考切片信号中概率满足第二预设条件的幅值；
确定所述幅值对应的慢时间方向上的慢时间采样点；
对所述慢时间采样点对应的快时间方向的切片信号的幅值进行矫正，得到矫正后雷达回波信号。
5.根据权利要求3所述的方法，其特征在于，所述对所述雷达回波信号进行信号增强处理，以增强所述雷达回波信号中的生命体征成分，包括：
采用巴特沃斯滤波方法，滤除所述雷达回波信号中的噪声成分；
和/或，采用信号自相关处理方法，提高所述雷达回波信号的信噪比。
6.一种生命体征信息提取装置，其特征在于，包括：
探测模块，用于采用超宽带雷达探测得到雷达回波信号；
确定模块，用于确定所述雷达回波信号中慢时间方向的各切片信号对应的信号能量；
提取模块，用于基于预设宽度的能量窗以及所述各切片信号对应的信号能量，从所述雷达回波信号中提取人体体征信号；
处理模块，用于对所述人体体征信号进行相干累积处理，并从处理后的人体体征信号中提取呼吸信号和心跳信号。
7.根据权利要求6所述的装置，其特征在于，所述提取模块，具体用于：
将所述各切片信号分组，得到多个切片信号组，其中，所述切片信号组包括所述雷达回波信号中慢时间方向的预设数量的相邻切片信号；
将多个所述切片信号组分别置于所述能量窗内，并确定所述能量窗内的多个所述切片信号组的信号能量；
从多个所述切片信号组中，提取对应信号能量满足第一预设条件的目标切片信号组；
将所述目标切片信号组中的各切片信号作为所述人体体征信号。
8.根据权利要求6或7所述的装置，其特征在于，还包括：
矫正模块，用于对所述雷达回波信号进行数值偏差矫正处理，以矫正所述雷达回波信号中因雷达不稳定性导致的数值偏差；
和/或，消除模块，用于采用线性趋势抑制法，对所述雷达回波信号进行背景波消除，以消除所述雷达回波信号中的恒定成分；
和/或，增强模块，用于对所述雷达回波信号进行信号增强处理，以增强所述雷达回波信号中的生命体征成分。
9.一种电子设备，其特征在于，包括：存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时，实现如权利要求1-5中任一项所述的方法。
10.一种存储有计算机指令的非瞬时计算机可读存储介质，其特征在于，所述计算机指令用于使所述计算机执行根据权利要求1-5中任一项所述的方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110974749.6A CN113786177B (zh) | 2021-08-24 | 2021-08-24 | 生命体征信息提取方法、装置及电子设备 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110974749.6A CN113786177B (zh) | 2021-08-24 | 2021-08-24 | 生命体征信息提取方法、装置及电子设备 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113786177A true CN113786177A (zh) | 2021-12-14 |
| CN113786177B CN113786177B (zh) | 2023-03-03 |
Family
ID=79182145
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110974749.6A Active CN113786177B (zh) | 2021-08-24 | 2021-08-24 | 生命体征信息提取方法、装置及电子设备 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113786177B (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114601432A (zh) * | 2022-01-28 | 2022-06-10 | 中国第一汽车股份有限公司 | 一种车载人体信息获取方法、装置及车辆 |
| CN114792109A (zh) * | 2022-05-16 | 2022-07-26 | 煤炭科学研究总院有限公司 | 生命体征信号提取方法、装置、电子设备及存储介质 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107167802A (zh) * | 2017-05-24 | 2017-09-15 | 北京大学 | 一种基于超宽带雷达的呼吸信号检测算法 |
| CN110507293A (zh) * | 2019-07-26 | 2019-11-29 | 中国电子科技集团公司第三十八研究所 | 一种超宽带穿墙雷达人体呼吸及心跳检测方法及系统 |
| CN110554381A (zh) * | 2019-08-30 | 2019-12-10 | 湖南正申科技有限公司 | 一种用于冲激脉冲式穿墙雷达的人体静止目标加速检测方法 |
| WO2020004721A1 (ko) * | 2018-06-27 | 2020-01-02 | 유메인주식회사 | 초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정 방법 |
| WO2020078249A1 (zh) * | 2018-10-15 | 2020-04-23 | 中国人民解放军第四军医大学 | 一种基于雷达信号的多目标生命探测方法及探测雷达 |
-
2021
- 2021-08-24 CN CN202110974749.6A patent/CN113786177B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107167802A (zh) * | 2017-05-24 | 2017-09-15 | 北京大学 | 一种基于超宽带雷达的呼吸信号检测算法 |
| WO2020004721A1 (ko) * | 2018-06-27 | 2020-01-02 | 유메인주식회사 | 초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정 방법 |
| WO2020078249A1 (zh) * | 2018-10-15 | 2020-04-23 | 中国人民解放军第四军医大学 | 一种基于雷达信号的多目标生命探测方法及探测雷达 |
| CN110507293A (zh) * | 2019-07-26 | 2019-11-29 | 中国电子科技集团公司第三十八研究所 | 一种超宽带穿墙雷达人体呼吸及心跳检测方法及系统 |
| CN110554381A (zh) * | 2019-08-30 | 2019-12-10 | 湖南正申科技有限公司 | 一种用于冲激脉冲式穿墙雷达的人体静止目标加速检测方法 |
Non-Patent Citations (2)
| Title |
|---|
| 聂镭: "复杂雷达回波特征提取及信号检测", 《中国优秀博硕士学位论文全文数据库(博士)信息科技辑》, no. 10, 15 October 2014 (2014-10-15), pages 68 - 108 * |
| 陈德肯: "基于雷达信号的非接触体征检测技术研究", 《中国优秀博硕士学位论文全文数据库(硕士)信息科技辑》, no. 9, 15 September 2019 (2019-09-15), pages 1 - 57 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114601432A (zh) * | 2022-01-28 | 2022-06-10 | 中国第一汽车股份有限公司 | 一种车载人体信息获取方法、装置及车辆 |
| CN114601432B (zh) * | 2022-01-28 | 2024-10-11 | 中国第一汽车股份有限公司 | 一种车载人体信息获取方法、装置及车辆 |
| CN114792109A (zh) * | 2022-05-16 | 2022-07-26 | 煤炭科学研究总院有限公司 | 生命体征信号提取方法、装置、电子设备及存储介质 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113786177B (zh) | 2023-03-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US20210321879A1 (en) | Method and apparatus for non-contact fast vital sign acquisition based on radar signal | |
| US20220378313A1 (en) | Multi-person vital signs monitoring using millimeter wave (mm-wave) signals | |
| CN113786177B (zh) | 生命体征信息提取方法、装置及电子设备 | |
| US20220142478A1 (en) | Radar cardiography: a precise cardiac data reconstruction method | |
| CN112336322B (zh) | 一种无接触式呼吸或心跳检测方法 | |
| CN107003387A (zh) | 用于运行机动车的雷达系统的方法和设备 | |
| US11860266B2 (en) | Detection system and method of detecting life using radar | |
| CN110837079B (zh) | 一种基于雷达的目标检测方法及装置 | |
| US20190320940A1 (en) | Apparatus and method for determining bio-information of target using impulse radar | |
| CN113702968A (zh) | 生命探测方法及其装置 | |
| US9964634B2 (en) | Imaging tissue motion estimation | |
| CN113786176B (zh) | 精确的毫米波雷达呼吸心跳测量方法、系统和存储介质 | |
| CN115969351A (zh) | 一种基于毫米波雷达实时计算呼吸的方法 | |
| CN115736872A (zh) | 一种基于毫米波雷达实时计算心率的方法 | |
| CN117838083A (zh) | 一种基于毫米波雷达的体征快速精确检测方法 | |
| CN112363160B (zh) | 一种基于宽带式信号的卧床跌落检测方法、介质、设备及装置 | |
| CN112617748A (zh) | 基于调频连续波雷达的体征数据处理方法及系统 | |
| CN110269642A (zh) | 基于分数阶傅里叶变换和小波变换的多普勒心率估计方法 | |
| CN117918812A (zh) | 一种基于雷达的非接触人体生命体征感知方法及系统 | |
| US11559215B2 (en) | Method for measuring high-accuracy realtime heart rate based on continuous-wave Doppler radar and radar system therefor | |
| Xian et al. | Improved Denoising Method for UWB Vital Signs Detection and Extraction | |
| CN107133624B (zh) | 一种目标检测方法及设备 | |
| US20240329198A1 (en) | Apparatus and method for determining a desired low-amplitude signal component of a radar receive signal of a radar sensor and radar system | |
| Elmakhzangy et al. | Improved Multi-Person Vital Signs Estimation Algorithm in Sitting and Standing Positions using MIMO FMCW Radar | |
| Xiang et al. | Heart rate estimation algorithm for radar signals based on wavelet threshold noise reduction and Resnet50 |
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 | ||
| CB02 | Change of applicant information | ||
| CB02 | Change of applicant information |
Address after: 100013 No. 5 Youth ditch East Road, Hepingli, Beijing, Chaoyang District Applicant after: General Coal Research Institute Co.,Ltd. Applicant after: LIAONING TECHNICAL University Address before: 100013 No. 5 Youth ditch East Road, Hepingli, Beijing, Chaoyang District Applicant before: CHINA COAL Research Institute Applicant before: LIAONING TECHNICAL University |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |