CN110811609A - Intelligent epileptic spike detection method based on fusion of adaptive template matching and machine learning algorithm - Google Patents

Intelligent epileptic spike detection method based on fusion of adaptive template matching and machine learning algorithm Download PDF

Info

Publication number
CN110811609A
CN110811609A CN201911033704.8A CN201911033704A CN110811609A CN 110811609 A CN110811609 A CN 110811609A CN 201911033704 A CN201911033704 A CN 201911033704A CN 110811609 A CN110811609 A CN 110811609A
Authority
CN
China
Prior art keywords
spike
electroencephalogram
template matching
machine learning
epileptic
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201911033704.8A
Other languages
Chinese (zh)
Other versions
CN110811609B (en
Inventor
王紫萌
吴端坡
冯维
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hangzhou Electronic Science and Technology University
Original Assignee
Hangzhou Electronic Science and Technology University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hangzhou Electronic Science and Technology University filed Critical Hangzhou Electronic Science and Technology University
Priority to CN201911033704.8A priority Critical patent/CN110811609B/en
Publication of CN110811609A publication Critical patent/CN110811609A/en
Application granted granted Critical
Publication of CN110811609B publication Critical patent/CN110811609B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/369Electroencephalography [EEG]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/725Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7264Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Biophysics (AREA)
  • General Health & Medical Sciences (AREA)
  • Veterinary Medicine (AREA)
  • Public Health (AREA)
  • Animal Behavior & Ethology (AREA)
  • Psychiatry (AREA)
  • Surgery (AREA)
  • Molecular Biology (AREA)
  • Artificial Intelligence (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Signal Processing (AREA)
  • Physiology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Computation (AREA)
  • Mathematical Physics (AREA)
  • Fuzzy Systems (AREA)
  • Psychology (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Abstract

The invention provides an intelligent epileptic spike detection method based on self-adaptive template matching and machine learning algorithm fusion, which comprises the following steps: (1) electroencephalogram (EEG) acquisition: selecting an experimental object, acquiring electroencephalogram data of an epileptic by using electroencephalogram acquisition equipment, and establishing an experimental database; (2) data preprocessing: performing band-pass filtering on the acquired original EEG data to obtain a standard EEG signal; (3) carrying out adaptive template matching spike detection: (4) the spike detection method based on machine learning comprises the following steps: firstly, segmenting a brain electrical signal into electroencephalogram segments with the length of 1s, then extracting time domain and frequency domain characteristics in each electroencephalogram segment, and constructing spike characteristic vectors; and training a random forest classification model by using the feature vectors to obtain a spike detection result based on machine learning. (5) And (3) fusing detection results: the detection methods of step S3 and step S4 are fused, and if detected as a spike by both S3 and S4, it is regarded as an epileptic spike.

Description

Intelligent epileptic spike detection method based on fusion of adaptive template matching and machine learning algorithm
Technical Field
The invention relates to the field of computers, in particular to an intelligent epileptic spike detection method based on self-adaptive template matching and machine learning algorithm fusion.
Background
Epilepsy is a chronic disease in which sudden abnormal discharges in cerebral neurons lead to transient cerebral dysfunction. Over six thousand five million people worldwide suffer from epilepsy, with about nine million people in epilepsy in china. Seizures are paroxysmal, repetitive and unpredictable and may occur at any age.
Electroencephalograms (EEG) are potential signals generated by the discharge of neurons in the brain, reflect the rhythmic activity of bioelectricity in the brain, and include a large amount of physiological and disease information. In clinical medicine, EEG signal processing can not only provide a diagnosis basis for some brain diseases, but also provide an effective treatment means for some brain diseases, and plays an important role in the detection of epilepsy.
Spike waves are typical epilepsy characteristic waveforms, are usually recorded in electroencephalograms, are sharp relative to background waveforms, have high amplitude and transient characteristics, and clinically, current epilepsy examination mainly identifies spike waves of electroencephalogram signals through human eye detection. At present, clinical examination of epilepsia electroencephalogram is mainly to identify spike waves in electroencephalogram signals through manual detection, but the efficiency is low, the subjectivity is strong, the accuracy of results cannot be guaranteed, and therefore the spike wave automatic detection technology receives more and more attention in recent years.
There are many methods for identifying spike waves, and it is common that wavelet analysis method performs wavelet decomposition on time-frequency characteristics of epileptic brain electrical signals, and takes wavelet coefficients as input signals of machine learning classifier and neural network to detect spike waves. However, because the characteristics of the mother wavelet and the spike wave are different, the background electroencephalogram inhibition is poor in the signals obtained by decomposition and reconstruction, and the spike wave extraction effect is not ideal. The detection method of the singular point of the signal of the wavelet transform modulus maximum is another common method for spike detection, but the method can only detect positive-phase spike and has higher false positive rate when detecting negative-phase spike. Morphological filtering is another approach to study spike extraction, which uses predefined structural elements to match signals based on their geometric features to extract signals with similar morphology. The method has the characteristics of easy algorithm gradual change, definite physical significance, practicality, effectiveness and the like, can decompose a signal containing complex components into parts with different physical significances, separate the signal from a background and keep the global or local main morphological characteristics of the signal, but a single opening-closing (OC) or closing-opening (CO) operation can cause a statistical bias phenomenon, so that the detected spike and an actual spike have certain deviation on the waveform and the position.
Disclosure of Invention
Aiming at the problems in the prior art, the invention provides an intelligent epileptic spike detection method based on the fusion of adaptive template matching and a machine learning method, so as to improve the identification rate of epileptic spikes.
In order to achieve the purpose, the invention is realized by the following scheme: an intelligent epileptic spike detection method based on self-adaptive template matching and machine learning algorithm fusion comprises the following steps:
step S1: acquiring electroencephalogram signals, selecting an experimental object, establishing an epileptic electroencephalogram database, and marking spike waves in each channel of the electroencephalogram signals;
step S2: preprocessing the electroencephalogram signals, and removing high-frequency components and artifacts by using a 5-order Butterworth band-pass filter.
Step S3: and (3) detecting the epileptic spike by adopting a self-adaptive template matching method.
Step S4: and detecting the epileptic spike by adopting a machine learning method.
Step S5: and fusing the detection results of the step S3 and the step S4 to obtain a final spike detection result.
According to an embodiment of the present invention, the sampling frequency in step S1 is 500Hz, and a large amount of electroencephalogram data is required to be taken as an experimental sample, and the experimental subject includes people of different sexes and different ages.
According to an embodiment of the invention, in the adaptive template matching spike detection process, morphological characteristics such as rising edge slope, falling edge slope, amplitude and duration of manually marked spikes are counted to establish a universal spike template.
According to an embodiment of the present invention, in the process of detecting epilepsy spike by using an adaptive template matching method, the method includes:
step S31, counting the characteristics of rising edge slope, falling edge slope, amplitude height, duration and the like of the spike waveform in the electroencephalogram data, and defining a universal template;
step S32, setting the window width to 300, and carrying out general template matching operation on the electroencephalogram signals according to the time sequence to obtain candidate spike signals;
step S33, performing K-means clustering on the candidate spikes, and dividing the candidate spikes into different classes according to different waveforms;
step S34, counting the number of candidate spike waves in each spike wave cluster, if the number is less than 5% of the total number of candidate spike waves, rejecting the class, and finally taking the centroid of the rest classes as a new template;
and step S35, performing new template matching by respectively using the centroid of each class as a template, and superposing the results to obtain a spike detection result.
According to an embodiment of the present invention, the K-means clustering includes:
step S331: randomly selecting k samples in a sample set as an initial clustering centroid;
step S332: and calculating the distance between each sample and the initial centroid, reclassifying according to the minimum distance, and classifying each candidate spike into the class of the centroid closest to the candidate spike to obtain a clustering result.
Step S333: and averaging the obtained samples in each class to be used as the centroid of the next clustering.
Step S334: and repeating the steps S332 and S333 until the centroid position is not changed any more, and finishing the clustering.
According to an embodiment of the present invention, in the process of detecting epilepsy spike by using a machine learning method, the method includes:
step S41: the electroencephalogram signal is divided into single-channel segments with the length of 1s, if the segments have spikes, the segments are marked as 1, and if the segments do not have spikes, the segments are marked as 2. And respectively extracting time domain characteristics and frequency domain characteristics of the electroencephalogram fragments to construct epilepsy spike characteristic vectors with strong robustness.
Step S42: and randomly dividing the feature vectors into a training set and a testing set, and training a plurality of decision trees in a random forest classifier by using a plurality of electroencephalogram signal samples in the training set to form a random forest model.
Step S43: and inputting the data in the test set into the trained random forest model to obtain a spike detection result based on a machine learning method.
According to an embodiment of the invention, in the training of the random forest model, the method comprises the following steps:
step S421: the new training set with the same number of samples as the training set is extracted and put back in the training sample set.
Step S422: randomly sampling without playback in the feature vector set to form a feature vector set to be selected;
step S423: and (4) according to the candidate feature training set obtained in the step (S422), calculating the optimal splitting mode of each node and splitting the node without pruning until the impurity degree of each leaf node reaches the specified requirement to form a decision tree.
Step S424: and repeating the steps S421 to S423 until all decision trees are generated and integrated to obtain the random forest model.
By adopting the technical scheme of the invention, the epileptic spike detection is carried out by fusing self-adaptive template matching and machine learning, so that the identification rate of the epileptic spike is greatly improved.
Drawings
Fig. 1 is a general flowchart of an intelligent epileptic spike detection method based on adaptive template matching and machine learning algorithm fusion according to the present invention.
FIG. 2 is a flow chart of adaptive template matching spike detection according to the present invention.
FIG. 3 is a K-means clustering flow chart of the present invention.
FIG. 4 is a flowchart illustrating a machine learning spike detection process according to the present invention.
FIG. 5 is a flow chart of the training of a random forest model for machine learning spike detection in accordance with the present invention.
Detailed Description
Electroencephalogram signals generally contain a lot of physiological information about human diseases, and especially play an important role in the aspect of detection of epilepsy. The electroencephalogram signals contain many epileptic characteristic waves, and spike waves are typical waveforms in the electroencephalogram signals. Therefore, spike detection is required to be carried out on the epileptic brain electrical signals for better research. The existing spike method is difficult to completely and accurately determine the spike position, so that the research on epileptic diseases is greatly influenced. In view of this, the present embodiment provides an intelligent epileptic spike detection method based on adaptive template matching and machine learning algorithm fusion.
In order to make the objects, implementations and innovations of the present invention more prominent, the present invention will be further described in detail with reference to the accompanying drawings and examples.
Fig. 1 is a general flowchart of an intelligent epileptic spike detection method based on adaptive template matching and machine learning algorithm fusion, including:
step S1: acquiring an electroencephalogram signal: selecting an experimental object, acquiring electroencephalogram data of an epileptic by using electroencephalogram acquisition equipment, and establishing an experimental database;
step S2: data preprocessing: performing Butterworth band-pass filtering on the acquired original EEG data to obtain a standard EEG signal;
step S3: self-adaptive template matching spike detection: firstly, defining a universal template according to the waveform characteristics of an epileptic spike, and carrying out universal template matching to obtain a candidate spike signal; then clustering the candidate spikes by using a K-means algorithm to obtain a plurality of classes; counting the number of candidate spikes in each class, and if the number of spikes is less than 5% of the total number of spikes, rejecting the class; and respectively using the screened class centers as new templates to perform self-adaptive template matching, and adding all matching results to obtain spike detection results.
Step S4: machine learning spike detection: firstly, segmenting a brain electrical signal into electroencephalogram segments with the length of 1s, then extracting time domain and frequency domain characteristics in each electroencephalogram segment, and constructing spike characteristic vectors; and training a random forest classification model by using the feature vectors to obtain a spike detection result based on machine learning.
Step S5: and (3) fusing detection results: the S3 spike detection method and the S4 spike detection method are fused, and if the spike is detected by both S3 and S4, the spike is regarded as an epileptic spike.
The intelligent epileptic spike detection method based on the fusion of adaptive template matching and machine learning algorithm proposed in the present embodiment is described in detail below with reference to fig. 1 to 5.
The intelligent epileptic spike detection method based on the fusion of the adaptive template matching and the machine learning algorithm provided by the embodiment starts in step S1, wherein a multi-lead electroencephalograph is used for collecting long-range monitoring electroencephalograms of a patient, the sampling frequency is 500Hz, the electrode distribution adopts the international 10-20 electroencephalogram collection standard, 19-channel electroencephalogram data are collected in total, and a large number of electroencephalograms of experimental bodies with different sexes and different ages are collected to obtain a plurality of electroencephalogram signal samples. Marking a plurality of electroencephalogram signal samples by a professional doctor, and marking the spike waveform in each channel of the electroencephalogram signal.
Then, step S2 is executed to perform preprocessing operation on the brain wave. A5-order Butterworth band-pass filter is adopted to filter frequency components above 32Hz and below 0.5Hz, and the interference of noise and artifacts is reduced.
And step S3, carrying out self-adaptive template matching on the electroencephalogram signals after the preprocessing operation to obtain a spike detection result. The adaptive template matching spike detection method will be described in detail below with reference to fig. 2.
Firstly, statistical analysis is carried out on spike waveforms marked in the electroencephalogram signal, and the average values of the rising edge slope, the falling edge slope, the peak value and the duration of all marked spike waveforms are respectively obtained and used as standards to establish a universal template (step S31). Then, setting the window width to 300, and performing a general template matching operation on the electroencephalogram signals in time sequence to obtain candidate spike signals (step S32). And performing K-means clustering on the candidate spikes, and classifying the candidate spikes into different classes according to different waveforms (step S33). And counting the number of candidate spikes in each spike cluster, if the number is less than 5% of the total candidate spike number, rejecting the class, and finally taking the centroid of the rest class as a new template (step S34). And (4) performing new template matching by using the centroid of each class as a template, and superposing the results to obtain a spike detection result (step S35).
As shown in fig. 3, the process of K-means clustering is as follows:
k samples are randomly selected among the n candidate spikes as an initial cluster centroid (step S331). And calculating the distance between each candidate spike and each initial centroid, reclassifying according to the minimum distance, and classifying each candidate spike into the class of the closest centroid to obtain a clustering result (step S332). The candidate spikes in each cluster are averaged as the centroid of the next cluster (step S333). And repeating the steps S332 and S333 until the position of the centroid is not changed any more or the clustering frequency reaches the requirement, finishing clustering and obtaining a clustering result (step S334). In this embodiment, the number k of the initial centroids is n, that is, each candidate spike is clustered as a centroid, and finally n clustering results are obtained.
And step S4, performing spike extraction on the preprocessed electroencephalogram signal by adopting a machine learning method. The machine learning spike detection method will be described in detail below with reference to fig. 3.
Firstly, segmenting the electroencephalogram signal of each channel into segments with the length of 1S, extracting a plurality of characteristic parameters of each segment, wherein the characteristic parameters comprise time domain characteristic parameters and frequency domain characteristic parameters, and constructing a characteristic vector corresponding to each electroencephalogram segment (step S41). The feature vectors are then divided into a training set and a test set, and the random forest classification model is trained using the data in the training set (step S42). The data in the test set is input into the random forest model, and the output result obtained after voting is the spike detection result after each decision tree, so that whether a spike exists in the segment can be detected (step S43).
The electroencephalogram segment obtained by segmentation in the step S41 is recorded as x (N), N is 1,2, …, N and N are the lengths of the electroencephalogram segments, the sampling frequency of the electroencephalogram signal is 500Hz, so that N is 500Hz, the rhythm wave is extracted through wavelet packet transformation before feature extraction, and because the spike wave frequency range is over 14Hz, the signal decomposition and reconstruction are carried out by using a db6 wavelet basis function to obtain β wave and gamma wave which are respectively recorded as x (N)1(n) and x2(n)。
The time domain characteristic parameters extracted in the step S41 comprise an original electroencephalogram signal x (n) and two rhythm wave signals x1(n) and x2Minimum, maximum, mean, standard deviation, kurtosis, skewness, and Hjo of (n)And (4) rth parameter. The minimum value Min and the maximum value Max are respectively the maximum value of the signal amplitude, the average value Mean is the electroencephalogram signal amplitude trend, and the formula is as follows:
the standard deviation SD reflects the difference between the amplitude and the average value of each sample point, and the formula is as follows:
Figure BDA0002250853620000072
wherein x (N) is electroencephalogram signal, N is the number of sampling points of x (N), and is the average value of the amplitudes of all the sampling points in x (N).
Kurtosis Kur represents the peak level of the data frequency distribution curve, and is given by the formula:
Figure BDA0002250853620000073
skewness Skaew represents the characteristic of amplitude asymmetry degree of the electroencephalogram signals, and the formula is as follows:
the Hjorth parameters include Hjorth mobility and Hjorth complexity:
hjorth mobility can be represented by the following equation:
the Hjorth complexity can be represented by the following equation:
Figure BDA0002250853620000083
wherein
Figure BDA0002250853620000084
dnfn=x(n)-x(n-1)。
The frequency domain characteristic parameter extracted in the step S41 comprises the energy E of two rhythm wavesiTwo rhythm wave energy to total signal energy ratio Ri
Extracting rhythm waves through wavelet packet transformation, and performing five-layer wavelet decomposition on the electroencephalogram signals by using a db6 wavelet function to obtain β waves and gamma waves which are respectively marked as x1(n) and x2(n)。
Two rhythm wave energy EiObtained from the following equation:
Figure BDA0002250853620000085
total energy E of signalallThe formula of (1) is as follows:
Figure BDA0002250853620000086
further, the energy ratio, R, of the rhythm wave can be calculatedi=Ei/Eall,i=1,2。
Step S42 is a training process of the random forest model, and the random forest classifier includes a plurality of decision trees, and its output class is determined by the maximum votes in the results of all the trees. And repeatedly and randomly selecting M samples from the M samples in the original training sample set by using a bootstrap resampling technology to generate a new training sample set, and then generating a random forest by using M individual decision tree classifiers. The essence of the random forest classifier is an improvement on a decision tree algorithm, a plurality of decision trees are combined together, and each tree is established by independent randomly extracted samples. Each tree in the forest has the same distribution, and the classification error depends on the classification capability of each tree and the relevance of each tree.
The data set is divided into a training set and a testing set, and the random forest model training process is described below with reference to fig. 5:
step S421: firstly, sampling with returning is carried out for M times from all the feature vector sets to form a feature set to be selected, and the number of samples in the feature set to be selected is the same as that of the samples in the original feature vector set.
Step S422: secondly, randomly selecting a certain number of feature vectors from the features to be selected, and selecting the optimal features.
Step S423: and (4) according to the candidate feature training set obtained in the step (S422), calculating the optimal splitting mode of each node and splitting the node without pruning until the impurity degree of each leaf node reaches the specified requirement to form a decision tree.
Step S424: and repeating the step S421 to the step S423 until all the decision trees stop growing, and generating a random forest.
And step S43, inputting the electroencephalogram data in the test set into a random forest model, obtaining a spike wave detection result after voting selection by a decision tree, determining an electroencephalogram segment where a spike wave is located, and further determining an electroencephalogram channel where the spike wave is located and a time point.
In step S5, the data in the test set is first subjected to adaptive template matching to obtain a spike detection result. And meanwhile, inputting the test set into a random forest model for classification to obtain a spike detection result. The results of the two methods are then compared in fusion and if detected as spikes by both methods at the same time, they are considered epileptic spikes.
The above description of the embodiments is only intended to facilitate the understanding of the method of the invention and its core idea. It should be noted that, for those skilled in the art, it is possible to make various improvements and modifications to the present invention without departing from the principle of the present invention, and those improvements and modifications also fall within the scope of the claims of the present invention.
The previous description of the disclosed embodiments is provided to enable any person skilled in the art to make or use the present invention. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the generic principles defined herein may be applied to other embodiments without departing from the spirit or scope of the invention. Thus, the present invention is not intended to be limited to the embodiments shown herein but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.

Claims (3)

1. The intelligent epilepsia spike detection method based on the fusion of adaptive template matching and a machine learning algorithm is characterized by comprising the following steps of:
step S1: collecting electroencephalogram signals; selecting an experimental object, acquiring electroencephalogram data of an epileptic by using electroencephalogram acquisition equipment, and establishing an experimental database;
step S2: preprocessing data; performing Butterworth band-pass filtering on the acquired original EEG data to obtain a standard EEG signal;
step S3: adaptive template matching spike detection; firstly, defining a universal template according to the waveform characteristics of an epileptic spike, and carrying out universal template matching to obtain a candidate spike signal; then clustering the candidate spikes by using a K-means algorithm to obtain a plurality of classes; counting the number of candidate spikes in each class, and if the number of spikes is less than 5% of the total number of spikes, rejecting the class; respectively using the screened class centers as new templates to perform self-adaptive template matching, and adding all matching results to obtain spike detection results;
step S4: machine learning spike detection; firstly, segmenting a brain electrical signal into electroencephalogram segments with the length of 1s, then extracting time domain and frequency domain characteristics in each electroencephalogram segment, and constructing spike characteristic vectors; training a random forest classification model by using the feature vectors to obtain a spike detection result based on machine learning;
step S5: fusing detection results; fusing the spike detection result of the step S3 and the spike detection result of the step S4, and if the spike signals are detected in the same segment in the steps S3 and S4, marking the final result as that the segment has the spike signals, and regarding the segment as an epileptic spike;
wherein the step S3 further includes:
step S31, counting the characteristics of rising edge slope, falling edge slope, amplitude height, duration and the like of the spike waveform in the electroencephalogram data, and defining a universal template;
step S32, setting the window width to 300, and carrying out general template matching operation on the electroencephalogram signals according to the time sequence to obtain candidate spike signals;
step S33, performing K-means clustering on the candidate spikes, and dividing the candidate spikes into different classes according to different waveforms;
step S34, counting the number of candidate spike waves in each spike wave cluster, if the number is less than 5% of the total number of candidate spike waves, rejecting the class, and finally taking the centroid of the rest classes as a new template;
step S35, new template matching is carried out by respectively using the mass center of each class as a template, and the results are superposed to obtain a spike detection result;
the step S4 further includes:
step S41, each channel of the electroencephalogram signal is divided into segments with the length of 1S, the time domain feature and the frequency domain feature of each segment are extracted, and the feature vector of each electroencephalogram segment is constructed;
step S42, dividing the feature vectors into a training set and a testing set, and training a random forest classification model by using data in the training set;
step S43, inputting the data in the test set into the random forest model, and obtaining an output result, which is a spike detection result, so as to detect whether there is a spike in this segment.
2. The intelligent epileptic spike detection method based on adaptive template matching fused with machine learning algorithm as claimed in claim 1, characterized in that in the data preprocessing, 5 th order IIR butterworth band pass filter with frequency range of 0.5-32Hz is used to remove noise and artifacts in EEG signals.
3. The intelligent epileptic spike detection method based on adaptive template matching and machine learning algorithm fusion as claimed in claim 1 or 2, characterized in that spike detection results are compared with detection results of adjacent channels, if the detection results have a "spike-to-spike" phenomenon in the adjacent channels, the spike is considered, otherwise, the result is discarded.
CN201911033704.8A 2019-10-28 2019-10-28 Epileptic spike intelligent detection device based on self-adaptive template matching and machine learning algorithm fusion Active CN110811609B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911033704.8A CN110811609B (en) 2019-10-28 2019-10-28 Epileptic spike intelligent detection device based on self-adaptive template matching and machine learning algorithm fusion

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911033704.8A CN110811609B (en) 2019-10-28 2019-10-28 Epileptic spike intelligent detection device based on self-adaptive template matching and machine learning algorithm fusion

Publications (2)

Publication Number Publication Date
CN110811609A true CN110811609A (en) 2020-02-21
CN110811609B CN110811609B (en) 2022-11-22

Family

ID=69551174

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911033704.8A Active CN110811609B (en) 2019-10-28 2019-10-28 Epileptic spike intelligent detection device based on self-adaptive template matching and machine learning algorithm fusion

Country Status (1)

Country Link
CN (1) CN110811609B (en)

Cited By (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111643076A (en) * 2020-05-13 2020-09-11 杭州电子科技大学 BECT spike intelligent detection method based on multi-channel electroencephalogram signals
CN112244871A (en) * 2020-09-25 2021-01-22 吉林大学 Amplitude integration electroencephalogram classification recognition system based on machine learning
CN112270314A (en) * 2020-12-22 2021-01-26 苏州国科康成医疗科技有限公司 Spike wave identification method and device, electronic equipment and computer readable storage medium
CN112336354A (en) * 2020-11-06 2021-02-09 山西三友和智慧信息技术股份有限公司 Epilepsy monitoring method based on EEG signal
CN112568868A (en) * 2020-10-16 2021-03-30 浙江赛微思生物科技有限公司 Automatic quantitative analysis method and device for electrophysiological signals of epilepsy model
CN112641449A (en) * 2020-12-18 2021-04-13 浙江大学 EEG signal-based rapid evaluation method for cranial nerve functional state detection
CN113197585A (en) * 2021-04-01 2021-08-03 燕山大学 Neuromuscular information interaction model construction and parameter identification optimization method
CN113397567A (en) * 2021-05-19 2021-09-17 中国航天科工集团第二研究院 Human behavior electroencephalogram signal classification method and system
CN113974653A (en) * 2021-11-30 2022-01-28 杭州妞诺霄云大数据科技有限公司 Optimized spike detection method and device based on Joyston index, storage medium and terminal
CN114027854A (en) * 2021-12-01 2022-02-11 杭州电子科技大学 BECT spike detection method based on optimal template matching and morphological feature extraction
CN114081508A (en) * 2021-10-28 2022-02-25 杭州电子科技大学 Spike detection method based on fusion of deep neural network and CCA (continuous channel assessment) features
CN114145755A (en) * 2021-12-21 2022-03-08 上海理工大学 Household epileptic seizure interactive intelligent monitoring system and method
CN114532994A (en) * 2022-03-23 2022-05-27 电子科技大学 Automatic detection method for unsupervised electroencephalogram high-frequency oscillation signals based on convolution variational self-encoder
CN115299946A (en) * 2022-08-25 2022-11-08 电子科技大学 Self-adaptive input signal channel screening circuit introducing detection result feedback
CN115670397A (en) * 2022-11-17 2023-02-03 北京中科心研科技有限公司 PPG artifact identification method and device, storage medium and electronic equipment
CN116584900A (en) * 2023-07-17 2023-08-15 博睿康医疗科技(上海)有限公司 Signal steady state detection method, detection model and epilepsy detection device
CN115299946B (en) * 2022-08-25 2024-05-28 电子科技大学 Self-adaptive input signal channel screening circuit introducing detection result feedback

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101259016A (en) * 2007-03-06 2008-09-10 李小俚 Method for real time automatically detecting epileptic character wave
CN106137185A (en) * 2016-06-21 2016-11-23 华南理工大学 A kind of epileptic chracter wave detecting method based on structure of transvers plate small echo
US20180140256A1 (en) * 2016-11-21 2018-05-24 International Business Machines Corporation Brain Wave Processing for Diagnosis of a Subject
CN108962391A (en) * 2018-05-02 2018-12-07 杭州电子科技大学 Epileptics prediction technique early period based on wavelet packet character and random forest
CN110338786A (en) * 2019-06-28 2019-10-18 北京师范大学 A kind of identification of epileptiform discharges and classification method, system, device and medium

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101259016A (en) * 2007-03-06 2008-09-10 李小俚 Method for real time automatically detecting epileptic character wave
CN106137185A (en) * 2016-06-21 2016-11-23 华南理工大学 A kind of epileptic chracter wave detecting method based on structure of transvers plate small echo
US20180140256A1 (en) * 2016-11-21 2018-05-24 International Business Machines Corporation Brain Wave Processing for Diagnosis of a Subject
CN108962391A (en) * 2018-05-02 2018-12-07 杭州电子科技大学 Epileptics prediction technique early period based on wavelet packet character and random forest
CN110338786A (en) * 2019-06-28 2019-10-18 北京师范大学 A kind of identification of epileptiform discharges and classification method, system, device and medium

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
DUANPO WU 等: "Automatic Epileptic Seizures Joint Detection Algorithm Based on Improved Multi-Domain Feature of cEEG and Spike Feature of aEEG", 《IEEE ACCESS》 *

Cited By (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111643076A (en) * 2020-05-13 2020-09-11 杭州电子科技大学 BECT spike intelligent detection method based on multi-channel electroencephalogram signals
CN112244871A (en) * 2020-09-25 2021-01-22 吉林大学 Amplitude integration electroencephalogram classification recognition system based on machine learning
CN112568868A (en) * 2020-10-16 2021-03-30 浙江赛微思生物科技有限公司 Automatic quantitative analysis method and device for electrophysiological signals of epilepsy model
CN112336354A (en) * 2020-11-06 2021-02-09 山西三友和智慧信息技术股份有限公司 Epilepsy monitoring method based on EEG signal
CN112641449A (en) * 2020-12-18 2021-04-13 浙江大学 EEG signal-based rapid evaluation method for cranial nerve functional state detection
CN112270314A (en) * 2020-12-22 2021-01-26 苏州国科康成医疗科技有限公司 Spike wave identification method and device, electronic equipment and computer readable storage medium
CN113197585B (en) * 2021-04-01 2022-02-18 燕山大学 Neuromuscular information interaction model construction and parameter identification optimization method
CN113197585A (en) * 2021-04-01 2021-08-03 燕山大学 Neuromuscular information interaction model construction and parameter identification optimization method
CN113397567A (en) * 2021-05-19 2021-09-17 中国航天科工集团第二研究院 Human behavior electroencephalogram signal classification method and system
CN114081508A (en) * 2021-10-28 2022-02-25 杭州电子科技大学 Spike detection method based on fusion of deep neural network and CCA (continuous channel assessment) features
CN114081508B (en) * 2021-10-28 2024-05-14 杭州电子科技大学 Spike detection method based on fusion of deep neural network and CCA (common cancer cell) characteristics
CN113974653A (en) * 2021-11-30 2022-01-28 杭州妞诺霄云大数据科技有限公司 Optimized spike detection method and device based on Joyston index, storage medium and terminal
CN113974653B (en) * 2021-11-30 2024-03-22 杭州妞诺霄云大数据科技有限公司 Method and device for detecting optimized spike based on about log index, storage medium and terminal
CN114027854A (en) * 2021-12-01 2022-02-11 杭州电子科技大学 BECT spike detection method based on optimal template matching and morphological feature extraction
CN114145755B (en) * 2021-12-21 2023-09-01 上海理工大学 Household epileptic seizure interactive intelligent monitoring system and method
CN114145755A (en) * 2021-12-21 2022-03-08 上海理工大学 Household epileptic seizure interactive intelligent monitoring system and method
CN114532994A (en) * 2022-03-23 2022-05-27 电子科技大学 Automatic detection method for unsupervised electroencephalogram high-frequency oscillation signals based on convolution variational self-encoder
CN115299946A (en) * 2022-08-25 2022-11-08 电子科技大学 Self-adaptive input signal channel screening circuit introducing detection result feedback
CN115299946B (en) * 2022-08-25 2024-05-28 电子科技大学 Self-adaptive input signal channel screening circuit introducing detection result feedback
CN115670397A (en) * 2022-11-17 2023-02-03 北京中科心研科技有限公司 PPG artifact identification method and device, storage medium and electronic equipment
CN116584900A (en) * 2023-07-17 2023-08-15 博睿康医疗科技(上海)有限公司 Signal steady state detection method, detection model and epilepsy detection device
CN116584900B (en) * 2023-07-17 2023-10-31 博睿康医疗科技(上海)有限公司 Signal steady state detection method, detection model and epilepsy detection device

Also Published As

Publication number Publication date
CN110811609B (en) 2022-11-22

Similar Documents

Publication Publication Date Title
CN110811609B (en) Epileptic spike intelligent detection device based on self-adaptive template matching and machine learning algorithm fusion
US20240023886A1 (en) Noninvasive method and system for sleep apnea detection
CN110338786B (en) Epileptic discharge identification and classification method, system, device and medium
CN109758145B (en) Automatic sleep staging method based on electroencephalogram causal relationship
CN112741638B (en) Medical diagnosis auxiliary system based on EEG signal
Kumari et al. Seizure detection in EEG using time frequency analysis and SVM
CN109602417A (en) Sleep stage method and system based on random forest
Boubchir et al. A review of feature extraction for EEG epileptic seizure detection and classification
CN110135285B (en) Electroencephalogram resting state identity authentication method and device using single-lead equipment
CN113057648A (en) ECG signal classification method based on composite LSTM structure
CN112869716B (en) Pulse feature identification system and method based on two-channel convolutional neural network
CN111643076A (en) BECT spike intelligent detection method based on multi-channel electroencephalogram signals
CN111832537B (en) Abnormal electrocardiosignal identification method and abnormal electrocardiosignal identification device
CN108420406A (en) Method based on pulse wave sleep stage
CN105574348B (en) A kind of heart beat cycle acquisition methods and device based on BCG signals
Mahesh et al. ECG arrhythmia classification based on logistic model tree
Wu et al. HFO detection in epilepsy: a stacked denoising autoencoder and sample weight adjusting factors-based method
CN111887811B (en) Brain abnormal discharge detection method and system based on electroencephalogram signal characteristics
Awang et al. Analysis of EEG signals by eigenvector methods
Biran et al. Automatic qrs detection and segmentation using short time fourier transform and feature fusion
Cheng et al. Multiview feature fusion representation for interictal epileptiform spikes detection
Liu et al. Diagnosis of AF based on time and frequency features by using a hierarchical classifier
Benlamri et al. An automated system for analysis and interpretation of epileptiform activity in the EEG
Sanghavi et al. Detection of atrial fibrillation in electrocardiogram signals using machine learning
Zhang et al. Evaluation of single-lead ECG signal quality with different states of motion

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