CN113786176B - 精确的毫米波雷达呼吸心跳测量方法、系统和存储介质 - Google Patents
精确的毫米波雷达呼吸心跳测量方法、系统和存储介质 Download PDFInfo
- Publication number
- CN113786176B CN113786176B CN202110939806.7A CN202110939806A CN113786176B CN 113786176 B CN113786176 B CN 113786176B CN 202110939806 A CN202110939806 A CN 202110939806A CN 113786176 B CN113786176 B CN 113786176B
- Authority
- CN
- China
- Prior art keywords
- signal
- signals
- heartbeat
- phase
- respiratory
- 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.)
- Active
Links
- 238000000691 measurement method Methods 0.000 title claims abstract description 8
- 230000000241 respiratory effect Effects 0.000 claims abstract description 100
- 238000005259 measurement Methods 0.000 claims abstract description 45
- 238000000034 method Methods 0.000 claims abstract description 36
- 238000005311 autocorrelation function Methods 0.000 claims description 32
- 230000029058 respiratory gaseous exchange Effects 0.000 claims description 28
- 238000001914 filtration Methods 0.000 claims description 17
- 238000001228 spectrum Methods 0.000 claims description 13
- 238000001514 detection method Methods 0.000 claims description 11
- 238000007781 pre-processing Methods 0.000 claims description 7
- 238000002592 echocardiography Methods 0.000 claims description 4
- 238000009531 respiratory rate measurement Methods 0.000 claims description 3
- 230000001131 transforming effect Effects 0.000 claims description 3
- 239000000284 extract Substances 0.000 claims 1
- 210000000038 chest Anatomy 0.000 description 6
- 210000001015 abdomen Anatomy 0.000 description 2
- 238000010009 beating Methods 0.000 description 2
- 230000004927 fusion Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 210000000115 thoracic cavity Anatomy 0.000 description 2
- 230000007547 defect Effects 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000001121 heart beat frequency Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000036391 respiratory frequency Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
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/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7221—Determining signal validity, reliability or quality
-
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Physiology (AREA)
- Animal Behavior & Ethology (AREA)
- Molecular Biology (AREA)
- Pathology (AREA)
- Signal Processing (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Radar, Positioning & Navigation (AREA)
- Surgery (AREA)
- Biophysics (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Remote Sensing (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Artificial Intelligence (AREA)
- Cardiology (AREA)
- Psychiatry (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Pulmonology (AREA)
- Power Engineering (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
本发明提供了一种精确的毫米波雷达呼吸心跳测量方法、系统和存储介质。首先根据回波信号幅度信息,得到目标初始位置信息A;接着提取位置A及其临近单元的相位信号并评估位置A及其临近单元的呼吸信号和心跳信号质量,最终根据呼吸信号和心跳信号质量确定呼吸心跳测量点。本发明综合利用了雷达回波的幅度信息和相位信息,从而得到了更加精确的呼吸心跳测量位置,同时提高了呼吸心跳估计精度。
Description
技术领域
本发明涉及呼吸心跳测量领域,具体涉及基于毫米波雷达的呼吸心跳测量。
背景技术
雷达发射电磁波到空间并接收回波,若空间中存在目标,目标表面会产生反射信号被雷达接收到。目标与雷达之间的距离体现在发射波形与接收波形之间的时延上。雷达按照发射信号的不同分为脉冲雷达和连续波雷达,本发明中讨论的是连续波雷达,更具体来说是线性调频连续波雷达(Chirp continuous wave radar)。在线性调频连续波雷达中,发射信号为Chirp信号,这是一种频率随着时间成线性变化的信号。在用线性调频连续波雷达测量目标距离时,目标与雷达的距离产生时延差,时延差导致发射信号与接收信号产生固定的频率差,通过测量该频率差进而推断出目标与雷达之间的距离信息。
人体在正常呼吸时,胸腔会随着呼吸过程产生相应的起伏,同样的,心脏的跳动亦会带来心脏外皮肤附近的振动。雷达在对人体进行探测时,人体胸腔的起伏和心脏的跳动会导致人体目标与雷达之间的距离发生微弱的变化。这种距离的微弱变化使得接收到的雷达信号产生相应的相位变化。通过测量该相位的变化,可估算出人体的呼吸状态与心跳状态。
在用雷达进行人体呼吸、心跳体征监测时,测量点的选取直接影响了呼吸、心跳信号的估计质量。通常的选取方式为人体目标位置回波信号幅度最大处,即:认为雷达回波信号幅度最大处为人体目标胸腔位置。然而,在实际中,由于雷达分辨率的限制及个体之间的差异等原因,雷达回波幅度最大处不一定是胸腔位置,也有可能是腹部或者其他位置,进而导致呼吸、心跳信号估计不准或失效。
现有技术CN112674738A公开了一种呼吸心跳信号的检测方法及装置。其中,该方法包括:获取雷达发射的电磁波信号遇到障碍物反射回来的回波信号;从回波信号中提取出目标回波信号,其中,目标回波信号为电磁波信号遇到人体的胸腔反射回来的信号;从目标回波信号中提取出人体的呼吸信号和心跳信号;依据呼吸信号确定人体的呼吸频率,并依据心跳信号确定人体的心跳频率。现有技术CN112168152A公开了一种呼吸心跳的检测方法、检测装置与计算机可读存储介质。该方法包括:获取原始雷达回波;对原始雷达回波进行预处理,得到呼吸心跳信号;对呼吸心跳信号进行滤波处理,得到呼吸信号和心跳信号;在时域确定呼吸信号的有效波峰数目和心跳信号的有效周期数目;根据呼吸信号的有效波峰数目和第一时间长度,确定呼吸信号的有效频率,第一时间长度为呼吸信号的时间长度;根据心跳信号的有效周期数目和第二时间长度,确定心跳信号的有效频率,第二时间长度为心跳信号的时间长度。
上述两项现有技术均采用直接从距离维信号中选择能量最高处作为人体目标回波信号的方式确定目标位置。接着,从该人体目标回波处,提取相位信息并进行后续呼吸、心跳信号的估计。然而,直接选取距离维信号中能量最高处作为人体目标回波信号是不合理的。这是因为,由于雷达分辨率的限制及个体之间的差异等原因,雷达回波幅度最大处不一定就是人体的胸腔位置(也有可能是腹部或者其他)。简言之,雷达回波幅度最大处的呼吸和心跳信号并不一定是最强的,导致呼吸、心跳信号估计不准或失效。此外,上述专利中,仅用距离维信息进行人体目标检测是有缺陷的,若在同个距离(不同方位)同时存在多个人,会导致多个人的呼吸、心跳相互叠加,造成后续呼吸、心跳信号估计错误。
因此,现有技术的上述方案中,存在着雷达回波幅度最大处与胸腔位置不匹配导致的呼吸心跳估计不准或失效的问题。
发明内容
为了克服上述现有技术的不足,本发明提供一种精确的毫米波雷达呼吸心跳测量方法、系统和存储介质。
本发明提供了一种精确的毫米波雷达呼吸心跳测量方法,其特征在于,包括:
步骤1,获取雷达发射的电磁波信号遇到障碍物后的反射回波信号;
步骤2,根据回波信号幅度信息,得到目标位置A;
步骤3,提取位置A及其临近单元的相位信号并做预处理;
步骤4,提取位置A及其临近单元的呼吸信号和心跳信号;
步骤5,评估位置A及其临近单元的呼吸信号和心跳信号质量,并根据结果确定最终的呼吸测量位置AB和心跳测量位置AH。
可选地,步骤1中,包括:
步骤1-1,对雷达原始回波信号进行距离维傅里叶变换;
步骤1-2,进行波束成形,得到探测空间中各个点的反射回波。
可选地,步骤1-1中,对采集到的每一根天线的雷达回波信号Secho进行傅里叶变换得到各个天线的距离维信号其中j为虚数符号,ΔR为不同距离的相位,t为快时间;
可选地,步骤1-2中,对雷达原始回波信号进行波束成形,包括:根据雷达天线的排布情况与雷达探测空间的几何关系,对各个雷达天线的回波信号进行加权,以增强感兴趣信号和抑制干扰信号;
具体地,对距离维信号SRange进行传统波束成形:计算空间目标(/>分别代表目标的水平、垂直入射角)到达雷达各个天线的时延/>其中(xi,zi)为第i根天线相对参考点的坐标;对不同天线接收的信号进行加权求和得到波束成形结果/>
可选地,步骤1中,雷达信号为线性调频连续波。
可选地,步骤2中,包括:选取回波信号中幅度最大值的位置作为目标位置A;
具体地,Smax=max(abs(SBF));
A=find(abs(SBF)==Smax);
其中,SBF为波束成形后的信号;abs(·)为取绝对值,max(·)表示求最大值,find(·)表示找满足条件的目标位置,A为找到的目标位置。
可选地,步骤3中,包括:以位置A为原点,提取位置A及其临近单元的相位信号,对相位信号进行解缠绕处理,得到处理后的相位信号;
具体地,PhA=phase(SBF(A));
PhNA1=phase(SBF(NA1));
PhNA2=phase(SBF(NA2));
……
PhNAn=phase(SBF(NAn));
其中,NA1、NA2、……、NAn为位置A临近单元的位置,共n个;PhA、PhNA1、PhNA2、……、PhNAn为求相位信号后的信号,共n+1个;SBF(·)为特定位置处的波束成形后的信号;phase(·)表示求相位信号;
按照时间顺序,对每一帧分别进行上述操作,得到每一个时刻的相位信息,并按照时间先后顺序排列得到一段时间内A及其周围临近单元NA处的相位信息,并对相位信号解缠绕操作:
PhaseA=unwrap([PhA(1),PhA(2),…PhA(t)]);
PhaseNA1=unwrap([PhNA1(1),PhNA1(2),…PhNA1(t)]);
PhaseNA2=unwrap([PhNA2(1),PhNA1(2),…PhNA2(t)]);
……
PhaseNAn=unwrap([PhNAn(1),PhNAn(2),…PhNAn(t)]);
其中,PhaseA、PhaseNA1、PhaseNA2、……、PhaseNAn为解缠绕后的信号,共n+1个;PhA(·)、PhNA1(·)、……、PhNAn(·)分别为信号的特定相位,每个信号都包含t个相位;unwrap(·)表示对相位信号解缠绕。
可选地,步骤4中,包括:
步骤4-1,用呼吸带通滤波器对处理后的位置A及其临近单元的相位信号进行带通滤波;
步骤4-2,用心跳带通滤波器对处理后的位置A及其临近单元的相位信号进行带通滤波;
步骤4-3,得到进行带通滤波后的信号;
具体地,呼吸信号表示为:
S_breath_A=filter(PhaseA,BP1);
S_breath_NA1=filter(PhaseNA1,BP1);
S_breath_NA2=filter(PhaseNA2,BP1);
……
S_breath_NAn=filter(PhaseNAn,BP1);
心跳信号表示为:
S_heart_A=filter(PhaseA,BP2);
S_heart_NA1=filter(PhaseNA1,BP2);
S_heart_NA2=filter(PhaseNA2,BP2);
……
S_heart_NAn=filter(PhaseNAn,BP2);
其中,filter(·)表示对信号进行滤波操作,BP1为通带0.1~0.5Hz的滤波器参数,BP2为通带0.8~2.5Hz的滤波器参数;S_breath_A、S_breath_NA1、S_breath_NA2、……S_breath_NAn为呼吸信号,共n+1个;S_heart_A、S_heart_NA1、S_heart_NA2、……S_heart_NAn为心跳信号,共n+1个。
可选地,步骤4-1中,呼吸带通滤波器通带范围为0.1-0.5Hz;可选地,步骤4-2中,心跳带通滤波器通带范围为0.8-2.5Hz。
可选地,步骤5中,包括:
步骤5-1,对位置A及其临近单元的呼吸信号质量进行评估;
步骤5-2,对位置A及其临近单元的心跳信号质量进行评估;
步骤5-3,根据结果确定最终的呼吸测量位置AB和心跳测量位置AH。
可选地,步骤5-1中,包括:
步骤5-1-1,对呼吸信号进行慢时间傅里叶变换,得到呼吸信号的慢时间频谱;
步骤5-1-2,保存频谱中呼吸信号频率范围内的最高频点处幅度值;其中,呼吸信号频率范围为0.1-0.5Hz;
步骤5-1-3,计算呼吸信号慢时间维自相关函数;
步骤5-1-4,保存自相关函数中呼吸信号周期范围内的峰值;其中,呼吸信号周期范围为2-10s;
步骤5-1-5,对呼吸信号范围内幅度值、自相关值进行加权,得到呼吸质量评分值;
可选地,步骤5-1中,对目标位置A处呼吸信号进行质量评估:计算A处呼吸信号的频谱及自相关函数ACF,
S_breath_A_Spec=abs(fft(S_breath_A))
S_breath_A_ACF=autocorr(S_breath_A)
其中,abs(·)为取绝对值,fft(·)表示求傅里叶变换,autocorr(·)表示求自相关函数;保留S_breath_A_Spec频率范围[0.1,0.5]Hz内幅度的最大值,记为BS1;保留S_breath_A_ACF周期范围[2,10]s内幅度的峰值,记为BS2;将BS1与BS2加权融合得到位置A处呼吸信号质量得分:ScoreBA=a*BS1+b*BS2,其中a,b为BS1、BS2的加权系数,可根据经验进行设置;
同样的,分别对NA1,NA2,…NAn处进行如上操作,得到NA1,NA2,…NAn处呼吸信号质量:ScoreBNA1,ScoreBNA2,…,ScoreBNAn;最后,从ScoreBA,ScoreBNA1,ScoreBNA2,…,ScoreBNAn中挑选得分最高值处位置作为呼吸测量点,并输出该位置的呼吸测量结果。
可选地,步骤5-2中,包括:
步骤5-2-1,对心跳信号进行慢时间傅里叶变换,得到心跳信号的慢时间频谱;
步骤5-2-2,保存频谱中心跳信号频率范围内最高频点处幅度值;其中,心跳信号频率范围为0.8-2.5Hz;
步骤5-2-3,计算心跳信号慢时间维自相关函数;
步骤5-2-4,保存自相关函数中心跳信号周期范围内峰值;其中,心跳信号周期范围为0.4-1.25s;
步骤5-2-5,对心跳信号范围内幅度值、自相关值进行加权,得到心跳质量评分值。
可选地,步骤5-2,包括:对目标位置A处心跳信号进行质量评估的方法为:计算A处心跳信号的频谱及自相关函数ACF,
S_heart_A_Spec=abs(fft(S_heart_A))
S_heart_A_ACF=autocorr(S_heart_A)
其中,abs(·)为取绝对值,fft(·)表示求傅里叶变换,autocorr(·)表示求自相关函数;保留S_heart_A_Spec频率范围[0.8,2.5]Hz内幅度的最大值,记为HS1;保留S_heart_A_ACF周期范围[0.4,1.25]s内幅度的峰值,记为HS2;将HS1与HS2加权融合得到位置A处呼吸信号质量得分:ScoreHA=c*HS1+d*HS2,其中c,d为HS1、HS2的加权系数,可根据经验进行设置;
同样的,分别对NA1,NA2,…NAn处进行如上操作,得到NA1,NA2,…NAn处心跳信号质量:ScoreHNA1,ScoreHNA2,…,ScoreHNAn。最后,从ScoreHA,ScoreHNA1,ScoreHNA2,…,ScoreHNAn中挑选得分最高值处位置作为心跳测量点,并输出该位置的心跳测量结果。
可选地,步骤5-3,包括:
步骤5-3-1,找出位置A及其临近单元的呼吸信号最高的质量评分值;
步骤5-3-2,找出位置A及其临近单元的心跳信号最高的质量评分值;
步骤5-3-3,将位置A及其临近单元呼吸信号及心跳信号得分值最大值处位置确定为呼吸和心跳测量位置。
本发明提供了一种精确的毫米波雷达呼吸心跳测量系统,其特征在于,包括:
反射回波信号获取模块,获取雷达发射的电磁波信号遇到障碍物后的反射回波信号;
目标位置获取模块,根据回波信号幅度信息,得到目标位置A;
相位信号获取模块,提取位置A及其临近单元的相位信号并做预处理;
呼吸信号和心跳信号获取模块,提取位置A及其临近单元的呼吸信号和心跳信号;
呼吸测量位置和心跳测量位置获取模块,评估位置A及其临近单元的呼吸信号和心跳信号质量,并根据结果确定最终的呼吸测量位置AB和心跳测量位置AH。
本发明提供了一种计算机可读存储介质,其特征在于,所述计算机可读存储介质包括存储的程序,其中,在所述程序运行时控制所述计算机可读存储介质所在设备执行上述精确的毫米波雷达呼吸心跳测量方法。
基于上述方案,本发明解决了现有技术中由于测量位置不准确,导致呼吸心跳估计失效或者错误的问题。与只用雷达回波的幅度信息进行目标定位的方法相比,本发明额外利用了相位信息进行位置精确定位,得到的呼吸心跳测量位置更加准确。本发明综合利用了雷达回波的幅度信息和相位信息,从而得到了更加精确的呼吸心跳测量位置,同时提高了呼吸心跳估计精度。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施方式,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为精确的呼吸、心跳检测位置定位流程图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合本发明具体实施例及相应的附图对本发明技术方案进行清楚、完整地描述。显然,所描述的实施例仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
相关术语解释:
雷达:利用电磁波探测目标的电子设备;
毫米波雷达:工作频率为30-300GHz(对应电磁波波长为10mm~1mm)的雷达;
FFT:快速傅里叶变换,一种将时域信号变换为频域信号的算法;
滤波:是一种将信号中特定波段频率滤除的操作;
波束成形:一种空域滤波方法,以增强感兴趣信号和抑制干扰信号为目的。
图1是根据本发明实施例的精确的呼吸、心跳检测位置定位流程图。该方法包括以下实施步骤:
步骤S101,获取雷达原始回波信号,对原始回波信号进行距离维傅里叶变换和波束成形,得到雷达探测空间的回波;
步骤S102,根据步骤S101的结果,选取雷达探测空间回波中幅度最大值处的位置作为目标位置信息并保存;
步骤S103,以步骤S102中目标位置信息为原点,分别提取它及其临近单元的相位信号,并对相位信号进行解缠绕处理;
步骤S104,对步骤S103中得到的相位信号分别进行带通滤波,得到目标位置及其临近单元的呼吸信号和心跳信号;
步骤S105,分别对步骤S104中目标位置及其临近单元的呼吸信号和心跳信号质量进行评估,根据评估质量得到呼吸心跳的准确监测点。
通过上述步骤,可以从雷达探测回波中准确地定位目标的呼吸心跳测量点,并得到更加精确的呼吸心跳估计频率(周期)。
在另一实施例中,本发明提供了一种精确的毫米波雷达呼吸心跳测量方法,其特征在于,包括:
步骤1,获取雷达发射的电磁波信号遇到障碍物后的反射回波信号;
步骤2,根据回波信号幅度信息,得到目标位置A;
步骤3,提取位置A及其临近单元的相位信号并做预处理;
步骤4,提取位置A及其临近单元的呼吸信号和心跳信号;
步骤5,评估位置A及其临近单元的呼吸信号和心跳信号质量,并根据结果确定最终的呼吸测量位置AB和心跳测量位置AH。
可选地,步骤1中,包括:
步骤1-1,对雷达原始回波信号进行距离维傅里叶变换;
步骤1-2,进行波束成形,得到探测空间中各个点的反射回波。
可选地,步骤1-1中,对采集到的每一根天线的雷达回波信号Secho进行傅里叶变换得到各个天线的距离维信号其中j为虚数符号,ΔR为不同距离的相位,t为快时间;
可选地,步骤1-2中,对雷达原始回波信号进行波束成形,包括:根据雷达天线的排布情况与雷达探测空间的几何关系,对各个雷达天线的回波信号进行加权,以增强感兴趣信号和抑制干扰信号;
具体地,对距离维信号SRange进行传统波束成形:计算空间目标(/>分别代表目标的水平、垂直入射角)到达雷达各个天线的时延/>其中(xi,zi)为第i根天线相对参考点的坐标;对不同天线接收的信号进行加权求和得到波束成形结果/>
可选地,步骤1中,雷达信号为线性调频连续波。
可选地,步骤2中,包括:选取回波信号中幅度最大值的位置作为目标位置A;
具体地,Smax=max(abs(SBF));
A=find(abs(SBF)==Smax);
其中,SBF为波束成形后的信号;abs(·)为取绝对值,max(·)表示求最大值,find(·)表示找满足条件的目标位置,A为找到的目标位置。
可选地,步骤3中,包括:以位置A为原点,提取位置A及其临近单元的相位信号,对相位信号进行解缠绕处理,得到处理后的相位信号;
具体地,PhA=phase(SBF(A));
PhNA1=phase(SBF(NA1));
PhNA2=phase(SBF(NA2));
……
PhNAn=phase(SBF(NAn));
其中,NA1、NA2、……、NAn为位置A临近单元的位置,共n个;PhA、PhNA1、PhNA2、……、PhNAn为求相位信号后的信号,共n+1个;SBF(·)为特定位置处的波束成形后的信号;phase(·)表示求相位信号;
按照时间顺序,对每一帧分别进行上述操作,得到每一个时刻的相位信息,并按照时间先后顺序排列得到一段时间内A及其周围临近单元NA处的相位信息,并对相位信号解缠绕操作:
PhaseA=unwrap([PhA(1),PhA(2),…PhA(t)]);
PhaseNA1=unwrap([PhNA1(1),PhNA1(2),…PhNA1(t)]);
PhaseNA2=unwrap([PhNA2(1),PhNA1(2),…PhNA2(t)]);
……
PhaseNAn=unwrap([PhNAn(1),PhNAn(2),…PhNAn(t)]);
其中,PhaseA、PhaseNA1、PhaseNA2、……、PhaseNAn为解缠绕后的信号,共n+1个;PhA(·)、PhNA1(·)、……、PhNAn(·)分别为信号的特定相位,每个信号都包含t个相位;unwrap(·)表示对相位信号解缠绕。
可选地,步骤4中,包括:
步骤4-1,用呼吸带通滤波器对处理后的位置A及其临近单元的相位信号进行带通滤波;
步骤4-2,用心跳带通滤波器对处理后的位置A及其临近单元的相位信号进行带通滤波;
步骤4-3,得到进行带通滤波后的信号;
具体地,呼吸信号表示为:
S_breath_A=filter(PhaseA,BP1);
S_breath_NA1=filter(PhaseNA1,BP1);
S_breath_NA2=filter(PhaseNA2,BP1);
……
S_breath_NAn=filter(PhaseNAn,BP1);
心跳信号表示为:
S_heart_A=filter(PhaseA,BP2);
S_heart_NA1=filter(PhaseNA1,BP2);
S_heart_NA2=filter(PhaseNA2,BP2);
……
S_heart_NAn=filter(PhaseNAn,BP2);
其中,filter(·)表示对信号进行滤波操作,BP1为通带0.1~0.5Hz的滤波器参数,BP2为通带0.8~2.5Hz的滤波器参数;S_breath_A、S_breath_NA1、S_breath_NA2、……S_breath_NAn为呼吸信号,共n+1个;S_heart_A、S_heart_NA1、S_heart_NA2、……S_heart_NAn为心跳信号,共n+1个。
可选地,步骤4-1中,呼吸带通滤波器通带范围为0.1-0.5Hz;可选地,步骤4-2中,心跳带通滤波器通带范围为0.8-2.5Hz。
可选地,步骤5中,包括:
步骤5-1,对位置A及其临近单元的呼吸信号质量进行评估;
步骤5-2,对位置A及其临近单元的心跳信号质量进行评估;
步骤5-3,根据结果确定最终的呼吸测量位置AB和心跳测量位置AH。
可选地,步骤5-1中,包括:
步骤5-1-1,对呼吸信号进行慢时间傅里叶变换,得到呼吸信号的慢时间频谱;
步骤5-1-2,保存频谱中呼吸信号频率范围内的最高频点处幅度值;其中,呼吸信号频率范围为0.1-0.5Hz;
步骤5-1-3,计算呼吸信号慢时间维自相关函数;
步骤5-1-4,保存自相关函数中呼吸信号周期范围内的峰值;其中,呼吸信号周期范围为2-10s;
步骤5-1-5,对呼吸信号范围内幅度值、自相关值进行加权,得到呼吸质量评分值;
可选地,步骤5-1中,对目标位置A处呼吸信号进行质量评估:计算A处呼吸信号的频谱及自相关函数ACF,
S_breath_A_Spec=abs(fft(S_breath_A))
S_breath_A_ACF=autocorr(S_breath_A)
其中,abs(·)为取绝对值,fft(·)表示求傅里叶变换,autocorr(·)表示求自相关函数;保留S_breath_A_Spec频率范围[0.1,0.5]Hz内幅度的最大值,记为BS1;保留S_breath_A_ACF周期范围[2,10]s内幅度的峰值,记为BS2;将BS1与BS2加权融合得到位置A处呼吸信号质量得分:ScoreBA=a*BS1+b*BS2,其中a,b为BS1、BS2的加权系数,可根据经验进行设置;
同样的,分别对NA1,NA2,…NAn处进行如上操作,得到NA1,NA2,…NAn处呼吸信号质量:ScoreBNA1,ScoreBNA2,…,ScoreBNAn;最后,从ScoreBA,ScoreBNA1,ScoreBNA2,…,ScoreBNAn中挑选得分最高值处位置作为呼吸测量点,并输出该位置的呼吸测量结果。
可选地,步骤5-2中,包括:
步骤5-2-1,对心跳信号进行慢时间傅里叶变换,得到心跳信号的慢时间频谱;
步骤5-2-2,保存频谱中心跳信号频率范围内最高频点处幅度值;其中,心跳信号频率范围为0.8-2.5Hz;
步骤5-2-3,计算心跳信号慢时间维自相关函数;
步骤5-2-4,保存自相关函数中心跳信号周期范围内峰值;其中,心跳信号周期范围为0.4-1.25s;
步骤5-2-5,对心跳信号范围内幅度值、自相关值进行加权,得到心跳质量评分值。
可选地,步骤5-2,包括:对目标位置A处心跳信号进行质量评估的方法为:计算A处心跳信号的频谱及自相关函数ACF,
S_heart_A_Spec=abs(fft(S_heart_A))
S_heart_A_ACF=autocorr(S_heart_A)
其中,abs(·)为取绝对值,fft(·)表示求傅里叶变换,autocorr(·)表示求自相关函数;保留S_heart_A_Spec频率范围[0.8,2.5]Hz内幅度的最大值,记为HS1;保留S_heart_A_ACF周期范围[0.4,1.25]s内幅度的峰值,记为HS2;将HS1与HS2加权融合得到位置A处呼吸信号质量得分:ScoreHA=c*HS1+d*HS2,其中c,d为HS1、HS2的加权系数,可根据经验进行设置;
同样的,分别对NA1,NA2,…NAn处进行如上操作,得到NA1,NA2,…NAn处心跳信号质量:ScoreHNA1,ScoreHNA2,…,ScoreHNAn。最后,从ScoreHA,ScoreHNA1,ScoreHNA2,…,ScoreHNAn中挑选得分最高值处位置作为心跳测量点,并输出该位置的心跳测量结果。
可选地,步骤5-3,包括:
步骤5-3-1,找出位置A及其临近单元的呼吸信号最高的质量评分值;
步骤5-3-2,找出位置A及其临近单元的心跳信号最高的质量评分值;
步骤5-3-3,将位置A及其临近单元呼吸信号及心跳信号得分值最大值处位置确定为呼吸和心跳测量位置。
在另一实施例中,本发明提供了一种精确的毫米波雷达呼吸心跳测量系统,其特征在于,包括:
反射回波信号获取模块,获取雷达发射的电磁波信号遇到障碍物后的反射回波信号;
目标位置获取模块,根据回波信号幅度信息,得到目标位置A;
相位信号获取模块,提取位置A及其临近单元的相位信号并做预处理;
呼吸信号和心跳信号获取模块,提取位置A及其临近单元的呼吸信号和心跳信号;
呼吸测量位置和心跳测量位置获取模块,评估位置A及其临近单元的呼吸信号和心跳信号质量,并根据结果确定最终的呼吸测量位置AB和心跳测量位置AH。
在另一实施例中,本发明提供了一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储有多个程序;在所述程序运行时控制所述计算机可读存储介质所在设备加载并执行上述精确的毫米波雷达呼吸心跳测量方法。
最后应说明的是:以上各实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述各实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。
Claims (8)
1.一种精确的毫米波雷达呼吸心跳测量方法,其特征在于,包括:
步骤1,获取雷达发射的电磁波信号遇到障碍物后的反射回波信号;
步骤2,根据回波信号幅度信息,得到目标位置A,具体包括:
选取回波信号中幅度最大值的位置作为目标位置A;
具体地,Smax=max(abs(SBF));
A=find(abs(SBF)==Smax);
其中,SBF为波束成形后的信号;abs(·)为取绝对值,max(·)表示求最大值,find(·)表示找满足条件的目标位置,A为找到的目标位置;
步骤3,提取位置A及其临近单元的相位信号并做预处理,具体包括:
以位置A为原点,提取位置A及其临近单元的相位信号,对相位信号进行解缠绕处理,得到处理后的相位信号;
具体地,PhA=phase(SBF(A));
PhNA1=phase(SBF(NA1));
PhNA2=phase(SBF(NA2));
……
PhNAn=phase(SBF(NAn));
其中,NA1、NA2、……、NAn为位置A临近单元的位置,共n个;PhA、PhNA1、PhNA2、……、PhNAn为求相位信号后的信号,共n+1个;SBF(·)为特定位置处的波束成形后的信号;phase(·)表示求相位信号;
按照时间顺序,对每一帧分别进行处理,得到每一个时刻的相位信息,并按照时间先后顺序排列得到一段时间内A及其周围临近单元NA处的相位信息,并对相位信号解缠绕操作:
PhaseA=unwrap([PhA(1),PhA(2),…PhA(t)]);
PhaseNA1=unwrap([PhNA1(1),PhNA1(2),…PhNA1(t)]);
PhaseNA2=unwrap([PhNA2(1),PhNA1(2),…PhNA2(t)]);
……
PhaseNAn=unwrap([PhNAn(1),PhNAn(2),…PhNAn(t)]);
其中,PhaseA、PhaseNA1、PhaseNA2、……、PhaseNAn为解缠绕后的信号,共n+1个;PhA(·)、PhNA1(·)、……、PhNAn(·)分别为信号的特定相位,每个信号都包含t个相位;unwrap(·)表示对相位信号解缠绕;
步骤4,提取位置A及其临近单元的呼吸信号和心跳信号;
步骤5,评估位置A及其临近单元的呼吸信号和心跳信号质量,并根据结果确定最终的呼吸测量位置AB和心跳测量位置AH。
2.根据权利要求1所述的方法,其特征在于,步骤1中,包括:
步骤1-1,对雷达原始回波信号进行距离维傅里叶变换;
对采集到的每一根天线的雷达回波信号Secho进行傅里叶变换得到各个天线的距离维信号其中j为虚数符号,ΔR为不同距离的相位,t为快时间;
步骤1-2,进行波束成形,得到探测空间中各个点的反射回波;
对雷达原始回波信号进行波束成形,包括:根据雷达天线的排布情况与雷达探测空间的几何关系,对各个雷达天线的回波信号进行加权,以增强感兴趣信号和抑制干扰信号;
具体地,对距离维信号SRange进行传统波束成形:计算空间目标到达雷达各个天线的时延/>其中(xi,zi)为第i根天线相对参考点的坐标;对不同天线接收的信号进行加权求和得到波束成形结果/>
其中,分别代表目标的水平、垂直入射角。
3.根据权利要求1所述的方法,其特征在于,步骤4中,包括:
步骤4-1,用呼吸带通滤波器对处理后的位置A及其临近单元的相位信号进行带通滤波;
步骤4-2,用心跳带通滤波器对处理后的位置A及其临近单元的相位信号进行带通滤波;
步骤4-3,得到进行带通滤波后的信号;
具体地,呼吸信号表示为:
S_breath_A=filter(PhaseA,BP1);
S_breath_NA1=filter(PhaseNA1,BP1);
S_breath_NA2=filter(PhaseNA2,BP1);
……
S_breath_NAn=filter(PhaseNAn,BP1);
心跳信号表示为:
S_heart_A=filter(PhaseA,BP2);
S_heart_NA1=filter(PhaseNA1,BP2);
S_heart_NA2=filter(PhaseNA2,BP2);
……
S_heart_NAn=filter(PhaseNAn,BP2);
其中,filter(·)表示对信号进行滤波操作,BP1为通带0.1~0.5Hz的滤波器参数,BP2为通带0.8~2.5Hz的滤波器参数;S_breath_A、S_breath_NA1、S_breath_NA2、……S_breath_NAn为呼吸信号,共n+1个;S_heart_A、S_heart_NA1、S_heart_NA2、……S_heart_NAn为心跳信号,共n+1个。
4.根据权利要求1所述的方法,其特征在于,步骤5中,包括:
步骤5-1,对位置A及其临近单元的呼吸信号质量进行评估;
步骤5-2,对位置A及其临近单元的心跳信号质量进行评估;
步骤5-3,根据结果确定最终的呼吸测量位置AB和心跳测量位置AH。
5.根据权利要求4所述的方法,其特征在于,步骤5-1中,包括:
步骤5-1-1,对呼吸信号进行慢时间傅里叶变换,得到呼吸信号的慢时间频谱;
步骤5-1-2,保存频谱中呼吸信号频率范围内的最高频点处幅度值;其中,呼吸信号频率范围为0.1-0.5Hz;
步骤5-1-3,计算呼吸信号慢时间维自相关函数;
步骤5-1-4,保存自相关函数中呼吸信号周期范围内的峰值;其中,呼吸信号周期范围为2-10s;
步骤5-1-5,对呼吸信号范围内幅度值、自相关值进行加权,得到呼吸质量评分值;
步骤5-2中,包括:
步骤5-2-1,对心跳信号进行慢时间傅里叶变换,得到心跳信号的慢时间频谱;
步骤5-2-2,保存频谱中心跳信号频率范围内最高频点处幅度值;其中,心跳信号频率范围为0.8-2.5Hz;
步骤5-2-3,计算心跳信号慢时间维自相关函数;
步骤5-2-4,保存自相关函数中心跳信号周期范围内峰值;其中,心跳信号周期范围为0.4-1.25s;
步骤5-2-5,对心跳信号范围内幅度值、自相关值进行加权,得到心跳质量评分值。
6.根据权利要求4所述的方法,其特征在于,步骤5-3,包括:
步骤5-3-1,找出位置A及其临近单元的呼吸信号最高的质量评分值;
步骤5-3-2,找出位置A及其临近单元的心跳信号最高的质量评分值;
步骤5-3-3,将位置A及其临近单元呼吸信号及心跳信号得分值最大值处位置确定为呼吸和心跳测量位置。
7.一种精确的毫米波雷达呼吸心跳测量系统,其特征在于,包括:
反射回波信号获取模块,获取雷达发射的电磁波信号遇到障碍物后的反射回波信号;
目标位置获取模块,根据回波信号幅度信息,得到目标位置A,具体包括:选取回波信号中幅度最大值的位置作为目标位置A;
具体地,Smax=max(abs(SBF));
A=find(abs(SBF)==Smax);
其中,SBF为波束成形后的信号;abs(·)为取绝对值,max(·)表示求最大值,find(·)表示找满足条件的目标位置,A为找到的目标位置;
相位信号获取模块,提取位置A及其临近单元的相位信号并做预处理,具体包括:
以位置A为原点,提取位置A及其临近单元的相位信号,对相位信号进行解缠绕处理,得到处理后的相位信号;
具体地,PhA=phase(SBF(A));
PhNA1=phase(SBF(NA1));
PhNA2=phase(SBF(NA2));
……
PhNAn=phase(SBF(NAn));
其中,NA1、NA2、……、NAn为位置A临近单元的位置,共n个;PhA、PhNA1、PhNA2、……、PhNAn为求相位信号后的信号,共n+1个;SBF(·)为特定位置处的波束成形后的信号;phase(·)表示求相位信号;
按照时间顺序,对每一帧分别进行处理,得到每一个时刻的相位信息,并按照时间先后顺序排列得到一段时间内A及其周围临近单元NA处的相位信息,并对相位信号解缠绕操作:
PhaseA=unwrap([PhA(1),PhA(2),…PhA(t)]);
PhaseNA1=unwrap([PhNA1(1),PhNA1(2),…PhNA1(t)]);
PhaseNA2=unwrap([PhNA2(1),PhNA1(2),…PhNA2(t)]);
……
PhaseNAn=unwrap([PhNAn(1),PhNAn(2),…PhNAn(t)]);
其中,PhaseA、PhaseNA1、PhaseNA2、……、PhaseNAn为解缠绕后的信号,共n+1个;PhA(·)、PhNA1(·)、……、PhNAn(·)分别为信号的特定相位,每个信号都包含t个相位;unwrap(·)表示对相位信号解缠绕;
呼吸信号和心跳信号获取模块,提取位置A及其临近单元的呼吸信号和心跳信号;
呼吸测量位置和心跳测量位置获取模块,评估位置A及其临近单元的呼吸信号和心跳信号质量,并根据结果确定最终的呼吸测量位置AB和心跳测量位置AH。
8.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质包括存储的程序,其中,在所述程序运行时控制所述计算机可读存储介质所在设备执行上述权利要求1-6任一项的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110939806.7A CN113786176B (zh) | 2021-08-17 | 2021-08-17 | 精确的毫米波雷达呼吸心跳测量方法、系统和存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110939806.7A CN113786176B (zh) | 2021-08-17 | 2021-08-17 | 精确的毫米波雷达呼吸心跳测量方法、系统和存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113786176A CN113786176A (zh) | 2021-12-14 |
CN113786176B true CN113786176B (zh) | 2024-03-15 |
Family
ID=78876006
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110939806.7A Active CN113786176B (zh) | 2021-08-17 | 2021-08-17 | 精确的毫米波雷达呼吸心跳测量方法、系统和存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113786176B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115025358B (zh) * | 2022-06-01 | 2024-03-08 | 广州小鹏汽车科技有限公司 | 呼吸引导方法、系统、车辆及存储介质 |
CN116577768B (zh) * | 2023-07-10 | 2023-09-22 | 厦门中为科学仪器有限公司 | 一种基于毫米波雷达的安检场景活体检测装置 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110507293A (zh) * | 2019-07-26 | 2019-11-29 | 中国电子科技集团公司第三十八研究所 | 一种超宽带穿墙雷达人体呼吸及心跳检测方法及系统 |
CN110879389A (zh) * | 2019-10-24 | 2020-03-13 | 中国人民解放军第四军医大学 | 基于多基地ir-uwb生物雷达信号的多人体目标识别定位方法 |
CN111568399A (zh) * | 2020-05-15 | 2020-08-25 | 中国人民解放军陆军军医大学 | 一种基于雷达的呼吸和心跳信号检测方法及系统 |
CN111680651A (zh) * | 2020-06-12 | 2020-09-18 | 武汉星巡智能科技有限公司 | 非接触式生命体征检测方法、装置、存储介质及系统 |
CN112168152A (zh) * | 2020-10-19 | 2021-01-05 | 北京清雷科技有限公司 | 呼吸心跳的检测方法、检测装置与计算机可读存储介质 |
CN112363139A (zh) * | 2020-11-02 | 2021-02-12 | 深圳大学 | 基于幅度特征的人体呼吸时长检测方法、装置和存储介质 |
CN112674738A (zh) * | 2020-12-07 | 2021-04-20 | 北京清雷科技有限公司 | 呼吸心跳信号的检测方法及装置 |
CN113064132A (zh) * | 2021-03-29 | 2021-07-02 | 中国人民解放军国防科技大学 | 一种基于连续信任函数的鲁棒雷达目标检测方法 |
-
2021
- 2021-08-17 CN CN202110939806.7A patent/CN113786176B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110507293A (zh) * | 2019-07-26 | 2019-11-29 | 中国电子科技集团公司第三十八研究所 | 一种超宽带穿墙雷达人体呼吸及心跳检测方法及系统 |
CN110879389A (zh) * | 2019-10-24 | 2020-03-13 | 中国人民解放军第四军医大学 | 基于多基地ir-uwb生物雷达信号的多人体目标识别定位方法 |
CN111568399A (zh) * | 2020-05-15 | 2020-08-25 | 中国人民解放军陆军军医大学 | 一种基于雷达的呼吸和心跳信号检测方法及系统 |
CN111680651A (zh) * | 2020-06-12 | 2020-09-18 | 武汉星巡智能科技有限公司 | 非接触式生命体征检测方法、装置、存储介质及系统 |
CN112168152A (zh) * | 2020-10-19 | 2021-01-05 | 北京清雷科技有限公司 | 呼吸心跳的检测方法、检测装置与计算机可读存储介质 |
CN112363139A (zh) * | 2020-11-02 | 2021-02-12 | 深圳大学 | 基于幅度特征的人体呼吸时长检测方法、装置和存储介质 |
CN112674738A (zh) * | 2020-12-07 | 2021-04-20 | 北京清雷科技有限公司 | 呼吸心跳信号的检测方法及装置 |
CN113064132A (zh) * | 2021-03-29 | 2021-07-02 | 中国人民解放军国防科技大学 | 一种基于连续信任函数的鲁棒雷达目标检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113786176A (zh) | 2021-12-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11850027B2 (en) | Multi-person vital signs monitoring using millimeter wave (mm-wave) signals | |
CN113786176B (zh) | 精确的毫米波雷达呼吸心跳测量方法、系统和存储介质 | |
CN113440120B (zh) | 一种基于毫米波雷达的人员呼吸心跳检测方法 | |
CN112754441B (zh) | 一种基于毫米波的非接触式心跳检测方法 | |
CN111856455A (zh) | 匹配不同雷达带宽的多目标心率呼吸测量方法及系统 | |
CN112336322A (zh) | 一种无接触式呼吸或心跳检测方法 | |
CN115460980A (zh) | 一种基于多普勒雷达的非接触式呼吸监测方法 | |
CN112220464A (zh) | 一种基于uwb雷达的人体呼吸心跳信号检测方法和系统 | |
CN115399747B (zh) | 一种基于毫米波雷达的多人生命体征监测方法 | |
Ling et al. | Non-contact heart rate monitoring based on millimeter wave radar | |
CN114259213A (zh) | 杂波背景下毫米波mimo雷达的邻近多目标生命体征检测方法 | |
CN114052709B (zh) | 一种稳健的毫米波雷达生命体征测量方法 | |
CN117838083A (zh) | 一种基于毫米波雷达的体征快速精确检测方法 | |
CN113702968A (zh) | 生命探测方法及其装置 | |
Liu et al. | Vital sign extraction in the presence of radar mutual interference | |
CN115736872A (zh) | 一种基于毫米波雷达实时计算心率的方法 | |
CN115736854A (zh) | 一种基于毫米波雷达的呼吸心跳监测系统 | |
CN114027805B (zh) | 一种基于微波雷达的生命体征测量方法及系统 | |
Jing et al. | Harmonics and intermodulation products-based fuzzy logic (HIPBFL) algorithm for vital sign frequency estimation using a UWB radar | |
JP2022185752A (ja) | 生体情報検出システム、プログラム、及び、生体情報検出方法 | |
CN113786177B (zh) | 生命体征信息提取方法、装置及电子设备 | |
CN115969351A (zh) | 一种基于毫米波雷达实时计算呼吸的方法 | |
KR101916591B1 (ko) | 레이더 신호의 주성분 분석을 이용한 생체 정보 결정 장치 및 방법 | |
KR102381262B1 (ko) | 레이더 시스템 및 그에서 수행되는 생체 신호 탐지방법 | |
CN114355329A (zh) | 一种调频连续波雷达生命体征距离检测的方法 |
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 |