CN108652613A - The method and device that signal time-frequency figure generates - Google Patents
The method and device that signal time-frequency figure generates Download PDFInfo
- Publication number
- CN108652613A CN108652613A CN201710204372.XA CN201710204372A CN108652613A CN 108652613 A CN108652613 A CN 108652613A CN 201710204372 A CN201710204372 A CN 201710204372A CN 108652613 A CN108652613 A CN 108652613A
- Authority
- CN
- China
- Prior art keywords
- interval series
- normal interval
- heart
- power spectrum
- time
- 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
Links
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
-
- 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
- A61B5/349—Detecting specific parameters of the electrocardiograph cycle
- A61B5/352—Detecting R peaks, e.g. for synchronising diagnostic apparatus; Estimating R-R interval
-
- 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
-
- 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/7235—Details of waveform analysis
- A61B5/725—Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
-
- 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/7253—Details of waveform analysis characterised by using transforms
Abstract
The present invention is suitable for processing of biomedical signals field, provides a kind of method and device that signal time-frequency figure generates.The method includes：Original electrocardiographicdigital data are obtained, phase signal between RR is determined according to the original electrocardiographicdigital data；The outlying interval sequence in phase signal between the RR is removed, normal interval series are obtained；The power spectrum of the normal interval series is calculated, and, after being segmented to the normal interval series, then calculate separately the power spectrum of the normal interval series after segmentation；After showing the power spectrum of the normal interval series and being segmented to the normal interval series, then calculate separately the power spectrum of the normal interval series after segmentation.Pass through the above method so that user can obtain more information within the same time, greatly improve the operation ease of user.
Description
Technical field
A kind of method that the embodiment of the present invention belongs to processing of biomedical signals field more particularly to signal time-frequency figure generates
And device.
Background technology
What the heart rate of normal person was not constant between, but can be influenced by factors such as autonomic nerves systems, pass through friendship
Sense and parasympathetic adjust its variation, while also having certain circadian rhythm.The change of this sinus rate at any time
Change is known as heart rate variability (Heart Rate Variability, HRV).
HRV be for characterizing sinus arrhythmia degree index, can be by instantaneous heart rate or gradually heartbeat heart rate
Variation (phase between RR) obtain.It is numerous studies have shown that heart rate increases when air-breathing, decreased heart rate when expiration；It is hindered using vagus nerve
Disconnected agent can be such that heart rate increases, and vagal tone enhancing can make decreased heart rate.It can be with HRV come Indirect evaluation autonomic nerves system
To the dominating role of heart.The spectrum component of HRV can be obtained by power Spectral Estimation, low frequency components generally believe at present
It is sympathetic nerve and vagus nerve is coefficient as a result, and radio-frequency component mainly reflects vagal activity.Clinically HRV
Analysis has non-invasive, simplicity, can be divided into and make a variation and make a variation when long in short-term, to a variety of diseases such as hypertension, myocardial infarction
Clinical diagnosis and prognosis evaluation are of great significance.
In general, the phase refers to the time interval between two QRS wave peaks in electrocardiographic wave between RR.The phase is at some between RR
It is a very important reference parameter during the diagnosis of the diseases such as common rhythm abnormality.It is drawn after phase processing between RR
Time-frequency figure can react the case where HRV to a certain extent, but the presentation mode of time-frequency figure is relatively simple at present, otherwise it is general
The information analyzed when only showing long or the information for only showing limited short-time analysis cannot be that user is same at an interface
When provide both integrated information, to be unfavorable for operator while obtain two different information.
Currently, the common linear analysis method of time-frequency figure has temporal analysis, frequency domain analysis and Time-Frequency Analysis Method.
Time Domain Analysis is simple, but specificity and sensitivity is not high；Frequency domain analysis describes the signal energy in a period of time
The case where changing with its frequency, cannot reflect the variation of real-time power；Although both methods physiological significance is more clear, all
The process of autonomic nerve variation cannot be described dynamically.Time-Frequency Analysis Method can simultaneously divide signal in temporal frequency domain
Analysis, with the no advantage of first two method, but its presentation mode is relatively simple, the letter analyzed when or generally only showing long
Breath or the information for only showing limited short-time analysis cannot be that user provides the comprehensive letter of the two simultaneously at an interface
Breath.
Therefore, it is desirable to provide a kind of new technical solution, to solve the above problems.
Invention content
An embodiment of the present invention provides the method and devices that a kind of signal time-frequency figure generates, it is intended to solve existing time-frequency figure
Analysis and the problem of short-time analysis when cannot show long simultaneously.
The embodiment of the present invention in a first aspect, provide a kind of method that signal time-frequency figure generates, the method includes：
Original electrocardiographicdigital data are obtained, phase signal between RR is determined according to the original electrocardiographicdigital data；
The outlying interval sequence in phase signal between the RR is removed, normal interval series are obtained；
The power spectrum of the normal interval series is calculated, and, after being segmented to the normal interval series, then distinguish
Calculate the power spectrum of the normal interval series after segmentation；
After showing the power spectrum of the normal interval series and being segmented to the normal interval series, then count respectively
The power spectrum of normal interval series after point counting section.
The second aspect of the embodiment of the present invention, provides a kind of device that signal time-frequency figure generates, and described device includes：
Original electrocardiographicdigital data capture unit determines RR for obtaining original electrocardiographicdigital data according to the original electrocardiographicdigital data
Between phase signal；
Normal interval series acquiring unit obtains normal for removing the outlying interval sequence between the RR in phase signal
Interval series；
Spectra calculation unit, the power spectrum for calculating the normal interval series, and, to it is described it is normal between phase sequence
After row are segmented, then calculate separately the power spectrum of the normal interval series after segmentation；
Time-frequency figure display unit, power spectrum for showing the normal interval series and to the normal interval series
After being segmented, then calculate separately the power spectrum of the normal interval series after segmentation.
In embodiments of the present invention, since long time-histories power spectrum mainly reflects the general trend in whole time, and in short-term
Journey power spectrum is higher to the quality requirement of signal, is mainly used for reflecting a certain local information, therefore, the same interface display
Information above so that user can obtain more information within the same time, greatly improve the operation ease of user.
Description of the drawings
Fig. 1 is the flow chart for the method that a kind of signal time-frequency figure that first embodiment of the invention provides generates；
Fig. 2 is a kind of detail flowchart for the step S12 that first embodiment of the invention provides；
Fig. 3 is that the electro-cardiologic signal waveforms that first embodiment of the invention provides are shown and editing interface schematic diagram；
Fig. 4 is another detail flowchart for the step S12 that first embodiment of the invention provides；
Fig. 5 is between a kind of original RR interval series randomly selected that first embodiment of the invention provides and removal exception RR
The comparative result figure of phase；
Fig. 6 is a kind of time-frequency image schematic diagram that first embodiment of the invention provides；
Fig. 7 is the structure chart for the device that a kind of signal time-frequency figure that second embodiment of the invention provides generates.
Specific implementation mode
In order to make the purpose , technical scheme and advantage of the present invention be clearer, with reference to the accompanying drawings and embodiments, right
The present invention is further elaborated.It should be appreciated that the specific embodiments described herein are merely illustrative of the present invention, and
It is not used in the restriction present invention.
In the embodiment of the present invention, original electrocardiographicdigital data are obtained, phase signal between RR is determined according to the original electrocardiographicdigital data, goes
Except the outlying interval sequence in phase signal between the RR, normal interval series are obtained, calculate the power of the normal interval series
Spectrum, and, after being segmented to the normal interval series, then the power spectrum of the normal interval series after segmentation is calculated separately,
After showing the power spectrum of the normal interval series and being segmented to the normal interval series, then after calculating separately segmentation
Normal interval series power spectrum.
In order to illustrate technical solutions according to the invention, illustrated below by specific embodiment.
Embodiment one：
Fig. 1 shows the flow chart for the method that a kind of signal time-frequency figure that first embodiment of the invention provides generates, and is described in detail
It is as follows：
Step S11 obtains original electrocardiographicdigital data, and phase signal between RR is determined according to the original electrocardiographicdigital data.
Optionally, the step S11 is specifically included：
The original electrocardiographicdigital data of A1, continuous acquisition preset duration.Due to being related to long time-histories time frequency analysis, and when long time-histories
Frequency analysis was advisable with 24 hours, it is therefore desirable to which the electrocardiogram (ECG) data of continuous acquisition long period (being greater than or equal to 24 hours), which is used as, to be divided
The original electrocardiographicdigital data of analysis.
The characteristic point of A2, the detection original electrocardiographicdigital data extract R crest locations and in the top of corresponding R crest locations
The heart belonging to it is marked to clap type, clapping type according to the normal heart determines phase signal between RR.
QRS wave detection is carried out to original electrocardiographicdigital data, extract R crest locations and is marked in the top of corresponding R crest locations
Remember that the heart belonging to it claps type, label normally is that ectopic beat or artifact are labeled as other types of letter, here
For the convenience of description, the abnormal heart is clapped the heart bat for being referred to as non-N labels.Position, which is clapped, with the latter heart subtracts previous heart bat position
It sets to obtain RR interval series RR1, RR2,…,RRk,…,RRn；The accuracy of R peak detections can influence subsequent analysis, so
When selecting electrocardiosignal (or electrocardiogram (ECG) data) as possible the electrocardiosignal of preferable that lead record of selection signal mass ratio as
Input.Currently, the method type of R waves detection is various, Technical comparing is ripe, and QRS wave detection method is unlimited here, for example, also may be used
To be detected using Wavelet Modulus Maxima Algorithm etc..
Step S12 removes the outlying interval sequence in phase signal between the RR, obtains normal interval series.
When acquiring electrocardiogram (ECG) data, gathered person needs the wearing acquisition instrument of long period, gatherer process center telecommunications
It number is also just inevitably influenced by some noises and artefact, certain subjects may be due to individual behavior or body itself
The reason of body, generates ectopic beat, therefore needs to pre-process RR interval series before power spectrum analysis, with removal
Outlying interval sequence between the RR in phase signal.The algorithm comparison of phase is simple between many removal exception RR at present, treatment effect
It is unsatisfactory.Although simple manual identification ectopic beat and artificial artefact, and abnormal point is removed to obtain the phase between normal NN
Method enables to subsequent analysis result reliability higher, but for the data of long time-histories, manual handle is comparatively laborious,
Therefore it needs to expend the more time, also, higher in abnormal point occurrence frequency, manual identification also can most probably leak
Fall the deletion of certain points.
Optionally, as shown in Fig. 2, more accurately normal interval series, the step S12 are specifically included in order to obtain：
The heart of B1, display ECG signal and label clap type.
B2, heart bat type modification instruction input by user is received, type is clapped according to the heart and changes instruction modification flase drop
The heart claps type.After reading electrocardiogram (electrocardiogram, ECG) signal, allow user to after QRS wave detects
The heart clap type and modify, user can substantially browse electrocardiosignal, and the heart of the apparent flase drop of manual modification is clapped, and is revised as
The suitable heart claps label, carries out next step again later, can improve the accuracy of phase between removal exception RR in this way.Of the invention real
It applies in example, user is allowed to clap type into edlin to the heart, than such as whether show the length (ms) and heart bat type of phase between each RR,
And it is that the correct heart claps type etc. that can change the flase drop heart to clap type, specific electro-cardiologic signal waveforms are shown with editing interface such as
Shown in Fig. 3.
B3, type is clapped between the modification heart after obtained RR phase signal carry out at equal intervals interpolation processing, go trend handle and
It is filtered, obtains normal interval series.Since the RR interval series after abnormal point deletion are a sampling heterogeneous
Sequence, it is therefore desirable to which interpolation processing at equal intervals is carried out to it.Specific method uses spline method.
Phase signal carries out trend and handles between the RR after interpolation at equal intervals, it is therefore an objective to eliminate some due to non-stationary
Interference caused by fluctuation, improves the accuracy of HRV spectra calculations.Typically using 3 order polynomials come to the phase between after interpolation
Sequence elimination trend term.
And the above-mentioned method being filtered mainly using wavelet transformation filters out a certain frequency, such as the frequency of 0.5Hz or more
Rate information.Specifically, signal is carried out by n-layer wavelet decomposition according to the frequency of sequence after interpolation, will be greater than the high fdrequency component of 0.5Hz
Coefficient sets to 0 then usage factor reproducing sequence and can be obtained new sequence.
Optionally, as shown in figure 4, more accurately normal interval series, the step S12 are specifically included in order to obtain：
B1 ', phase signal between the corresponding RR of improper heart bat type is deleted.Phase between RR and the heart are clapped into type as inputting,
The heart that non-N labels are deleted on the basis of QRS wave detection and partially manual modification claps the corresponding previous phase between the latter RR,
Most of abnormal heart can be clapped and be deleted, but the dystopy pace-making heart that can only delete detection in this way is clapped, it is also possible to existed since R waves leak
The abnormal point that inspection and R waves are examined more.
B2 ', the RR after the improper heart claps type will be deleted between phase signal carry out not overlapping segmentation, calculate between each section of RR
The median of time value.Here not overlapping segmentation refers to that there is no the parts of overlapping between two adjacent sectionals.By phase sequence between RR
The not overlapping segmentation of row progress, each segment length is N, it is assumed that per one piece of data comprising a points of m (generally taking 3~14), if between RR
Phase sequence length is not the integral multiple of m, then final stage data length is overall length divided by the remainder of m, is divided into n sections altogether, n is more than
1。
B3 ', judge between each RR in each section time value with the difference of the median of corresponding section whether preset respectively
In difference range (such as in preset percentage range).
Between B4 ', each RR in each section time value with the difference of the median of corresponding section in preset difference range
When, retain time value between the RR, otherwise, delete time value between the RR, to obtain normal interval series.
Assuming that using the data of kth (k is less than or equal to n) section as input, the median of the segment data is calculated
MedRR；Time value between each RR in this section of cycle criterion, if the difference of time value and MedRR are in certain difference model between RR
In enclosing (when the difference range is percentage range, which can usually be set as 20%~30%), then retain the point
RR between time value otherwise delete the phase between the RR.The step is repeated, until each section of RR interval data is disposed, to obtain
Normal RR interval series.
Fig. 5 is shown the original RR interval series randomly selected from MIT-BIH databases and is implemented using the present invention
The comparative result figure of phase between secondary discrimination step removal exception RR between the removal exception RR that example provides in phase algorithm, can from Fig. 5
Know, by secondary discrimination, the obtained normal interval series of RR are more acurrate.
Step S13 calculates the power spectrum of the normal interval series, and, the normal interval series are segmented
Afterwards, then calculate separately segmentation after normal interval series power spectrum.
The calculating of power spectrum includes mainly two-part content.First part is frequency division when being applied to whole RR interval datas
The method of analysis calculates power spectrum, and second part is the phase segmentation between RR, and calculates the power spectrum data per one piece of data, the present invention
The calculation that this two parts content is selected in embodiment is different, but the not calculation of power limitation spectrum, can also select
With same computational methods, as long as it is the same in last display mode to ensure.
The common Time-Frequency Analysis Method of first part has Short Time Fourier Transform, wavelet transformation and non-linear time frequency analysis
Method etc., the embodiment of the present invention is unlimited to determine Time-frequency Analysis method, and a kind of appropriate method can be selected according to the advantage and disadvantage of method
Spectra calculation is carried out, after the completion of calculating, result is deposited in a time-frequency matrix.
The RR interval series after preprocessing process are subjected to not overlapping segmentation in second part, namely per one piece of data
It is not be overlapped between the section of front, and be continuous between section and section；The phase is divided s sections between whole RR, and each segment length is preferably 5
Minute (can certainly select more than 5 minutes) is denoted as the first duration, calculates every one piece of data difference power spectrum in s sections,
And each section of power spectrum data is stored in a matrix of consequence, it is s sections total.
Step S14, after showing the power spectrum of the normal interval series and being segmented to the normal interval series,
The power spectrum of the normal interval series after segmentation is calculated separately again.
In the step, when power spectrum is the power spectrum of normal interval series, obtained time-frequency figure is long time-histories power spectrum
Corresponding time-frequency figure；When power spectrum be the normal interval series are segmented after, then calculate separately segmentation after it is normal between
When the power spectrum of phase sequence, obtain being the corresponding spectrogram of journey power spectrum in short-term.
Optionally, the step S14 is specifically included：Determine the corresponding color of different power spectrum, and according to determining face
Color marks corresponding power spectrum in the time-frequency figure of display.Specifically, by the time-frequency matrix obtained in step S13 by power spectral value
Size is mapped with a preset color matrix, is the power spectral value calculated in time-frequency matrix, with different face
Color table shows the size of power spectral value, you can generates the time-frequency figure of phase signal between RR.
Optionally, include in step S14：
The corresponding time-frequency figure of power spectrum of C1, the display normal interval series.
C2, time-frequency figure click commands are received.
At the time of the trigger point that C3, basis send out the time-frequency figure click commands corresponds to, the moment corresponding RR is searched
Between phase segments.
C4, generate phase tendency chart between the corresponding spectrogram of the segments and RR, the spectrogram be to it is described it is normal between
After phase sequence is segmented, then calculate separately the power spectrum of the normal interval series after segmentation.
Specifically, when the user clicks after time-frequency figure, at the time of correspondence according to the trigger point of click, first duration is shown
Phase tendency chart between corresponding spectrogram and RR, phase tendency chart is distributed in same interface with time-frequency figure between spectrogram and RR.Optionally,
When showing for the first time, can display segment (such as 5 minutes) interval series power spectrum and RR between phase tendency chart, it is normal obtaining
After the power spectrum of interval series, then show the power spectrum of the normal interval series (at this point, wouldn't show that the phase becomes between spectrogram and RR
Gesture figure).
Optionally, in order to enable user obtains more detailed information, before or after the C3, including：
Show the ecg wave form information from the moment.Specifically, phase tendency chart is same between time-frequency figure, spectrogram, RR
The electrocardio heartbeat data for also having acquisition of one interface display, click time-frequency figure, at the time of correspondence according to the trigger point of click, in the heart
At the time of electrograph region can also navigate to specific, the person of being conveniently operated is checked.
As shown in fig. 6, the image of the top is the time-frequency figure of generation in Fig. 6, abscissa x is expressed as time, ordinate y
It is expressed as frequency, indicates the size of power spectral value with color on corresponding position (x, y).In the power spectrum curve drawn out
On, the size of current time-frequency figure power can be intuitively reacted by the difference of shade, in the centre of Fig. 6, including display
Spectrogram and RR between phase tendency chart the two figures.In the bottom of Fig. 6, it is shown that ecg wave form information (or signal).
Certainly, if phase tendency chart, ecg wave form information between the time-frequency figure at a certain moment, spectrogram, RR currently has been displayed,
When triggering again or clicking other moment on time-frequency figure, phase segments between other moment corresponding RR can be found first,
Then phase tendency chart between the corresponding spectrogram of the segments and RR is generated, the figure of the preceding paragraph is override；Other moment pair
Phase tendency chart will be shown in the lower section of time-frequency figure, while the electrocardio of the subject from the moment between the HRV spectrograms answered and RR
Shape information is shown in the bottom at interface.
Power spectrum is divided into the corresponding spectrogram of short time-history analysis and the time-frequency figure of long time-history analysis, long time-history analysis main
Reflect the general trend in whole time, and short time-history analysis is higher to the quality requirement of signal, is mainly used for reflecting a certain
The emphasis of the information of part, the two is different；And ecg wave form is the most basic information of collected subject electrocardiosignal, with
The essential information of upper information and subject in the same interface display, when can convenient for users to or doctor read and analysis,
Second is that certain foundation can be provided for clinic.
In first embodiment of the invention, original electrocardiographicdigital data are obtained, determine that the phase believes between RR according to the original electrocardiographicdigital data
Number, the outlying interval sequence in phase signal between the RR is removed, normal interval series are obtained, calculates the normal interval series
Power spectrum, and, after being segmented to the normal interval series, then calculate separately the power of the normal interval series after segmentation
Spectrum after showing the power spectrum of the normal interval series and being segmented to the normal interval series, then calculates separately point
The power spectrum of normal interval series after section, i.e. the time-frequency figure including long time-histories power spectrum and the in short-term corresponding frequency spectrum of journey power spectrum
Figure.Since long time-histories power spectrum mainly reflects the general trend in whole time, and journey power spectrum wants the quality of signal in short-term
It asks higher, is mainly used for reflecting a certain local information, therefore, the same interface display information above so that user is same
One time interior energy enough obtains more information, greatly improves the operation ease of user.
It should be understood that in embodiments of the present invention, size of the sequence numbers of the above procedures is not meant to the elder generation of execution sequence
Afterwards, the execution sequence of each process should be determined by its function and internal logic, the implementation process structure without coping with the embodiment of the present invention
At any restriction.
Embodiment two：
Fig. 7 shows the structure chart for the device that a kind of signal time-frequency figure that second embodiment of the invention provides generates, in order to
Convenient for explanation, illustrate only and the relevant part of the embodiment of the present invention.
The signal time-frequency figure generate device include：Original electrocardiographicdigital data capture unit 71, normal interval series obtain single
Member 72, spectra calculation unit 73, time-frequency figure display unit 74.Wherein：
Original electrocardiographicdigital data capture unit 71 is determined for obtaining original electrocardiographicdigital data according to the original electrocardiographicdigital data
Phase signal between RR.
Optionally, the original electrocardiographicdigital data capture unit 71 includes：
Original electrocardiographicdigital data continuous acquisition module is used for the original electrocardiographicdigital data of continuous acquisition preset duration.Continuous acquisition
Duration should be greater than or be equal to 24 hours.
Characteristic point detection module, the characteristic point for detecting the original electrocardiographicdigital data extract R crest locations and in correspondence
The heart of the overlay mark of R crest locations belonging to it claps type, and clapping type according to the normal heart determines phase signal between RR.Specifically,
It is N that the normal heart, which is clapped type mark, and ectopic beat or artifact are labeled as other types of letter, here for the convenience of description,
The abnormal heart is clapped the heart bat for being referred to as non-N labels.
Normal interval series acquiring unit 72 obtains just for removing the outlying interval sequence between the RR in phase signal
Normal interval series.
Optionally, the normal interval series acquiring unit 72 includes：
The heart claps type mark module, for showing that the heart of ECG signal and label claps type.
The heart claps type and changes command reception module, type modification instruction is clapped for receiving the heart input by user, according to described
The heart that the heart claps type modification instruction modification flase drop claps type.
RR interval series preprocessing modules are inserted at equal intervals for phase signal clapping the RR obtained after type to the modification heart
Value processing goes trend to handle and be filtered, and obtains normal interval series.Wherein, batten can be used in interpolation processing at equal intervals
Interpolation method；It is typically to eliminate trend term to the interval series after interpolation using 3 order polynomials to go trend processing；And above-mentioned filter
Wave processing mainly filters out a certain frequency using the method for wavelet transformation, such as the frequency information of 0.5Hz or more.Specifically, according to
Signal is carried out n-layer wavelet decomposition by the frequency of sequence after interpolation, and the high fdrequency component coefficient that will be greater than 0.5Hz sets to 0 and then utilize system
Number reproducing sequence can be obtained new sequence.
Optionally, the normal interval series acquiring unit 72 includes：
Phase removing module between improper RR claps phase signal between the corresponding RR of type for deleting the improper heart.
The median computing module of time value between RR is carried out for phase signal between deleting the RR after the improper heart claps type
Not overlapping segmentation, calculates the median of time value between each section of RR.RR interval series are subjected to not overlapping segmentation, each segment length is
N, it is assumed that per one piece of data comprising a points of m (generally taking 3~14), if RR interval series length is not the integral multiple of m, finally
One piece of data length is overall length divided by the remainder of m, is divided into n sections altogether, n is more than 1.
The median comparison module of time value between RR, for judging time value and corresponding section between each RR in each section respectively
Whether the difference of median is in preset difference range.Wherein, when preset difference range is preset percentage range,
The preset percentage can be set as 20%~30%.
Time value retains selecting module between RR, the difference for time value and the median of corresponding section between each RR in each section
When value is in preset difference range, retains time value between the RR, otherwise, delete time value between the RR, with obtain it is normal between the phase
Sequence.
Spectra calculation unit 73, the power spectrum for calculating the normal interval series, and, to it is described it is normal between the phase
After sequence is segmented, then calculate separately the power spectrum of the normal interval series after segmentation.
The calculating of power spectrum includes mainly two-part content.First part is frequency division when being applied to whole RR interval datas
The method of analysis calculates power spectrum, and second part is the phase segmentation between RR, and calculates the power spectrum data per one piece of data.
The common Time-Frequency Analysis Method of first part has Short Time Fourier Transform, wavelet transformation and non-linear time frequency analysis
Method etc., the embodiment of the present invention is unlimited to determine Time-frequency Analysis method, and a kind of appropriate method can be selected according to the advantage and disadvantage of method
Spectra calculation is carried out, after the completion of calculating, result is deposited in a time-frequency matrix.
The RR interval series after preprocessing process are subjected to not overlapping segmentation in second part, namely per one piece of data
It is not be overlapped between the section of front, and be continuous between section and section；The phase is divided s sections between whole RR, and each segment length is preferably 5
Minute (can certainly select more than 5 minutes) is denoted as the first duration, calculates every one piece of data difference power spectrum in s sections,
And each section of power spectrum data is stored in a matrix of consequence, it is s sections total.
Time-frequency figure display unit 74, power spectrum for showing the normal interval series and to it is described it is normal between phase sequence
After row are segmented, then calculate separately the power spectrum of the normal interval series after segmentation.
Optionally, the time-frequency figure display unit 74 includes：
The corresponding color determination module of power spectrum, for determining the corresponding color of different power spectrum；
Power spectrum mark module, for marking corresponding power spectrum in the time-frequency figure of display according to determining color.
Optionally, the time-frequency figure display unit 74 includes：
Time-frequency figure display module, for showing the corresponding time-frequency figure of the power spectrum of the normal interval series.
Time-frequency figure click commands receiving module, for receiving time-frequency figure click commands.
Phase segments searching module between RR, for according to send out the time-frequency figure click commands trigger point it is corresponding when
It carves, searches phase segments between the moment corresponding RR.
Spectrogram generation module, for generating phase tendency chart between the corresponding spectrogram of the segments and RR, the frequency spectrum
Figure is after being segmented to the normal interval series, then calculates separately the power spectrum of the normal interval series after segmentation.
Specifically, when the user clicks after time-frequency figure, at the time of correspondence according to the trigger point of click, first duration is shown
Phase tendency chart between corresponding spectrogram and RR, phase tendency chart is distributed in same interface with time-frequency figure between spectrogram and RR.
Optionally, the device of signal time-frequency figure generation includes：
Electrocardio wave information display unit, for showing the ecg wave form information from the moment.Specifically, with time-frequency
The electrocardio heartbeat data for also having acquisition of the same interface display of phase tendency chart, click time-frequency figure, according to point between figure, spectrogram, RR
At the time of the trigger point hit corresponds to, at the time of electrocardio graph region can also navigate to specific, the person of being conveniently operated is checked.
In second embodiment of the invention, since long time-histories power spectrum mainly reflects the general trend in whole time, and it is short
Time-histories power spectrum is higher to the quality requirement of signal, is mainly used for reflecting a certain local information, and therefore, the same interface is aobvious
Show information above so that user can obtain more information within the same time, greatly improve the operation ease of user.
Those of ordinary skill in the art may realize that lists described in conjunction with the examples disclosed in the embodiments of the present disclosure
Member and algorithm steps can be realized with the combination of electronic hardware or computer software and electronic hardware.These functions are actually
It is implemented in hardware or software, depends on the specific application and design constraint of technical solution.Professional technician
Each specific application can be used different methods to achieve the described function, but this realization is it is not considered that exceed
The scope of the present invention.
It is apparent to those skilled in the art that for convenience and simplicity of description, the system of foregoing description,
The specific work process of device and unit, can refer to corresponding processes in the foregoing method embodiment, and details are not described herein.
In several embodiments provided herein, it should be understood that disclosed systems, devices and methods, it can be with
It realizes by another way.For example, the apparatus embodiments described above are merely exemplary, for example, the unit
It divides, only a kind of division of logic function, formula that in actual implementation, there may be another division manner, such as multiple units or component
It can be combined or can be integrated into another system, or some features can be ignored or not executed.Another point, it is shown or
The mutual coupling, direct-coupling or communication connection discussed can be the indirect coupling by some interfaces, device or unit
It closes or communicates to connect, can be electrical, machinery or other forms.
The unit illustrated as separating component may or may not be physically separated, aobvious as unit
The component shown may or may not be physical unit, you can be located at a place, or may be distributed over multiple
In network element.Some or all of unit therein can be selected according to the actual needs to realize the mesh of this embodiment scheme
's.
In addition, each functional unit in each embodiment of the present invention can be integrated in a processing unit, it can also
It is that each unit physically exists alone, it can also be during two or more units be integrated in one unit.
It, can be with if the function is realized in the form of SFU software functional unit and when sold or used as an independent product
It is stored in a computer read/write memory medium.Based on this understanding, technical scheme of the present invention is substantially in other words
The part of the part that contributes to existing technology or the technical solution can be expressed in the form of software products, the meter
Calculation machine software product is stored in a storage medium, including some instructions are used so that a computer equipment (can be
People's computer, server or network equipment etc.) it performs all or part of the steps of the method described in the various embodiments of the present invention.
And storage medium above-mentioned includes：USB flash disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), arbitrary access are deposited
The various media that can store program code such as reservoir (RAM, Random Access Memory), magnetic disc or CD.
The above description is merely a specific embodiment, but scope of protection of the present invention is not limited thereto, any
Those familiar with the art in the technical scope disclosed by the present invention, can easily think of the change or the replacement, and should all contain
Lid is within protection scope of the present invention.Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.
Claims (11)
1. a kind of method that signal time-frequency figure generates, which is characterized in that the method includes：
Original electrocardiographicdigital data are obtained, phase signal between RR is determined according to the original electrocardiographicdigital data；
The outlying interval sequence in phase signal between the RR is removed, normal interval series are obtained；
The power spectrum of the normal interval series is calculated, and, after being segmented to the normal interval series, then calculate separately
The power spectrum of normal interval series after segmentation；
After showing the power spectrum of the normal interval series and being segmented to the normal interval series, then calculate separately point
The power spectrum of normal interval series after section.
2. according to the method described in claim 1, it is characterized in that, the acquisition original electrocardiographicdigital data, according to the original heart
Electric data determine phase signal between RR, specifically include：
The original electrocardiographicdigital data of continuous acquisition preset duration；
The characteristic point of the original electrocardiographicdigital data is detected, extracts R crest locations and in its institute of overlay mark of corresponding R crest locations
The heart of category claps type, and clapping type according to the normal heart determines phase signal between RR.
3. according to the method described in claim 2, it is characterized in that, outlying interval sequence between the removal RR in phase signal
Row, obtain normal interval series, specifically include：
Show that the heart of ECG signal and label claps type；
It receives the heart input by user and claps type modification instruction, the heart that type modification instruction modification flase drop is clapped according to the heart claps class
Type；
After clapping type to the modification heart obtained RR phase signal carry out at equal intervals interpolation processing, go trend handle and filtering at
Reason, obtains normal interval series.
4. according to the method described in claim 2, it is characterized in that, outlying interval sequence between the removal RR in phase signal
Row, obtain normal interval series, specifically include：
It deletes the improper heart and claps phase signal between the corresponding RR of type；
Phase signal carries out not overlapping segmentation between the RR after the improper heart claps type will be deleted, between the RR of each section of calculating in time value
Digit；
Judge between each RR in each section time value with the difference of the median of corresponding section whether in preset difference range respectively
It is interior；
Between each RR in each section time value with the difference of the median of corresponding section in preset difference range when, retain institute
Between stating RR otherwise time value deletes time value between the RR, to obtain normal interval series.
5. method according to any one of claims 1 to 4, which is characterized in that the display normal interval series
Power spectrum and after being segmented to the normal interval series, then calculate separately the power of the normal interval series after segmentation
Spectrum, including：
Show the corresponding time-frequency figure of the power spectrum of the normal interval series；
Receive time-frequency figure click commands；
At the time of correspondence according to the trigger point for sending out the time-frequency figure click commands, phase segmentation between the moment corresponding RR is searched
Number；
Generate phase tendency chart between the corresponding spectrogram of the segments and RR, the spectrogram be to the normal interval series into
After row segmentation, then calculate separately the power spectrum of the normal interval series after segmentation.
6. according to the method described in claim 5, it is characterized in that, generating the corresponding spectrogram of the segments and RR described
Between before or after phase tendency chart, including：
Show the ecg wave form information from the moment.
7. according to claim 1-4 any one of them methods, which is characterized in that described to show the corresponding time-frequency of the power spectrum
Figure, including：
Determine the corresponding color of different power spectrum；
According to determining color corresponding power spectrum is marked in the time-frequency figure of display.
8. the device that a kind of signal time-frequency figure generates, which is characterized in that described device includes：
Original electrocardiographicdigital data capture unit determines the phase between RR for obtaining original electrocardiographicdigital data according to the original electrocardiographicdigital data
Signal；
Normal interval series acquiring unit, for removing the outlying interval sequence between the RR in phase signal, obtain it is normal between the phase
Sequence；
Spectra calculation unit, the power spectrum for calculating the normal interval series, and, to the normal interval series into
After row segmentation, then calculate separately the power spectrum of the normal interval series after segmentation；
Time-frequency figure display unit, power spectrum for showing the normal interval series and carries out the normal interval series
After segmentation, then calculate separately the power spectrum of the normal interval series after segmentation.
9. device according to claim 8, which is characterized in that the original electrocardiographicdigital data capture unit includes：
Original electrocardiographicdigital data continuous acquisition module is used for the original electrocardiographicdigital data of continuous acquisition preset duration；
Characteristic point detection module, the characteristic point for detecting the original electrocardiographicdigital data extract R crest locations and in corresponding R waves
The heart of the overlay mark of peak position belonging to it claps type, and clapping type according to the normal heart determines phase signal between RR.
10. device according to claim 9, which is characterized in that the normal interval series acquiring unit includes：
The heart claps type mark module, for showing that the heart of ECG signal and label claps type；
The heart claps type and changes command reception module, claps type modification instruction for receiving the heart input by user, is clapped according to the heart
The heart that type changes instruction modification flase drop claps type；
RR interval series preprocessing modules are carried out for phase signal clapping the RR obtained after type to the modification heart at interpolation at equal intervals
It manages, trend is gone to handle and be filtered, obtain normal interval series.
11. device according to claim 9, which is characterized in that the normal interval series acquiring unit includes：
Phase removing module between improper RR claps phase signal between the corresponding RR of type for deleting the improper heart；
The median computing module of time value between RR is not handed over for phase signal between deleting the RR after the improper heart claps type
Folded segmentation, calculates the median of time value between each section of RR；
The median comparison module of time value between RR, the middle position for judging time value and corresponding section between each RR in each section respectively
Whether several differences is in preset percentage；
Time value retains selecting module between RR, exists with the difference of the median of corresponding section for time value between each RR in each section
When in preset percentage, retain time value between the RR, otherwise, delete time value between the RR, to obtain normal interval series.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710204372.XA CN108652613B (en) | 2017-03-30 | 2017-03-30 | Method and device for generating signal time-frequency diagram |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710204372.XA CN108652613B (en) | 2017-03-30 | 2017-03-30 | Method and device for generating signal time-frequency diagram |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108652613A true CN108652613A (en) | 2018-10-16 |
| CN108652613B CN108652613B (en) | 2020-11-03 |
Family
ID=63786446
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710204372.XA Active CN108652613B (en) | 2017-03-30 | 2017-03-30 | Method and device for generating signal time-frequency diagram |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108652613B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109953752A (en) * | 2019-03-05 | 2019-07-02 | 上海朗朗信息科技有限公司 | HRV analysis method, device, equipment and storage medium under arrhythmia conditions |
| CN111685754A (en) * | 2020-06-19 | 2020-09-22 | 苏州大学 | Heart rate calculation method and system for wearable ECG acquisition equipment |
| CN114098756A (en) * | 2021-11-08 | 2022-03-01 | 杭州菲诗奥医疗科技有限公司 | Cardiopulmonary coupling analysis method based on single-channel ECG signal |
| CN114431843A (en) * | 2022-01-25 | 2022-05-06 | 华中科技大学 | Portable intelligent medical measuring equipment |
Citations (20)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1113739A (en) * | 1994-06-23 | 1995-12-27 | 王湘生 | Function testing system for vegetative nervous system and its method |
| US5718235A (en) * | 1992-10-06 | 1998-02-17 | Gw Scientific, Inc. | Detection of abnormal and induction of normal heart rate variability |
| CN1162195C (en) * | 1998-02-26 | 2004-08-18 | 欧姆龙株式会社 | Exercise machine, physical strength evaluating method, and pulsation meter |
| US20050251056A1 (en) * | 2004-05-10 | 2005-11-10 | Medpond, Llc | Method and apparatus for processing respiration data and assessing autonomic function |
| CN101953682A (en) * | 2010-10-15 | 2011-01-26 | 张辉 | Heartbeat detection method based on cuff device |
| CN101969841A (en) * | 2008-03-14 | 2011-02-09 | 皇家飞利浦电子股份有限公司 | Modifying a psychophysiological state of a subject |
| CN102693546A (en) * | 2011-03-23 | 2012-09-26 | 苏州大学 | Signal time frequency characteristic image generation method and device thereof |
| CN102088904B (en) * | 2008-07-11 | 2013-03-20 | 皇家飞利浦电子股份有限公司 | Ballistocardiogram analysis method and device |
| CN103006210A (en) * | 2013-01-11 | 2013-04-03 | 山东师范大学 | Sinus heart rate turbulence trend detection method based on piecewise linearization |
| CN103110417A (en) * | 2013-02-28 | 2013-05-22 | 华东师范大学 | Automatic electrocardiogram recognition system |
| CN103479349A (en) * | 2013-09-25 | 2014-01-01 | 深圳市理邦精密仪器股份有限公司 | Electrocardiosignal data acquisition and processing method and system |
| CN103792502A (en) * | 2012-10-26 | 2014-05-14 | 美国西门子医疗解决公司 | Automatic system for timing in imaging |
| CN103961089A (en) * | 2014-05-27 | 2014-08-06 | 山东师范大学 | Sinus heart rate turbulence tendency detecting method based on segmented straight line fitting |
| CN104127194A (en) * | 2014-07-14 | 2014-11-05 | 华南理工大学 | Depression evaluating system and method based on heart rate variability analytical method |
| CN104161509A (en) * | 2014-08-08 | 2014-11-26 | 申岱 | Heart rate variability analyzing method based on amplitude spectrum and instruments |
| CN104545887A (en) * | 2014-12-24 | 2015-04-29 | 深圳先进技术研究院 | Method and device for identifying artifact electrocardiograph waveforms |
| WO2015136567A1 (en) * | 2014-03-12 | 2015-09-17 | Universita' Politecnica Delle Marche | Method for the filtering of electrocardiographic signals or the like, usable in particular for the monitoring of fetal heart function |
| CN105686801A (en) * | 2015-12-25 | 2016-06-22 | 四川大学 | An atrial fibrillation spontaneous termination predication method based on inherent time scale decomposing |
| CN105796096A (en) * | 2014-12-31 | 2016-07-27 | 深圳先进技术研究院 | Heart rate variability analysis method, heart rate variability analysis system and terminal |
| CN105979863A (en) * | 2013-10-30 | 2016-09-28 | Apn健康有限责任公司 | Heartbeat detection and categorization |
-
2017
- 2017-03-30 CN CN201710204372.XA patent/CN108652613B/en active Active
Patent Citations (20)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5718235A (en) * | 1992-10-06 | 1998-02-17 | Gw Scientific, Inc. | Detection of abnormal and induction of normal heart rate variability |
| CN1113739A (en) * | 1994-06-23 | 1995-12-27 | 王湘生 | Function testing system for vegetative nervous system and its method |
| CN1162195C (en) * | 1998-02-26 | 2004-08-18 | 欧姆龙株式会社 | Exercise machine, physical strength evaluating method, and pulsation meter |
| US20050251056A1 (en) * | 2004-05-10 | 2005-11-10 | Medpond, Llc | Method and apparatus for processing respiration data and assessing autonomic function |
| CN101969841A (en) * | 2008-03-14 | 2011-02-09 | 皇家飞利浦电子股份有限公司 | Modifying a psychophysiological state of a subject |
| CN102088904B (en) * | 2008-07-11 | 2013-03-20 | 皇家飞利浦电子股份有限公司 | Ballistocardiogram analysis method and device |
| CN101953682A (en) * | 2010-10-15 | 2011-01-26 | 张辉 | Heartbeat detection method based on cuff device |
| CN102693546A (en) * | 2011-03-23 | 2012-09-26 | 苏州大学 | Signal time frequency characteristic image generation method and device thereof |
| CN103792502A (en) * | 2012-10-26 | 2014-05-14 | 美国西门子医疗解决公司 | Automatic system for timing in imaging |
| CN103006210A (en) * | 2013-01-11 | 2013-04-03 | 山东师范大学 | Sinus heart rate turbulence trend detection method based on piecewise linearization |
| CN103110417A (en) * | 2013-02-28 | 2013-05-22 | 华东师范大学 | Automatic electrocardiogram recognition system |
| CN103479349A (en) * | 2013-09-25 | 2014-01-01 | 深圳市理邦精密仪器股份有限公司 | Electrocardiosignal data acquisition and processing method and system |
| CN105979863A (en) * | 2013-10-30 | 2016-09-28 | Apn健康有限责任公司 | Heartbeat detection and categorization |
| WO2015136567A1 (en) * | 2014-03-12 | 2015-09-17 | Universita' Politecnica Delle Marche | Method for the filtering of electrocardiographic signals or the like, usable in particular for the monitoring of fetal heart function |
| CN103961089A (en) * | 2014-05-27 | 2014-08-06 | 山东师范大学 | Sinus heart rate turbulence tendency detecting method based on segmented straight line fitting |
| CN104127194A (en) * | 2014-07-14 | 2014-11-05 | 华南理工大学 | Depression evaluating system and method based on heart rate variability analytical method |
| CN104161509A (en) * | 2014-08-08 | 2014-11-26 | 申岱 | Heart rate variability analyzing method based on amplitude spectrum and instruments |
| CN104545887A (en) * | 2014-12-24 | 2015-04-29 | 深圳先进技术研究院 | Method and device for identifying artifact electrocardiograph waveforms |
| CN105796096A (en) * | 2014-12-31 | 2016-07-27 | 深圳先进技术研究院 | Heart rate variability analysis method, heart rate variability analysis system and terminal |
| CN105686801A (en) * | 2015-12-25 | 2016-06-22 | 四川大学 | An atrial fibrillation spontaneous termination predication method based on inherent time scale decomposing |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109953752A (en) * | 2019-03-05 | 2019-07-02 | 上海朗朗信息科技有限公司 | HRV analysis method, device, equipment and storage medium under arrhythmia conditions |
| CN111685754A (en) * | 2020-06-19 | 2020-09-22 | 苏州大学 | Heart rate calculation method and system for wearable ECG acquisition equipment |
| CN111685754B (en) * | 2020-06-19 | 2022-05-10 | 苏州大学 | Heart rate calculation method and system for wearable ECG acquisition equipment |
| CN114098756A (en) * | 2021-11-08 | 2022-03-01 | 杭州菲诗奥医疗科技有限公司 | Cardiopulmonary coupling analysis method based on single-channel ECG signal |
| CN114431843A (en) * | 2022-01-25 | 2022-05-06 | 华中科技大学 | Portable intelligent medical measuring equipment |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108652613B (en) | 2020-11-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11147500B2 (en) | Electrocardiogram processing system for delineation and classification | |
| CN109171712A (en) | Auricular fibrillation recognition methods, device, equipment and computer readable storage medium | |
| CN108652613A (en) | The method and device that signal time-frequency figure generates | |
| de Chazal et al. | Automatic classification of ECG beats using waveform shape and heart beat interval features | |
| AU2021250923B2 (en) | Computer-implemented method of handling electrocardiogram data | |
| US20220031223A1 (en) | Electrocardiogram processing system for delineation and classification | |
| CN109998529A (en) | P wave, T wave detecting method in a kind of ECG signal based on Gaussian function fitting | |
| CN105902263A (en) | Method and device for assessing irregularity degree at RR intervals of electrocardiogram waveform | |
| JP5449207B2 (en) | System for morphological feature analysis of physiological data | |
| US11672464B2 (en) | Electrocardiogram processing system for delineation and classification | |
| CN1180746C (en) | Cardioelectric spectrometer | |
| Sanamdikar et al. | Extraction of different features of ECG signal for detection of cardiac arrhythmias by using wavelet transformation Db 6 | |
| JP2009082664A (en) | System of analysis, display, storage, filing and management for comprehensive processing of electrocardiograph data | |
| CN112438735B (en) | Electrocardiogram P wave detection method, analysis device and storage medium | |
| WO2010000009A1 (en) | Improved detection of cardiac dysfunction | |
| JP2001224565A (en) | Statistical mapping for physiological state of mammalian heart | |
| Sanamdikar et al. | Analysis of several characteristics of ECG signal for cardiac arrhythmia detection | |
| Suarez Leon et al. | PyECG: A software tool for the analysis of the QT interval in the electrocardiogram | |
| US20230277107A1 (en) | Systems and methods for visual tracing of abnormality in diagnostic and monitoring clinical applications | |
| Konstantin et al. | Noise-resilient Automatic Interpretation of Holter ECG Recordings | |
| Everss-Villalba et al. | Clinical severity of noise in ECG | |
| CN117224133A (en) | Electrocardiogram blood pressure linkage display method, device, equipment and readable storage medium | |
| Nafees et al. | Classification of Atrio-Ventricular Reentrant Tachycardia using Intracardiac Signals | |
| Singh et al. | Signal Quality Evaluation and Processing for QRS Detection in ECG based Smart Healthcare Systems | |
| Egorov et al. | Noise-Resilient Automatic Interpretation of Holter ECG Recordings |
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 |