CN114002655A - 超宽带连续波生物雷达数据处理方法及系统 - Google Patents
超宽带连续波生物雷达数据处理方法及系统 Download PDFInfo
- Publication number
- CN114002655A CN114002655A CN202111269767.0A CN202111269767A CN114002655A CN 114002655 A CN114002655 A CN 114002655A CN 202111269767 A CN202111269767 A CN 202111269767A CN 114002655 A CN114002655 A CN 114002655A
- Authority
- CN
- China
- Prior art keywords
- matrix
- range
- doppler spectrum
- time
- distance
- 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
- 238000003672 processing method Methods 0.000 title claims description 12
- 239000011159 matrix material Substances 0.000 claims abstract description 151
- 238000001228 spectrum Methods 0.000 claims abstract description 112
- 238000000034 method Methods 0.000 claims abstract description 30
- 238000012545 processing Methods 0.000 claims abstract description 29
- 230000006399 behavior Effects 0.000 claims abstract description 14
- 238000004458 analytical method Methods 0.000 claims abstract description 12
- 238000001914 filtration Methods 0.000 claims abstract description 7
- 239000013598 vector Substances 0.000 claims description 21
- 230000001629 suppression Effects 0.000 claims description 16
- 230000009466 transformation Effects 0.000 claims description 12
- 230000003044 adaptive effect Effects 0.000 claims description 7
- 238000013507 mapping Methods 0.000 claims description 3
- 238000012163 sequencing technique Methods 0.000 claims description 3
- 239000000126 substance Substances 0.000 claims description 3
- 238000013135 deep learning Methods 0.000 abstract description 12
- 230000009471 action Effects 0.000 description 7
- 238000004891 communication Methods 0.000 description 5
- 238000010801 machine learning Methods 0.000 description 5
- 230000036544 posture Effects 0.000 description 5
- 101100269850 Caenorhabditis elegans mask-1 gene Proteins 0.000 description 4
- 101000845188 Homo sapiens Tetratricopeptide repeat protein 4 Proteins 0.000 description 3
- 102100031279 Tetratricopeptide repeat protein 4 Human genes 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 238000010187 selection method Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000002790 cross-validation Methods 0.000 description 2
- 238000007405 data analysis Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000000018 DNA microarray Methods 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 238000004377 microelectronic Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- 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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
-
- 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
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明提出了一种超宽带连续波生物雷达数据处理方法及系统。该方法为：持续接收超宽带连续波生物雷达回波信号；对N个周期的雷达回波信号进行混频、滤波和快速傅里叶变换，得到N×M维的距离时间矩阵RT；在距离时间矩阵RT的列方向进行动目标显示处理,得到动目标显示处理后的距离时间矩阵RTMTI；根据距离时间矩阵RT或动目标显示处理后的距离时间矩阵RTMTI得到距离多普勒谱矩阵；选取距离多普勒谱矩阵中的最优距离单元，对其对应的时序信号进行时频分析，得到时间多普勒谱；将该时间多普勒谱用于人体行为识别。该方法能够在保证较高识别率的前提下大大减小深度学习网络模型参数，能获得更优秀的行为识别性能。
Description
技术领域
本发明涉及信号处理领域，具体涉及一种超宽带连续波生物雷达数据处理方法及系统。
背景技术
人体姿态识别技术对于描述人体姿态、预测人体行为至关重要，近年来广泛应用于智能控制、医疗康复、辅助教练、体感游戏、安防监控等领域。相比于可穿戴式惯性传感器，视频、压力、射频类等非入侵式环境传感器不需要用户佩戴，也不存在电池供电受限等问题；相比于视频、音频等感知设备，以雷达为代表的射频类传感器能够更好地保护用户隐私，能更好地适应各种光照、噪音环境，有些频段甚至可以穿墙探测。
目前用于人体姿态识别的雷达多采用连续波多普勒体制，为窄带雷达，不具备高的距离分辨率。针对这种体制雷达的数据分析方法主要是时频分析方法，如短时傅里叶变换，分析雷达回波信号中各种频率成分随时间的变化规律，从中提取特征指标后输入机器学习模型，或将时频分析方法得到的时间多普勒谱直接送入深度学习网络，以区分不同人体姿态动作类型。
随着微电子和芯片技术的发展，近年来国外涌现出多款超宽带连续波生物芯片雷达，这种体制的雷达同时具备高的距离分辨率和速度分辨率，能够为人体姿态识别提供更加多维的信息，从而提升识别效率。在得到雷达回波数据多维度的信息后，可从中提取特征指标后输入机器学习模型，或直接将多维度图谱送入深度学习网络，以区分不同人体姿态动作类型。
在雷达数据分析方法上不同于窄带连续波生物雷达，超宽带连续波生物雷达的主要区别在于具有随时间变化的多个距离分辨单元。在提取特征时，多采用滑窗重叠取样分析特征指标。由于采用滑窗重叠取样计算，时频分析方法的计算量相对傅里叶变换较大，因此通常情况下(即线上实时处理计算资源有限的情况下)，不可能对每个距离单元的时序信号做时频分析。因此，如何挑选合适的距离单元对应的时序信号进一步时频分析，是超宽带连续波生物雷达数据处理需要解决的关键技术问题。
发明内容
为了克服上述现有技术中存在的缺陷，本发明的目的是提供一种超宽带连续波生物雷达数据处理方法及系统。
为了实现本发明的上述目的，本发明提供了一种超宽带连续波生物雷达数据处理方法,包括以下步骤：
持续接收超宽带连续波生物雷达回波信号；
对N个周期的雷达回波信号进行混频、滤波和快速傅里叶变换，得到N×M维的距离时间矩阵RT，M为超宽带连续波生物雷达的距离分辨单元数量；
在距离时间矩阵RT的慢时间方向进行动目标显示处理,得到动目标显示处理后的距离时间矩阵RTMTI；
根据距离时间矩阵RT或动目标显示处理后的距离时间矩阵RTMTI得到距离多普勒谱矩阵；
选取距离多普勒谱矩阵中的最优距离单元，对其对应的时序信号进行时频分析，得到时间多普勒谱；
将该时间多普勒谱用于分类识别任务。
距离多普勒谱矩阵获得步骤为：对距离时间矩阵RT的所有行向量进行高通滤波，然后沿时间维做傅里叶变换，或直接对动目标显示处理后的距离时间矩阵RTMTI沿时间维做傅里叶变换，并将傅里叶变换后的零频点移到频谱的中间，得到距离多普勒矩阵RD,对距离多普勒矩阵RD进行归一化，找出距离多普勒矩阵RD的最大值，将距离多普勒矩阵RD的所有元素除以该最大值；对归一化的距离多普勒矩阵RD取绝对值和对数变换，将得到的矩阵定义为距离多普勒谱矩阵。
该方法能够在保证较高识别率的前提下大大减小深度学习网络模型参数，能获得更优秀的分类识别性能，特别适用于人体行为识别。
本方法的优选方案：得到距离多普勒谱矩阵后，先对距离多普勒谱矩阵进行自适应门限动态杂波抑制，然后再选取最优距离单元。
距离多普勒谱矩阵进行自适应门限动态杂波抑制步骤为：
A、计算距离多普勒谱矩阵的中值，记为b；定义变换因子d，并赋初值d＝1；
B、计算门限Th3＝b×d，将距离多普勒谱矩阵中小于门限Th3的值设为0，大于或等于门限的值设为1，得到距离多普勒谱矩阵掩模二值矩阵，记为Mask2；
C、计算距离多普勒谱矩阵掩模矩阵Mask2的四个角矩阵的非零元素的个数，分别记为Cn5、Cn6、Cn7、Cn8，计算四个角矩阵的非零元素个数的和，记为CNS2，其中CNS2＝Cn5+Cn6+Cn7+Cn8；
D、在距离多普勒谱矩阵掩模矩阵Mask2上和下分别取两个条带，计算距离多普勒谱矩阵掩模矩阵Mask2上下条带内非零元素的个数，分别记为Cnup、Cndown，选取二者中的较小值记为CNM；
E、当CNS2≥Th4且CNM≥Th5时，令d＝d-△d，更新d值，△d是步长，其中Th4为Mask2的四个角矩阵的非零元素个数和的门限值；Th5为Mask2上下条带内非零元素的个数门限值,为正整数；
F、重复执行步骤B至E，直到CNS2<Th4并且CNM<Th5终止，得到最终门限Th3＝b×d；将距离多普勒谱矩阵中大于或等于最终门限Th3的值保持原值，小于最终门限Th3的值等于Th3。
Mask2的四个角矩阵的大小均为
Mask2上下条带的大小均为
其中，
CM≤k2M，k2是小于1的正实数，Nfft和M分别为距离多普勒谱矩阵的行和列。Mask2的非零元素的数量门限值
q3是小于1的正实数。
该距离多普勒谱矩阵的自适应门限动态杂波抑制可以突显目标，抑制环境干扰，有益于后期特征的提取，减少网络参数、降低网络复杂度。
本方法的优选方案：最优距离单元的选取方法为：
对距离多普勒谱矩阵中元素大小进行排序，形成向量X；
设置百分比p，将向量X的变换范围映射到区间[0,100]，寻找向量X中p对应百分位数的元素，记为x；
将距离多普勒谱矩阵中小于等于x的元素维持值不变，大于x的元素赋值为x，得到更新的距离多普勒谱矩阵；
对更新后的距离多普勒谱矩阵按照多普勒频率轴方向求和，得到向量SRD，其大小为1×M，M为距离单元的个数；
计算向量SRD最大值对应的距离单元，将其作为最优距离单元。
这里选择多普勒频率轴方向的最大的距离单元作为最优距离单元，可以抑制环境中运动目标多径干扰的伪影带入的虚假分量，而且兼顾到稳定的多普勒最大值，从而可以比较科学、完整地包含动作行为频率随时间变化的信息。
本申请还提出了一种超宽带连续波生物雷达数据处理系统，包括处理器和存储器，所述处理器和存储器之间通信连接，相互通信，所述存储器用于存放至少一可执行指令，所述可执行指令使所述处理器执行如上述的超宽带连续波生物雷达数据处理方法的操作。
本发明的有益效果是：本发明能够在保证较高识别率的前提下大大减小深度学习网络模型参数；其采用最优距离单元选取方法得到的时间多普勒谱比用现有方法选取距离单元得到的时间多普勒谱在同样深度学习网络模型框架下能获得更优秀的行为识别性能，特别适用于人体行为识别。
本发明的附加方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解，其中：
图1是超宽带连续波生物雷达数据处理方法的流程图。
具体实施方式
下面详细描述本发明的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，仅用于解释本发明，而不能理解为对本发明的限制。
在本发明的描述中，除非另有规定和限定，需要说明的是，术语“安装”、“相连”、“连接”应做广义理解，例如，可以是机械连接或电连接，也可以是两个元件内部的连通，可以是直接相连，也可以通过中间媒介间接相连，对于本领域的普通技术人员而言，可以根据具体情况理解上述术语的具体含义。
如图1所示，本发明提供了一种超宽带连续波生物雷达数据处理方法，该方法具体介绍如下：
超宽带连续波生物雷达可发射和接收步进频或调频连续波信号，因此持续接收超宽带连续波生物雷达回波信号。
对每个周期的雷达回波信号进行混频、滤波和快速傅里叶变换，即进行下变频、距离压缩处理，得到包含M个距离单元的一维行向量；对N个周期信号进行该处理并按照列方向进行排列，得到一个N×M维的距离时间矩阵RT，M为超宽带连续波生物雷达的距离分辨单元数量。
如果直接对距离时间矩阵RT进行取绝对值和对数变换，得到距离时间谱矩阵，此时得到的谱矩阵包含大量静态环境杂波，无法直接用于动作行为识别。因此，本实施例先在距离时间矩阵RT的慢时间方向(列方向)进行动目标显示处理(MTI),得到动目标显示处理后的距离时间矩阵RTMTI。动目标显示处理(MTI)最简单的一种执行方法就是对N×M维距离时间矩阵RT沿着慢时间方向(列方向)做差分，本实施例中的动目标显示处理(MTI)可采用该方法。
根据距离时间矩阵RT或动目标显示处理后的距离时间矩阵RTMTI得到距离多普勒谱矩阵。
本实施例中，距离多普勒谱矩阵获得步骤为：对距离时间矩阵RT或动目标显示处理后的距离时间矩阵RTMTI的所有行向量进行高通滤波，然后沿时间维做傅里叶变换，并将零频点移到频谱的中间，得到距离多普勒矩阵RD,对距离多普勒矩阵RD进行归一化，找出距离多普勒矩阵RD的最大值，将距离多普勒矩阵RD的所有元素除以该最大值；对归一化的距离多普勒矩阵RD取绝对值和对数变换，将得到的矩阵定义为距离多普勒谱矩阵。
对距离多普勒谱矩阵进行自适应门限动态杂波抑制。
本实施例中，距离多普勒谱矩阵自适应门限动态杂波抑制步骤为：
A、计算距离多普勒谱矩阵的中值，记为b；定义变换因子d，并赋初值d＝1。
B、计算门限Th3＝b×d，将距离多普勒谱矩阵中小于门限Th3的值设为0，大于或等于门限的值设为1，得到距离多普勒谱矩阵掩模二值矩阵，记为Mask2。
C、计算距离多普勒谱矩阵掩模矩阵Mask2的四个角矩阵的非零元素的个数，分别记为Cn5、Cn6、Cn7、Cn8，计算Mask2的四个角矩阵的非零元素个数的和，记为CNS2，其中CNS2＝Cn5+Cn6+Cn7+Cn8。这里四个角矩阵的范围应根据距离多普勒谱矩阵的大小来选取，本实施例中，Mask2的四个角矩阵的大小均为
其中，
CM≤k2M，k2是小于1的正实数，可应根据雷达系统设备的噪声系数情况具体调整，本实施例中优选但不限于为0.02，即
以及CM≤M/50，Nfft和M分别为距离多普勒谱矩阵的行数和列数。
该步骤中，由于经过频谱变换和零频点中间移动，距离多普勒谱矩阵中的目标相对集中在零频(矩阵的中心点)附近，因此四个角矩阵没有目标出现，因此省去对目标位置的判断步骤。
D、在距离多普勒谱矩阵掩模矩阵Mask2上和下分别取两个条带，大小均为
计算距离多普勒谱矩阵掩模矩阵Mask2上下条带内非零元素的个数，分别记为Cnup、Cndown，选取二者中的较小值记为CNM，有CNM＝min(Cnup,Cndown)。
E、设置Mask2的四个角矩阵的非零元素个数和的门限值Th4与Mask2上下条带内非零元素的个数门限值Th5，Th5可以根据频谱的具体噪声情况进行微调，本实施例中优选但不限于为1，当CNS2≥Th4且CNM≥Th5时，令d＝d-△d，更新d值，△d是步长，其中Mask2的非零元素的数量门限值
q3是小于1的正实数，可以根据频谱的具体噪声情况进行微调，本实施例中q3优选但不限于为0.05,Mask2上下条带内非零元素的个数门限值Th5为正整数。
F、重复执行步骤B至E，直到CNS2<Th4并且CNM<Th5终止，得到最终门限Th3＝b×d；将距离多普勒谱矩阵中大于或等于最终门限Th3的值保持原值，小于最终门限Th3的值等于Th3，至此，完成距离多普勒谱矩阵的自适应门限动态杂波抑制。
选取距离多普勒谱矩阵中的最优距离单元，最好是基于自适应杂波抑制后的距离多普勒谱矩阵选取最优距离单元，对其对应的时序信号进行时频分析，如短时傅里叶变换，得到时间多普勒谱；将该时间多普勒谱用于人体行为识别。
具体地，最优距离单元的选取方法为：
对距离多普勒谱矩阵中元素大小进行排序，形成向量X。
设置百分比p，将向量X的变换范围映射到区间[0,100]，寻找向量X中p对应百分位数的元素，记为x，本p为小于100的正实数，可根据距离多普勒谱矩阵的噪声背景调整，本实施例中优选为98。
将距离多普勒谱矩阵中小于等于x的元素维持值不变，大于x的元素赋值为x，得到更新的距离多普勒谱矩阵。
对更新后的距离多普勒谱矩阵按照多普勒频率轴方向(列方向)求和，得到向量SRD，其大小为1×M，M为距离单元的个数。
计算向量SRD最大值对应的距离单元，将其作为最优距离单元。
这里在最优距离单元选取的时候，之所以选择和最大而不是距离多普勒上多普勒值最大的距离单元，是因为这样可以抑制环境中运动目标多径干扰的伪影带入的虚假分量，而且兼顾到稳定的多普勒最大值，从而可以比较科学、完整地包含动作行为频率随时间变化的信息。
优选的，在得到时间多普勒谱后，可以参照距离多普勒谱矩阵自适应门限动态杂波抑制方法，对时间多普勒谱进行动态杂波抑制，然后再用于人体行为识别。
本实施例中，还可根据动目标显示处理后的距离时间矩阵RTMTI得到距离时间谱矩阵，然后对对距离时间谱矩阵进行自适应门限动态杂波抑制。
距离时间谱矩阵获得步骤为：对距离时间矩阵RTMTI进行归一化，找出距离时间矩阵RTMTI中的最大值，然后将距离时间矩阵RTMTI的所有元素除以该最大值；对归一化的距离时间矩阵RTMTI取绝对值和对数变换，得到距离时间谱矩阵。
距离时间谱矩阵自适应门限动态杂波抑制步骤为：
S1，计算距离时间谱矩阵的中值，记为a；定义变换因子c，并赋初值c＝1。
S2，计算门限Th1＝a×c，将距离时间谱矩阵中小于门限Th1的值设为0，大于或等于门限的值设为1，得到距离时间谱矩阵掩模二值矩阵，记为Mask1。
S3，计算距离时间谱矩阵掩模矩阵Mask1的四个角矩阵的非零元素的个数，分别记为Cn1、Cn2、Cn3、Cn4，四个角的顺序不分先后。这里四个角矩阵的范围应根据距离时间谱矩阵的大小来选取，本实施例中，Mask1的四个角矩阵的大小均为CN×CM,CN≤k1N，CM≤k1M，N和M分别为距离时间谱矩阵的行数和列数，这里系数k1是小于1的正实数，可根据雷达系统设备的噪声系数情况具体调整，通常为1/50。
S4，排除Cn1、Cn2、Cn3、Cn4中存在目标的位置：分别判断四个角矩阵的非零元素的个数Cn1、Cn2、Cn3、Cn4是否大于或等于设定值,如果是，则认为大于或等于设定值的非零元素的个数所对应的角矩阵处为包含目标，不对该角矩阵的元素个数进行求和,本实施例中，设定值优选但不限于为q1×CN×CM，q1为小于1的正实数，本实施例中优选为9/10，即如果是目标，非零元素的个数则会占据90％的角矩阵位置；如果存在小于设定值的非零元素的个数，则认为该非零元素的个数所对应的角矩阵为非目标处角矩阵，计算所有非目标处角矩阵非零元素个数的和，记为CNS1，当CNS1≥Th2时，令c＝c-△c，更新c值，其中△c是步长，Th2是Mask1的非目标处角矩阵非零元素的数量门限值，本实施例中，非零元素的数量门限值Th2＝q2CN×CM,q2是小于1的正实数，本实施例中，q2优选为1/10。这里所提到的目标为检测时需要识别提取的有效信号。
S5，重复步骤S2至S4，直到CNS1<Th2终止，得到最终门限Th1＝a×c；将距离时间谱矩阵中大于或等于最终门限Th1的值保持原值，小于最终门限Th1的值等于Th1。至此，完成距离时间谱矩阵的自适应门限动态杂波抑制。
采用经上述自适应门限动态杂波抑制后的距离时间谱矩阵、距离多普勒谱矩阵进行人体行为识别，可以突显目标，抑制环境和系统抖动引入的干扰，同样有益于后期特征的提取，减少网络参数、降低网络复杂度。
本申请还提出了一种超宽带连续波生物雷达数据处理系统，包括处理器和存储器，所述处理器和存储器之间通信连接，相互通信，所述存储器用于存放至少一可执行指令，所述可执行指令使所述处理器执行如上述的超宽带连续波生物雷达数据处理方法的操作。
本申请所侧重于机器学习(包含深度学习)网络模型前端的雷达数据处理方法，不对机器学习(包含深度学习)网络模型进行描述，具体实施时，采用现有的机器学习(包含深度学习)网络模型即可。
为了验证提出方法的有效性和可靠性，对本申请所提出的方法进行了大量的实验。首先用超宽带生物雷达采集了36人的3240次行为动作数据，将分析得到的距离时间谱、距离多普勒谱和时间多普勒谱的数据集，通过上下翻转和左右翻转等数据扩充方法后，送入搭建的深度学习网络进行5折交叉验证分析，网络中所有的随机种子数进行固定，计算5折交叉验证的平均识别性能，评价指标包括准确率、精确度、灵敏度、特异性和F1-Score。大量实验结果证实：用本方法得到的距离时间谱比原始距离时间谱的对于动作行为的识别性能要好；用本方法得到的距离时间谱、距离多普勒谱和时间多普勒谱能够在保证较高识别率的前提下大大减小深度学习网络模型参数；用最优距离单元选取方法得到的时间多普勒谱比用现有方法选取距离单元得到的时间多普勒谱在同样深度学习网络模型框架下能获得更优秀的行为识别性能。
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
尽管已经示出和描述了本发明的实施例，本领域的普通技术人员可以理解：在不脱离本发明的原理和宗旨的情况下可以对这些实施例进行多种变化、修改、替换和变型，本发明的范围由权利要求及其等同物限定。
Claims (9)
1.一种超宽带连续波生物雷达数据处理方法,其特征在于，包括以下步骤：
持续接收超宽带连续波生物雷达回波信号；
对N个周期的雷达回波信号进行混频、滤波和快速傅里叶变换，得到N×M维的距离时间矩阵RT，M为超宽带连续波生物雷达的距离分辨单元数量；
在距离时间矩阵RT的慢时间方向进行动目标显示处理,得到动目标显示处理后的距离时间矩阵RTMTI；
根据距离时间矩阵RT或动目标显示处理后的距离时间矩阵RTMTI得到距离多普勒谱矩阵；
选取距离多普勒谱矩阵中的最优距离单元，对其对应的时序信号进行时频分析，得到时间多普勒谱；
将该时间多普勒谱用于分类识别任务。
2.根据权利要求1所述的超宽带连续波生物雷达数据处理方法,其特征在于，距离多普勒谱矩阵获得步骤为：对距离时间矩阵RT的所有行向量进行高通滤波，然后沿时间维做傅里叶变换，或直接对动目标显示处理后的距离时间矩阵RTMTI沿时间维做傅里叶变换，并将傅里叶变换后的零频点移到频谱的中间，得到距离多普勒矩阵RD,对距离多普勒矩阵RD进行归一化，找出距离多普勒矩阵RD的最大值，将距离多普勒矩阵RD的所有元素除以该最大值；对归一化的距离多普勒矩阵RD取绝对值和对数变换，将得到的矩阵定义为距离多普勒谱矩阵。
3.根据权利要求1所述的超宽带连续波生物雷达数据处理方法,其特征在于，得到距离多普勒谱矩阵后，先对距离多普勒谱矩阵进行自适应门限动态杂波抑制，然后再选取最优距离单元。
4.根据权利要求3所述的超宽带连续波生物雷达数据处理方法,其特征在于，距离多普勒谱矩阵进行自适应门限动态杂波抑制步骤为：
A、计算距离多普勒谱矩阵的中值，记为b；定义变换因子d，并赋初值d＝1；
B、计算门限Th3＝b×d，将距离多普勒谱矩阵中小于门限Th3的值设为0，大于或等于门限的值设为1，得到距离多普勒谱矩阵掩模二值矩阵，记为Mask2；
C、计算距离多普勒谱矩阵掩模矩阵Mask2的四个角矩阵的非零元素的个数，分别记为Cn5、Cn6、Cn7、Cn8，计算四个角矩阵的非零元素个数的和，记为CNS2，其中CNS2＝Cn5+Cn6+Cn7+Cn8；
D、在距离多普勒谱矩阵掩模矩阵Mask2上和下分别取两个条带，计算距离多普勒谱矩阵掩模矩阵Mask2上下条带内非零元素的个数，分别记为Cnup、Cndown，选取二者中的较小值记为CNM；
E、当CNS2≥Th4且CNM≥Th5时，令d＝d-△d，更新d值，△d是步长，其中Th4为Mask2的四个角矩阵的非零元素个数和的门限值；Th5为Mask2上下条带内非零元素的个数门限值,为正整数；
F、重复执行步骤B至E，直到CNS2<Th4并且CNM<Th5终止，得到最终门限Th3＝b×d；将距离多普勒谱矩阵中大于或等于最终门限Th3的值保持原值，小于最终门限Th3的值等于Th3。
5.根据权利要求1至4任一项所述的超宽带连续波生物雷达数据处理方法,其特征在于，最优距离单元的选取方法为：
对距离多普勒谱矩阵中元素大小进行排序，形成向量X；
设置百分比p，将向量X的变换范围映射到区间[0,100]，寻找向量X中p对应百分位数的元素，记为x；
将距离多普勒谱矩阵中小于或等于x的元素维持值不变，大于x的元素赋值为x，得到更新的距离多普勒谱矩阵；
对更新后的距离多普勒谱矩阵按照多普勒频率轴方向求和，得到向量SRD，其大小为1×M，M为距离单元的个数；
计算向量SRD最大值对应的距离单元，将其作为最优距离单元。
6.根据权利要求4所述的超宽带连续波生物雷达数据处理方法,其特征在于，Mask2的四个角矩阵的大小均为
Mask2上下条带的大小均为
其中，
CM≤k2M，k2是小于1的正实数，Nfft和M分别为距离多普勒谱矩阵的行和列。
7.根据权利要求4所述的超宽带连续波生物雷达数据处理方法,其特征在于，Mask2的非零元素的数量门限值
q3是小于1的正实数。
8.根据权利要求1所述的超宽带连续波生物雷达数据处理方法,其特征在于，对时间多普勒谱矩阵进行动态杂波抑制，然后再用于人体行为识别。
9.一种超宽带连续波生物雷达数据处理系统，其特征在于，包括处理器和存储器，所述处理器和存储器之间通信连接，相互通信，所述存储器用于存放至少一可执行指令，所述可执行指令使所述处理器执行如权利要求1至8任一项所述的超宽带连续波生物雷达数据处理方法的操作。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111269767.0A CN114002655B (zh) | 2021-10-29 | 2021-10-29 | 超宽带连续波生物雷达数据处理方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111269767.0A CN114002655B (zh) | 2021-10-29 | 2021-10-29 | 超宽带连续波生物雷达数据处理方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114002655A true CN114002655A (zh) | 2022-02-01 |
| CN114002655B CN114002655B (zh) | 2024-06-25 |
Family
ID=79925038
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111269767.0A Active CN114002655B (zh) | 2021-10-29 | 2021-10-29 | 超宽带连续波生物雷达数据处理方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114002655B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114580535A (zh) * | 2022-03-04 | 2022-06-03 | 中国人民解放军空军军医大学 | 基于注意力多视角学习网络的多基雷达人体行为融合识别方法、装置及介质 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107607917A (zh) * | 2017-08-04 | 2018-01-19 | 西安电子科技大学 | 一种基于两级多普勒相关性判别的海杂波抑制方法 |
| CN110286368A (zh) * | 2019-07-10 | 2019-09-27 | 北京理工大学 | 一种基于超宽带雷达的老人跌倒检测方法 |
| CN111568399A (zh) * | 2020-05-15 | 2020-08-25 | 中国人民解放军陆军军医大学 | 一种基于雷达的呼吸和心跳信号检测方法及系统 |
| CN112363133A (zh) * | 2020-10-12 | 2021-02-12 | 山东省科学院自动化研究所 | 一种基于多普勒雷达的目标检测方法、介质及装置 |
| CN112816942A (zh) * | 2019-11-15 | 2021-05-18 | 安讯士有限公司 | 方法、计算机程序产品、装置和调频连续波雷达系统 |
-
2021
- 2021-10-29 CN CN202111269767.0A patent/CN114002655B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107607917A (zh) * | 2017-08-04 | 2018-01-19 | 西安电子科技大学 | 一种基于两级多普勒相关性判别的海杂波抑制方法 |
| CN110286368A (zh) * | 2019-07-10 | 2019-09-27 | 北京理工大学 | 一种基于超宽带雷达的老人跌倒检测方法 |
| CN112816942A (zh) * | 2019-11-15 | 2021-05-18 | 安讯士有限公司 | 方法、计算机程序产品、装置和调频连续波雷达系统 |
| CN111568399A (zh) * | 2020-05-15 | 2020-08-25 | 中国人民解放军陆军军医大学 | 一种基于雷达的呼吸和心跳信号检测方法及系统 |
| CN112363133A (zh) * | 2020-10-12 | 2021-02-12 | 山东省科学院自动化研究所 | 一种基于多普勒雷达的目标检测方法、介质及装置 |
Non-Patent Citations (1)
| Title |
|---|
| 祁富贵 等;: "SFCW生物雷达人体细粒度运动信号微多普勒特征增强方法研究", 《中国医疗设备》, vol. 31, no. 02, 25 February 2016 (2016-02-25), pages 39 - 43 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114580535A (zh) * | 2022-03-04 | 2022-06-03 | 中国人民解放军空军军医大学 | 基于注意力多视角学习网络的多基雷达人体行为融合识别方法、装置及介质 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114002655B (zh) | 2024-06-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110286368B (zh) | 一种基于超宽带雷达的老人跌倒检测方法 | |
| Dekker et al. | Gesture recognition with a low power FMCW radar and a deep convolutional neural network | |
| CN111695420A (zh) | 一种手势识别方法以及相关装置 | |
| CN111157988B (zh) | 一种基于rdtm与atm融合的手势雷达信号处理方法 | |
| CN107561509B (zh) | 一种机载毫米波雷达电力线检测方法 | |
| CN111142102B (zh) | 一种呼吸数据计算方法以及相关设备 | |
| CN111178331A (zh) | 雷达图像识别系统、方法、设备及计算机可读存储介质 | |
| KR102221839B1 (ko) | 잡음 환경 내 실시간 심장 박동수 탐지 장치 및 그 방법 | |
| CN112990026A (zh) | 基于对抗训练的无线信号感知模型构建、感知方法及系统 | |
| CN113447905A (zh) | 双毫米波雷达人体跌倒检测装置及检测方法 | |
| CN111901028B (zh) | 一种基于多天线上csi的人体行为识别方法 | |
| CN115877376A (zh) | 基于多头自注意力机制的毫米波雷达手势识别方法及识别系统 | |
| CN114002655A (zh) | 超宽带连续波生物雷达数据处理方法及系统 | |
| CN113064483A (zh) | 一种手势识别方法以及相关装置 | |
| Zhang et al. | Temporal-range-doppler features interpretation and recognition of hand gestures using MmW FMCW radar sensors | |
| CN115469303A (zh) | 一种检测人体姿态和生命体征的认知生物雷达方法与装置 | |
| Krishnamoorthy et al. | Temporal and spectral analysis of internal carotid artery Doppler signal for normal and abnormal flow detection | |
| Yang et al. | Extraction and denoising of human signature on radio frequency spectrums | |
| CN114002656B (zh) | 超宽带连续波生物雷达数据杂波抑制方法及系统 | |
| CN112380903A (zh) | 一种基于WiFi-CSI信号增强的人体活动识别方法 | |
| CN114002656A (zh) | 超宽带连续波生物雷达数据杂波抑制方法及系统 | |
| CN115937977A (zh) | 一种基于多维特征融合的少样本人体动作识别方法 | |
| Song et al. | High-accuracy gesture recognition using mm-wave radar based on convolutional block attention module | |
| Li et al. | Objective evaluation of clutter suppression for micro-Doppler spectrograms of hand gesture/sign language based on pseudo-reference image | |
| CN115345216A (zh) | 一种融合先验信息的fmcw雷达干扰消除方法 |
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 |