CN116250844B - Electrocardiosignal noise reduction optimization method and system based on condition generation countermeasure network - Google Patents
Electrocardiosignal noise reduction optimization method and system based on condition generation countermeasure network Download PDFInfo
- Publication number
- CN116250844B CN116250844B CN202310225141.2A CN202310225141A CN116250844B CN 116250844 B CN116250844 B CN 116250844B CN 202310225141 A CN202310225141 A CN 202310225141A CN 116250844 B CN116250844 B CN 116250844B
- Authority
- CN
- China
- Prior art keywords
- noise reduction
- model
- electrocardiosignal
- deep neural
- signal
- 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.)
- Active
Links
- 230000009467 reduction Effects 0.000 title claims abstract description 90
- 238000000034 method Methods 0.000 title claims abstract description 55
- 238000005457 optimization Methods 0.000 title claims abstract description 43
- 238000013528 artificial neural network Methods 0.000 claims abstract description 61
- 238000004364 calculation method Methods 0.000 claims abstract description 22
- 238000010606 normalization Methods 0.000 claims abstract description 9
- 238000012545 processing Methods 0.000 claims abstract description 8
- 239000012634 fragment Substances 0.000 claims abstract description 5
- 230000006870 function Effects 0.000 claims description 35
- 210000002569 neuron Anatomy 0.000 claims description 24
- 230000008569 process Effects 0.000 claims description 23
- 230000015654 memory Effects 0.000 claims description 17
- 238000012360 testing method Methods 0.000 claims description 14
- 238000012549 training Methods 0.000 claims description 13
- 230000000694 effects Effects 0.000 claims description 11
- 238000004422 calculation algorithm Methods 0.000 description 8
- 238000005070 sampling Methods 0.000 description 7
- 230000002457 bidirectional effect Effects 0.000 description 6
- 230000006403 short-term memory Effects 0.000 description 6
- 238000013527 convolutional neural network Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 238000007635 classification algorithm Methods 0.000 description 3
- 238000013135 deep learning Methods 0.000 description 3
- 238000003745 diagnosis Methods 0.000 description 3
- 208000019622 heart disease Diseases 0.000 description 3
- 101001121408 Homo sapiens L-amino-acid oxidase Proteins 0.000 description 2
- 102100026388 L-amino-acid oxidase Human genes 0.000 description 2
- 230000004913 activation Effects 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000012821 model calculation Methods 0.000 description 2
- 230000011218 segmentation Effects 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 101000827703 Homo sapiens Polyphosphoinositide phosphatase Proteins 0.000 description 1
- 102100023591 Polyphosphoinositide phosphatase Human genes 0.000 description 1
- 101100012902 Saccharomyces cerevisiae (strain ATCC 204508 / S288c) FIG2 gene Proteins 0.000 description 1
- 101100233916 Saccharomyces cerevisiae (strain ATCC 204508 / S288c) KAR5 gene Proteins 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 206010003119 arrhythmia Diseases 0.000 description 1
- 230000006793 arrhythmia Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 210000003205 muscle Anatomy 0.000 description 1
- 230000037081 physical activity Effects 0.000 description 1
- 230000000306 recurrent effect Effects 0.000 description 1
- 238000007634 remodeling Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 238000012706 support-vector machine Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
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]
-
- 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
-
- 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/7264—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Signal Processing (AREA)
- Veterinary Medicine (AREA)
- Heart & Thoracic Surgery (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Pathology (AREA)
- Public Health (AREA)
- Medical Informatics (AREA)
- Artificial Intelligence (AREA)
- Cardiology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physiology (AREA)
- Psychiatry (AREA)
- Mathematical Physics (AREA)
- Evolutionary Computation (AREA)
- Theoretical Computer Science (AREA)
- Power Engineering (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Software Systems (AREA)
- Fuzzy Systems (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
Description
技术领域Technical Field
本发明涉及心电信号降噪技术领域,尤其涉及一种基于条件生成对抗网络的心电信号降噪优化方法及系统。The present invention relates to the technical field of electrocardiogram signal noise reduction, and in particular to an electrocardiogram signal noise reduction optimization method and system based on a conditional generative adversarial network.
背景技术Background technique
本部分的陈述仅仅是提供了与本发明相关的背景技术信息,不必然构成在先技术。The statements in this section merely provide background information related to the present invention and do not necessarily constitute prior art.
心电图是一种通过采集体表固定位置电生理信号来反映心脏活动情况的技术,目前已经被广泛用于心脏疾病的诊断,是心脏疾病医学诊断的金标准之一。但心电图信号(下称心电信号)采集过程易受到身体活动、电路噪声等各种噪声的干扰,进而影响心电信号用于心脏疾病的诊断。Electrocardiogram (ECG) is a technology that reflects heart activity by collecting electrophysiological signals at fixed locations on the body surface. It has been widely used in the diagnosis of heart disease and is one of the gold standards for medical diagnosis of heart disease. However, the acquisition process of ECG signals (hereinafter referred to as ECG signals) is easily interfered by various noises such as body activity and circuit noise, which in turn affects the use of ECG signals for the diagnosis of heart disease.
传统的心电信号降噪方法,如:滤波器、傅里叶变换、小波分解等,一般都是针对某类单一噪声,且噪声时频特征应与心电信号本身有明确的差异。随着深度学习的发展,针对多种噪声的深度学习降噪算法越来越多,如:降噪自编码器、卷积自编码器等方法,其中,基于全连接网络的对抗性降噪、基于卷积自编码器的生成对抗降噪等方法,一定程度上解决了上述算法中存在的心电信号需要在降噪前按照心拍进行分割,不能处理混合噪声等问题。Traditional ECG signal denoising methods, such as filters, Fourier transforms, wavelet decomposition, etc., are generally targeted at a single type of noise, and the time-frequency characteristics of the noise should be clearly different from the ECG signal itself. With the development of deep learning, there are more and more deep learning denoising algorithms for various noises, such as denoising autoencoders, convolutional autoencoders, etc. Among them, adversarial denoising based on fully connected networks and generative adversarial denoising based on convolutional autoencoders have, to a certain extent, solved the problems in the above algorithms, such as the need to segment ECG signals according to heart beats before denoising and the inability to handle mixed noise.
但是,上述算法还存在模型学习数据分布的潜力挖掘不够,模型缺乏合理、系统的优化,模型计算复杂度较高,难以在远程医疗设备上部署的限制。However, the above algorithms still have limitations: the model learning data distribution potential is not fully explored, the model lacks reasonable and systematic optimization, the model calculation complexity is high, and it is difficult to deploy on telemedicine equipment.
比如:for example:
现有技术公开了基于生成对抗网络进行心电信号降噪的方法,其中的生成器网络多采用全连接、卷积等网络结构,但无论采用卷积神经网络(CNN)或双向长短时记忆网络(BiLSTM)进行建模,由于缺少高效降噪网络模型建模的优化方法,导致这些现有技术所建立模型的结构较复杂、规模较大、计算复杂度较高,对部署环境要求苛刻。The prior art discloses methods for denoising ECG signals based on generative adversarial networks, wherein the generator networks mostly adopt fully connected, convolutional and other network structures. However, regardless of whether a convolutional neural network (CNN) or a bidirectional long short-term memory network (BiLSTM) is used for modeling, due to the lack of an optimization method for modeling an efficient denoising network model, the models established by these prior arts have a complex structure, a large scale, and a high computational complexity, and have stringent requirements on the deployment environment.
发明内容Summary of the invention
为了解决上述问题,本发明提出了一种基于条件生成对抗网络的心电信号降噪优化方法及系统,基于深度神经网络和条件生成对抗网络构建心电信号降噪模型架构,同时提出了一种兼顾降噪性能和计算开销的指标,对该降噪模型的内部结构、学习框架、数据分割、模型复杂度等方面进行了系统化优化设计,使其在花费较低的计算开销的同时获得更优的降噪效果。In order to solve the above problems, the present invention proposes an ECG signal denoising optimization method and system based on conditional generative adversarial network, constructs an ECG signal denoising model architecture based on deep neural network and conditional generative adversarial network, and proposes an indicator that takes into account both denoising performance and computational overhead. The internal structure, learning framework, data segmentation, model complexity and other aspects of the denoising model are systematically optimized and designed, so that better denoising effect can be achieved with lower computational overhead.
在一些实施方式中,采用如下技术方案:In some embodiments, the following technical solutions are adopted:
一种基于条件生成对抗网络的心电信号降噪优化方法,包括:An electrocardiogram signal denoising optimization method based on a conditional generative adversarial network comprises:
获取心电信号,将所述心电信号按照样本长度T进行分割,对分割后的数据片段进行最大最小值归一化处理;Acquire an electrocardiogram signal, segment the electrocardiogram signal according to a sample length T, and perform maximum and minimum value normalization processing on the segmented data segments;
将处理后的数据片段输入至训练好的心电信号降噪模型中,得到降噪后的心电信号;Input the processed data fragments into the trained ECG signal denoising model to obtain the denoised ECG signal;
其中,心电信号降噪模型基于深度神经网络和条件生成对抗网络构建,包括生成器和判别器;对于心电信号降噪模型中的每一种建模参量,按一定步长分别建立建模参量不同数值下(通常≥3种情况)的降噪模型以获得与其对应降噪性能与计算复杂度比率指标(The ratio of Signal-to-NoiseRatio to Computational Complexity,SNR-CC),采用拟合函数预测每一种建模参量的对应最优值,通过验证,以确定模型的降噪性能和计算开销达到最优;Among them, the ECG signal denoising model is constructed based on deep neural network and conditional generative adversarial network, including generator and discriminator; for each modeling parameter in the ECG signal denoising model, denoising models with different values of modeling parameters (usually ≥3 cases) are established at a certain step size to obtain the corresponding denoising performance and computational complexity ratio index (The ratio of Signal-to-NoiseRatio to Computational Complexity, SNR-CC), and the fitting function is used to predict the corresponding optimal value of each modeling parameter, and through verification, it is determined that the denoising performance and computational overhead of the model are optimal;
所述建模参量包括:样本长度T、生成器中编码器和解码器的深度神经网络层数N及每一层的神经元数量Lx和判别器的深度神经网络层数M及每一层的神经元数量Rx中的至少一种或多种。The modeling parameters include: at least one or more of the sample length T, the number of deep neural network layers N of the encoder and decoder in the generator and the number of neurons Lx in each layer, and the number of deep neural network layers M of the discriminator and the number of neurons Rx in each layer.
作为进一步的方案,所述生成器是由深度神经网络组成的降噪自编码器(DAE)所构成的,DAE中包含了N层深度神经网络组成的编码器以及N层深度神经网络组成的解码器,训练过程中的判别器由M层深度神经网络组成的二分类深度神经网络构成;As a further solution, the generator is composed of a denoising autoencoder (DAE) composed of a deep neural network, the DAE includes an encoder composed of N layers of deep neural networks and a decoder composed of N layers of deep neural networks, and the discriminator in the training process is composed of a binary classification deep neural network composed of M layers of deep neural networks;
生成器的输入为长度T的心电信号数据片段,输出为降噪后的信号;判别器在心电信号降噪模型训练时与生成器联合进行抗博弈学习。The input of the generator is an ECG signal data segment of length T, and the output is a denoised signal; the discriminator performs anti-game learning together with the generator during the training of the ECG signal denoising model.
作为进一步的方案,所述生成器的损失函数在条件生成对抗网络的损失函数的基础上,增加了降噪后的心电信号与原始心电信号之间的差距ldist,以及降噪后的心电信号与原始心电信号之间的最大局部误差lmax。As a further solution, the loss function of the generator adds the difference l dist between the denoised ECG signal and the original ECG signal, and the maximum local error l max between the denoised ECG signal and the original ECG signal, based on the loss function of the conditional generative adversarial network.
作为进一步的方案,心电信号降噪模型的降噪性能与计算复杂度比率指标具体为:模型测试时降噪后信号的平均信噪比,与模型处理单个数据样本所需时间的比值。As a further solution, the denoising performance and computational complexity ratio indicator of the ECG signal denoising model is specifically: the ratio of the average signal-to-noise ratio of the denoised signal during model testing to the time required for the model to process a single data sample.
作为进一步的方案,采用拟合函数预测样本长度T的最优值的具体过程如下:As a further solution, the specific process of using the fitting function to predict the optimal value of the sample length T is as follows:
按照样本长度T对心电信号进行分割,并进行数据归一化,获得对应样本长度T的数据集;The ECG signal is segmented according to the sample length T, and the data is normalized to obtain a data set corresponding to the sample length T;
设定多个不同样本长度T的取值,分别得到每一个样本长度T对应的数据集;Set multiple different sample lengths T to obtain data sets corresponding to each sample length T;
固定其他建模参量不变,根据设定的T值和对应的数据集,分别进行降噪模型的建模和训练;Other modeling parameters are fixed unchanged, and the denoising model is modeled and trained according to the set T value and the corresponding data set;
对于训练后的心电信号降噪模型,分别计算模型的降噪性能与计算复杂度比率指标;并依据不同样本长度T和对应的SNR-CC数据进行曲线拟合得到拟合函数f(T),再对拟合后的函数求导得到f′(T),并计算得出f′(T)=0时,T的取值T0。采用T0重新建立数据库,对降噪模型进行训练、测试和验证。若模型降噪后的信号平均信噪比不小于设定的期望值,则当前样本长度T0即为T的最优值。For the trained ECG signal denoising model, the denoising performance and computational complexity ratio of the model are calculated respectively; and the fitting function f(T) is obtained by curve fitting based on different sample lengths T and corresponding SNR-CC data, and then the fitting function is derived to obtain f′(T), and the value of T when f′(T)= 0 is calculated. The database is re-established using T 0 , and the denoising model is trained, tested and verified. If the average signal-to-noise ratio of the signal after denoising by the model is not less than the set expected value, the current sample length T 0 is the optimal value of T.
同时,建模参量优选时还可联合其他可能会影响模型部署的指标一起进行模型的优化,如:模型参量的占用内存大小等。结合上述建模参量优选过程,对模型参量的占用内存大小加以限制,可得到相应的分段拟合函数,再对其进行分析和预测,以确定在分段内的最优建模参量值。At the same time, when optimizing modeling parameters, other indicators that may affect model deployment can also be combined to optimize the model, such as the memory size occupied by model parameters, etc. Combined with the above modeling parameter optimization process, the memory size occupied by model parameters is restricted, and the corresponding piecewise fitting function can be obtained, which is then analyzed and predicted to determine the optimal modeling parameter value within the segment.
对于生成器中编码器和解码器的深度神经网络层数N,以及生成器中每层深度神经网络神经元数量Lx,分别采用与样本长度T值相同的优化方法,固定非待优化建模参量,选取相应待优化建模参量的最优值。For the number of deep neural network layers N of the encoder and decoder in the generator, and the number of neurons L x in each layer of the deep neural network in the generator, the same optimization method as the sample length T value is used, the non-optimized modeling parameters are fixed, and the optimal values of the corresponding modeling parameters to be optimized are selected.
对于判别器中深度神经网络的层数M,以及判别器中每层深度神经网络神经元数量Rx,分别采用与样本长度T值相同的优化方法,选取相应的最优值;For the number of layers M of the deep neural network in the discriminator and the number of neurons R x in each layer of the deep neural network in the discriminator, the same optimization method as the sample length T is used to select the corresponding optimal value;
在判别器中深度神经网络的优化参量优化完成后,联合已优化的生成器参量一起再构建降噪模型,进一步验证已优化参量的降噪效果,直到模型降噪性能与计算复杂度比率指标在最大值位置附近趋于平稳。After the optimization parameters of the deep neural network in the discriminator are optimized, the denoising model is rebuilt together with the optimized generator parameters to further verify the denoising effect of the optimized parameters until the model denoising performance and computational complexity ratio index stabilizes near the maximum value position.
选取不同种类的深度神经网络和条件生成对抗网络构建对应心电信号降噪模型;分别采用上述的优化方法进行对应降噪模型建模参量的优化,对比优化后对应模型的降噪性能以及降噪性能与计算复杂度比率指标,选取性能最优的降噪模型,作为最终的心电信号降噪模型。Different types of deep neural networks and conditional generative adversarial networks are selected to construct corresponding ECG signal denoising models. The above-mentioned optimization methods are used to optimize the modeling parameters of the corresponding denoising models, and the denoising performance of the optimized corresponding models and the ratio of denoising performance to computational complexity are compared. The denoising model with the best performance is selected as the final ECG signal denoising model.
在另一些实施方式中,采用如下技术方案:In other embodiments, the following technical solutions are adopted:
一种基于条件生成对抗网络的心电信号降噪优化系统,包括:An electrocardiogram signal denoising optimization system based on conditional generative adversarial network, comprising:
数据获取模块,用于获取心电信号,将所述心电信号按照样本长度T进行分割,对分割后的数据片段进行最大最小值归一化处理;A data acquisition module is used to acquire an ECG signal, segment the ECG signal according to a sample length T, and perform maximum and minimum value normalization processing on the segmented data segments;
降噪模块,用于将处理后的数据片段输入至训练好的心电信号降噪模型中,得到降噪后的心电信号;A denoising module, used to input the processed data segments into a trained ECG signal denoising model to obtain a denoised ECG signal;
其中,心电信号降噪模型基于深度神经网络和条件生成对抗网络构建,包括生成器和判别器;对于心电信号降噪模型中的每一种建模参量,按设定步长选取不同数值,分别建立建模参量在不同数值下的降噪模型并获得与其对应降噪性能与计算复杂度比率指标,采用拟合函数预测每一种建模参量的最优值,以使得模型的降噪性能和计算开销达到最优;Among them, the ECG signal denoising model is constructed based on deep neural networks and conditional generative adversarial networks, including generators and discriminators. For each modeling parameter in the ECG signal denoising model, different values are selected according to the set step size, and denoising models with different modeling parameters are established respectively, and the corresponding denoising performance and computational complexity ratio indicators are obtained. The fitting function is used to predict the optimal value of each modeling parameter, so that the denoising performance and computational overhead of the model are optimized.
所述建模参量包括:样本长度T、生成器中编码器和解码器的深度神经网络层数N及每一层的神经元数量Lx和判别器的深度神经网络层数M及每一层的神经元数量Rx中的至少一种或多种。The modeling parameters include: at least one or more of the sample length T, the number of deep neural network layers N of the encoder and decoder in the generator and the number of neurons Lx in each layer, and the number of deep neural network layers M of the discriminator and the number of neurons Rx in each layer.
在另一些实施方式中,采用如下技术方案:In other embodiments, the following technical solutions are adopted:
一种终端设备,其包括处理器和存储器,处理器用于实现指令;存储器用于存储多条指令,所述指令适于由处理器加载并执行上述的基于条件生成对抗网络的心电信号降噪优化方法。A terminal device comprises a processor and a memory, wherein the processor is used to implement instructions; the memory is used to store multiple instructions, wherein the instructions are suitable for being loaded by the processor and executing the above-mentioned electrocardiogram signal denoising optimization method based on conditional generative adversarial network.
与现有技术相比,本发明的有益效果是:Compared with the prior art, the present invention has the following beneficial effects:
综上,建立并提出一种目前缺少的能够同时兼顾计算复杂度和降噪性能的指标及根据该指标对模型的建模进行充分、合理、系统的优化方法,能够解决上述几种现有的降噪方法实施时所出现的建模时间长、降噪模型计算复杂度较高,运算开销大等问题。In summary, an indicator that is currently lacking and can take into account both computational complexity and denoising performance is established and proposed, as well as a method for fully, reasonably and systematically optimizing the model based on the indicator. This can solve the problems of long modeling time, high computational complexity of the denoising model, and high computational overhead that arise when implementing the above-mentioned existing denoising methods.
(1)本发明将深度神经网络与条件生成对抗网络进行结合,来实现心电信号的降噪优化;其中深度神经网络可以学习样本中的真实数据的分布,从而剔除噪声;条件生成对抗网络是一种改进的生成对抗网络,它将带噪信号作为条件与原本的输入信号(降噪后的信号或原始信号)相结合作为判别器的输入,这使得生成器和判别器能更好的进行对抗博弈,进而控制生成器的学习方向更接近于样本中的真实数据分布。(1) The present invention combines a deep neural network with a conditional generative adversarial network to achieve noise reduction optimization of ECG signals; wherein the deep neural network can learn the distribution of real data in the sample, thereby removing noise; the conditional generative adversarial network is an improved generative adversarial network, which combines the noisy signal as a condition with the original input signal (the denoised signal or the original signal) as the input of the discriminator, which enables the generator and the discriminator to better engage in adversarial game, thereby controlling the learning direction of the generator to be closer to the real data distribution in the sample.
(2)本发明提出降噪性能与计算复杂度比率指标,用来评价模型降噪性能和模型计算开销,并利用所提指标与相应建模参量实验数据进行曲线拟合,获得拟合函数并求导,进而预测建模参量的最优取值,再通过实验验证优选出的各个待优化建模参量的取值。该优化方法可利用少量的实验,实现模型参数的组合优化,使得模型在达到相同降噪性能的前提下,实现尽可能小的模型计算开销。(2) The present invention proposes a noise reduction performance and computational complexity ratio index to evaluate the model noise reduction performance and model computational overhead, and uses the proposed index to perform curve fitting with the experimental data of the corresponding modeling parameters to obtain the fitting function and derive it, thereby predicting the optimal value of the modeling parameter, and then verifying the value of each modeling parameter to be optimized through experiments. This optimization method can use a small amount of experiments to achieve combined optimization of model parameters, so that the model can achieve the smallest possible model computational overhead while achieving the same noise reduction performance.
(3)本发明降噪模型生成器的损失函数在常规CGAN网络的损失函数的基础上,还添加了降噪信号与纯净心电信号之间的距离以及最大局部误差。其中降噪信号与纯净心电信号之间的距离体现了降噪后的信号与纯净心电信号之间的整体差异,能够控制生成器生成数据的大方向;而降噪信号与纯净心电信号之间的最大局部误差体现了降噪后信号与纯净信号之间待改进的局部差异,即生成器生成数据时需要改进的信号细节。当生成器的损失函数中添加了以上两个部分就能够在整体方向和局部细节上对生成器的数据生成过程进行把控,进而获得更好的降噪效果。(3) The loss function of the noise reduction model generator of the present invention is based on the loss function of the conventional CGAN network, and also adds the distance between the noise reduction signal and the pure ECG signal and the maximum local error. The distance between the noise reduction signal and the pure ECG signal reflects the overall difference between the noise reduction signal and the pure ECG signal, and can control the general direction of the data generated by the generator; while the maximum local error between the noise reduction signal and the pure ECG signal reflects the local difference between the noise reduction signal and the pure signal that needs to be improved, that is, the signal details that need to be improved when the generator generates data. When the above two parts are added to the loss function of the generator, the data generation process of the generator can be controlled in terms of the overall direction and local details, thereby obtaining a better noise reduction effect.
本发明的其他特征和附加方面的优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本方面的实践了解到。Other features and advantages of additional aspects of the present invention will be given in part in the following description, and in part will become obvious from the following description, or will be learned through the practice of the present invention.
附图说明BRIEF DESCRIPTION OF THE DRAWINGS
图1为本发明实施例中的心电信号降噪模型优化过程示意图;FIG1 is a schematic diagram of an optimization process of an ECG signal denoising model in an embodiment of the present invention;
图2为本发明实施例中的心电信号降噪模型整体结构示意图;FIG2 is a schematic diagram of the overall structure of an ECG signal denoising model according to an embodiment of the present invention;
图3为本发明实施例中的建模参量T与SNR-CC的拟合函数曲线示意图;FIG3 is a schematic diagram of a fitting function curve of a modeling parameter T and SNR-CC in an embodiment of the present invention;
图4为本发明实施例中的基于条件生成对抗网络的心电信号降噪过程示意图。FIG4 is a schematic diagram of an electrocardiogram signal denoising process based on a conditional generative adversarial network in an embodiment of the present invention.
具体实施方式Detailed ways
应该指出,以下详细说明都是例示性的,旨在对本发明提供进一步的说明。除非另有指明,本发明使用的所有技术和科学术语具有与本发明所属技术领域的普通技术人员通常理解的相同含义。It should be noted that the following detailed descriptions are illustrative and intended to provide further explanation of the present invention. Unless otherwise specified, all technical and scientific terms used in the present invention have the same meanings as those commonly understood by those skilled in the art to which the present invention belongs.
需要注意的是,这里所使用的术语仅是为了描述具体实施方式,而非意图限制根据本发明的示例性实施方式。如在这里所使用的,除非上下文另外明确指出,否则单数形式也意图包括复数形式,此外,还应当理解的是,当在本说明书中使用术语“包含”和/或“包括”时,其指明存在特征、步骤、操作、器件、组件和/或它们的组合。It should be noted that the terms used herein are only for describing specific embodiments and are not intended to limit exemplary embodiments according to the present invention. As used herein, unless the context clearly indicates otherwise, the singular form is also intended to include the plural form. In addition, it should be understood that when the terms "comprising" and/or "including" are used in this specification, it indicates the presence of features, steps, operations, devices, components and/or combinations thereof.
实施例一Embodiment 1
在一个或多个实施方式中,公开了一种基于条件生成对抗网络的心电信号降噪优化方法,具体包括如下:In one or more embodiments, a method for optimizing electrocardiogram signal noise reduction based on a conditional generative adversarial network is disclosed, which specifically includes the following:
(1)获取心电信号,将心电信号按照样本长度T进行分割,对分割后的数据片段进行最大最小值归一化处理;(1) obtaining an ECG signal, segmenting the ECG signal according to a sample length T, and performing maximum and minimum value normalization processing on the segmented data segments;
(2)将处理后的数据片段输入至训练好的心电信号降噪模型中,得到降噪后的心电信号;(2) inputting the processed data segments into the trained ECG signal denoising model to obtain the denoised ECG signal;
其中,心电信号降噪模型基于深度神经网络和条件生成对抗网络构建;通过计算心电信号降噪模型的降噪性能与计算复杂度比率指标,通过曲线拟合预测最优的建模参量值,通过验证,以确定模型的降噪性能和计算开销达到最优。其中,建模参量可以是样本长度T、编码器和解码器的深度神经网络层数N及每一层的神经元数量Lx、判别器的深度神经网络层数M及每一层的神经元数量Rx中的一种或多种。The ECG signal denoising model is constructed based on a deep neural network and a conditional generative adversarial network. The denoising performance and computational complexity ratio of the ECG signal denoising model is calculated, and the optimal modeling parameter value is predicted by curve fitting. The modeling parameter is verified to determine whether the denoising performance and computational overhead of the model are optimal. The modeling parameter may be one or more of the sample length T, the number of deep neural network layers N of the encoder and decoder and the number of neurons Lx in each layer, the number of deep neural network layers M of the discriminator and the number of neurons Rx in each layer.
本实施例中,结合图2,心电信号降噪模型基于深度神经网络(如:双向长短时记忆网络、卷积神经网络、循环神经网络等)和条件生成对抗网络组合构建;心电信号降噪模型由生成器和判别器两部分构成,其中,生成器可由深度神经网络组成的降噪自编码器(DAE)所构成的,DAE中包含了N层深度神经网络组成的编码器以及N层深度神经网络组成的解码器,判别器由M层深度神经网络组成的二分类深度神经网络所构成。In this embodiment, in combination with Figure 2, the ECG signal denoising model is constructed based on a combination of a deep neural network (such as a bidirectional long short-term memory network, a convolutional neural network, a recurrent neural network, etc.) and a conditional generative adversarial network; the ECG signal denoising model consists of two parts: a generator and a discriminator. The generator can be composed of a denoising autoencoder (DAE) composed of a deep neural network, and the DAE includes an encoder composed of N layers of deep neural networks and a decoder composed of N layers of deep neural networks. The discriminator is composed of a binary classification deep neural network composed of M layers of deep neural networks.
生成器的输入为长度T的带噪信号输出为降噪后的信号/>判别器的输入为带噪信号与降噪信号的组合或带噪信号与原始信号的组合,输出为真或假。The input of the generator is a noisy signal of length T The output is the signal after noise reduction/> The input of the discriminator is a combination of the noisy signal and the denoised signal or a combination of the noisy signal and the original signal, and the output is true or false.
下面结合图1对心电信号降噪模型的训练和优化过程进行说明。The training and optimization process of the ECG signal denoising model is described below in conjunction with FIG1 .
a)信号加噪。a) Signal noise addition.
在公开或特定的心电数据库中选择较干净的心电信号作为原始信号(xn),从心电噪声数据库中选择多种常见的心电噪声作为噪声信号(en),原始信号和噪声信号的采样频率均为f。根据不同输入信噪比(Signal to Noise Ratio,SNR),将多种单一噪声和多种混合噪声分别与原始信号进行叠加处理,得到带噪信号 Select a relatively clean ECG signal from a public or specific ECG database as the original signal (x n ), and select a variety of common ECG noises from the ECG noise database as the noise signal (e n ). The sampling frequency of the original signal and the noise signal is f. According to different input signal-to-noise ratios (SNR), multiple single noises and multiple mixed noises are superimposed on the original signal to obtain a noisy signal.
b)模型待优化建模参量集设定。b) Setting of modeling parameter set to be optimized.
心电信号降噪模型的待优化建模参量集C包括:生成器建模参量(样本长度T、编码器和解码器深度神经网络层数N、每一层的神经元数量Lx)、判别器建模参量(判别器深度神经网络层数M、每一层的神经元数量Rx)。The set of modeling parameters C to be optimized for the ECG signal denoising model includes: generator modeling parameters (sample length T, number of encoder and decoder deep neural network layers N, number of neurons in each layer L x ), and discriminator modeling parameters (number of discriminator deep neural network layers M, number of neurons in each layer R x ).
首先,应根据以下规律预设上述降噪模型的待优化建模参量。First, the modeling parameters to be optimized of the above-mentioned denoising model should be preset according to the following rules.
样本长度T既决定了模型输入的信息量的大小,也决定了生成器和判别器的起始层神经元数量(即模型起始尺寸)。通常心率范围为50-150次/分钟,即单次心跳的周期为0.4-1.2s,在采样率为f时,为保证每个样本段中至少包含一个心跳周期,通常T应在≥1.2f的范围内选择多个值进行优选。The sample length T determines both the amount of information input to the model and the number of neurons in the starting layer of the generator and discriminator (i.e., the starting size of the model). The heart rate usually ranges from 50 to 150 beats per minute, i.e., the cycle of a single heartbeat is 0.4 to 1.2 seconds. When the sampling rate is f, in order to ensure that each sample segment contains at least one heartbeat cycle, T should usually be selected from multiple values within the range of ≥1.2f for optimization.
为了使降噪模型的生成器和判别器能够在训练时更好的对抗博弈,以达到更好的降噪效果。一般来说生成器与判别器所用的神经网络结构应该是相似的,且其中M应该大于等于N。即,应首先根据组成生成器的深度神经网络特点选定一个合适的N,然后通过选定M为大于等于N的值,并依据所选深度神经网络类型设置相应网络层的神经元数量(Lx、Rx)。In order to make the generator and discriminator of the denoising model better able to compete against each other during training, so as to achieve better denoising effect. Generally speaking, the neural network structures used by the generator and the discriminator should be similar, and M should be greater than or equal to N. That is, a suitable N should be selected first according to the characteristics of the deep neural network that constitutes the generator, and then M should be selected as a value greater than or equal to N, and the number of neurons (L x , R x ) of the corresponding network layer should be set according to the selected deep neural network type.
c)数据分割及归一化:按照样本长度(T)对带噪信号和原始信号(xn)进行分割,得到相应长度的心电信号数据段,并分别将其最大最小值归一化,进而组成相应的数据集(DT)。数据集(DT)中,占总数据C%的数据用于训练,1-C%的数据用于测试。c) Data segmentation and normalization: The noisy signal is segmented according to the sample length (T). The original signal (x n ) is segmented to obtain ECG signal data segments of corresponding lengths, and the maximum and minimum values are normalized to form the corresponding data set (D T ). In the data set (D T ), C% of the total data is used for training, and 1-C% of the data is used for testing.
其中,数据进行最大最小归一化的计算公式如下:Among them, the calculation formula for maximum and minimum normalization of data is as follows:
其中,xn为待归一化信号第n个采样点的值,xmin是待归一化信号采样点中的最小值,xmax是待归一化信号采样点中的最大值。Wherein, xn is the value of the nth sampling point of the signal to be normalized, xmin is the minimum value among the sampling points of the signal to be normalized, and xmax is the maximum value among the sampling points of the signal to be normalized.
d)模型训练。d) Model training.
利用准备好的数据集(DT),通过对抗博弈学习,训练构建好的心电信号降噪模型。当判别器不能分辩输入的是否为原始信号或降噪信号时,认为对抗学习达到纳什均衡,此时模型收敛。Using the prepared data set (D T ), the constructed ECG signal denoising model is trained through adversarial game learning. When the discriminator cannot distinguish whether the input is the original signal or the denoised signal, it is considered that the adversarial learning has reached a Nash equilibrium and the model has converged.
e)降噪模型测试:使用训练好的生成器,输入数据集中测试带噪信号,得到降噪后的信号,并计算其信噪比(SNR)、均方根误差(Root mean square error,RMSE)等降噪性能指标用于降噪效果评估。随后,通过统计此时降噪模型的参数量、参数所占内存、计算量等参数用于评估该模型的计算消耗(开销)。通过综合对比不同待优化建模参量取值下的模型降噪性能指标和模型计算开销,优选待优化建模参量的取值。e) Noise reduction model test: Use the trained generator to input the noisy signal in the data set to obtain the denoised signal, and calculate its signal-to-noise ratio (SNR), root mean square error (RMSE) and other noise reduction performance indicators for noise reduction effect evaluation. Subsequently, the number of parameters of the noise reduction model at this time, the memory occupied by the parameters, the amount of calculation and other parameters are statistically used to evaluate the computational consumption (overhead) of the model. By comprehensively comparing the model noise reduction performance indicators and model calculation overhead under different values of the modeling parameters to be optimized, the values of the modeling parameters to be optimized are preferred.
为了方便综合对比上述降噪性能和计算复杂度指标,本实施例中提出了降噪性能与计算复杂度比率(The ratio of Signal-to-Noise Ratio to ComputationalComplexity,SNR-CC)这一指标,其计算方法如下:In order to facilitate the comprehensive comparison of the above-mentioned noise reduction performance and computational complexity indicators, the ratio of noise reduction performance to computational complexity (SNR-CC) is proposed in this embodiment, and its calculation method is as follows:
其中,O为SNR-CC的值,代表了模型降噪性能与计算复杂度之间的比值,该指标越大,表示在降噪信号达到相同信噪比时所消耗的计算时间(计算开销)越小。SNRAverage是模型测试时降噪后信号的平均信噪比。ts是训练好的降噪模型(使用其中的生成器部分)计算单个样本所需要的时间,也在一定程度上反映了降噪模型的参数量、参数所占内存、计算量等指标的综合影响。当SNRAverage大于或等于期望值S0时,降噪模型的SNR-CC越大越好。其中,S0为所优化模型降噪性能即平均SNR的期望值,可通过参考其他先进降噪算法性能指标得到。Among them, O is the value of SNR-CC, which represents the ratio between the model denoising performance and the computational complexity. The larger the index is, the smaller the computational time (computational overhead) consumed when the denoised signal reaches the same signal-to-noise ratio. SNR Average is the average signal-to-noise ratio of the denoised signal during model testing. t s is the time required for the trained denoising model (using the generator part) to calculate a single sample, which also reflects the comprehensive influence of the parameters of the denoising model, the memory occupied by the parameters, the amount of computation, and other indicators to a certain extent. When SNR Average is greater than or equal to the expected value S 0 , the SNR-CC of the denoising model is as large as possible. Among them, S 0 is the expected value of the denoising performance of the optimized model, that is, the average SNR, which can be obtained by referring to the performance indicators of other advanced denoising algorithms.
优选算法模型架构时,首先根据建模参量p和SNR-CC对应数据进行曲线拟合,得到拟合函数f(p),图3给出了样本长度T的不同取值与其对应的SNR-CC取值的拟合函数示意图。对拟合函数求导得到f′(p)。然后通过计算得到当f′(p)=0时,p的取值p0。常见的三阶多项式拟合函数及其导数如下所示:When optimizing the algorithm model architecture, firstly, curve fitting is performed based on the modeling parameter p and the corresponding SNR-CC data to obtain the fitting function f(p). Figure 3 shows a schematic diagram of the fitting function for different values of sample length T and its corresponding SNR-CC value. The fitting function is derived to obtain f′(p). Then, when f′(p)=0, the value of p is obtained by calculation p 0. The common third-order polynomial fitting function and its derivative are shown below:
O=f(p)=ap3+bp2+cp+d(3)O=f(p)=ap 3 +bp 2 +cp+d(3)
O′=f′(p)=3ap2+2bp+c(4)O′=f′(p)=3ap 2 +2bp+c(4)
随后,根据建模参量p0重新进行建模并进行训练和测试,计算出此时的信噪比SNR值与此参量对应的SNR-CC值O0,验证此时模型的降噪性能和计算损耗。即,当f′(p)=0且SNRAverage≥S0时,认为是模型当前情况下该建模参量的优选取值。Then, the model is re-modeled and trained and tested according to the modeling parameter p 0 , and the signal-to-noise ratio SNR value at this time and the SNR-CC value O 0 corresponding to this parameter are calculated to verify the noise reduction performance and calculation loss of the model at this time. That is, when f′(p)=0 and SNR Average ≥S 0 , it is considered to be the optimal value of the modeling parameter under the current model situation.
此时,如果对降噪模型参数所占用的内存大小有所要求,就应根据其具体限制将所得拟合函数f(p)变换成为一个分段函数。At this point, if there are requirements on the memory size occupied by the denoising model parameters, the obtained fitting function f(p) should be transformed into a piecewise function according to its specific restrictions.
O=f(p)=ap3+bp2+cp+d(SNRAverage≥S0且MC≤MC0)(5)O = f(p) = ap 3 + bp 2 + cp + d (SNR Average ≥ S 0 and MC ≤ MC 0 ) (5)
其中MC为模型参数占用内存大小,MC0为所期望的模型参数占用内存的最大值。当SNRAverage≥S0且MC≤MC0时,使f′(p)能取到的最小值的pmin被认为是模型当前情况下该建模参量的优选取值。Where MC is the memory size occupied by the model parameters, and MC 0 is the maximum value of the memory occupied by the model parameters. When SNR Average ≥ S 0 and MC ≤ MC 0 , p min, which is the minimum value that f′(p) can take, is considered to be the optimal value of the modeling parameter under the current model situation.
此外,还应采用在待优化建模参量取优选值时所构建的降噪模型,对测试数据进行降噪处理,然后利用基于心电信号特征的心电分类算法对降噪前、降噪后、原始的心电信号进行分类性能评估,用于评估降噪模型保留具有医学价值信息的能力。In addition, the denoising model constructed when the modeling parameters to be optimized take the optimal values should be used to denoise the test data, and then the ECG classification algorithm based on the characteristics of the ECG signal should be used to evaluate the classification performance of the original ECG signals before and after denoising, in order to evaluate the ability of the denoising model to retain information of medical value.
f)生成器及判别器模型优化。f) Generator and discriminator model optimization.
降噪处理时只使用训练好的降噪模型中的生成器部分,因此,T、N是模型优化过程中两个较重要的参数。但由于M能够影响判别器性能进而响应整个训练过程中的博弈情况,所以M值也必须要进行优选。During denoising, only the generator part of the trained denoising model is used. Therefore, T and N are two important parameters in the model optimization process. However, since M can affect the performance of the discriminator and thus respond to the game situation during the entire training process, the M value must also be optimized.
经过分析,T、N、M的取值、生成器每个深度神经网络层的神经元数量(Lx)及生成器每个深度神经网络层的神经元数量(Rx),均会极大的影响模型的降噪效果。因此,需对上述模型建模参量进行优化。After analysis, the values of T, N, M, the number of neurons in each deep neural network layer of the generator (L x ) and the number of neurons in each deep neural network layer of the generator (R x ) will greatly affect the denoising effect of the model. Therefore, it is necessary to optimize the modeling parameters of the above models.
模型建模参量的一般优化过程如下:固定当前非优化建模参量不变,按一定步长分别建立待优化建模参量不同数值下(通常≥3种情况)的降噪模型,按照上述步骤b)-e),获得对应降噪性能与计算复杂度比率指标,利用拟合函数预测待优化建模参量的最优值,进而以最优值构建训练降噪模型,并测试降噪模型的综合性能加以验证。然后,按照顺序依次优化剩余的建模参量取值。The general optimization process of modeling parameters is as follows: fix the current non-optimized modeling parameters unchanged, establish denoising models with different values of the modeling parameters to be optimized (usually ≥3 cases) according to a certain step size, obtain the corresponding denoising performance and computational complexity ratio index according to the above steps b)-e), use the fitting function to predict the optimal value of the modeling parameter to be optimized, and then build a training denoising model with the optimal value, and test the comprehensive performance of the denoising model for verification. Then, optimize the remaining modeling parameter values in order.
由于模型最终通过生成器来实现降噪处理,因此应首先优化生成器模型。具体可采用如下方案:首先进行建模参量T的取值优选,并在优选完成后,采用相应的数据集(DT),通过循环上述的优选步骤,按照N、Lx顺序,在一定范围内依次优选生成器模型(即,N、Lx的取值)。Since the model ultimately uses the generator to achieve noise reduction, the generator model should be optimized first. Specifically, the following scheme can be adopted: first, the value of the modeling parameter T is optimized, and after the optimization is completed, the corresponding data set (D T ) is used to cycle the above optimization steps and optimize the generator model (i.e., the value of N and L x ) in sequence within a certain range.
然后,按照同样的方法再对判别器中待优化建模参量(M、Rx)的取值进行优化。由于N、M值能分别影响对抗博弈中生成器和判别器的性能,因此,两者还存在一定的互相影响关系。因此,首先在生成器优化中确定N的优选值,然后在判别器中深度神经网络的优化参量优化完成后,联合已优化的生成器参量一起再构建降噪模型,进一步验证已优化参量的降噪效果,直到模型降噪性能与计算复杂度比率指标在高位(最大值位置附近)趋于平稳。Then, the values of the modeling parameters (M, R x ) to be optimized in the discriminator are optimized in the same way. Since the values of N and M can affect the performance of the generator and the discriminator in the adversarial game respectively, there is a certain mutual influence between the two. Therefore, first determine the preferred value of N in the generator optimization, and then after the optimization parameters of the deep neural network in the discriminator are optimized, the denoising model is rebuilt together with the optimized generator parameters to further verify the denoising effect of the optimized parameters until the model denoising performance and computational complexity ratio index tends to be stable at a high level (near the maximum value position).
最终,得到优选的T、N、M、Lx、Rx的取值以及训练好的降噪算法优化模型。Finally, the optimal values of T, N, M, L x , and R x and the trained denoising algorithm optimization model are obtained.
心电信号降噪模型经过训练和优化后,在实施降噪处理时,通过该优化模型中的生成器部分进行降噪处理;具体过程如下:After the ECG signal denoising model is trained and optimized, when implementing the denoising process, the denoising process is performed through the generator part of the optimized model; the specific process is as follows:
将实际采集到心电信号数据按照优选的T值进行固定长度分割,然后对分割后的数据片段进行最大最小值归一化得到待处理的带噪数据。随后将其输入到训练好的降噪算法优化模型的生成器部分中进行计算,进而得到降噪后的信号。The actual collected ECG signal data is segmented into fixed lengths according to the preferred T value, and then the segmented data segments are normalized to the maximum and minimum values to obtain the noisy data to be processed. Subsequently, it is input into the generator part of the trained denoising algorithm optimization model for calculation, and then the denoised signal is obtained.
实施例二Embodiment 2
本实施例的心电降噪模型由双向长短时记忆网络(bidirectional long short-term memory,BiLSTM)和条件生成对抗网络CGAN构建,模型采用CGAN的学习架构,由双向长短时记忆网络组成的生成器和判别器构成。The ECG denoising model of this embodiment is constructed by a bidirectional long short-term memory network (BiLSTM) and a conditional generative adversarial network (CGAN). The model adopts the learning architecture of CGAN and consists of a generator and a discriminator composed of a bidirectional long short-term memory network.
模型初始结构如图2所示,其中,生成器是由2层(N)BiLSTM组成的编码器、2层(N)BiLSTM组成的解码器以及一层全连接层组成的降噪自编码器所构成的,判别器是由3层(M)BiLSTM、一层全连接层以及激活函数Sigmoid组成的二分类深度学习网络所构成的。生成器和判别器输入样本长度为T,生成器中编码器和解码器的神经元数量均为Lx={400,200},判别器中神经元数量为Rx={400,200,100}。激活函数Sigmoid的计算公式如下:The initial structure of the model is shown in Figure 2, where the generator is composed of an encoder consisting of 2 layers (N) BiLSTM, a decoder consisting of 2 layers (N) BiLSTM, and a denoising autoencoder consisting of a fully connected layer. The discriminator is composed of a binary classification deep learning network consisting of 3 layers (M) BiLSTM, a fully connected layer, and an activation function Sigmoid. The input sample length of the generator and discriminator is T. The number of neurons in the encoder and decoder in the generator is L x = {400, 200}, and the number of neurons in the discriminator is R x = {400, 200, 100}. The calculation formula of the activation function Sigmoid is as follows:
生成器的输入为带噪信号,输出为降噪后的信号;判别器的输入为带噪信号与降噪后的信号的组合或带噪信号与原始信号的组合,输出为真或假。The input of the generator is a noisy signal, and the output is a denoised signal; the input of the discriminator is a combination of a noisy signal and a denoised signal or a combination of a noisy signal and an original signal, and the output is true or false.
本实施例中,心电降噪模型生成器的损失函数在常规CGAN网络的损失函数的基础上,还添加了降噪后的信号与纯净心电信号之间的差距ldist以及最大局部误差lmax。In this embodiment, the loss function of the ECG denoising model generator is based on the loss function of the conventional CGAN network, and further adds the difference l dist between the denoised signal and the pure ECG signal and the maximum local error l max .
生成器的损失函数具体为:The loss function of the generator is specifically:
其中,为带噪信号,服从带噪信号分布/> 是生成器输入为带噪信号/>时的输出,/>是判别器输入为/>和/>时的输出,λ1和λ2是分别是ldist和lmax的权重系数,分别取0.7和0.2。in, is a noisy signal, obeying the noisy signal distribution/> The generator input is a noisy signal/> Output when, /> The discriminator input is /> and/> is the output when , λ 1 and λ 2 are the weight coefficients of l dist and l max , which are 0.7 and 0.2 respectively.
降噪后的信号与纯净心电信号之间的差距ldist的计算公式如下:The calculation formula of the difference l dist between the denoised signal and the pure ECG signal is as follows:
其中,为降噪后的信号,xn为原始信号,T为样本长度。in, is the signal after noise reduction, xn is the original signal, and T is the sample length.
降噪后的信号与纯净心电信号之间最大局部误差lmax的计算公式如下:The calculation formula of the maximum local error l max between the denoised signal and the pure ECG signal is as follows:
此时模型的整体损失函数为:At this time, the overall loss function of the model is:
其中,x为原始信号,服从原始信号分布pdata(x);为带噪信号,服从带噪信号分布 Where x is the original signal, which obeys the original signal distribution p data (x); is a noisy signal, which obeys the distribution of noisy signals
模型采用Pytorch编程实现,采用RMSProp(Root Mean Square Prop均方根传递)优化器进行参数更新,学习率设置为0.0001。模型训练和测试运行在Dell T640服务器上(NVIDIA GTX 3090 24GB)。模型训练及优化流程,与实施例一中相同,具体实施过程如下:The model is implemented using Pytorch programming, and the RMSProp (Root Mean Square Prop) optimizer is used for parameter update, and the learning rate is set to 0.0001. Model training and testing are run on a Dell T640 server (NVIDIA GTX 3090 24GB). The model training and optimization process is the same as in Example 1, and the specific implementation process is as follows:
a)在MIT-BIH心率失常数据库中选择标号为103、105、111、116、122、205、213、219、223、230的心电信号作为原始信号(xn)。从心电噪声数据库中选择肌肉伪影、电极运动伪影、基线漂移三种常见的心电噪声作为噪声信号(en)。根据将输入SNR分别设置为0dB、1dB、2dB、3dB、4dB、5dB,采用三种单一噪声、三种两两叠加的混合噪声以及一种三种噪声叠加的混合噪声分别与原始信号进行叠加处理,得到带有不同输入SNR及不同种类噪声的带噪信号上述心电信号的采样率为360Hz,因此样本长度T应≥432个采样点。a) Select ECG signals numbered 103, 105, 111, 116, 122, 205, 213, 219, 223, and 230 from the MIT-BIH arrhythmia database as the original signals (x n ). Select three common ECG noises, muscle artifacts, electrode motion artifacts, and baseline drift, from the ECG noise database as noise signals (e n ). According to setting the input SNR to 0dB, 1dB, 2dB, 3dB, 4dB, and 5dB, three single noises, three mixed noises superimposed on each other, and one mixed noise superimposed on three noises are respectively superimposed on the original signal to obtain noisy signals with different input SNRs and different types of noises. The sampling rate of the above ECG signal is 360 Hz, so the sample length T should be ≥ 432 sampling points.
b)依据优选策略,首先进行预先设定待优化建模参量集C。根据BiLSTM网络的数据特征提取能力,暂定T=512、N=2、M=3、生成器的Lx={400,200}、判别器Rx={400,200,100}。b) According to the optimization strategy, the modeling parameter set C to be optimized is first pre-set. According to the data feature extraction capability of the BiLSTM network, T=512, N=2, M=3, Lx ={400, 200} of the generator, and Rx ={400, 200, 100} of the discriminator are tentatively set.
c)将上述带有不同输入SNR和噪声的带噪信号和原始信号按照对应的T值(此时为512),进行分割,得到不同长度的心电信号数据段,并分别将其最大最小值归一化,进而组成相应的数据集(D512)。数据集其中80%数据用于训练,20%的数据用于测试。c) The noisy signals and original signals with different input SNRs and noises are segmented according to the corresponding T value (512 in this case) to obtain ECG signal data segments of different lengths, and their maximum and minimum values are normalized respectively to form a corresponding data set (D 512 ). 80% of the data in the data set is used for training and 20% of the data is used for testing.
d)对本实施例的心电降噪模型进行训练,通过对抗博弈学习过程,采用不同的数据集进行模型的训练,不断优化生成器和判别器。当判别器不能分辩输入的是否为原始信号或降噪信号时,或者生成的降噪后信号的SNR足够高或RMSE足够低时,模型收敛,训练完毕。d) The ECG denoising model of this embodiment is trained by using different data sets through the adversarial game learning process to continuously optimize the generator and the discriminator. When the discriminator cannot distinguish whether the input is the original signal or the denoised signal, or when the SNR of the generated denoised signal is high enough or the RMSE is low enough, the model converges and the training is completed.
e)此时,经过统计和计算以该建模参量(T=512)构建的降噪模型性能如下:SNR=28.71dB,SNR-CC=97.32dB/ms。设定模型期望降噪性能指标S0=28dB,此时模型降噪性能SNR>S0。e) At this time, after statistics and calculation, the performance of the noise reduction model constructed with the modeling parameters (T=512) is as follows: SNR=28.71dB, SNR-CC=97.32dB/ms. The model expected noise reduction performance index S 0 =28dB is set, and the model noise reduction performance SNR>S 0 .
f)生成器模型优化:首先,以T值的优选为例,其过程如下:设定多个样本长度(T)值(如:512、1024、2048、4096等),获得多种样本长度的数据集(DT),固定其他建模参数构建降噪模型,采用数据集(DT)训练该降噪模型。再利用训练好的生成器,输入相应的数据集(DT)中用于测试带噪信号,得到降噪后的信号,计算其信噪比(SNR)、SNR-CC降噪性能和计算开销指标,并依据T值及SNR-CC值进行了曲线拟合,拟合结果如图3所示,当f′(T)=0时,拟合结果T0=2127,为便于计算,这里可取T0=2048。采用T0重新建模后,综合对比各项指标,模型降噪性能SNR>S0达到了较好的水平(SNR=36.21dB),SNR-CC=109.74dB/ms,验证了其计算开销较小。随后,采用基于心电信号特征的心电分类算法对降噪前、降噪后、原始的心电信号进行分类性能对比,评估降噪模型保留具有医学价值信息的能力。f) Generator model optimization: First, taking the optimization of T value as an example, the process is as follows: set multiple sample length (T) values (such as 512, 1024, 2048, 4096, etc.), obtain data sets ( DT ) of various sample lengths, fix other modeling parameters to build a denoising model, and use the data set ( DT ) to train the denoising model. Then use the trained generator to input the corresponding data set ( DT ) for testing the noisy signal to obtain the denoised signal, calculate its signal-to-noise ratio (SNR), SNR-CC denoising performance and computational overhead indicators, and perform curve fitting based on T value and SNR-CC value. The fitting results are shown in Figure 3. When f′(T)=0, the fitting result T 0 =2127. For the convenience of calculation, T 0 =2048 can be taken here. After remodeling with T 0 , comprehensive comparison of various indicators showed that the model denoising performance SNR>S 0 reached a good level (SNR=36.21dB), SNR-CC=109.74dB/ms, which verified that its computational overhead was small. Subsequently, an ECG classification algorithm based on ECG signal features was used to compare the classification performance of the original ECG signals before and after denoising to evaluate the ability of the denoising model to retain medically valuable information.
重复上述步骤,通过综合对比降噪性能指标、保留医学信息的能力、模型的计算开销,即计算和分析SNR-CC、建模参量与SNR-CC拟合函数及其导数,进而预测和优化生成器中的其他建模参量(N和Lx)。先后依次确定N和Lx的优化取值。Repeat the above steps, and through comprehensive comparison of noise reduction performance indicators, the ability to retain medical information, and the computational overhead of the model, that is, calculating and analyzing SNR-CC, modeling parameters, and SNR-CC fitting functions and their derivatives, predict and optimize other modeling parameters (N and L x ) in the generator. Determine the optimal values of N and L x in turn.
随后,重复上述步骤,再通过固定生成器的待优化建模参量的取值,依次对M和Rx进行优选,即优化判别器模型。其中,确定M的优选值后,再重复上述步骤进行N值的再次优化,以确定N的有效性。Then, repeat the above steps, and then optimize M and R x in turn by fixing the values of the modeling parameters to be optimized of the generator, that is, optimize the discriminator model. After determining the optimal value of M, repeat the above steps to optimize the value of N again to determine the effectiveness of N.
结果显示,当N的取值范围为{1,2,3,4}和M的取值范围为{3,4,5,6}时,得出当N=3、M=4、生成器的Lx={800,400,200}、判别器Rx={800,400,200,100}时,其降噪后的心电信号平均SNR=45.12dB,模型SNR-CC=131.54dB/ms,此时该模型计算消耗比其他近似降噪效果的算法要低,达到计算消耗与降噪效果较好的平衡,即此时模型的SNR-CC值可达高位平稳。此外,还应采用基于心电信号特征的心电分类算法(支持向量机)对降噪前、降噪后和原始的心电信号进行了异常心电的四分类性能评估,采用所提模型降噪后的信号的分类准确率接近原始信号的分类准确率(≥95%),证实了降噪后信号保留了大量的有医学价值信息。The results show that when the value range of N is {1, 2, 3, 4} and the value range of M is {3, 4, 5, 6}, it is concluded that when N=3, M=4, the generator's Lx ={800, 400, 200}, and the discriminator's Rx ={800, 400, 200, 100}, the average SNR of the denoised ECG signal is 45.12dB, and the model SNR-CC is 131.54dB/ms. At this time, the computational consumption of the model is lower than that of other algorithms with approximate denoising effects, achieving a good balance between computational consumption and denoising effect, that is, the SNR-CC value of the model can reach a high and stable level. In addition, an ECG classification algorithm (support vector machine) based on ECG signal features was used to evaluate the four-category performance of abnormal ECG before, after and after denoising. The classification accuracy of the denoised signal using the proposed model was close to that of the original signal (≥95%), which confirmed that the denoised signal retained a large amount of medically valuable information.
结合图4,对于训练后的心电降噪模型,采集得到心电信号数据,将其按照T=2048的固定长度进行分割,然后对分割后的数据片段进行最大最小值归一化得到待处理的带噪心电信号片段。随后将其输入到训练优化完成的心电降噪模型中的生成器中进行降噪处理,就可得到降噪后的信号。Combined with Figure 4, for the trained ECG denoising model, the ECG signal data is collected and segmented according to a fixed length of T = 2048, and then the segmented data segments are normalized to obtain the noisy ECG signal segments to be processed. Then, the segments are input into the generator in the trained and optimized ECG denoising model for denoising, and the denoised signal can be obtained.
实施例三Embodiment 3
本实施例中,分别构建两种降噪模型:In this embodiment, two noise reduction models are constructed respectively:
第①种:基于卷积神经网络与条件生成对抗网络构建的降噪模型;The first type: a denoising model based on convolutional neural networks and conditional generative adversarial networks;
第②种:基于双向长短时记忆网络与条件生成对抗网络构建的降噪模型。The second type: a denoising model based on a bidirectional long short-term memory network and a conditional generative adversarial network.
首先,基于相同的硬件计算平台或系统,采用实施例一中所提方法对两种降噪模型分别进行优化,分别得到不同深度神经网络组成的优化后的降噪模型。然后,对比优化后的两种降噪模型的降噪性能及SNR-CC等指标,确定两者中更优的降噪模型。First, based on the same hardware computing platform or system, the two denoising models are optimized using the method proposed in Example 1 to obtain optimized denoising models composed of different deep neural networks. Then, the denoising performance and SNR-CC and other indicators of the two optimized denoising models are compared to determine the better denoising model.
第①种降噪模型,经过优化后在其最佳建模参量下构建的模型,经过测试,其SNR=41.03dB,SNR-CC=2.28dB/ms。The first noise reduction model is constructed under the optimal modeling parameters after optimization. After testing, its SNR=41.03dB and SNR-CC=2.28dB/ms.
第②种降噪模型,经过优化后在其最佳建模参量下构建的模型,经过测试,其SNR=45.12dB,模型SNR-CC=131.54dB/ms。The second noise reduction model is constructed under the optimal modeling parameters after optimization. After testing, its SNR=45.12dB and model SNR-CC=131.54dB/ms.
由此可见,基于双向长短时记忆网络与条件生成对抗网络构建的降噪模型,其降噪性能、降噪性价比均优于基于卷积神经网络与条件生成对抗网络构建的降噪模型。It can be seen that the denoising model constructed based on the bidirectional long short-term memory network and the conditional generative adversarial network has better denoising performance and denoising cost-effectiveness than the denoising model constructed based on the convolutional neural network and the conditional generative adversarial network.
实施例四Embodiment 4
在一个或多个实施方式中,公开了一种基于条件生成对抗网络的心电信号降噪优化系统,包括:In one or more embodiments, a system for optimizing electrocardiogram signal noise reduction based on a conditional generative adversarial network is disclosed, comprising:
数据获取模块,用于获取心电信号,将所述心电信号按照长度T进行分割,对分割后的数据片段进行最大最小值归一化处理;A data acquisition module is used to acquire an ECG signal, segment the ECG signal according to a length T, and perform maximum and minimum value normalization processing on the segmented data segments;
降噪模块,用于将处理后的数据片段输入至训练好的心电信号降噪模型中,得到降噪后的心电信号;A denoising module, used to input the processed data segments into a trained ECG signal denoising model to obtain a denoised ECG signal;
其中,心电信号降噪模型基于深度神经网络和条件生成对抗网络构建,包括生成器和判别器;对于心电信号降噪模型中的每一种建模参量,按设定步长选取不同数值,分别建立建模参量在不同数值下的降噪模型并获得与其对应降噪性能与计算复杂度比率指标,采用拟合函数预测每一种建模参量的最优值,以使得模型的降噪性能和计算开销达到最优;Among them, the ECG signal denoising model is constructed based on deep neural networks and conditional generative adversarial networks, including generators and discriminators. For each modeling parameter in the ECG signal denoising model, different values are selected according to the set step size, and denoising models with different modeling parameters are established respectively, and the corresponding denoising performance and computational complexity ratio indicators are obtained. The fitting function is used to predict the optimal value of each modeling parameter, so that the denoising performance and computational overhead of the model are optimized.
所述建模参量包括:样本长度T、生成器中编码器和解码器的深度神经网络层数N及每一层的神经元数量Lx和判别器的深度神经网络层数M及每一层的神经元数量Rx中的至少一种或多种。The modeling parameters include: at least one or more of the sample length T, the number of deep neural network layers N of the encoder and decoder in the generator and the number of neurons Lx in each layer, and the number of deep neural network layers M of the discriminator and the number of neurons Rx in each layer.
需要说明的是,上述各模块的具体实现方式以及在实施例一中进行了详细的说明,此处不再详述。It should be noted that the specific implementation methods of the above modules have been described in detail in Example 1 and will not be described in detail here.
实施例五Embodiment 5
在一个或多个实施方式中,公开了一种终端设备,包括服务器,所述服务器包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现实施例一中的基于条件生成对抗网络的心电信号降噪优化方法。为了简洁,在此不再赘述。In one or more embodiments, a terminal device is disclosed, including a server, the server including a memory, a processor, and a computer program stored in the memory and executable on the processor, and the processor implements the electrocardiogram signal denoising optimization method based on a conditional generative adversarial network in Embodiment 1 when executing the program. For the sake of brevity, it will not be described in detail here.
应理解,本实施例中,处理器可以是中央处理单元CPU,处理器还可以是其他通用处理器、数字信号处理器DSP、专用集成电路ASIC,现成可编程门阵列FPGA或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。It should be understood that in this embodiment, the processor may be a central processing unit CPU, and the processor may also be other general-purpose processors, digital signal processors DSP, application-specific integrated circuits ASIC, off-the-shelf programmable gate arrays FPGA or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. The general-purpose processor may be a microprocessor or the processor may also be any conventional processor, etc.
存储器可以包括只读存储器和随机存取存储器,并向处理器提供指令和数据、存储器的一部分还可以包括非易失性随机存储器。例如,存储器还可以存储设备类型的信息。The memory may include a read-only memory and a random access memory, and provide instructions and data to the processor. A portion of the memory may also include a non-volatile random access memory. For example, the memory may also store information about the device type.
在实现过程中,上述方法的各步骤可以通过处理器中的硬件的集成逻辑电路或者软件形式的指令完成。In the implementation process, each step of the above method can be completed by an integrated logic circuit of hardware in a processor or an instruction in the form of software.
上述虽然结合附图对本发明的具体实施方式进行了描述,但并非对本发明保护范围的限制,所属领域技术人员应该明白,在本发明的技术方案的基础上,本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。Although the above describes the specific implementation mode of the present invention in conjunction with the accompanying drawings, it is not intended to limit the scope of protection of the present invention. Technical personnel in the relevant field should understand that various modifications or variations that can be made by technical personnel in the field without creative work on the basis of the technical solution of the present invention are still within the scope of protection of the present invention.
Claims (9)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310225141.2A CN116250844B (en) | 2023-03-03 | 2023-03-03 | Electrocardiosignal noise reduction optimization method and system based on condition generation countermeasure network |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310225141.2A CN116250844B (en) | 2023-03-03 | 2023-03-03 | Electrocardiosignal noise reduction optimization method and system based on condition generation countermeasure network |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116250844A CN116250844A (en) | 2023-06-13 |
CN116250844B true CN116250844B (en) | 2024-04-26 |
Family
ID=86682353
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310225141.2A Active CN116250844B (en) | 2023-03-03 | 2023-03-03 | Electrocardiosignal noise reduction optimization method and system based on condition generation countermeasure network |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116250844B (en) |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110879254A (en) * | 2018-09-05 | 2020-03-13 | 哈尔滨工业大学 | An Acoustic Emission Signal Detection Method for Rail Cracks Based on Improved Least Square Generative Adversarial Network |
CN111067507A (en) * | 2019-12-26 | 2020-04-28 | 常熟理工学院 | ECG Signal Denoising Method Based on Generative Adversarial Network and Policy Gradient |
CN111564160A (en) * | 2020-04-21 | 2020-08-21 | 重庆邮电大学 | Voice noise reduction method based on AEWGAN |
WO2020248492A1 (en) * | 2019-06-14 | 2020-12-17 | 平安科技(深圳)有限公司 | Method and device for denoising oct image based on cyclic generative adversarial network |
CN112487914A (en) * | 2020-11-25 | 2021-03-12 | 山东省人工智能研究院 | ECG noise reduction method based on deep convolution generation countermeasure network |
CN113963213A (en) * | 2021-10-27 | 2022-01-21 | 上海交通大学 | Method and system for removing antagonistic noise aiming at antagonistic sample of deep neural network |
CN114129171A (en) * | 2021-12-01 | 2022-03-04 | 山东省人工智能研究院 | Electrocardiosignal noise reduction method based on improved residual error dense network |
CN114418886A (en) * | 2022-01-19 | 2022-04-29 | 电子科技大学 | Robustness denoising method based on deep convolution self-encoder |
CN114428234A (en) * | 2021-12-23 | 2022-05-03 | 西安电子科技大学 | Recognition method of radar high-resolution range image noise reduction based on GAN and self-attention |
CN114966860A (en) * | 2022-05-13 | 2022-08-30 | 重庆科技学院 | Seismic data denoising method based on convolutional neural network |
CN115442191A (en) * | 2022-11-08 | 2022-12-06 | 中国人民解放军战略支援部队航天工程大学 | Communication signal noise reduction method and system based on relative average generative adversarial network |
CN115429285A (en) * | 2022-09-26 | 2022-12-06 | 哈尔滨理工大学 | Electrocardiosignal denoising method based on generation countermeasure network |
WO2023279366A1 (en) * | 2021-07-09 | 2023-01-12 | Oppo广东移动通信有限公司 | Noise reduction method based on transfer learning, terminal device, network device and storage medium |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10531415B2 (en) * | 2018-03-02 | 2020-01-07 | DeepSig Inc. | Learning communication systems using channel approximation |
-
2023
- 2023-03-03 CN CN202310225141.2A patent/CN116250844B/en active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110879254A (en) * | 2018-09-05 | 2020-03-13 | 哈尔滨工业大学 | An Acoustic Emission Signal Detection Method for Rail Cracks Based on Improved Least Square Generative Adversarial Network |
WO2020248492A1 (en) * | 2019-06-14 | 2020-12-17 | 平安科技(深圳)有限公司 | Method and device for denoising oct image based on cyclic generative adversarial network |
CN111067507A (en) * | 2019-12-26 | 2020-04-28 | 常熟理工学院 | ECG Signal Denoising Method Based on Generative Adversarial Network and Policy Gradient |
CN111564160A (en) * | 2020-04-21 | 2020-08-21 | 重庆邮电大学 | Voice noise reduction method based on AEWGAN |
CN112487914A (en) * | 2020-11-25 | 2021-03-12 | 山东省人工智能研究院 | ECG noise reduction method based on deep convolution generation countermeasure network |
WO2023279366A1 (en) * | 2021-07-09 | 2023-01-12 | Oppo广东移动通信有限公司 | Noise reduction method based on transfer learning, terminal device, network device and storage medium |
CN113963213A (en) * | 2021-10-27 | 2022-01-21 | 上海交通大学 | Method and system for removing antagonistic noise aiming at antagonistic sample of deep neural network |
CN114129171A (en) * | 2021-12-01 | 2022-03-04 | 山东省人工智能研究院 | Electrocardiosignal noise reduction method based on improved residual error dense network |
CN114428234A (en) * | 2021-12-23 | 2022-05-03 | 西安电子科技大学 | Recognition method of radar high-resolution range image noise reduction based on GAN and self-attention |
CN114418886A (en) * | 2022-01-19 | 2022-04-29 | 电子科技大学 | Robustness denoising method based on deep convolution self-encoder |
CN114966860A (en) * | 2022-05-13 | 2022-08-30 | 重庆科技学院 | Seismic data denoising method based on convolutional neural network |
CN115429285A (en) * | 2022-09-26 | 2022-12-06 | 哈尔滨理工大学 | Electrocardiosignal denoising method based on generation countermeasure network |
CN115442191A (en) * | 2022-11-08 | 2022-12-06 | 中国人民解放军战略支援部队航天工程大学 | Communication signal noise reduction method and system based on relative average generative adversarial network |
Non-Patent Citations (2)
Title |
---|
An ECG Signal Denoising Method Using Conditional Generative Adversarial Net;Xiaoyu Wang,Bingchu Chen,等;《IEEE JOURNAL OF BIOMEDICAL AND HEALTH INFORMATICS》;20220731;第26卷(第7期);全文 * |
基于生成对抗网络的心电信号降噪算法及性能分析;王晓宇;《中国优秀硕士学位论文全文数据库(医药卫生科技辑)》;20230215;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN116250844A (en) | 2023-06-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111163690B (en) | Arrhythmia detection method, device, electronic equipment and computer storage medium | |
CN109948647B (en) | Electrocardiogram classification method and system based on depth residual error network | |
CN110193127B (en) | Music sleep assisting method and device, computer equipment and storage medium | |
CN109431492B (en) | Simulation and reconstruction method of ECG lead data based on neural network algorithm | |
CN111067507B (en) | Electrocardiosignal denoising method based on generation of countermeasure network and strategy gradient | |
CN115363586B (en) | A psychological stress level assessment system and method based on pulse wave signal | |
CN116662742B (en) | Brain electrolysis code method based on hidden Markov model and mask empirical mode decomposition | |
CN115470832A (en) | A blockchain-based ECG signal data processing method | |
CN115281685A (en) | Sleep stage identification method, device and computer-readable storage medium based on abnormal detection | |
CN116172522A (en) | An Anesthesia Depth Monitoring Method Based on Neural Network | |
CN116250844B (en) | Electrocardiosignal noise reduction optimization method and system based on condition generation countermeasure network | |
CN118762826B (en) | Method for predicting coronary artery stenosis and FFR based on big data of coronary artery diagnosis and treatment | |
CN118366225B (en) | Hand rehabilitation exercise monitoring method and system based on image processing | |
CN118918638A (en) | Burn patient rehabilitation training intelligent guidance system based on deep learning | |
CN118557187A (en) | Blood glucose concentration prediction method based on ELM-BiLSTM-attention | |
EP3770921A1 (en) | Method and system for pressure autoregulation based synthesizing of photoplethysmogram signal | |
Yu et al. | Work-in-progress: On the feasibility of lightweight scheme of real-time atrial fibrillation detection using deep learning | |
CN115546156A (en) | Thyroid ultrasound abnormal region capturing method and system and electronic equipment | |
CN112434630B (en) | Training method of continuous motion information prediction model and training device and equipment thereof | |
CN113538451A (en) | Magnetic resonance image segmentation method and device for deep venous thrombosis, electronic equipment and storage medium | |
Grzywalski et al. | Interactive Lungs Auscultation with Reinforcement Learning Agent | |
KR102751384B1 (en) | Fundus image processing apparatus and method for acquiring a normal fundus image based on the target image of the fundus | |
Van Gorp et al. | Aleatoric Uncertainty Estimation of Overnight Sleep Statistics Through Posterior Sampling Using Conditional Normalizing Flows | |
Duan et al. | Uncertainty Detection and Reduction in Neural Decoding of EEG Signals | |
CN116913513A (en) | Nine-channel heart sound and Bi-LSTM based non-invasive coronary heart disease prediction method and device |
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 |