申请/专利权人：江苏华康信息技术有限公司;南京贺普检测仪器有限公司

申请日：2018-03-23

公开（公告）日：2021-01-12

公开（公告）号：CN108814590B

主分类号：A61B5/366(20210101)

分类号：A61B5/366(20210101);A61B5/349(20210101);A61B5/24(20210101)

优先权：

专利状态码：有效-授权

法律状态：2021.01.12#授权;2018.12.11#实质审查的生效;2018.11.16#公开

摘要：本发明公开了一种心电QRS波群的检测方法，包括使用极值法寻找QRS波中的波峰和波谷，再通过找到的极值点，结合等电位段电位值、幅值信息来检出QRS波群中的各个子波形，确定QRS波群的时间位置，从而可清楚地显示QRS波群形态及计算出心率；QRS波群的精确检测是心电图自动诊断的基础。因为只有在确定QRS波群的检测后才有可能计算心率即每分钟心跳的次数、心率变异性及其他心电图各波段时间间隔测量和幅度测量。

主权项：1.一种心电QRS波群的检测方法，其特征在于，包括如下步骤：A、获取采样率为fs的原始心电波形幅值数据S，并将心电波形S分成以N秒为时间单位的处理段信号X，N为8～120间的任意一值；B、将信号X经过包络处理得到包络信号Xe，然后在包络信号Xe上利用阈值法定位QRS波群的位置得到集合Bxe；在包络信号Xe上利用面积法找出每个QRS波群的起始点QRSon和终末点QRSoff的预估位置；C、在集合Bxe中选取第n个QRS波群位置Bxen，求出该QRS波群的RR间期RRxen＝Bxen-Bxen-1；从信号X中截取一片段Xios，截取的起始位置为Bxen-1+RRxen×k3，终止位置为Bxen-RRxen×0.5-k3，其中0.3＜k3＜0.5，即第n-1个心跳的T波之后到第n个心跳的P波之前的一段数据；作出片段Xios中幅值范围的直方图Gp，得到Xios数据取值的统计分布；在直方图Gp中找出Xios幅值分布个数最多的子区间b，该子区间b代表的电位即为第n次心跳之前的等电位段电位vios；D、从信号X中第n个QRS波群位置Bxen周围前后截取一数据段qn，在数据段qn中找到所有极大值和极小值点的集合为极值点P和QRS波群中绝对幅值最大的波QM；E、极值点P中满足最小QRS可识别波条件的为主极点Dn，再从主极点Dn中筛选出特征点和关键特征点其中最小QRS可识别波条件为：幅度大于ρmin微伏且持续时间大于dmin毫秒，其中20uV＜ρmin＜80uV，6ms＜dmin＜16ms；F、将特征点和关键特征点中所有位置小于QRSon和大于QRSoff的特征点删除，形成新的代表QRS波群的特征点Rn和关键特征点G、采用新的Rn或者并结合步骤C中的等电位段的电势值vios，确定QRS波群中的子波形Q、R、S、r＇、s＇、r＂、s＂位置；H、重复C到H步，直至检出集合Bxe中所有QRS的特征点，用Bn＝{Bnj}表示；I、重复A到I步，直到识别完全部原始心电波形幅值数据S。

