CN114869296A - 定位qrs-t波参考水平的方法、系统、设备及介质 - Google Patents

定位qrs-t波参考水平的方法、系统、设备及介质 Download PDF

Info

Publication number
CN114869296A
CN114869296A CN202210797009.4A CN202210797009A CN114869296A CN 114869296 A CN114869296 A CN 114869296A CN 202210797009 A CN202210797009 A CN 202210797009A CN 114869296 A CN114869296 A CN 114869296A
Authority
CN
China
Prior art keywords
section
qrs
band
pass filtering
linear regression
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
CN202210797009.4A
Other languages
English (en)
Other versions
CN114869296B (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.)
Changsha Biovision Software Technology Co ltd
Original Assignee
Changsha Biovision Software Technology Co ltd
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 Changsha Biovision Software Technology Co ltd filed Critical Changsha Biovision Software Technology Co ltd
Priority to CN202210797009.4A priority Critical patent/CN114869296B/zh
Publication of CN114869296A publication Critical patent/CN114869296A/zh
Application granted granted Critical
Publication of CN114869296B publication Critical patent/CN114869296B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316—Modalities, i.e. specific diagnostic methods
    • A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/346—Analysis of electrocardiograms
    • A61B5/349—Detecting specific parameters of the electrocardiograph cycle
    • A61B5/366—Detecting abnormal QRS complex, e.g. widening
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316—Modalities, i.e. specific diagnostic methods
    • A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/346—Analysis of electrocardiograms
    • A61B5/349—Detecting specific parameters of the electrocardiograph cycle
    • A61B5/355—Detecting T-waves
    • 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/7225—Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
    • 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)
  • Cardiology (AREA)
  • Molecular Biology (AREA)
  • General Health & Medical Sciences (AREA)
  • Biophysics (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Pathology (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Signal Processing (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physiology (AREA)
  • Psychiatry (AREA)
  • Power Engineering (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Abstract

本发明公开了一种定位QRS‑T波参考水平的方法、系统、设备及介质，方法包括：获取输入信号；根据输入信号的PR段时长和PR段采样数判断输入信号的PR段是否满足一元线性回归条件，若输入信号的PR段满足一元线性回归条件，将输入信号进行滤波，得到带通滤波信号；将带通滤波信号的PR段进行一元线性回归，得到一元线性回归方程；根据一元线性回归方程判断带通滤波信号的PR段是否存在特殊波形；若带通滤波信号的PR段不存在特殊波形，判断带通滤波信号的PR段是否倾斜，若带通滤波信号的PR段不存在倾斜，根据一元线性回归方程定位QRS‑T波参考水平。多轮决策判断PR段是否适用，并提高定位QRS‑T波参考水平的精度。

Description

定位QRS-T波参考水平的方法、系统、设备及介质
技术领域
本发明涉及心电图定位测量技术领域，特别涉及一种定位QRS-T波参考水平的方法、系统、设备及介质。
背景技术
心血管疾病发病率越来越高，已成为目前对人类生命构成最大的危害疾病之一。而心电图是检测、诊断和预防这类疾病的主要有效方法。传统的PC监护仪价格昂贵，体积庞大，不便移动且主要集中在医院，而无法实时监护患者的病情，给病人和医生带来很大不便。针对这些弊端，也随着嵌入式和网络通讯技术的飞速发展，克服传统的移动心电监测仪应运而生，它具有低成本、体积小、可靠性高、操作简单等优点，适用于个人，家庭，中小型医院和社区医疗单位。移动心电监护设备能够快速对动态心电信号进行实时检测和心脏疾病的自动诊断，以便能为用户提供紧急救护、疾病预警、医学咨询和指导等多种服务，必然，选择一种准确快速适用的心率检测算法，意义重大。
目前现有技术大多都是针对QRS波群进行检测计算，该方法运算难度大，处理大批量数据时效率不高，无法针对QRS波群的特定段进行检测分析，能针对QRS波群的特定段进行检测分析的方法又无法针对PR段倾斜、PR段过短等不利于直接定位QRS-T波参考水平的情况，调整定位方式。
发明内容
本发明旨在至少解决现有技术中存在的技术问题。为此，本发明提出一种定位QRS-T波参考水平的方法、系统、设备及介质，能够通过多轮决策来判断PR段是否适用于定位QRS-T波参考水平，并得到QRS-T波参考水平。
第一方面，本发明的实施例提供了一种定位QRS-T波参考水平的方法，包括以下步骤：
获取输入信号；
根据所述输入信号的PR段时长和PR段采样数判断所述输入信号的PR段是否满足一元线性回归条件，若所述输入信号的PR段满足一元线性回归条件，将所述输入信号进行滤波，得到带通滤波信号；
将所述带通滤波信号的PR段进行一元线性回归，得到一元线性回归方程；
根据所述一元线性回归方程判断所述带通滤波信号的PR段是否存在特殊波形；若所述带通滤波信号的PR段不存在特殊波形，判断所述带通滤波信号的PR段是否倾斜，若所述带通滤波信号的PR段不存在倾斜，根据所述一元线性回归方程定位QRS-T波参考水平。
根据本发明的实施例，至少具有如下技术效果：
本发明通过输入信号的PR段时长和PR段采样数，判断PR段是否满足一元线性回归条件，判断能否使用PR段进行一元线性回归来计算平均电势，通过带通滤波，排除心电信号的基线漂移、肌电干扰、工频噪声等影响，使得进行一元线性回归时数据更加精确并增强了PR段拟合效果，使用一元线性回归能够拟合PR段直线，作为决策依据判断是否存在特殊波形或者倾斜，判断拟合PR段能否用于平均电势的统计，满足条件时将一元线性回归方程用于平均电势的统计，提高了QRS-T波参考水平的定位精度，整体上看进行多轮决策，针对不同的PR段形态，使用相应的QRS-T波参考水平的计算方式，使QRS-T波参考水平的定位精度更好。
根据本发明的一些实施例，所述根据所述输入信号的PR段时长和PR段采样数判断所述输入信号的PR段是否满足一元线性回归条件，包括以下步骤：
判断所述PR段时长是否小于最短时长阈值或所述PR段采样数是否小于最少采样数阈值；
所述PR段时长不小于最短时长阈值且所述PR段采样数不小于最少采样数阈值，所述输入信号的PR段满足一元线性回归条件；所述PR段时长小于最短时长阈值或所述PR段采样数小于最少采样数阈值，所述输入信号的PR段不满足一元线性回归条件，令QRS-T起始电势为QRS-T波参考水平。
根据本发明的一些实施例，所述将所述输入信号进行滤波，包括以下步骤：
对所述输入信号进行0.05Hz的高通滤波，得到高通滤波信号；
对所述高通滤波信号进行35Hz的低通滤波，得到带通滤波信号。
根据本发明的一些实施例，所述根据所述一元线性回归方程判断所述带通滤波信号的PR段是否存在特殊波形，包括以下步骤：
计算所述一元线性回归方程的回归平方和以及总平分和；
根据所述回归平方和以及所述总平分和计算判定系数；
所述判定系数若小于判定系数阈值，则存在特殊波形。
根据本发明的一些实施例，所述判断所述带通滤波信号的PR段是否存在倾斜，包括以下步骤：
根据所述一元回归线性方程获取所述带通滤波信号的PR段的拟合斜率；
判断所述拟合斜率是否属于正常斜率范围，若所述拟合斜率属于正常斜率范围，则所述带通滤波信号的PR段不存在倾斜，若所述拟合斜率不属于正常斜率范围，则所述带通滤波信号的PR段存在倾斜。
根据本发明的一些实施例，若所述带通滤波信号的PR段存在特殊波形或存在倾斜，所述定位QRS-T波参考水平的方法还包括以下步骤：
令QRS-T起始电势为QRS-T波参考水平。
根据本发明的一些实施例，所述根据所述一元线性回归方程定位QRS-T波参考水平，包括以下步骤：
计算所述一元线性回归方程在所述带通滤波信号的PR段上的残差数组；
将所述残差数组取绝对值，并增序排序得到残差绝对值数组；
对于所述残差绝对值数组的前
Figure 875317DEST_PATH_IMAGE002
个元素，获取并累加所述元素对应采样点的电势
Figure 74348DEST_PATH_IMAGE004
得到
Figure 713140DEST_PATH_IMAGE006
，所述
Figure 419934DEST_PATH_IMAGE008
表示所述判定系数，所述
Figure 703147DEST_PATH_IMAGE010
表示所述最少采样数阈值；
根据所述
Figure 119085DEST_PATH_IMAGE012
和所述
Figure 854960DEST_PATH_IMAGE014
，计算拟合平均电势
Figure 601330DEST_PATH_IMAGE016
；
令所述拟合平均电势
Figure 863684DEST_PATH_IMAGE016
为QRS-T波参考水平。
第二方面，本发明的实施例提供了一种定位QRS-T波参考水平的系统，包括：
信号获取模块，获取输入信号；
带通滤波模块，根据所述输入信号的PR段时长和PR段采样数判断所述输入信号的PR段是否满足一元线性回归条件，若所述输入信号的PR段满足一元线性回归条件，将所述输入信号进行滤波，得到带通滤波信号；
一元线性回归模块，将所述带通滤波信号的PR段进行一元线性回归，得到一元线性回归方程；
参考水平判断模块，根据所述一元线性回归方程判断所述带通滤波信号的PR段是否存在特殊波形；若所述带通滤波信号的PR段不存在特殊波形，判断所述带通滤波信号的PR段是否倾斜，若所述带通滤波信号的PR段不存在倾斜，根据所述一元线性回归方程定位QRS-T波参考水平。
第三方面，本发明的实施例提供了一种电子设备，包括至少一个控制处理器和用于与所述至少一个控制处理器通信连接的存储器；所述存储器存储有可被所述至少一个控制处理器执行的指令，所述指令被所述至少一个控制处理器执行，以使所述至少一个控制处理器能够执行第一方面所述的定位QRS-T波参考水平的方法。
第四方面，本发明的实施例提供了一种计算机存储介质，所述计算机可读存储介质存储有计算机可执行指令，所述计算机可执行指令用于使计算机执行第一方面所述的定位QRS-T波参考水平的方法。
需要注意的是，本发明的第二方面至第四方面与现有技术之间的有益效果与上述的定位QRS-T波参考水平的方法的有益效果相同，此处不再细述。
本发明的附加方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解，其中：
图1是本发明一个实施例提供的一种定位QRS-T波参考水平的方法流程图；
图2是本发明一个实施例提供的判断输入信号的PR段是否满足一元线性回归条件的流程图；
图3是本发明一个实施例提供的输入信号进行滤波的流程图；
图4是本发明一个实施例提供的判断输入信号的PR段是否存在特殊波形的流程图；
图5是本发明一个实施例提供的判断输入信号的PR段是否倾斜的流程图；
图6是本发明一个实施例提供的根据一元线性回归方程定位QRS-T波参考水平的流程图；
图7是本发明一个实施例提供的一种定位QRS-T波参考水平的系统结构图；
图8是本发明一个实施例提供的一种电子设备结构图；
图9是本发明一个实施例提供的一种定位QRS-T波参考水平的方法示意图。
具体实施方式
下面详细描述本发明的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，仅用于解释本发明，而不能理解为对本发明的限制。
在本发明的描述中，需要理解的是，术语“中心”、“纵向”、“横向”、“长度”、“宽度”、“厚度”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”、“轴向”、“径向”、“周向”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。此外，限定有“第一”、“第二”的特征可以明示或者隐含地包括一个或者更多个该特征。在本发明的描述中，除非另有说明，“多个”的含义是两个或两个以上。
在本发明的描述中，需要说明的是，除非另有明确的规定和限定，术语“安装”、“相连”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通。对于本领域的普通技术人员而言，可以具体情况理解上述术语在本发明中的具体含义。
QRS-T波参考水平是定位测量的重要中间量，是测量QRS波群振幅、T波振幅等诊断依据的前提条件，为了解决无法针对PR段倾斜、PR段过短等不利于直接定位QRS-T波参考水平的情况调整定位方式的问题，本发明提出一种定位QRS-T波参考水平的方法、系统、设备及介质，能够通过多轮决策来判断PR段是否适用于定位QRS-T波参考水平，并得到QRS-T波参考水平。
参照图1，在本发明的一些实施例中，提供了一种定位QRS-T波参考水平的方法，包括以下步骤：
步骤S101、获取输入信号。
步骤S200、根据输入信号的PR段时长和PR段采样数判断输入信号的PR段是否满足一元线性回归条件，若输入信号的PR段满足一元线性回归条件，将输入信号进行滤波，得到带通滤波信号。
步骤S300、将带通滤波信号的PR段进行一元线性回归，得到一元线性回归方程。
步骤S400、根据一元线性回归方程判断带通滤波信号的PR段是否存在特殊波形；若带通滤波信号的PR段不存在特殊波形，判断带通滤波信号的PR段是否倾斜，若带通滤波信号的PR段不存在倾斜，根据一元线性回归方程定位QRS-T波参考水平。
根据本发明实施例的一种定位QRS-T波参考水平的方法，通过输入信号的PR段时长和PR段采样数，判断PR段是否满足一元线性回归条件，判断能否使用PR段进行一元线性回归来计算平均电势，通过带通滤波，排除心电信号的基线漂移、肌电干扰、工频噪声等影响，使得进行一元线性回归时数据更加精确并增强了PR段拟合效果，使用一元线性回归能够拟合PR段直线，作为决策依据判断是否存在特殊波形或者倾斜，判断拟合PR段能否用于平均电势的统计，满足条件时将一元线性回归方程用于平均电势的统计，提高了QRS-T波参考水平的定位精度，整体上看进行多轮决策，针对不同的PR段形态，使用相应的QRS-T波参考水平的计算方式，使QRS-T波参考水平的定位精度更好。
参照图2，在本发明的一些实施例中，根据输入信号的PR段时长和PR段采样数判断输入信号的PR段是否满足一元线性回归条件，包括以下步骤：
步骤S201、判断PR段时长是否小于最短时长阈值或PR段采样数是否小于最少采样数阈值。
步骤S202、PR段时长不小于最短时长阈值且PR段采样数不小于最少采样数阈值，输入信号的PR段满足一元线性回归条件；PR段时长小于最短时长阈值或PR段采样数小于最少采样数阈值，输入信号的PR段不满足一元线性回归条件，令QRS-T起始电势为QRS-T波参考水平。
通过判断输入信号的PR段时长和PR段采样数是否满足一元线性回归条件进行筛选，过滤了不满足一元线性回归条件的输入信号，经过这轮决策能减轻运算负担，提高定位性能。
参照图3，在本发明的一些实施例中，将输入信号进行滤波，包括以下步骤：
步骤S203、对输入信号进行0.05Hz的高通滤波，得到高通滤波信号。
步骤S204、对高通滤波信号进行35Hz的低通滤波，得到带通滤波信号。
通过对输入信号进行带通滤波，高通滤波过滤基线漂移干扰，低通滤波过滤肌电干扰和工频噪声干扰，降低了基线漂移、肌电干扰、工频噪声的影响，使得后续的测算更加准确。
参照图4，在本发明的一些实施例中，根据一元线性回归方程判断带通滤波信号的PR段是否存在特殊波形，包括以下步骤：
步骤S401、计算一元线性回归方程的回归平方和以及总平分和。
步骤S402、根据回归平方和以及总平分和计算判定系数。
步骤S403、判定系数若小于判定系数阈值，则存在特殊波形。
通过一元线性回归方程的回归平方和以及总平方和计算判定系数，通过判定系数与判定系数阈值对比判定是否存在特殊波形，能够将存在特殊波形的输入信号过滤，防止出现误判的情况。
参照图5，在本发明的一些实施例中，判断带通滤波信号的PR段是否存在倾斜，包括以下步骤：
步骤S404、根据一元回归线性方程获取带通滤波信号的PR段的拟合斜率。
步骤S405、判断拟合斜率是否属于正常斜率范围，若拟合斜率属于正常斜率范围，则带通滤波信号的PR段不存在倾斜，若拟合斜率不属于正常斜率范围，则带通滤波信号的PR段存在倾斜。
通过对拟合斜率判断带通滤波信号的PR段存不存在倾斜，过滤存在倾斜的带通滤波信号，保证最后的输入信号是能够进行定位QRS-T波参考水平的，防止出行误判的情况。
在本发明的一些实施例中，若带通滤波信号的PR段存在特殊波形或存在倾斜，定位QRS-T波参考水平的方法还包括以下步骤：
步骤406、令QRS-T起始电势为QRS-T波参考水平。
若判定带通滤波信号的PR段存在特殊波形或存在倾斜，将不会进行后续操作，将QRS-T起始电势作为QRS-T波参考水平，保证了后续进行定位QRS-T波参考水平不会出现误判的情况。
参照图6，在本发明的一些实施例中，步骤407、根据一元线性回归方程定位QRS-T波参考水平，包括以下步骤：
步骤4071、计算一元线性回归方程在带通滤波信号的PR段上的残差数组。
步骤4072、将残差数组取绝对值，并增序排序得到残差绝对值数组。
步骤4073、对于残差绝对值数组的前
Figure 325890DEST_PATH_IMAGE018
个元素，获取并累加元素对应采样点的电势
Figure 925892DEST_PATH_IMAGE020
得到
Figure 131745DEST_PATH_IMAGE022
，
Figure 248606DEST_PATH_IMAGE024
表示判定系数，
Figure 147292DEST_PATH_IMAGE026
表示最少采样数阈值。
步骤4074、根据
Figure 733125DEST_PATH_IMAGE028
和
Figure 477090DEST_PATH_IMAGE030
，计算拟合平均电势
Figure 714036DEST_PATH_IMAGE032
。
步骤4075、令拟合平均电势
Figure 829629DEST_PATH_IMAGE034
为QRS-T波参考水平。
通过计算带通滤波信号的PR段上的残差数组并增序排序，选取前
Figure 152026DEST_PATH_IMAGE018
个元素，将在回归分析中置信度比较低的数据进行筛选，将置信度较低的数据不参与回归线拟合，使得拟合平均电势更加准确，定位QRS-T波参考水平更加精准。
参照图9，为了便于本领域人员理解，本发明的具体实施例提供了一种定位QRS-T波参考水平的方法，包括以下步骤：
第一步、根据PR段时长与采样数，判断输入信号的PR段是否满足一元线性回归条件。若PR段时长
Figure DEST_PATH_IMAGE036
小于最短时长阈值
Figure DEST_PATH_IMAGE038
或PR段采样数
Figure DEST_PATH_IMAGE040
少于最少采样数阈值
Figure DEST_PATH_IMAGE042
，则判断为在心电图中PR段不明显、应以QRS起始电势为QRS-T波参考水平，并跳转至第七步；否则，继续执行第二步。
第二步、将输入信号进行带通滤波，降低基线漂移、肌电干扰、工频噪声的影响，提高测算的准确性。
获取输入信号的滤波参数，包括高通滤波频率、低通滤波频率、陷波滤波频率；
对输入信号进行0.05Hz的高通滤波，过滤基线漂移干扰，得到高通滤波信号；
对得到的高通滤波信号进行35Hz的低通滤波，过滤肌电干扰和工频噪声干扰，得到带通滤波信号。
第三步、对带通滤波信号的PR段进行一元线性回归，得到滤波后PR段信号的一元线性回归方程。
第四步、评估一元线性回归方程。
计算一元线性回归方程的回归平方和SSR、总平方和SST；
根据回归平方和SSR、总平方和SST，计算判定系数
Figure DEST_PATH_IMAGE044
；
若
Figure 409002DEST_PATH_IMAGE044
小于判定系数阈值
Figure DEST_PATH_IMAGE046
，则判断为拟合效果未达到最低标准，可能存在特殊波形，应以QRS-T起始电势为QRS-T波参考水平，并跳转至第七步；否则，继续执行第五步。
第五步、判断PR段是否存在上斜或下斜。
通过得到的一元线性回归方程，获取PR段拟合斜率
Figure DEST_PATH_IMAGE048
；
若
Figure 438137DEST_PATH_IMAGE048
不属于PR段正常斜率范围
Figure DEST_PATH_IMAGE050
，则判断为PR段存在上斜或下斜，不能作为QRS起始部，应以QRS起始电势为QRS-T波参考水平，并跳转至第七步；否则，继续执行第六步。
第六步、根据一元线性回归方程，统计并计算QRS-T波参考水平。
计算一元线性回归方程在滤波后PR段上的残差数组
Figure DEST_PATH_IMAGE052
；
将残差数组
Figure 881888DEST_PATH_IMAGE052
取绝对值，并增序排序；
遍历排序后的残差绝对值数组。对于前
Figure 35789DEST_PATH_IMAGE018
个元素，获取并累加其对应采样点的滤波后电势
Figure 26617DEST_PATH_IMAGE054
，累加结果存于
Figure 847942DEST_PATH_IMAGE056
；
根据
Figure 383966DEST_PATH_IMAGE018
和
Figure 290742DEST_PATH_IMAGE056
，计算拟合平均电势
Figure 258829DEST_PATH_IMAGE058
；
令拟合平均电势
Figure 324874DEST_PATH_IMAGE058
为QRS-T波参考水平
Figure 907165DEST_PATH_IMAGE060
，并跳转至第八步。
第七步、获取QRS起点电势
Figure 943647DEST_PATH_IMAGE062
，令
Figure 371218DEST_PATH_IMAGE062
为QRS-T波参考水平
Figure 26190DEST_PATH_IMAGE060
，并跳转至第八步。
第八步、返回QRS-T波参考水平
Figure 779382DEST_PATH_IMAGE060
。
根据本发明的实施例，至少具有如下技术效果：
本发明通过输入信号的PR段时长和PR段采样数，判断PR段是否满足一元线性回归条件，判断能否使用PR段进行一元线性回归来计算平均电势，通过带通滤波，排除心电信号的基线漂移、肌电干扰、工频噪声等影响，使得进行一元线性回归时数据更加精确并增强了PR段拟合效果，使用一元线性回归能够拟合PR段直线，作为决策依据判断是否存在特殊波形或者倾斜，判断拟合PR段能否用于平均电势的统计，满足条件时将一元线性回归方程用于平均电势的统计，提高了QRS-T波参考水平的定位精度，整体上看进行多轮决策，针对不同的PR段形态，使用相应的QRS-T波参考水平的计算方式，使QRS-T波参考水平的定位精度更好。
参照图7，本发明的一个实施例，提供了一种定位QRS-T波参考水平的系统1000，包括信号获取模块1001、带通滤波模块1002、一元线性回归模块1003、参考水平判断模块1004，其中：
信号获取模块1001，获取输入信号；
带通滤波模块1002，根据输入信号的PR段时长和PR段采样数判断输入信号的PR段是否满足一元线性回归条件，若输入信号的PR段满足一元线性回归条件，将输入信号进行滤波，得到带通滤波信号；
一元线性回归模块1003，将带通滤波信号的PR段进行一元线性回归，得到一元线性回归方程；
参考水平判断模块1004，根据一元线性回归方程判断带通滤波信号的PR段是否存在特殊波形；若带通滤波信号的PR段不存在特殊波形，判断带通滤波信号的PR段是否倾斜，若带通滤波信号的PR段不存在倾斜，根据一元线性回归方程定位QRS-T波参考水平。
需要说明的是，由于本实施例中的一种定位QRS-T波参考水平的系统与上述的一种定位QRS-T波参考水平的方法基于相同的发明构思，因此，方法实施例中的相应内容同样适用于本装置实施例，此处不再详述。
参照图8，本发明的另一个实施例，还提供了一种电子设备，该电子设备6000可以是任意类型的智能终端，例如手机、平板电脑、个人计算机等。
具体的，电子设备6000包括：一个或多个控制处理器6001和存储器6002，图8中以一个控制处理器6001与一个存储器6002为例，控制处理器6001和存储器6002可以通过总线或者其他方式连接，图8中以通过总线连接为例。
存储器6002作为一种非暂态计算机可读存储介质，可用于存储非暂态软件程序、非暂态性计算机可执行程序以及模块，如本发明实施例中的一种电子设备对应的程序指令/模块；
控制处理器6001通过运行存储在存储器6002中的非暂态软件程序、指令以及模块，从而执行一种定位QRS-T波参考水平的方法的各种功能应用以及数据处理，即实现上述方法实施例的一种定位QRS-T波参考水平的方法。
存储器6002可以包括存储程序区和存储数据区，其中，存储程序区可存储操作系统、至少一个功能所需要的应用程序；存储数据区可存储定位QRS-T波参考水平的方法的使用所创建的数据等。此外，存储器6002可以包括高速随机存取存储器，还可以包括非暂态存储器，例如至少一个磁盘存储器件、闪存器件、或其他非暂态固态存储器件。在一些实施方式中，存储器6002可选包括相对于控制处理器6001远程设置的存储器，这些远程存储器可以通过网络连接至该电子设备6000。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。
在一个或者多个模块存储在存储器6002中，当被该一个或者多个控制处理器6001执行时，执行上述方法实施例中的定位QRS-T波参考水平的方法，例如执行以上描述的图1至图6的方法步骤。
存储器作为一种非暂态计算机可读存储介质，可用于存储非暂态软件程序以及非暂态性计算机可执行程序。此外，存储器可以包括高速随机存取存储器，还可以包括非暂态存储器，例如至少一个磁盘存储器件、闪存器件、或其他非暂态固态存储器件。在一些实施方式中，存储器可选包括相对于处理器远程设置的存储器，这些远程存储器可以通过网络连接至该处理器。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。
需要说明的是，由于本实施例中的一种电子设备与上述的一种定位QRS-T波参考水平的方法基于相同的发明构思，因此，方法实施例中的相应内容同样适用于本装置实施例，此处不再详述。
本发明的一个实施例，还提供一种计算机可读存储介质，存储有计算机可执行指令，计算机可执行指令用于执行：如上述实施例的定位QRS-T波参考水平的方法。
需要说明的是，由于本实施例中的一种计算机可读存储介质与上述的一种定位QRS-T波参考水平的方法基于相同的发明构思，因此，方法实施例中的相应内容同样适用于本装置实施例，此处不再详述。
本领域普通技术人员可以理解，上文中所公开方法中的全部或某些步骤、系统可以被实施为软件、固件、硬件及其适当的组合。某些物理组件或所有物理组件可以被实施为由处理器，如中央处理器、数字信号处理器或微处理器执行的软件，或者被实施为硬件，或者被实施为集成电路，如专用集成电路。这样的软件可以分布在计算机可读介质上，计算机可读介质可以包括计算机存储介质（或非暂时性介质）和通信介质（或暂时性介质）。如本领域普通技术人员公知的，术语计算机存储介质包括在用于存储数据（诸如计算机可读指令、数据结构、程序模块或其他数据）的任何方法或技术中实施的易失性和非易失性、可移除和不可移除介质。计算机存储介质包括但不限于RAM、ROM、EEPROM、闪存或其他存储器技术、CD-ROM、数字多功能盘（DVD）或其他光盘存储、磁盒、磁带、磁盘存储或其他磁存储装置、或者可以用于存储期望的数据并且可以被计算机访问的任何其他的介质。此外，本领域普通技术人员公知的是，通信介质通常包含计算机可读指令、数据结构、程序模块或者诸如载波或其他传输机制之类的调制数据信号中的其他数据，并且可包括任何数据递送介质。
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示意性实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
尽管已经示出和描述了本发明的实施例，本领域的普通技术人员可以理解：在不脱离本发明的原理和宗旨的情况下可以对这些实施例进行多种变化、修改、替换和变型，本发明的范围由权利要求及其等同物限定。

Claims (10)

1.一种定位QRS-T波参考水平的方法，其特征在于，包括以下步骤：
获取输入信号；
根据所述输入信号的PR段时长和PR段采样数判断所述输入信号的PR段是否满足一元线性回归条件，若所述输入信号的PR段满足一元线性回归条件，将所述输入信号进行滤波，得到带通滤波信号；
将所述带通滤波信号的PR段进行一元线性回归，得到一元线性回归方程；
根据所述一元线性回归方程判断所述带通滤波信号的PR段是否存在特殊波形；若所述带通滤波信号的PR段不存在特殊波形，判断所述带通滤波信号的PR段是否倾斜，若所述带通滤波信号的PR段不存在倾斜，根据所述一元线性回归方程定位QRS-T波参考水平。
2.根据权利要求1所述的定位QRS-T波参考水平的方法，其特征在于，所述根据所述输入信号的PR段时长和PR段采样数判断所述输入信号的PR段是否满足一元线性回归条件，包括以下步骤：
判断所述PR段时长是否小于最短时长阈值或所述PR段采样数是否小于最少采样数阈值；
所述PR段时长不小于最短时长阈值且所述PR段采样数不小于最少采样数阈值，所述输入信号的PR段满足一元线性回归条件；所述PR段时长小于最短时长阈值或所述PR段采样数小于最少采样数阈值，所述输入信号的PR段不满足一元线性回归条件，令QRS-T起始电势为QRS-T波参考水平。
3.根据权利要求2所述的定位QRS-T波参考水平的方法，其特征在于，所述将所述输入信号进行滤波，包括以下步骤：
对所述输入信号进行0.05Hz的高通滤波，得到高通滤波信号；
对所述高通滤波信号进行35Hz的低通滤波，得到带通滤波信号。
4.根据权利要求3所述的定位QRS-T波参考水平的方法，其特征在于，所述根据所述一元线性回归方程判断所述带通滤波信号的PR段是否存在特殊波形，包括以下步骤：
计算所述一元线性回归方程的回归平方和以及总平分和；
根据所述回归平方和以及所述总平分和计算判定系数；
所述判定系数若小于判定系数阈值，则存在特殊波形。
5.根据权利要求4所述的定位QRS-T波参考水平的方法，其特征在于，所述判断所述带通滤波信号的PR段是否存在倾斜，包括以下步骤：
根据所述一元回归线性方程获取所述带通滤波信号的PR段的拟合斜率；
判断所述拟合斜率是否属于正常斜率范围，若所述拟合斜率属于正常斜率范围，则所述带通滤波信号的PR段不存在倾斜，若所述拟合斜率不属于正常斜率范围，则所述带通滤波信号的PR段存在倾斜。
6.根据权利要求5所述的定位QRS-T波参考水平的方法，其特征在于，若所述带通滤波信号的PR段存在特殊波形或存在倾斜，所述定位QRS-T波参考水平的方法还包括以下步骤：
令QRS-T起始电势为QRS-T波参考水平。
7.根据权利要求6所述的定位QRS-T波参考水平的方法，其特征在于，所述根据所述一元线性回归方程定位QRS-T波参考水平，包括以下步骤：
计算所述一元线性回归方程在所述带通滤波信号的PR段上的残差数组；
将所述残差数组取绝对值，并增序排序得到残差绝对值数组；
对于所述残差绝对值数组的前
Figure DEST_PATH_IMAGE002
个元素，获取并累加所述元素对应采样点的电势
Figure DEST_PATH_IMAGE004
得到
Figure DEST_PATH_IMAGE006
，所述
Figure DEST_PATH_IMAGE008
表示所述判定系数，所述
Figure DEST_PATH_IMAGE010
表示所述最少采样数阈值；
根据所述
Figure DEST_PATH_IMAGE012
和所述
Figure DEST_PATH_IMAGE014
，计算拟合平均电势
Figure DEST_PATH_IMAGE016
；
令所述拟合平均电势
Figure 599255DEST_PATH_IMAGE016
为QRS-T波参考水平。
8.一种定位QRS-T波参考水平的系统，其特征在于，包括：
信号获取模块，获取输入信号；
带通滤波模块，根据所述输入信号的PR段时长和PR段采样数判断所述输入信号的PR段是否满足一元线性回归条件，若所述输入信号的PR段满足一元线性回归条件，将所述输入信号进行滤波，得到带通滤波信号；
一元线性回归模块，将所述带通滤波信号的PR段进行一元线性回归，得到一元线性回归方程；
参考水平判断模块，根据所述一元线性回归方程判断所述带通滤波信号的PR段是否存在特殊波形；若所述带通滤波信号的PR段不存在特殊波形，判断所述带通滤波信号的PR段是否倾斜，若所述带通滤波信号的PR段不存在倾斜，根据所述一元线性回归方程定位QRS-T波参考水平。
9.一种电子设备，其特征在于：包括至少一个控制处理器和用于与所述至少一个控制处理器通信连接的存储器；所述存储器存储有可被所述至少一个控制处理器执行的指令，所述指令被所述至少一个控制处理器执行，以使所述至少一个控制处理器能够执行如权利要求1至7任一项所述的定位QRS-T波参考水平的方法。
10.一种计算机可读存储介质，其特征在于：所述计算机可读存储介质存储有计算机可执行指令，所述计算机可执行指令用于使计算机执行如权利要求1至7任一项所述的定位QRS-T波参考水平的方法。
CN202210797009.4A 2022-07-08 2022-07-08 定位qrs-t波参考水平的方法、系统、设备及介质 Active CN114869296B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210797009.4A CN114869296B (zh) 2022-07-08 2022-07-08 定位qrs-t波参考水平的方法、系统、设备及介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210797009.4A CN114869296B (zh) 2022-07-08 2022-07-08 定位qrs-t波参考水平的方法、系统、设备及介质

Publications (2)

Publication Number Publication Date
CN114869296A true CN114869296A (zh) 2022-08-09
CN114869296B CN114869296B (zh) 2022-10-04

Family

ID=82683115

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210797009.4A Active CN114869296B (zh) 2022-07-08 2022-07-08 定位qrs-t波参考水平的方法、系统、设备及介质

Country Status (1)

Country Link
CN (1) CN114869296B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090234239A1 (en) * 2005-05-17 2009-09-17 Benjamin Shani Analysis of electrocardiogram signals
US20110040200A1 (en) * 2009-08-12 2011-02-17 Antony Louis Piriyakumar Douglas System and method for analyzing an electrocardiogram signal
CN102793539A (zh) * 2012-08-31 2012-11-28 深圳市理邦精密仪器股份有限公司 心电向量图检测分析方法及系统
CN108618777A (zh) * 2017-03-20 2018-10-09 深圳市理邦精密仪器股份有限公司 心电信号中st段信号的处理方法及装置、心电监测设备
CN109199366A (zh) * 2018-10-12 2019-01-15 山东理工大学 基于变分模态分解的心电信号qrs波检测方法
JP2020162894A (ja) * 2019-03-29 2020-10-08 国立大学法人信州大学 心電図解析装置，方法およびプログラム

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090234239A1 (en) * 2005-05-17 2009-09-17 Benjamin Shani Analysis of electrocardiogram signals
US20110040200A1 (en) * 2009-08-12 2011-02-17 Antony Louis Piriyakumar Douglas System and method for analyzing an electrocardiogram signal
CN102793539A (zh) * 2012-08-31 2012-11-28 深圳市理邦精密仪器股份有限公司 心电向量图检测分析方法及系统
CN108618777A (zh) * 2017-03-20 2018-10-09 深圳市理邦精密仪器股份有限公司 心电信号中st段信号的处理方法及装置、心电监测设备
CN109199366A (zh) * 2018-10-12 2019-01-15 山东理工大学 基于变分模态分解的心电信号qrs波检测方法
JP2020162894A (ja) * 2019-03-29 2020-10-08 国立大学法人信州大学 心電図解析装置，方法およびプログラム

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
陈清启: "ST段改变与测量方法", 《临床心电学杂志》 *

Also Published As

Publication number Publication date
CN114869296B (zh) 2022-10-04

Similar Documents

Publication Publication Date Title
CN207015288U (zh) 一种操作车辆的系统
US9042973B2 (en) Apparatus and method for measuring physiological signal quality
EP2945537B1 (en) Detection of sleep apnea using respiratory signals
JP2008520384A (ja) 限られたプロセッサ資源で呼吸数を実時間判定するための方法及びシステム
US20130331723A1 (en) Respiration monitoring method and system
CN109009084B (zh) 多导联心电信号的qrs波群校验方法、装置及设备、介质
CN112869752B (zh) 一种心电信号采集装置及质量等级评估和qrs波检测方法
CN104462744A (zh) 适用于心血管远程监测系统的数据质量控制方法
CN111370124A (zh) 一种基于面手部识别和大数据的健康分析系统及方法
US8712518B2 (en) Determining phase-specific parameters of a physiological variable
CN117457173B (zh) 一种用于心内科的可穿戴设备心律失常监测系统
CN114869296B (zh) 定位qrs-t波参考水平的方法、系统、设备及介质
KR101029386B1 (ko) 실시간 심전도 모니터링 시스템을 위한 큐알에스 검출 방법
CN110840401A (zh) 一种心率值测定方法、装置、设备及可读存储介质
John et al. An evaluation of ECG data fusion algorithms for wearable IoT sensors
CN117238510A (zh) 一种基于深度学习的脓毒症预测方法及系统
WO2019001040A1 (zh) 用于亚裔人群的肺癌检测方法及系统
Domazet et al. Improving the QRS detection for one-channel ECG sensor
KR101319624B1 (ko) 심전도 신호로부터 알-파 검출 방법 및 이를 이용한 심박변이도 측정 방법
CN114159075B (zh) Qrs波优化装置、系统及存储介质
CN112741623B (zh) 一种基于心电的糖化血红蛋白检测系统和方法
JP2010213809A (ja) 生体信号分析装置
CN111345814B (zh) 一种心电信号中心拍的分析方法、装置、设备和存储介质
CN113520368A (zh) 一种咳嗽监测的方法、系统和存储设备
Kim et al. Detection of R-peaks in ECG signal by adaptive linear neuron (ADALINE) artificial Neural Network

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