基于时频域熵特征的子宫肌电信号分析及早产预测方法技术

一种基于时频域熵特征的子宫肌电信号分析及早产预测方法，设以频率f

【技术实现步骤摘要】
基于时频域熵特征的子宫肌电信号分析及早产预测方法
本专利技术涉及一种基于时频域熵特征的子宫肌电信号(EHG)分析及早产预测方法
技术介绍
早产(PretermBirth，怀孕小于37周的新生儿)是全球范围的公共健康问题之一，即使在欧洲和其他一些发达国家，早产约占新生儿7％左右。早产不仅会给家庭带来严重的情感与经济负担，也是社会医疗支付的重要组成部分。虽然围产期和新生儿的护理水平已经有了很大的提高，但早产儿仍然具有很高的发病率和死亡率。据统80％的新生儿死亡病例与早产直接或间接相关。为此，早产诊断与预防方法一直是临床医学的研究热点之一。引起早产的直接原因是子宫过早地产生强有力的收缩行为。宫缩(uterinecontraction)是整个孕期中的普遍现象，也是孕妇在孕期求医主要原因。然而医生却很难正确判断这些宫缩是否真正会导致分娩的产生：在收治的被诊断为有早产风险的孕妇中，有近50％没有出现真正的分娩症状，直接的结果是被误诊为早产的胎儿接受了不必要的医学干预；而临床诊断为没有早产风险的孕妇中，却有近20％出现了早产。有效地诊断早产并进行相应治疗不仅能够降低早产儿的死亡率，而还有助于加快相关药物的研发。现有诊断早产的常规方法包括：子宫分娩力计(Tocodynamometry),宫内压力计(Intra-uterinepressurecatheter)，胎儿纤连蛋白(fetalfibronectin),宫颈长度测量(cervicallengthmeasurement)等，但这些方法都存在一定缺陷，或无法获得可靠的预测效果。不管是早产还是非早产，宫缩都是伴随着子宫肌细胞中动作电位(ActionPotential)的产生与传播而产生的。随着怀孕过程的进行，肌细胞的可激发性(excitability)变强，细胞间的电耦合作用也加大，使得动作电位在整个子宫内传播，从而产生强有力的收缩力将胎儿压出体外。子宫壁上的动作电位及其传播可以通过置于孕妇腹部的电极进行采集，得到子宫肌电信号(electrohysterogram,EHG)。子宫肌电信号(EHG)分析，是一种潜在有效的早产预测方法。准确的早产预测方法，可以使有可能发生早产的孕妇得到及时的医学治疗，避免给整个家庭带来巨大的经济及精神负担，具有良好的应用前景。早产预测大都采用基于机器学习的分类方法，因此样本数据获取是决定性因素之一。医疗行业的数据分析是前景广阔的研究和发展领域。医疗数据分析的核心是数据。为了推进全世界范围内基于医疗数据的研究和疾病诊断方法研究，已经建立起了开源的生理信号数据库PhysioBank。其中关于孕期子宫EHG信号的数据库有两个：TPEHG数据库(THETERM-PRETERMEHGDATABASE)和冰16电极子宫肌电信号数据库(ICELANDIC16-ELECTRODEELECTROHYSTEROGRAMDATABASE)。TPEHG数据库包含了300位孕妇在怀孕中期(20-30周)记录的子宫EHG信号，其中38位最后出现早产，剩下的262位正常分娩。区别于TPEHG数据，冰岛16电极EHG数据库记录45位孕妇在整个孕期多次采集的子宫肌电信号。这些开源数据促进了基于子宫EHG信号的早产预测方法研究，并取了一系列成果。除了数据样本之外，如何从样本数据中选取有效特征是决定早产预测效果的另一关键因素。EHG信号特征可以分为三类：时间特征(Timing-related),幅度特征(Amplitude-related)和频率特征(Frequency-related)。考虑到子宫收缩力与爆发性动作电位(Bursts)密切相关，时间特征主要指爆发性动作电位持续时间(interburstduration),周期，标准差等。基于这些特征的早产预测结果差异性大：有文献指出即使分娩前后7天的EHG信号也没有表现出明显的差别；而又有文献则指出临产前期爆发性动作电位持续时间及频率要明显增加。幅度特征通常以EHG信号的RMS或者其功率的峰值强度来表示，表征的是采集到的子宫EHG信号的强弱。由于信号采集过程中易受到外部因素影响(脂肪厚度、皮肤电导率等)，幅度特征的早产症状表征功能不强。常用的频率特征有：峰值频率(peakfrequency)，中值频率(medianfrequency)和均值功率频率(meanpowerfrequency)等。由于子宫EHG信号是非稳态的，单纯的频谱分析并不能反映怀孕不断发展的过程。为了表征这种随时间不断变化的性质，结合时间频率特性的特征提取(小波分析)方法受到重视，并被广泛应用。为了进一步提高预测结果的正确性，新的特征不断被加入，如近似熵(ApproximateEntropy)，采样熵(SampleEntropy)等，然而这些特征并没有很明显的区分性。研究人员发现，EHG信号的信噪比低，且其中存在很多不同频率的噪声，例如孕妇呼吸的噪声主要存在于0.2HZ-0.34HZ。EHG信号中存在的噪声，使特征的有效性减弱。Mischi等人对熵特征进行改进，提出了一种基于向量夹角信息的原始距离度量的修改方法，旨在限制大幅度波动和尖峰的影响，对EHG信号中存在的噪声进行处理，提高了熵特征的有效性。但是，EHG信号的幅度信息对于早产分类十分重要，这种改进方法并没有考虑这一点。且单纯的时域分析无法对EHG信号中不同频率的噪声进行很好的处理，对早产分类的效果提升有限。相比于整个怀孕周期，EHG信号的采集时间极短，计算得到的熵特征无法表征怀孕不断发展的过程。为了取得好的早产预测效果，我们对从EHG信号中提取的熵特征进行改进。
技术实现思路
为了克服现有技术的不足，为解决从EHG信号提取有效熵特征用于早产预测的问题，本专利技术提出了基于时频域熵特征的子宫肌电信号(EHG)分析及早产预测方法，使用短时傅里叶变换将EHG信号分解成不同频率成分的时序信号，计算不同频率成分的熵特征，并利用信号的记录时间进一步对熵特征进行改进，接着应用主成分分析技术对所有时频域熵特征进行分析，选择前20个主成分作为新提出的熵特征，从而达到去噪的目的。最后基于新提出的熵特征，训练机器学习模型，生成一个高准确率的早产预测系统。为了解决上述技术问题，本专利技术采用如下的技术方案：一种基于时频域熵特征的子宫肌电信号分析及早产预测方法，包括如下步骤：步骤1，设以频率fs采集子宫肌电信号(EHG)，得到时序信号Xi(i＝1,2,3,…,N)；步骤2，通过短时傅里叶变换将时序信号Xi分解成不同频率分量yω(t)；步骤3，针对不同频率ωx的分量信号yω(t)计算相应频率下的近似熵特征和样本熵特征步骤4，利用信号的记录时间tr分别对每个和进行改进得到和步骤5，对得到的所有特征应用主成分分析技术，选择最优数量的主成分得到特征，对得到的所有特征应用主成分分析技术，选择最优数量的主成分得到特征，去除EHG信号中噪声的影响；步骤6，基于新特征和选择合适的机器学习模型，使用样本来训练模型，提高模型早产预测的准确率；步骤7，使用训练好的机器学习模型，本文档来自技高网...

【技术保护点】
1.一种基于时频域熵特征的子宫肌电信号分析及早产预测方法，其特征在于，所述方法包括如下步骤：/n步骤1，设以频率f

【技术特征摘要】
1.一种基于时频域熵特征的子宫肌电信号分析及早产预测方法，其特征在于，所述方法包括如下步骤：
步骤1，设以频率fs采集子宫肌电信号EHG，得到时序信号Xi(i＝1,2,3,…,N)；
步骤2，通过短时傅里叶变换将时序信号Xi分解成不同频率分量yω(t)；
步骤3，针对不同频率ωx的分量信号yω(t)计算相应频率下的近似熵特征和样本熵特征
步骤4，利用信号的记录时间tr分别对每个和进行改进得到和
步骤5，对得到的所有特征应用主成分分析技术，选择最优数量的主成分得到特征，对得到的所有特征应用主成分分析技术，选择最优数量的主成分得到特征，去除EHG信号中噪声的影响；
步骤6，基于新特征和选择合适的机器学习模型，使用样本来训练模型，提高模型早产预测的准确率；
步骤7，使用训练好的机器学习模型，构建一个高准确率的早产预测系统。


2.如权利要求1所述的基于时频域熵特征的子宫肌电信号分析及早产预测方法，其特征在于，所述步骤2中，设以采集频率fs获得的子宫EHG时序信号Xi(i＝1,2,3,…,N)，引入短时傅里叶变换，将原始EHG信号分解成不同频率时序信号，即：
yω(t)＝STFT(Xi,fs,T0)(1)<...

【专利技术属性】
技术研发人员：楼航晓，许金山，胡凤丹，陈镇钦，汪梦婷，
申请(专利权)人：浙江工业大学，
类型：发明
国别省市：浙江;33