全文数据：一种心电QRS波群的检测方法及其心电分析方法技术领域[0001]本发明涉及一种心电信号中QRS波群的识别，尤其涉及一种心电QRS波群的检测方法及其心电分析方法。背景技术[0002]在心电检测领域测量心电图时，心电信号质量的好坏是心电有效分析的前提，信号质量差的心电信号会对心电分析起到误导作用，在单导联便携式心电采集中尤为突出，因为单导联便携式心电采集存在突然的波形抖动，运动伪波以及容易引入较大工频干扰等噪声。[0003]为了实现心电的正确的分析诊断，心电波形的质量必须要好，过多的噪声以及大幅的抖动会导致心电自动分析判断的准确度变差，出现误检漏检的情况。为了兼顾单导联便携式心电采集能快速测量的优点和心电自动分析的功能，有必要对心电波形的质量做实时分析，对波形质量差的地方不予以自动分析判断，避免出现误检测的情况。[0004]同时心电图的自动分析与诊断对于心脏功能检测、心血管疾病的诊断和预防有着重要的意义，尤其是QRS波群的精确检测，尤其是宽度、幅度和形态，是心电图自动诊断的基础。因为只有在确定QRS波群的检测后才有可能计算心率、心率变异性及其他心电图各波段时间间隔测量和幅度测量。如QRS波群形态正常，宽度在正常范围（60-100毫秒）内，则表示心室的激动是由房室交界部或更高的部位传下的，称为室上性QRS波群，属于正常现象；否则就表示心室是由房室交界部以下的异位节律点所激动，称为室性QRS波群，属于异常现象，或是存在有室内传导障碍等。发明内容[0005]发明目的：本发明的第一目的是提供一种可精准找出QRS波群位置并准确勾勒出QRS波群形体特征的心电QRS波群的检测方法。[0006]本发明的第二目的是提供一种基于该心电QRS波群的检测方法的心电分析方法。[0007]提供基于心电QRS波群检测方法的一种心电分析方法。[0008]技术方案:为实现以上目的，本发明公开了一种心电QRS波群的检测方法，包括如下步骤：[0009]A、获取采样率为fs的原始心电波形幅值数据S，并将心电波形S分成以~秒为时间单位的处理段信号X，N为8〜120间的任意一值；[0010]B、将信号X经过包络处理得到包络信号Xe，然后在包络信号Xe上利用阈值法定位QRS波群的位置得到集合Bxe;在包络信号Xe上利用面积法找出每个QRS波群的起始点QRSoJP终末点QRSciff的预估位置；[0011]C、在集合Bxe中选取第η个QRS波群位置Bxeη，求出该QRS波群的RR间期RRxeη=Bxen-Bxeη-1;从信号X中截取一片段Xicis，截取的起始位置为Bxen-1+RRxenXk3，终止位置为^11-1«^11\0.5;3，其中0.3〈1«〈0.5，即第11-1个心跳的1'波之后到第11个心跳的P波之前的一段数据;作出片段Xlcis中幅值范围的直方图Gp，得到Xlcis数据取值的统计分布;在直方图Gp中找出Xitjs幅值分布个数最多的子区间b，该子区间b代表的电位即为第η次心跳之前的等电位段电位Vicis;[0012]D、从信号X中第η个QRS波群位置Bxeη周围截取一数据段qn，在数据段找到所有极大值和极小值点的集合为极值点P和QRS波群中绝对幅值最大的波Qm;[0013]E、极值点P中满足最小QRS可识别波条件的为主极点Dn，再从主极点Dn中筛选出特征点氧1和关键特征点IfsV[0014]F、将特征点|„和关键特征点Afs中所有位置小于QRSon和大于QRSofT的特征点删除，形成新的代表QRS波群的特征点Rn和关键特征点Kps;[0015]G、采用新的1或者并结合步骤C中的等电位段的电势值Vicis，确定QRS波群中的子波形Q、R、S、r\s\r"、s"位置；[001ό]H、重复C到H步，直至检出集合Bxe中所有QRS的特征点，用Bn={Bnj}表示；[0017]I、重复A到I步，直到识别完全部原始心电波形幅值数据S。[0018]其中步骤B中信号X的处理方法、利用阀值法和面积法的定位方法包括如下步骤：[0019]Bl、将信号X经过一个带通滤波器，其频率范围为5Hz〜25Hz，得到凸显出QRS波群特征且滤除干扰波的信号Xfiltered;对信号Xfiltered做微分后再平方得到信号Xdiff，对Xdiff进行积分得到包络信号Xe;[0020]B2、将包络信号Xe3按0.2N秒的长度分段得到5个包络片段Xei，其中i=1一5，在每个包络片段内找到最大值XEmaxi，然后求出这5个最大值的中值，并用中值计算出阈值其中0.2〈kK0.8;再在包络信号Xe3中寻找一个大于阈值TH的区间[Icive^I1Ciwer]，该区间即为QRS波群的位置所在区域;然后确定QRS波群的坐标位置Bxeη=argmaxXeI_r，Ilower，η代表第η个QRS波，其中arg为自变量argument的英文缩写，argmax即代表取式子达到最大值时的变量取值，Bxe3η是I_r〜Ii_2间的使Xe^大的一个值;依次利用阀值在Xe上找出所有QRS波群的位置，记为集合Bxe;[0021]B3、在包络信号Xe上，截取第η个QRS波包络信号：[0022]Β4、利用公式中第t个点的前向面是积分索引，0.08〈1〈0.15为求面积时的时间跨度;得到夂=^中所有数据的前向面积i，求出使前向面积蓋最大的t值，即为第η个QRS波的预估终末点QRSoff;[0023]Β5、利用公式求出;iefS中第t个点的后向面积五〇，1是积分索引，0.08〈1〈0.15为求面积时的时间跨度;得到14@中所有数据的后向面积毛求出使后向面积其最大的t值，即为第η个QRS波的预估起始点QRSon;[0024]Β6、重复M〜Β6步，直到确定所有由Bxe定位的QRS波群的预估起止点。[0025]优选的，所述步骤D中获取极值点P和Qm的方法包括如下步骤：[0026]Dl、信号X在Bxeη处往前ω_=〇·07〜0·15fs、往后ω+=〇·1〜〇·23fs个点，截取信号X中的QRS波群的数据段qn=qi，…，q」，…，qu，其中ω=ω_+ω+表示qn的数据长度，qj表示qn中第j个数据；[0027]D2、在数据段qn中找到所有极大值Pk和极小值点Lo,极大值Pk和极小值点Lo合称为极值点P;然后再找出极大值Pk中的最大值Pmax=Vmax，Imax和极小值Lo中的最小值Lmin=Vmin，Imin，I、V分别表不点在qn中的位置和数据幅值大小；[0028]D3、定义QRS波群中最明显的波顶点为QRS波群中绝对幅值最大的波Qm，如果以_-ViosIk4XIVmin-VicisI则Qm点为Lmin点，否则QM点为PmaxA，其中2〈k4〈10〇[0029]再者，所述步骤E中最小QRS可识别波条件为：幅度大于Pmin微伏且持续时间大于dmin_秒，其中20uV〈Pmin〈80uV，6ms〈dmin〈16ms。[0030]优选的，所述步骤E中筛选特征点和关键特征点Pfs的具体方法包括如下步骤：[0031]E1、从极值点P中选择一个极值点W，用以下公式在处理段qn中搜索极值点W的左右支撑区间，矿：[0034]其中j和qj为极值点Pj在qn中的位置和幅值，τ为QRS最大的有意义的生理时间宽度，取80〜160ms;i、a、b和k均是为了求解支撑区间的辅助变量，没有特别含义;qx是qn中第X个数据，与W—个意义；[0035]E2、若极值点Pj的左右支撑区间f，矿不存在，或者矿-f则确认特征点ρ」为QRS中的特征点，如果不满足则确认特征点Pj不是QRS中的关键特征点，其中3〇°〈€〈0〈65°;对〇11中的所有极点进行判断，最后得到QRS中的特征点A.=[0038]E5、如及n中的极值点幅值满足条件：Ρρ.ί?5·}，Ριιάη〈ρ〇[«〈150μν，其中表不两点间的绝对幅度差，是步骤Ε4中找到的关键点，Pmin是最小可识别波幅值阈值，PQRS是为了抗干扰而设置的比大的幅度阈值，即为关键特征点及fS。[0039]进一步，所述步骤G中确定QRS波群中的子波形位置的方法包括如下步骤：[0040]G1、当01α1，再与每分割段的包络差ei作比较;当Si匕Xδ时，说明此分隔段的包络有突变，即为包络检测不合格；[0069]S204、在处理段X内，获取每一分割段Xi的幅值方差信，其中，m为^中数据的序号;再求取处理段X内波形片段幅值方差值平均值乘以比例系数k2,其中2k2l，再与各分割段X1的幅值方差值E1作比较，当：吋，说明此分割段的方差有突变，即为方差检测不合格；[0070]S205、在处理段X内，利用快速傅里叶变换，将心电片段信号Xi的时域信号转化为频域信号，即为功率谱信号；对1〜5Hz的幅值进行积分得到功率其中Re代表求复数的实部，Im代表求复数的虚部，j指某一频段的范围，fij即为功率谱信号;求取处理段X内波形的平均功率值P，乘以比例系数k3，其中2k31，再与各分割段Xi的功率Pi作比较；当d,X7时，说明此处理段的1〜5Hz功率有突变，即为功率谱检测不合格；[0071]S206、在处理段X内，得到步骤S5中的心电片段信号Xi的功率谱信号后，对5〜40Hz的幅值进行积分得到功率，其中Re代表求复数的实部，Im代表求复数的虚部，Psi代表心电信号的功率;再对40〜IOOHz的幅值进行积分得到功率其中Re代表求复数的实部，Im代表求复数的虚部，Pm代表高频噪声信号的功率;计算两者的信噪比kSNRl=PSlPm，当信噪比小于阈值Gsrn时，其中,大于2,说明此处理段信噪比过小，噪声所占比重过大，即信噪比检测不合格；[0072]S207、根据参数包络差、方差值、信号在1〜5Hz之间的功率和信噪比的合格情况对心电波形进行质量等级划分，其中心电波形的质量等级划分方法为：当信噪比参数不合格时，该处理段心电波形的质量等级系数直接评定为3,即波形存在严重噪声;在信噪比参数合格的前提下，根据包络差、方差值、信号在1〜5Hz之间的功率的合格情况来评定心电波形质量，当三个参数全部合格时，该处理段心电波形质量等级系数评定为0,即波形良好；当三个参数中有两个或一个合格时，心电波形质量等级系数评定为1，即波形较差；当三个参数全部不合格时，心电波形质量等级系数评定为2,即波形差；心电波形的质量分成4个等级，分别用〇〜3表示质量等级系数;质量等级系数0表示心电波形良好，质量等级系数1表示心电波形较差，质量等级系数2表示心电波形差，质量等级系数3表示心电波形存在严重噪声。即从波形质量的程度“〇”表示良好，“Γ表示较差，“2”表示差，“3”表示严重噪声。[0073]有益效果:与现有技术相比，本发明具有以下显著优点：[0074]1本发明通过步骤B中对信号进行包络处理，凸显QRS波群位置信息，削弱杂波对其影响，有效提高QRS波群的精准位置的准确度；[0075]2本发明通过步骤C中统计法得到等电位段电位，再利用幅度、斜率和极值点信息识别出QRS波群中各个关键点，各关键点、等电位段电位和预估的QRS波群的起始点和终末点，最终确定QRS波群中各子波的位置，并准确勾勒出QRS波群的形态特征；同时本发明可求出各个子波的幅度与宽度，进而可以为心电QRS波群形态特征分析提供详细的特征参数；[0076]3本发明利用心电波形幅值数据计算心电波形的包络差、方差值、信号在1〜5Hz之间的功率和信噪比四个参数，分析四个参数值是否合格，根据参数包络差、方差值、信号在1〜5Hz之间的功率和信噪比的合格情况对心电波形进行质量等级划分，为心电信号自动分析筛选合格的心电数据，确保待检测心电信号的可用性;从而进一步保障后续精准检测到QRS波群的位置；[0077]4本发明利用心电信号质量辨识方法对待分析的心电波形进行质量评估，对心电波形质量等级好的心电特征进行心电参数的计算和异常结果判断，对于心电波形质量等级差的心电信号发出提示，提高心电信号分析的可信度以及准确度;本发明最后结合波形质量分析法对心电波质量的甄别，有效地分析心电波形。附图说明[0078]图1为本发明的流程框图；[0079]图2为本发明步骤B中信号X的处理方法的流程框图；[0080]图3为本发明步骤B中信号X处理后的示意图；[0081]图4为本发明步骤M中面积法的示意图；[0082]图5为本发明步骤B中QRS波群的预估位置；[0083]图6为本发明步骤C中的直方图Gp的示意图；[0084]图7为本发明步骤E中极值点的左右区间示意图；[0085]图8为本发明步骤G5中QRS波的形态示意图；[0086]图9为本发明中QRS命名规则示意图一；[0087]图10为本发明中QRS命名规则示意图二；[0088]图11为本发明中QRS命名规则示意图三；[0089]图12为本发明中QRS波的具体形态示意图；[0090]图13为本发明的心电分析方法的流程框图；[0091]图14为本发明中心电信号质量辨识方法的流程框图；[0092]图15为本发明中心电波形中各波的示意图；[0093]图16为本发明实施例中的实验结果示意图一；[0094]图17为本发明实施例中的实验结果示意图二。具体实施方式[0095]下面结合附图对本发明的技术方案作进一步说明。[0096]QRS波群的精确检测是心电图自动诊断的基础，因为只有在确定QRS波群的检测后才有可能计算心率心率为一分钟内出现的QRS波群的个数，即每分钟心跳的次数）、心率变异性及其他心电图各波段时间间隔测量和幅度测量。[0097]如图1所示，本发明一种心电QRS波群的检测方法，包括如下步骤：[0098]A、获取采样率为fs的原始心电波形幅值数据S，并将心电波形S分成以~秒为时间单位的处理段信号X，N为8〜120间的任意一值；[0099]B、将信号X经过包络处理得到包络信号Xe，然后在包络信号Xe上利用阈值法定位QRS波群的位置得到集合Bxe;在包络信号Xe上利用面积法找出每个QRS波群的起始点QRScir^P终末点QRSciff的预估位置；[0100]如图2所示，其中步骤B中信号X的处理方法、利用阀值法和面积法的定位方法包括如下步骤：[0101]BI、将信号X经过一个带通滤波器，其频率范围为5Hz〜25Hz，得到凸显出QRS波群特征且滤除干扰波的信号Xfiltered;对信号Xfiltered做微分后再平方得到信号Xdiff，对Xdiff进行积分得到包络信号Xe，如图3所示；[0102]B2、将包络信号Xe按0.2N秒的长度分段得到5个包络片段Xei，其中ί=1···5,在每个包络片段内找到最大值XEmaxi，然后求出这5个最大值的中值恐并用中值计算出阈值TH=Is1XHwukc，其中〇.2〈1^〈0.8;再在包络信号Xe3中寻找一个大于阈值TH的区间[Icive^I1Ciwer]，该区间即为QRS波群的位置所在区域;然后确定QRS波群的坐标位置Bxeη=argmaxXeI_r，Ilower，η代表第η个QRS波，其中arg为自变量argument的英文缩写，argmax即代表取式子达到最大值时的变量取值，Bxe3η是I_r〜Ii_2间的使Xe^大的一个值;依次利用阀值在Xe上找出所有QRS波群的位置，记为集合Bxe;[0103]B3、在包络信号Xe上，截取第η个QRS波包络信号Bxeη+0.2Xfs；[0104]B4、利用公式'中第t个点的前向面积是积分索引，0.08fs〈W〈0.15fs为求面积时的时间跨度;得到1；^^中所有数据的前向面积;S，求出使前向面积S最大的t值，即为第η个QRS波的预估终末点QRSoff，如图4所示；[0105]B5、利用公式中第t个点的后向面积Atj是积分索引，0.08〈1〈0.15为求面积时的时间跨度;得到心^中所有数据的后向面积I求出使后向面积S最大的t值，即为第η个QRS波的预估起始点QRSon;[0106]Β6、重复M〜Β6步，直到确定所有由Bxe定位的QRS波群的预估起止点，如图5所示；[0107]C、在集合Bxe中选取第η个QRS波群位置Bxeη，求出该QRS波群的RR间期RRxeη=Bxen-Bxeη-1;从信号X中截取一片段Xicis，截取的起始位置为Bxen-1+RRxenXk3，终止位置为^1〇-1«^11\0.5-1«，其中0.3〈1«〈0.5，即第11-1个心跳的1'波之后到第11个心跳的P波之前的一段数据;作出片段Xlcis中幅值范围的直方图Gp，得到Xlcis数据取值的统计分布，如图6所不;在直方图Gp中找出Xios幅值分布个数最多的子区间b，该子区间b代表的电位即为第η次心跳之前的等电位段电位Vlcis;[0108]D、从信号X中第η个QRS波群位置Bxeη周围截取一数据段qn，在数据段^中找到所有极大值和极小值点的集合为极值点P和QRS波群中绝对幅值最大的波Qm;其中获取极值点P和Qm的方法包括如下步骤：[0109]Dl、信号X在Bxeη处往前ω_=〇·07〜0·15fs、往后ω+=〇·1〜〇·23fs个点，截取信号X中的QRS波群的数据段qn=qi，…，q」，…，qu，其中ω=ω_+ω+表示qn的数据长度，qj表示qn中第j个数据；[0110]D2、在数据段qn中找到所有极大值Pk和极小值点Lo,极大值Pk和极小值点Lo合称为极值点P;然后再找出极大值Pk中的最大值Pmax=Vmax，Imax和极小值Lo中的最小值Lmin=Vmin，Imin，I、V分别表不点在qn中的位置和数据幅值大小；[0111]D3、定义QRS波群中最明显的波顶点为QRS波群中绝对幅值最大的波Qm，如果以_-ViosIk4XIVmin-VicisI则Qm点为Lmin点，否则QM点为PmaxA，其中2〈k4〈10;[0112]E、极值点P中满足最小QRS可识别波条件的为主极点Dn，再从主极点Dn中筛选出特征点仏和关键特征点最小QRS可识别波条件为：幅度大于Pmin微伏且持续时间大于CUin毫秒，其中2〇1^〈0*〈8〇1^，61118〈1_〈161118;筛选特征点馬1和关键特征点及^的具体方法包括如下步骤：[0113]E1、从极值点P中选择一个极值点W，用以下公式在处理段qn中搜索极值点w的左右支撑区间f，Zfi[0116]其中j和qj为极值点？」在^中的位置和幅值，τ为QRS最大的有意义的生理时间宽度，取80〜160ms;i、a、b和k均是为了求解支撑区间的辅助变量，没有特别含义;qx是qn中第X个数据，与^一个意义，如图7所示；[0117]E2、若极值点Pj的左右支撑区间f，矿不存在，或者矿-tana°则确认特征点pj为QRS中的特征点，如果不满足则确认特征点Pj不是QRS中的关键特征点，其中3〇°〈€〈0〈65°;对〇11中的所有极点进行判断，最后得到QRS中的特征点[0120]E5、如及n中的极值点幅值满足条件%Rs}，.Pmin〈PQRs〈150yV，其中表不两点间的绝对幅度差，是步骤E4中找到的关键点，ρ_是最小可识别波幅值阈值，pQRs是为了抗干扰而设置的比大的幅度阈值，即为关键特征点%[0121]F、将特征点和关键特征点中所有位置小于QRSon和大于QRSoff的特征点删除，形成新的代表QRS波群的特征点Rn和关键特征点;[0122]G、采用新的匕或者并结合步骤C中的等电位段的电势值Vlcis，确定QRS波群中的子波形Q、R、S、r\s\r"、s"位置；[0123]其中确定QRS波群中的子波形位置的方法包括如下步骤：[0124]G1、当0hl，再与每分割段的包络差ei作比较;当ei:C1Xδ时，说明此分隔段的包络有突变，即为包络检测不合格；[0147]S204、在处理段X内，获取每一分割段Xi的幅值方差值其中，!!!为^中数据的序号;再求取处理段X内波形片段幅值方差值平均值民，乘以比例系数k2,其中2k2l，再与各分割段X1的幅值方差值E1作比较，当：时，说明此分割段的方差有突变，即为方差检测不合格；[0148]S205、在处理段X内，利用快速傅里叶变换，将心电片段信号Xi的时域信号转化为频域信号，即为功率谱信号；对1〜5Hz的幅值进行积分得到功率，其中Re代表求复数的实部，Im代表求复数的虚部，j指某一频段的范围,fdj即为功率谱信号;求取处理段X内波形的平均功率值P:，乘以比例系数k3,其中2k3l，再与各分割段Xi的功率Pi作比较；当]时，说明此处理段的1〜5Hz功率有突变，即为功率谱检测不合格；[0149]S206、在处理段X内，得到步骤S5中的心电片段信号Xi的功率谱信号后，对5〜40Hz的幅值进行积分得到功率其中Re代表求复数的实部，Im代表求复数的虚部，Psi代表心电信号的功率;再对40〜IOOHz的幅值进行积分得到功率其中Re代表求复数的实部，Im代表求复数的虚部，Pm代表高频噪声信号的功率;计算两者的信噪比kSNRl=PSlPm，当信噪比小于阈值Rsrn时，其中:¾大于2,说明此处理段信噪比过小，噪声所占比重过大，即信噪比检测不合格；[0150]S207、根据参数包络差、方差值、信号在1〜5Hz之间的功率和信噪比的合格情况对心电波形进行质量等级划分，其中心电波形的质量等级划分方法为：当信噪比参数不合格时，该处理段心电波形的质量等级系数直接评定为3,即波形存在严重噪声;在信噪比参数合格的前提下，根据包络差、方差值、信号在1〜5Hz之间的功率的合格情况来评定心电波形质量，当三个参数全部合格时，该处理段心电波形质量等级系数评定为0,即波形良好；当三个参数中有两个或一个合格时，心电波形质量等级系数评定为1，即波形较差；当三个参数全部不合格时，心电波形质量等级系数评定为2,即波形差；心电波形的质量分成4个等级，分别用〇〜3表示质量等级系数;质量等级系数0表示心电波形良好，质量等级系数1表示心电波形较差，质量等级系数2表示心电波形差，质量等级系数3表示心电波形存在严重噪声。即从波形质量的程度“〇”表示良好，“Γ表示较差，“2”表示差，“3”表示严重噪声；[0151]S3、当心电信号质量等级系数为3且持续时间大于0.3N秒时，提示“信号质量噪声过大，无法分析”返回至步骤A重新获取待分析的心电波形；否则进行心电参数的计算和异常结果判断;其中心电参数的计算和异常结果判断的步骤为：[0152]S301、依据心电QRS波群的检测方法找出QRS波群中各个子波形的起止时间点、幅度参数和代表QRS波群形态的Bn;然后再找出心电特征包括P、T波的起止时间点及其幅度参数，构成了一次心跳；[0153]S302、取心电波形质量等级系数小于3的心电特征用于计算心电参数RR间期和心率，用心电波形质量等级为〇的心电特征计算心电参数PR间期、QRS宽度、QT间期、ST段高度、P波平均电轴、QRS平均电轴、T波平均电轴、Q波宽度和R波高度；[0154]S303、利用心电特征和上述心电参数进行心电异常判断得出心电图异常结果；[0155]S4、当心电信号质量等级系数为1或2且持续时间大于0.3N秒时，提示“信号质量较差”并输出计算得到的心电参数和心电图异常结果;否则认为心电波形质量较好，直接输出计算得到的心电参数和心电图异常结果。[0156]实施例[0157]如图15所示，本发明一种基于心电QRS波群的检测方法的心电分析方法，包括如下步骤：[0158]S1、获取10秒的待分析的心电波形；[0159]S2、依据心电信号质量辨识方法对待分析的心电波形进行质量评估，得到10秒的心电波形质量等级，其中心电波形质量评估以Is为单位，共得到10个数据分别代表第1秒到第10秒的质量等级系数；[0160]如图13所示，其中心电信号质量辨识方法包括如下步骤：[0161]S201、获取采样率为fs的原始的心电波形幅值数据，并将心电波形幅值数据分成以10秒为时间单位的处理段X;[0162]S202、在处理段X内，将心电信号以1秒为单位进行片段分割，分割段为X1，其中i=1…10，xi的数据长度为M=IXfs=fs;[0163]S203、以处理段X为单位，取出每一心电波形片段Xi的幅值最大值mmajP最小值mmin构成包络点，通过比较各个包络点得到分割段X1的包络差61=!1^^1^111，再求取处理段X内波形包络差的平均值：，乘以比例系数1^，其中2hl，再与每分割段的包络差ei作比较;当ei:C1X即寸，说明此分隔段的包络有突变，即为包络检测不合格；[0164]S204、在处理段X内，获取每一分割段Xi的幅值方差值其中m为Xl*数据的序号;再求取处理段X内波形片段幅值方差值平均值馬，乘以比例系数k2,其中2k2l，再与各分割段X1的幅值方差值E1作比较，当吋，说明此分割段的方差有突变，即为方差检测不合格；[0165]S205、在处理段X内，利用快速傅里叶变换，将心电片段信号Xi的时域信号转化为频域信号，即为功率谱信号；对1〜5Hz的幅值进行积分得到功率其中Re代表求复数的实部，Im代表求复数的虚部，j指某一频段的范围，fij即为功率谱信号;求取处理段X内波形的平均功率值，乘以比例系数k3，其中2k31，再与各分割段Xi的功率Pi作比较;当时，说明此处理段的1〜5Hz功率有突变，即为功率谱检测不合格；[0166]S206、在处理段X内，得到步骤S5中的心电片段信号Xi的功率谱信号后，对5〜40Hz的幅值进行积分得到功率，其中Re代表求复数的实部，Im代表求复数的虚部，Psi代表心电信号的功率;再对40〜IOOHz的幅值进行积分得到功率其中Re代表求复数的实部，Im代表求复数的虚部，Pm代表高频噪声信号的功率;计算两者的信噪比kSNRl=PSlPm，当信噪比小于阈值ftSRN时，其中:kSRN大于2,说明此处理段信噪比过小，噪声所占比重过大，即信噪比检测不合格；[0167]S207、根据参数包络差、方差值、信号在1〜5Hz之间的功率和信噪比的合格情况对心电波形进行质量等级划分，其中心电波形的质量等级划分方法为：当信噪比参数不合格时，该处理段心电波形的质量等级系数直接评定为3,即波形存在严重噪声;在信噪比参数合格的前提下，根据包络差、方差值、信号在1〜5Hz之间的功率的合格情况来评定心电波形质量，当三个参数全部合格时，该处理段心电波形质量等级系数评定为0,即波形良好；当三个参数中有两个或一个合格时，心电波形质量等级系数评定为1，即波形较差；当三个参数全部不合格时，心电波形质量等级系数评定为2,即波形差；心电波形的质量分成4个等级，分别用〇〜3表示质量等级系数;质量等级系数0表示心电波形良好，质量等级系数1表示心电波形较差，质量等级系数2表示心电波形差，质量等级系数3表示心电波形存在严重噪声。即从波形质量的程度“〇”表示良好，“Γ表示较差，“2”表示差，“3”表示严重噪声；[0168]S3、当心电信号质量等级系数为3且持续时间大于3秒时，提示“信号质量噪声过大，无法分析”返回至步骤A重新获取待分析的心电波形；否则进行心电参数的计算和异常结果判断;其中心电参数的计算和异常结果判断的步骤为：[0169]S301、依据心电QRS波群的检测方法找出QRS波群中各个子波形的起止时间点、幅度参数和代表QRS波群形态的Bn;然后再找出心电特征包括P、T波的起止时间点及其幅度参数，构成了一次心跳；[0170]S302、取心电波形质量等级系数小于3的心电特征用于计算心电参数RR间期和心率，用心电波形质量等级为〇的心电特征计算心电参数PR间期、QRS宽度、QT间期、ST段高度、P波平均电轴、QRS平均电轴、T波平均电轴、Q波宽度和R波高度；[0171]S303、利用心电特征和上述心电参数进行心电异常判断得出心电图异常结果；[0172]S4、当心电信号质量等级系数为1或2且持续时间大于3秒时，提示“信号质量较差”并输出计算得到的心电参数和心电图异常结果；否则认为心电波形质量较好，直接输出计算得到的心电参数和心电图异常结果。[0173]如图14所示，本发明中所提及的心电波形P波为心房激动波，也称为心房的除极波;QRS波群为心室激动的综合波群，也称为心室的除极波，其第一个向下的波为Q波、任何向上的波都叫R波、在R波以后的任何向下的波都叫S波;T波称为心室的复极波，而心房的复极波很小，且埋于心室的除极波中，不易辨认，未给于专门命名。[0174]如图16和图17所示，本发明以10为时间单位进行心电分析的实验数据图和结果，图中方波表示质量辩识结果的标识，标识方波中：〇表示波形质量好和良好;1、2表示波形质量较差和差;3表示波形质量严重差。

