

本发明涉及心血管疾病诊断技术领域,具体涉及一种基于心电心音的儿童心血管疾病诊断方法。The present invention relates to the technical field of cardiovascular disease diagnosis, and in particular to a method for diagnosing cardiovascular disease in children based on electrocardiogram and heart sounds.
据统计,目前1000名新生儿中,大概有8人患有先天性心脏病,其中心室间隔缺损(Ventricular Septal Defect,VSD)是最为常见的先天性心脏病。根据心室间隔缺损患者的情况,有时必须及早进行手术,因此先天性心脏病这类心血管疾病的精确诊断,对儿童具有重大意义。According to statistics, currently, about 8 out of 1,000 newborns have congenital heart disease, of which ventricular septal defect (Ventricular Septal Defect, VSD) is the most common congenital heart disease. Depending on the condition of patients with ventricular septal defects, early surgery is sometimes necessary, so accurate diagnosis of cardiovascular diseases such as congenital heart disease is of great significance to children.
在心血管疾病的诊断方法中,听诊古时候就有,即根据心音判别心脏的健康状况,但是听诊对于医生来说,需要多年的经验积累。近年来随着计算机技术的不断发展,可以利用计算机从心音信号中抽取与心脏病或心脏功能相关的特征参数。因此,采用计算机进行心血管疾病的判别和心脏功能的评价成为可能。Among the diagnostic methods for cardiovascular diseases, auscultation has been around since ancient times, which is to judge the health status of the heart based on heart sounds. However, auscultation requires many years of experience for doctors. In recent years, with the continuous development of computer technology, computers can be used to extract characteristic parameters related to heart disease or cardiac function from heart sound signals. Therefore, it becomes possible to use computers to identify cardiovascular diseases and evaluate cardiac function.
目前,大多从心音信号中提取特征波形的峰值的时间幅度、时间间隔作为评价指标,进行血管疾病的判断,并没有考虑到心电信号和心音信号的结合,且评价的维度也局限于时域,没有考虑到频域、能量域等维度,因此,本申请提出了一种基于心电心音的儿童心血管疾病诊断方法。At present, most methods extract the time amplitude and time interval of the peak value of the characteristic waveform from the heart sound signal as evaluation indicators to judge vascular diseases. The combination of the ECG signal and the heart sound signal is not considered, and the evaluation dimension is also limited to the time domain. , does not consider dimensions such as frequency domain and energy domain. Therefore, this application proposes a method for diagnosing cardiovascular diseases in children based on ECG heart sounds.
发明内容Contents of the invention
针对现有技术中的缺陷,本发明提供了一种基于心电心音的儿童心血管疾病诊断方法;所述方法包括以下步骤:In view of the deficiencies in the prior art, the present invention provides a method for diagnosing cardiovascular diseases in children based on electrocardiogram and heart sounds; the method includes the following steps:
S1:获取待诊断患儿的心电信号和心音信号;S1: Obtain the ECG signal and heart sound signal of the child to be diagnosed;
S2:对所述心电信号和心音信号进行处理,得到多维心音特征;S2: Process the electrocardiogram signal and heart sound signal to obtain multi-dimensional heart sound characteristics;
S3:将所述多维心音特征输入诊断模型,得到诊断结果,其中,所述诊断结果为心血管疾病类型。S3: Input the multi-dimensional heart sound features into the diagnostic model to obtain a diagnostic result, where the diagnostic result is the type of cardiovascular disease.
优选地,S2中,对心电信号和心音信号进行处理,得到心电心音特征的Preferably, in S2, the ECG signal and the heart sound signal are processed to obtain the ECG and heart sound characteristics.
方法包括:Methods include:
S21:对心电信号和心音信号进行预处理;S21: Preprocess ECG signals and heart sound signals;
S22:根据预处理后的心电信号和心音信号,获取初始心音特征;S22: Obtain the initial heart sound characteristics based on the preprocessed ECG signal and heart sound signal;
S23:通过时域、频域、能量域分析方法对所述初始心音特征进行分析,得到多维心音特征,所述多维心音特征包括时域心音特征、频域心音特征和能量域心音特征。S23: Analyze the initial heart sound features through time domain, frequency domain, and energy domain analysis methods to obtain multidimensional heart sound features. The multidimensional heart sound features include time domain heart sound features, frequency domain heart sound features, and energy domain heart sound features.
优选地,S21中,对心电信号和心音信号进行预处理的方法包括:Preferably, in S21, the method for preprocessing the ECG signal and the heart sound signal includes:
S211:对心电信号和心音信号进行截取;S211: Intercept ECG signals and heart sound signals;
S212:对采样后的心电信号和心音信号进行滤波降噪处理;S212: Perform filtering and noise reduction processing on the sampled ECG signals and heart sound signals;
S213:对滤波降噪后的心电信号和心音信号进行归一化处理。S213: Normalize the ECG signal and heart sound signal after filtering and noise reduction.
优选地,S22中,根据预处理后的心电信号和心音信号,获取初始心音特征的方法包括:Preferably, in S22, the method for obtaining initial heart sound features based on the preprocessed ECG signal and heart sound signal includes:
S221:识别预处理后的心电信号的QRS波;S221: Identify the QRS wave of the preprocessed ECG signal;
S222:根据QRS波的Q波,从预处理后的心音信号中识别出第一心音的持续和幅值信息,根据所述QRS波的T波,从预处理后的心音信号中识别出第二心音的持续时间和幅值信息;S222: According to the Q wave of the QRS wave, identify the duration and amplitude information of the first heart sound from the preprocessed heart sound signal, and according to the T wave of the QRS wave, identify the second heart sound signal from the preprocessed heart sound signal. Duration and amplitude information of heart sounds;
S223:根据第一心音和第二心音的持续时间和幅值信息,获取收缩期的持 续时间和幅值信息、舒张期的持续时间和幅值信息;S223: According to the duration and amplitude information of the first heart sound and the second heart sound, obtain the duration and amplitude information of the systolic phase and the duration and amplitude information of the diastolic phase;
S224:将所述第一心音、第二心音、收缩期和舒张期的持续时间和幅值信息作为初始心音特征。S224: Use the duration and amplitude information of the first heart sound, the second heart sound, systole and diastole as initial heart sound features.
优选地,S23中,通过时域分析方法对所述初始心音特征进行分析,得到时域心音特征的方法包括:分别计算第一心音、收缩期、第二心音和舒张期两两之间的幅值比;分别计算相邻第一心音、收缩期、第二心音和舒张期之间的时间差;计算EMAT%;计算各个幅值比、各个时间差、EMAT%在整个心动周期的均值,得到11个时域心音特征。Preferably, in S23, the initial heart sound characteristics are analyzed through a time domain analysis method. The method of obtaining the time domain heart sound characteristics includes: respectively calculating the first heart sound, systole, second heart sound and diastole. Amplitude ratio; calculate the time difference between adjacent first heart sound, systole, second heart sound and diastole respectively; calculate EMAT%; calculate the average of each amplitude ratio, each time difference and EMAT% in the entire cardiac cycle, and get 11 time domain heart sound features.
优选地,S23中,通过频域分析方法对所述初始心音特征进行分析,得到时域心音特征的方法包括:对所述初始心音特征进行傅里叶变换,得到对应的频谱序列,其中,所述频谱序列包括第一心音频谱序列、第二心音频谱序列、收缩期频谱序列和舒张期频谱序列;分别提取第一心音频谱序列、第二心音频谱序列、收缩期频谱序列和舒张期频谱序列中频率为多个指定频率的数值,求取各个数值在所有心动周期的均值,并形成新的频谱序列;得到4个频域心音特征。Preferably, in S23, the initial heart sound features are analyzed through a frequency domain analysis method, and the method for obtaining the time domain heart sound features includes: performing Fourier transform on the initial heart sound features to obtain a corresponding spectrum sequence, wherein: The spectrum sequence includes a first heart rate spectrum sequence, a second heart rate spectrum sequence, a systolic spectrum sequence, and a diastolic spectrum sequence; the first heart rate spectrum sequence, the second heart rate spectrum sequence, the systolic spectrum sequence, and the diastolic spectrum are extracted respectively. The frequencies in the sequence are multiple values of specified frequencies, and the average value of each value in all cardiac cycles is calculated and a new spectrum sequence is formed; 4 frequency domain heart sound characteristics are obtained.
优选地,所述第一心音频谱序列、第二心音频谱序列、收缩期、舒张期对应的多个指定频率为60HZ、70HZ、80HZ、90HZ、100HZ、110HZ、120HZ、130HZ、80HZ、150HZ。Preferably, the plurality of designated frequencies corresponding to the first heartbeat spectrum sequence, the second heartbeat spectrum sequence, the systolic phase, and the diastolic phase are 60HZ, 70HZ, 80HZ, 90HZ, 100HZ, 110HZ, 120HZ, 130HZ, 80HZ, and 150HZ.
优选地,S23中,通过频域分析方法对所述初始心音特征进行分析,得到时域心音特征的方法包括:分别计算第一心音、收缩期、第二心音和舒张期两两之间的能量比,并求各个能量比在所有心动周期的均值,得到6个能量域心音特征。Preferably, in S23, the initial heart sound characteristics are analyzed through a frequency domain analysis method, and the method for obtaining the time domain heart sound characteristics includes: separately calculating the first heart sound, systolic period, second heart sound and diastolic period. Energy ratio, and find the average value of each energy ratio in all cardiac cycles to obtain 6 energy domain heart sound characteristics.
优选地,S3中,所述诊断模型包括输入层、隐藏层、输出层和softmax函数,其中,所述输入层由21个神经元组成,所述21个神经元对应多维心音 特征;所述隐藏层分别有7和8个神经元;所述输出层由个8神经元组成,对应心血管疾病的类型;最后使用一个softmax函数,用于解决多分类问题而创建。Preferably, in S3, the diagnostic model includes an input layer, a hidden layer, an output layer and a softmax function, wherein the input layer consists of 21 neurons, and the 21 neurons correspond to multi-dimensional heart sound features; the hidden The layers have 7 and 8 neurons respectively; the output layer consists of 8 neurons, corresponding to the types of cardiovascular diseases; finally, a softmax function is used to create a multi-classification problem.
本发明的有益效果体现在:本发明实施例提供了一种基于心电心音的儿童心血管疾病诊断方法,结合了心音信号和心电信号,充分考虑了心血管疾病在心音信号上的体现,提取了包括时域、频域和能量域心音特征的多维心音特征,通过构建诊断模型和多维心音特征,实现了心血管疾病的精准预测。另外,本发明还提供了一种多维心音特征的提取方法,为诊断模型进行心血管疾病精准预测提供了基础支持。本发明具体实用高效等优点。The beneficial effects of the present invention are reflected in: embodiments of the present invention provide a method for diagnosing cardiovascular diseases in children based on electrocardiogram and heart sounds, which combines heart sound signals and electrocardiogram signals and fully considers the manifestation of cardiovascular diseases in heart sound signals. Multi-dimensional heart sound features including time domain, frequency domain and energy domain heart sound features are extracted, and accurate prediction of cardiovascular diseases is achieved by building a diagnostic model and multi-dimensional heart sound features. In addition, the present invention also provides a method for extracting multi-dimensional heart sound features, which provides basic support for the diagnostic model to accurately predict cardiovascular diseases. The invention has the advantages of practicality, high efficiency and the like.
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案,下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍。在所有附图中,类似的元件或部分一般由类似的附图标记标识。附图中,各元件或部分并不一定按照实际的比例绘制。In order to more clearly explain the specific embodiments of the present invention or the technical solutions in the prior art, the drawings that need to be used in the description of the specific implementations or the prior art will be briefly introduced below. Throughout the drawings, similar elements or portions are generally identified by similar reference numerals. In the drawings, elements or parts are not necessarily drawn to actual scale.
图1为本发明实施例所提供的一种基于心电心音的儿童心血管疾病诊断方法的流程图;Figure 1 is a flow chart of a method for diagnosing cardiovascular diseases in children based on electrocardiogram and heart sounds provided by an embodiment of the present invention;
图2为本发明实施例所提供的诊断模型的构建图。Figure 2 is a construction diagram of a diagnostic model provided by an embodiment of the present invention.
下面将结合附图对本发明技术方案的实施例进行详细的描述。以下实施例仅用于更加清楚地说明本发明的技术方案,因此只作为示例,而不能以此来限制本发明的保护范围。The embodiments of the technical solution of the present invention will be described in detail below with reference to the accompanying drawings. The following examples are only used to illustrate the technical solutions of the present invention more clearly, and are therefore only examples and cannot be used to limit the scope of the present invention.
需要注意的是,除非另有说明,本申请使用的技术术语或者科学术语应当为本发明所属领域技术人员所理解的通常意义。It should be noted that, unless otherwise stated, the technical terms or scientific terms used in this application should have the usual meanings understood by those skilled in the art to which this invention belongs.
如图1所示,本发明实施例提供了一种基于心电心音的儿童心血管疾病诊 断方法;所述方法包括以下步骤:As shown in Figure 1, an embodiment of the present invention provides a method for diagnosing cardiovascular disease in children based on electrocardiogram and heart sounds; the method includes the following steps:
步骤一:获取待诊断患儿的心电信号和心音信号;Step 1: Obtain the ECG signal and heart sound signal of the child to be diagnosed;
需要说明的,获取待诊断患儿的心电信号和心音信号的方法为按照常规心电图操作方法将心电图十二导联放好,在此基础上同步放置高精心音探头,同步记录心音的特征,得到待诊断患儿的心电信号和心音信号。It should be noted that the method of obtaining the ECG signal and heart sound signal of the child to be diagnosed is to place the twelve leads of the ECG according to the conventional ECG operation method, and on this basis, simultaneously place the high-precision heart sound probe and simultaneously record the characteristics of the heart sound. Obtain the ECG signal and heart sound signal of the child to be diagnosed.
步骤二:对所述心电信号和心音信号进行处理,得到多维心音特征;Step 2: Process the electrocardiogram signal and heart sound signal to obtain multi-dimensional heart sound characteristics;
在本发明实施例中,对心电信号和心音信号进行处理,得到心电心音特征的方法包括:In the embodiment of the present invention, the method of processing the ECG signal and the heart sound signal to obtain the ECG and heart sound characteristics includes:
A:对心电信号和心音信号进行预处理;其中,对心电信号和心音信号进行预处理的方法包括:对心电信号和心音信号进行截取;对采样后的心电信号和心音信号进行滤波降噪处理;对滤波降噪后的心电信号和心音信号进行归一化处理;A: Preprocess the ECG signal and heart sound signal; among them, the method of preprocessing the ECG signal and heart sound signal includes: intercepting the ECG signal and heart sound signal; performing sampled ECG signal and heart sound signal. Filtering and noise reduction processing; normalizing the ECG signal and heart sound signal after filtering and noise reduction;
具体的,对心电信号和心音信号进行截取时,需要保证被截取的心电信号和心音信号至少包括11个心动周期,以满足诊断模型的输入序列的长度需求。滤波降噪处理采用无限冲激相应数字带组滤波器进行滤波,该滤波器为现有技术,其工作过程在此不做赘述。Specifically, when intercepting ECG signals and heart sound signals, it is necessary to ensure that the intercepted ECG signals and heart sound signals include at least 11 cardiac cycles to meet the length requirements of the input sequence of the diagnostic model. The filtering and noise reduction process uses an infinite impulse corresponding digital band group filter for filtering. This filter is an existing technology, and its working process will not be described in detail here.
B:根据预处理后的心电信号和心音信号,获取初始心音特征;方法具体为识别预处理后的心电信号的QRS波;根据QRS波的Q波,从预处理后的心音信号中识别出第一心音的持续和幅值信息,根据所述QRS波的T波,从预处理后的心音信号中识别出第二心音的持续时间和幅值信息;根据第一心音和第二心音的持续时间和幅值信息,获取收缩期的持续时间和幅值信息、舒张期的持续时间和幅值信息;将所述第一心音、第二心音、收缩期和舒张期的持续时间和幅值信息作为初始心音特征;B: Obtain the initial heart sound characteristics based on the preprocessed ECG signal and heart sound signal; the method is specifically to identify the QRS wave of the preprocessed ECG signal; based on the Q wave of the QRS wave, identify the preprocessed heart sound signal The duration and amplitude information of the first heart sound is obtained, and the duration and amplitude information of the second heart sound are identified from the preprocessed heart sound signal according to the T wave of the QRS wave; according to the first heart sound and the second heart sound The duration and amplitude information of the heart sound is obtained, and the duration and amplitude information of the systolic phase and the duration and amplitude information of the diastolic phase are obtained; the duration of the first heart sound, the second heart sound, the systolic phase and the diastolic phase are and amplitude information as initial heart sound features;
C:通过时域、频域、能量域分析方法对所述初始心音特征进行分析,得到 多维心音特征,所述多维心音特征包括时域心音特征、频域心音特征和能量域心音特征。C: Analyze the initial heart sound features through time domain, frequency domain, and energy domain analysis methods to obtain multidimensional heart sound features. The multidimensional heart sound features include time domain heart sound features, frequency domain heart sound features, and energy domain heart sound features.
需要说明的,通过时域分析方法对所述初始心音特征进行分析,得到时域心音特征的方法包括:分别计算第一心音、收缩期、第二心音和舒张期两两之间的幅值比;分别计算相邻第一心音、收缩期、第二心音和舒张期之间的时间差;分别计算每个心动周期的EMAT%;计算各个幅值比、各个时间差、各个EMAT%;得到11个时域心音特征。It should be noted that the initial heart sound characteristics are analyzed through the time domain analysis method, and the method of obtaining the time domain heart sound characteristics includes: separately calculating the amplitudes between the first heart sound, systole, second heart sound and diastole. Ratio; calculate the time difference between adjacent first heart sound, systole, second heart sound and diastole respectively; calculate EMAT% of each cardiac cycle respectively; calculate each amplitude ratio, each time difference and each EMAT%; get 11 time domain heart sound characteristics.
具体的,计算每个心动周期中上述所述的各个幅值比和各个时间差和各个EMAT%构建6个幅值比序列、4个时间差序列、1个EMAT%序列,得到11个时域心音特征。Specifically, each amplitude ratio, each time difference and each EMAT% mentioned above in each cardiac cycle are calculated to construct 6 amplitude ratio sequences, 4 time difference sequences, and 1 EMAT% sequence, and 11 time domain heart sound features are obtained. .
其中,EMAT是从心电图中QRS波起始至第一心音开始(二尖瓣关闭)的时间,包括电-机械延迟时间和二尖瓣关闭前的左室收缩期,是左室收缩前期的一部分,EMAT%是指EMAT在RR间期所占的比例,经研究表明,EMAT%对心血管疾病的具有显著影响,因此,本发明实施例将EMAT%序列作为一个时域心音特征。Among them, EMAT is the time from the onset of the QRS wave in the electrocardiogram to the onset of the first heart sound (mitral valve closure), including the electro-mechanical delay time and the left ventricular systolic period before mitral valve closure. It is the presystolic period of the left ventricle. In part, EMAT% refers to the proportion of EMAT in the RR interval. Studies have shown that EMAT% has a significant impact on cardiovascular disease. Therefore, the embodiment of the present invention uses the EMAT% sequence as a time-domain heart sound feature.
为保证序列长度的一致性,若序列长度是预设序列长度的倍数,则根据倍数对序列进行分段,求取各个分段中各个数值的平均值作为新的序列,若序列长度不是预设序列的倍数,则从序列的第一数值进行删除,直到序列长度为预设序列长度的倍数。具体的,本发明实施例所提供的序列长度为10,之前步骤截取的心动周期为12,则构建的时间差序列长度为11,此时删除时间长度序列的第一个数值,得到新的时间差序列,新的时间差序列为预设序列长度的1倍,则将1作为分段参数将时间差序列分为10段,求取10段的平均值,得到长度为10的时间差序列。In order to ensure the consistency of the sequence length, if the sequence length is a multiple of the preset sequence length, the sequence will be segmented according to the multiple, and the average value of each value in each segment will be calculated as the new sequence. If the sequence length is not the preset sequence length, Multiples of the sequence, delete the first value of the sequence until the sequence length is a multiple of the preset sequence length. Specifically, the sequence length provided by the embodiment of the present invention is 10, and the cardiac cycle intercepted in the previous step is 12, then the length of the constructed time difference sequence is 11. At this time, the first value of the time length sequence is deleted to obtain a new time difference sequence. , the new time difference sequence is 1 times the length of the preset sequence, then use 1 as the segmentation parameter to divide the time difference sequence into 10 segments, calculate the average of the 10 segments, and obtain a time difference sequence with a length of 10.
需要说明的,通过频域分析方法对所述初始心音特征进行分析,得到时域 心音特征的方法包括:对所述初始心音特征进行傅里叶变换,得到对应的频谱序列,其中,所述频谱序列包括第一心音频谱序列、第二心音频谱序列、收缩期频谱序列和舒张期频谱序列;分别提取第一心音频谱序列、第二心音频谱序列、收缩期频谱序列和舒张期频谱序列中频率为多个指定频率的数值,求取各个数值在所有心动周期的均值,并形成新的频谱序列;得到4个频域心音特征。It should be noted that the method of analyzing the initial heart sound features through the frequency domain analysis method to obtain the time domain heart sound features includes: performing Fourier transform on the initial heart sound features to obtain the corresponding spectrum sequence, wherein the spectrum The sequence includes the first heart rate spectrum sequence, the second heart rate spectrum sequence, the systolic spectrum sequence and the diastolic spectrum sequence; respectively extract the first heart rate spectrum sequence, the second heart rate spectrum sequence, the systolic spectrum sequence and the diastolic spectrum sequence. The frequency is a number of values at specified frequencies. The average value of each value in all cardiac cycles is calculated and a new spectrum sequence is formed; 4 frequency domain heart sound features are obtained.
由于第一心音和第二心音的频谱分布是50~100HZ,若超过100HZ就判定为存在杂音,并且预设的序列长度为10,因此,在本发明实施例中选定的5个指定频率的正常组和5个指定频率的异常组构建频谱心音特征,所述第一心音频谱序列、第二心音频谱序列、收缩期、舒张期对应的多个指定频率为60HZ、70HZ、80HZ、90HZ、100HZ、110HZ、120HZ、130HZ、80HZ、150HZ,并从而得到4个长度为10的频域心音特征。Since the spectrum distribution of the first heart sound and the second heart sound is 50-100HZ, if it exceeds 100HZ, it is determined that there is a murmur, and the preset sequence length is 10. Therefore, the five designated frequencies selected in the embodiment of the present invention The normal group and the abnormal group of 5 specified frequencies construct spectral heart sound characteristics. The multiple specified frequencies corresponding to the first heart sound spectrum sequence, the second heart sound spectrum sequence, the systolic period, and the diastolic period are 60HZ, 70HZ, 80HZ, and 90HZ. , 100HZ, 110HZ, 120HZ, 130HZ, 80HZ, 150HZ, and thus obtain 4 frequency domain heart sound features with a length of 10.
需要说明的,通过频域分析方法对所述初始心音特征进行分析,得到时域心音特征的方法包括:分别计算第一心音、收缩期、第二心音和舒张期两两之间的能量比,并求各个能量比在所有心动周期的均值,得到6个能量域心音特征。为保证能量域心音特征的长度满足要求,对序列长度进行处理的方法与上述获取用于输对时域心音特征序列长度进行处理的方法相同,在此不做赘述。It should be noted that the method of analyzing the initial heart sound characteristics through the frequency domain analysis method and obtaining the time domain heart sound characteristics includes: separately calculating the energy ratio between the first heart sound, systole, second heart sound and diastole. , and calculate the average value of each energy ratio in all cardiac cycles to obtain 6 energy domain heart sound characteristics. In order to ensure that the length of the energy domain heart sound feature meets the requirements, the method for processing the sequence length is the same as the above-mentioned method for obtaining the sequence length for transmitting the time domain heart sound feature, and will not be described again here.
步骤三:将所述多维心音特征输入诊断模型,得到诊断结果,其中,所述诊断结果为心血管疾病类型。Step 3: Input the multi-dimensional heart sound features into the diagnostic model to obtain a diagnostic result, where the diagnostic result is the type of cardiovascular disease.
如图2所示,在本发明实施例中,所述诊断模型包括输入层、隐藏层、输出层和softmax函数,其中,所述输入层由21个神经元组成,所述21个神经元对应多维心音特征;所述隐藏层分别有7和8个神经元,所述隐藏层为2层双向长短期记忆网络模型;所述输出层由个8个神经元组成,对应心血管疾病的类型;最后使用一个softmax函数,用于解决多分类问题而创建,其输出结果为多个心血管疾病的概率,其概率为0-1之间的数值,且得到的概率之和为1。As shown in Figure 2, in the embodiment of the present invention, the diagnostic model includes an input layer, a hidden layer, an output layer and a softmax function, wherein the input layer is composed of 21 neurons, and the 21 neurons correspond to Multi-dimensional heart sound features; the hidden layer has 7 and 8 neurons respectively, and the hidden layer is a 2-layer bidirectional long short-term memory network model; the output layer is composed of 8 neurons, corresponding to the type of cardiovascular disease; Finally, a softmax function is used, which is created to solve multi-classification problems. The output result is the probability of multiple cardiovascular diseases. The probability is a value between 0-1, and the sum of the obtained probabilities is 1.
在其他实施例中,所述输出层的神经元可以是8个也可以是50个或者100个,可以根据现有的常见心血管疾病进行设置,若增加输出层神经元就需要增加诊断模型训练时的样本,以保证诊断模型诊断的精确性。In other embodiments, the number of neurons in the output layer can be 8, 50, or 100, and can be set according to existing common cardiovascular diseases. If the output layer neurons are added, diagnostic model training needs to be added. time samples to ensure the accuracy of diagnosis by the diagnostic model.
在其他实施例中,所述诊断模型可以是全连接神经网络模型,2层全连接层作为全连接神经网络模型的隐藏层,S型生长曲线层作为全连接神经网络模型的输出层;当所述诊断模型为全连接神经网络模型是,其输出结果为多个类型的心血管疾病的概率,其概率为0-1的数值,且输出的概率的和不一定为1。In other embodiments, the diagnostic model may be a fully connected neural network model, with two fully connected layers as hidden layers of the fully connected neural network model, and an S-shaped growth curve layer as the output layer of the fully connected neural network model; when the The diagnostic model is a fully connected neural network model, and its output results are the probabilities of multiple types of cardiovascular diseases. The probabilities are 0-1 values, and the sum of the output probabilities is not necessarily 1.
综上,本发明实施例提供了一种基于心电心音的儿童心血管疾病诊断方法,结合了心音信号和心电信号,充分考虑了心血管疾病在心音信号上的体现,提取了包括时域、频域和能量域心音特征的多维心音特征,通过构建诊断模型和多维心音特征,实现了心血管疾病的精准预测。另外,本发明还提供了一种多维心音特征的提取方法,为诊断模型进行心血管疾病精准预测提供了基础支持。本发明具体实用高效等优点。In summary, embodiments of the present invention provide a method for diagnosing cardiovascular diseases in children based on ECG heart sounds, which combines heart sound signals and ECG signals, fully considers the manifestation of cardiovascular diseases in heart sound signals, and extracts time domain , multi-dimensional heart sound features of frequency domain and energy domain heart sound features, and by building diagnostic models and multi-dimensional heart sound features, accurate prediction of cardiovascular diseases is achieved. In addition, the present invention also provides a method for extracting multi-dimensional heart sound features, which provides basic support for the diagnostic model to accurately predict cardiovascular diseases. The invention has the advantages of practicality, high efficiency and the like.
最后应说明的是:以上各实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述各实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的范围,其均应涵盖在本发明的权利要求和说明书的范围当中。Finally, it should be noted that the above embodiments are only used to illustrate the technical solution of the present invention, but not to limit it. Although the present invention has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that: The technical solutions described in the foregoing embodiments can still be modified, or some or all of the technical features can be equivalently replaced; and these modifications or substitutions do not deviate from the essence of the corresponding technical solutions from the technical solutions of the embodiments of the present invention. scope, they should be covered by the claims and the scope of the description of the present invention.
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| PCT/CN2022/103931WO2024007152A1 (en) | 2022-07-05 | 2022-07-05 | Method for diagnosing pediatric cardiovascular diseases based on electrocardiographic and phonocardiographic signals |
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| PCT/CN2022/103931WO2024007152A1 (en) | 2022-07-05 | 2022-07-05 | Method for diagnosing pediatric cardiovascular diseases based on electrocardiographic and phonocardiographic signals |
| Publication Number | Publication Date |
|---|---|
| WO2024007152A1true WO2024007152A1 (en) | 2024-01-11 |
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| PCT/CN2022/103931CeasedWO2024007152A1 (en) | 2022-07-05 | 2022-07-05 | Method for diagnosing pediatric cardiovascular diseases based on electrocardiographic and phonocardiographic signals |
| Country | Link |
|---|---|
| WO (1) | WO2024007152A1 (en) |
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| TWI880187B (en)* | 2023-03-24 | 2025-04-11 | 萬潤科技股份有限公司 | Conveying devices, conveying equipment and conveying systems |
| CN120613109A (en)* | 2025-08-06 | 2025-09-09 | 北京麦迪克斯科技有限公司 | A heart sound and electrocardiogram acquisition and analysis system |
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20050222515A1 (en)* | 2004-02-23 | 2005-10-06 | Biosignetics Corporation | Cardiovascular sound signature: method, process and format |
| US20090192401A1 (en)* | 2008-01-25 | 2009-07-30 | Sourabh Ravindran | Method and system for heart sound identification |
| US20130226019A1 (en)* | 2010-08-25 | 2013-08-29 | Diacoustic Medical Devices (Pty) Ltd | System and method for classifying a heart sound |
| US20140276132A1 (en)* | 2013-03-15 | 2014-09-18 | Andreas J. Schriefl | Automated Diagnosis-Assisting Medical Devices Utilizing Rate/Frequency Estimation |
| CN111329508A (en)* | 2020-03-02 | 2020-06-26 | 浙江大学 | Heart murmur intelligent analysis method for precordial disease screening |
| CN112971839A (en)* | 2021-02-05 | 2021-06-18 | 云南大学 | Heart sound classification method based on feedforward convolutional neural network |
| CN115040135A (en)* | 2022-06-09 | 2022-09-13 | 张福伟 | Children cardiovascular disease diagnosis method based on electrocardio and heart sounds |
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20050222515A1 (en)* | 2004-02-23 | 2005-10-06 | Biosignetics Corporation | Cardiovascular sound signature: method, process and format |
| US20090192401A1 (en)* | 2008-01-25 | 2009-07-30 | Sourabh Ravindran | Method and system for heart sound identification |
| US20130226019A1 (en)* | 2010-08-25 | 2013-08-29 | Diacoustic Medical Devices (Pty) Ltd | System and method for classifying a heart sound |
| US20140276132A1 (en)* | 2013-03-15 | 2014-09-18 | Andreas J. Schriefl | Automated Diagnosis-Assisting Medical Devices Utilizing Rate/Frequency Estimation |
| CN111329508A (en)* | 2020-03-02 | 2020-06-26 | 浙江大学 | Heart murmur intelligent analysis method for precordial disease screening |
| CN112971839A (en)* | 2021-02-05 | 2021-06-18 | 云南大学 | Heart sound classification method based on feedforward convolutional neural network |
| CN115040135A (en)* | 2022-06-09 | 2022-09-13 | 张福伟 | Children cardiovascular disease diagnosis method based on electrocardio and heart sounds |
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| TWI880187B (en)* | 2023-03-24 | 2025-04-11 | 萬潤科技股份有限公司 | Conveying devices, conveying equipment and conveying systems |
| CN120613109A (en)* | 2025-08-06 | 2025-09-09 | 北京麦迪克斯科技有限公司 | A heart sound and electrocardiogram acquisition and analysis system |
| Publication | Publication Date | Title |
|---|---|---|
| CN111759345B (en) | Heart valve abnormality analysis method, system and device based on convolutional neural network | |
| Leung et al. | Classification of heart sounds using time-frequency method and artificial neural networks | |
| CN110537910A (en) | Non-invasive screening system for coronary heart disease based on joint analysis of ECG and heart sound signals | |
| Khan et al. | Artificial neural network-based cardiovascular disease prediction using spectral features | |
| JPH021217A (en) | Apparatus and method for detecting heart disease | |
| Javed et al. | A signal processing module for the analysis of heart sounds and heart murmurs | |
| Choi et al. | Selection of wavelet packet measures for insufficiency murmur identification | |
| Sedighian et al. | Pediatric heart sound segmentation using Hidden Markov Model | |
| CN112971839B (en) | Heart sound classification method based on feedforward convolution neural network | |
| Sinha et al. | Backpropagation artificial neural network classifier to detect changes in heart sound due to mitral valve regurgitation | |
| WO2024007152A1 (en) | Method for diagnosing pediatric cardiovascular diseases based on electrocardiographic and phonocardiographic signals | |
| CN113080991A (en) | Method, system, diagnosis device and storage medium for predicting and diagnosing heart failure based on CNN model and LSTM model | |
| Nizam et al. | Hilbert-envelope features for cardiac disease classification from noisy phonocardiograms | |
| Deperlioglu | Segmentation of heart sounds by re-sampled signal energy method | |
| CN115040135A (en) | Children cardiovascular disease diagnosis method based on electrocardio and heart sounds | |
| CN115281689A (en) | Heart disease comprehensive diagnosis platform based on artificial intelligence and computer vision | |
| CN111329508A (en) | Heart murmur intelligent analysis method for precordial disease screening | |
| Akbari et al. | Digital Subtraction Phonocardiography (DSP) applied to the detection and characterization of heart murmurs | |
| CN117349600A (en) | Heart sound and heart electricity combined diagnosis method and system based on dual-mode dual input | |
| Golpaygani et al. | Detection and identification of S1 and S2 heart sounds using wavelet decomposition method | |
| Ari et al. | On a robust algorithm for heart sound segmentation | |
| Yue et al. | Time-frequency analysis of heart sounds in telemedicine consulting system for auscultation | |
| Nigam et al. | A simplicity-based fuzzy clustering approach for detection and extraction of murmurs from the phonocardiogram | |
| CN109498055B (en) | Heart sound signal acquisition and analysis system and signal analysis method | |
| Botha et al. | Autonomous auscultation of the human heart employing a precordial electro-phonocardiogram and ensemble empirical mode decomposition |
| Date | Code | Title | Description |
|---|---|---|---|
| 121 | Ep: the epo has been informed by wipo that ep was designated in this application | Ref document number:22949740 Country of ref document:EP Kind code of ref document:A1 | |
| NENP | Non-entry into the national phase | Ref country code:DE | |
| 122 | Ep: pct application non-entry in european phase | Ref document number:22949740 Country of ref document:EP Kind code of ref document:A1 |