权利要求：I.一种心电QRS波群的检测方法，其特征在于，包括如下步骤：A、获取采样率为fs的原始心电波形幅值数据S，并将心电波形S分成以~秒为时间单位的处理段信号X，N为8〜120间的任意一值；B、将信号X经过包络处理得到包络信号Xe，然后在包络信号Xe上利用阈值法定位QRS波群的位置得到集合Bxe;在包络信号Xe上利用面积法找出每个QRS波群的起始点QRScin和终末点QRSciff的预估位置；C、在集合Bxe中选取第η个QRS波群位置Bxeη，求出该QRS波群的RR间期RRxen=Bxen-Bxeη-1;从信号X中截取一片段Xicis，截取的起始位置为Bxen-1+RRxenXk3，终止位置为^611-1«^11\0.5-1«，其中0.3k4XIVmin-VicisI则Qm点为Lmin点，否则Qm点为Pmx点，其中2hl，再与每分割段的包络差作比较;兰eiC1XB时，说明此分隔段的包络有突变，即为包络检测不合格；5204、在处理段X内，获取每一分割段Xi的幅值方差值，其中，111为中数据的序号;再求取处理段X内波形片段幅值方差值平均值筠，乘以比例系数k2，其中2k21，再与各分割段Xi的幅值方差值Ei作比较，当Epk2X筠时，说明此分割段的方差有突变，即为方差检测不合格；5205、在处理段X内，利用快速傅里叶变换，将心电片段信号X1的时域信号转化为频域信号，即为功率谱信号；对1〜5Hz的幅值进行积分得到功率，其中Re代表求复数的实部，Im代表求复数的虚部，j指某一频段的范围，fdj即为功率谱信号;求取处理段X内波形的平均功率值P，乘以比例系数k3，其中2k31，再与各分割段Xi的功率Pi作比较；当P,k;xP时，说明此处理段的1〜5Hz功率有突变，即为功率谱检测不合格；5206、在处理段X内，得到步骤S5中的心电片段信号Xi的功率谱信号后，对5〜40Hz的幅值进行积分得到功率.其中Re代表求复数的实部，Im代表求复数的虚部，Psi代表心电信号的功率;再对40〜IOOHz的幅值进行积分得到功率其中Re代表求复数的实部，Im代表求复数的虚部，Pm代表高频噪声信号的功率;计算两者的信噪比ksNRl=PSlPm，当信噪比小于阈值£SRN时，其中Ksrn〈于2,说明此处理段信噪比过小，噪声所占比重过大，即信噪比检测不合格；S207、根据参数包络差、方差值、信号在1〜5Hz之间的功率和信噪比的合格情况对心电波形进行质量等级划分，其中心电波形的质量等级划分方法为：当信噪比参数不合格时，该处理段心电波形的质量等级系数直接评定为3,即波形存在严重噪声;在信噪比参数合格的前提下，根据包络差、方差值、信号在1〜5Hz之间的功率的合格情况来评定心电波形质量，当三个参数全部合格时，该处理段心电波形质量等级系数评定为〇,即波形良好；当三个参数中有两个或一个合格时，心电波形质量等级系数评定为1，即波形较差；当三个参数全部不合格时，心电波形质量等级系数评定为2,即波形差；心电波形的质量分成4个等级，分别用0〜3表示质量等级系数;质量等级系数0表示心电波形良好，质量等级系数1表示心电波形较差，质量等级系数2表示心电波形差，质量等级系数3表示心电波形存在严重噪声。即从波形质量的程度“〇”表示良好，“Γ表示较差，“2”表示差，“3”表示严重噪声。

百度查询： 江苏华康信息技术有限公司;南京贺普检测仪器有限公司 一种心电QRS波群的检测方法及其心电分析方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD