CN108836317A - Dynamic electrocardiogram time scatter plot conversed analysis method - Google Patents
Dynamic electrocardiogram time scatter plot conversed analysis method Download PDFInfo
- Publication number
- CN108836317A CN108836317A CN201810680078.0A CN201810680078A CN108836317A CN 108836317 A CN108836317 A CN 108836317A CN 201810680078 A CN201810680078 A CN 201810680078A CN 108836317 A CN108836317 A CN 108836317A
- Authority
- CN
- China
- Prior art keywords
- scatter plot
- electrocardiogram
- time
- wave
- time scatter
- 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.)
- Pending
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]
- 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/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Cardiology (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Biophysics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
The present invention relates to a kind of dynamic electrocardiogram time scatter plot conversed analysis methods, include the following steps：S1：The ambulatory ecg signal that acquisition comes is carried out to automatically analyze identification R wave position；S2-1：According to the R wave position that step S1 is identified, using the interval of adjacent R wave as ordinate, the interval of adjacent R wave is also used as abscissa, draws time scatter plot；S2-2：The range of choosing will be enclosed on time scatter plot, carries out reverse operation；S2-3：The multiple electrocardiogram segments found by time scatter plot reverse method are shown, selects any electrocardiogram segment to amplify, checks for comparison.Dynamic electrocardiogram time scatter plot conversed analysis method provided by the invention can effectively auxiliary operation person filter out different interphase ranges and the ecg wave form of period, to carry out batch edit operation.
Description
Technical field
The present invention relates to a kind of dynamic electrocardiogram time scatter plot conversed analysis methods.
Background technique
Scatter plot is a kind of important means of nonlinear analysis electrocardiosignal, is with dynamic ECG especially for analysis
The long route cardiogram figure of representative is not replaced.First electrocardio scatter plot science summit on May 10th, 2011 and electrocardio
Scatterplot graphics composition is vertical, and meeting is in Beijing holds, from Beijing and participant expert in all parts of the country around Non-linear Chaos Theory, non-
How linear method is applied to cardiac diagnosis, how to recognize and use electrocardio scatter plot, electrocardio scatter plot to dynamic RR interphase
Physiological and pathological information, ventricular fibrillation and the spin wave theory of expression and the offer of sequence variation are discussed.Since then, it is dissipated based on electrocardio
The ambulatory ECG analysis of point diagram is paid more and more attention.Common Lorenz scatter plot uses the company in 24 hours long route cardiogram figures
Continuous RR interphase (Rn, Rn+1) does scatter plot as scatterplot coordinate in rectangular coordinate system.It can reflect the variation of adjacent R R interphase,
The nonlinear characteristic of heart rate variability is disclosed, dynamic ECG can be qualitatively judged.But Lorenz scatter plot lacks
Timing information in former electrocardiogram cannot reflect the time trend of HRV, especially be easy to mix the scatter plot of different time
Together.In order to more effectively be classified and retained the timing information that electrocardiosignal has to the not normal rhythm of the heart, when occurring
Between scatter plot, but how to solve to time scatter plot each point carry out ecg wave form control be anxious to be resolved ask
Topic.
The content of present invention
The present invention is directed to overcome above-mentioned the deficiencies in the prior art, a kind of dynamic electrocardiogram time scatter plot conversed analysis side is provided
Each of time scatter plot scatterplot or the scatterplot of selection area can inversely be found the electrocardio wave corresponding to it by method, the method
Shape.
In order to achieve the above object, the technical solution adopted by the present invention is as follows：
A kind of dynamic electrocardiogram time scatter plot conversed analysis method, includes the following steps：
S1：The ambulatory ecg signal that acquisition comes is carried out to automatically analyze identification R wave position；
S2-1：According to the R wave position that step S1 is identified, using the interval of adjacent R wave as ordinate, the interval of adjacent R wave
Also it is used as abscissa, draws time scatter plot；
S2-2：The range of choosing will be enclosed on time scatter plot, carries out reverse operation；
S2-3：The multiple electrocardiogram segments found by time scatter plot reverse method are shown, any heart is selected
The amplification of electrograph segment, is checked for comparison.
Preferably, the step S1 includes the following steps：
S1-1：12 are carried out to patient and leads dynamic electrocardiogram (ECG) data acquisition；
S1-2：12 next to acquisition lead electrocardiogram (ECG) data and are removed interference, then carry out R wave positioning again.
Preferably, in the step S2-1, with (Rn, Rn) as the scatterplot coordinate in rectangular coordinate system, wherein Rn is represented
N-th of RR interphase length, ordinate unit millisecond, abscissa unit minute, and mark time scatterplot, the normal sinus of ventricular premature beat
The time scatterplot of property heartbeat.
Preferably, in the step S2-2, range x1 to the x2 for selecting the x coordinate of range to represent the time is enclosed, encloses and y is selected to sit
Mark represents the range y1 to y2 of RR interphase, and the period x1 for filtering out the generation of RR interphase is determined according to selected time range<x<
X2 determines the RR interphase y1 of screening further according to selected RR interphase range<y<y2.
Method can effectively auxiliary operation person will be different for dynamic electrocardiogram time scatter plot conversed analysis provided by the invention
The ecg wave form of interphase range and period filter out, to carry out batch edit operation.
Detailed description of the invention
Below in conjunction with the embodiment that attached drawing provides, the present invention will be described in detail, wherein：
Fig. 1 is flow diagram of the embodiment of the present invention；
Fig. 2 is electrocardio of embodiment of the present invention time scatter plot schematic diagram；
Fig. 3 is electrocardio of embodiment of the present invention scatter plot selected areas and is found by time scatter plot conversed analysis method
Corresponding multiple electrocardiogram segments；
Electrocardiogram enlarged drawing corresponding to the electrocardiogram segment that Fig. 4 finds for the embodiment of the present invention.
Specific embodiment
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete
Site preparation description, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.It is based on
Embodiment in the present invention, it is obtained by those of ordinary skill in the art without making creative efforts every other
Embodiment shall fall within the protection scope of the present invention.
Referring to Fig. 1, technical solution provided by the invention：
Step S1：The ambulatory ecg signal that acquisition comes is carried out to automatically analyze identification R wave position：This step includes：
S1-1：The patient for looking for a ventricular arrhythmia more carries out 12 and leads dynamic electrocardiogram (ECG) data acquisition, acquisition time
It is 24 hours, pays attention to carrying out skin cleaning treatment, the data interference for guaranteeing that acquisition comes is smaller.
S1-2：12 lead electrocardiogram (ECG) data to what acquisition came and be removed interference by modes such as time domain filterings, then again into
The positioning of row R wave.
Step S2：The use of time scatter plot conversed analysis technology：
S2-1：The R wave position identified first according to step S1, using the interval of adjacent R wave as ordinate, adjacent R wave
Interval is also used as abscissa, draws time scatter plot as shown in Figure 2：Sat with (Rn, Rn) as the scatterplot in rectangular coordinate system
Mark, wherein Rn represents n-th of RR interphase length, ordinate unit ms, abscissa unit minute.And ventricular premature beat is represented with green
Time scatterplot, black represents the time scatterplot of normal sinus heartbeat.
S2-2:Then it encloses the range of choosing on time scatter plot with mouse as shown in Figure 2, carries out reverse operation：Enclose choosing
The x coordinate of range represents range x1 to the x2 of time, encloses the range y1 for selecting y-coordinate to represent RR interphase to y2, according to selected
Time range determines the period x1 for filtering out the generation of RR interphase<x<X2 determines screening further according to selected RR interphase range
RR interphase y1<y<y2.
S2-3:The multiple electrocardiogram segments found by time scatter plot reverse method are shown as shown in Figure 3, are selected
Enlarged drawing as shown in Figure 4 can be shown by selecting any electrocardiogram segment, be checked for comparison.
Above-described embodiment is used for illustrative purposes only, and is not limitation of the present invention, in relation to the general of technical field
Logical technical staff without departing from the spirit and scope of the present invention can be with various changes can be made and modification, therefore institute
There is equivalent technical solution also to should belong to scope of the invention, scope of patent protection of the invention should be limited by each claim.
Claims (4)
1. a kind of dynamic electrocardiogram time scatter plot conversed analysis method, it is characterised in that include the following steps：
S1：The ambulatory ecg signal that acquisition comes is carried out to automatically analyze identification R wave position；
S2-1：According to the R wave position that step S1 is identified, using the interval of adjacent R wave as ordinate, the interval of adjacent R wave is also made
For abscissa, time scatter plot is drawn；
S2-2：The range of choosing will be enclosed on time scatter plot, carries out reverse operation；
S2-3：The multiple electrocardiogram segments found by time scatter plot reverse method are shown, any electrocardiogram is selected
Segment amplification, is checked for comparison.
2. a kind of dynamic electrocardiogram time scatter plot conversed analysis method according to claim 1, it is characterised in that the step
Rapid S1 includes the following steps：
S1-1：12 are carried out to patient and leads dynamic electrocardiogram (ECG) data acquisition；
S1-2：12 next to acquisition lead electrocardiogram (ECG) data and are removed interference, then carry out R wave positioning again.
3. a kind of dynamic electrocardiogram time scatter plot conversed analysis method according to claim 1, it is characterised in that：
In the step S2-1, with (Rn, Rn) as the scatterplot coordinate in rectangular coordinate system, wherein Rn represents n-th of RR interphase
Length, ordinate unit millisecond, abscissa unit minute, and mark the time scatterplot of ventricular premature beat, normal sinus heartbeat when
Between scatterplot.
4. a kind of dynamic electrocardiogram time scatter plot conversed analysis method according to claim 1, it is characterised in that：
In the step S2-2, range x1 to the x2 for selecting the x coordinate of range to represent the time is enclosed, is enclosed between selecting y-coordinate to represent RR
The range y1 to y2 of phase determines the period x1 for filtering out the generation of RR interphase according to selected time range<x<X2, further according to institute
The RR interphase range of choosing determines the RR interphase y1 of screening<y<y2.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810680078.0A CN108836317A (en) | 2018-06-27 | 2018-06-27 | Dynamic electrocardiogram time scatter plot conversed analysis method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810680078.0A CN108836317A (en) | 2018-06-27 | 2018-06-27 | Dynamic electrocardiogram time scatter plot conversed analysis method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN108836317A true CN108836317A (en) | 2018-11-20 |
Family
ID=64199801
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810680078.0A Pending CN108836317A (en) | 2018-06-27 | 2018-06-27 | Dynamic electrocardiogram time scatter plot conversed analysis method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108836317A (en) |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109567794A (en) * | 2018-12-06 | 2019-04-05 | 厦门纳龙科技有限公司 | A kind of antialiasing method of abnormal P wave and device based on the superposition of P wave |
| CN109745037A (en) * | 2019-02-14 | 2019-05-14 | 广州视源电子科技股份有限公司 | Image processing method, device, equipment and storage medium |
| CN109770893A (en) * | 2019-03-08 | 2019-05-21 | 东南大学 | The method and device of atrial fibrillation position are quickly positioned in a kind of Holter analysis system |
| CN110009711A (en) * | 2019-05-20 | 2019-07-12 | 厦门纳龙科技有限公司 | A kind of heartbeat assemblage characteristic scatterplot drawing generating method, terminal device and storage medium |
| CN110403598A (en) * | 2019-08-05 | 2019-11-05 | 迪姆软件(北京)有限公司 | The processing method and processing device of three-dimensional electrocardio scatter plot |
| CN111543979A (en) * | 2020-05-13 | 2020-08-18 | 许祥林 | Method for outputting vector cardiogram through conventional leads |
| CN111685758A (en) * | 2020-06-23 | 2020-09-22 | 深圳星康医疗科技有限公司 | Novel display mode of electrocardiographic wave pattern strip |
| CN111920399A (en) * | 2020-09-04 | 2020-11-13 | 迪姆软件(北京)有限公司 | Analysis method and device for heart rate variability |
| CN112120689A (en) * | 2020-08-13 | 2020-12-25 | 深圳星康医疗科技有限公司 | Method and device for quickly marking long-range electrocardiogram data |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103479349A (en) * | 2013-09-25 | 2014-01-01 | 深圳市理邦精密仪器股份有限公司 | Electrocardiosignal data acquisition and processing method and system |
-
2018
- 2018-06-27 CN CN201810680078.0A patent/CN108836317A/en active Pending
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103479349A (en) * | 2013-09-25 | 2014-01-01 | 深圳市理邦精密仪器股份有限公司 | Electrocardiosignal data acquisition and processing method and system |
Non-Patent Citations (1)
| Title |
|---|
| 向晋涛等: "时间RR间期散点图及其逆向技术", 《中国心脏起搏与心电生理杂志》 * |
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109567794A (en) * | 2018-12-06 | 2019-04-05 | 厦门纳龙科技有限公司 | A kind of antialiasing method of abnormal P wave and device based on the superposition of P wave |
| CN109567794B (en) * | 2018-12-06 | 2022-03-25 | 厦门纳龙健康科技股份有限公司 | Abnormal P wave anti-aliasing method and device based on P wave superposition |
| CN109745037A (en) * | 2019-02-14 | 2019-05-14 | 广州视源电子科技股份有限公司 | Image processing method, device, equipment and storage medium |
| CN109770893A (en) * | 2019-03-08 | 2019-05-21 | 东南大学 | The method and device of atrial fibrillation position are quickly positioned in a kind of Holter analysis system |
| CN110009711B (en) * | 2019-05-20 | 2021-11-16 | 厦门纳龙科技有限公司 | Heart beat combination characteristic scatter diagram generation method, terminal device and storage medium |
| CN110009711A (en) * | 2019-05-20 | 2019-07-12 | 厦门纳龙科技有限公司 | A kind of heartbeat assemblage characteristic scatterplot drawing generating method, terminal device and storage medium |
| CN110403598A (en) * | 2019-08-05 | 2019-11-05 | 迪姆软件(北京)有限公司 | The processing method and processing device of three-dimensional electrocardio scatter plot |
| CN111543979A (en) * | 2020-05-13 | 2020-08-18 | 许祥林 | Method for outputting vector cardiogram through conventional leads |
| CN111543979B (en) * | 2020-05-13 | 2024-02-06 | 许祥林 | Method for outputting electrocardiographic vector diagram by conventional leads |
| CN111685758A (en) * | 2020-06-23 | 2020-09-22 | 深圳星康医疗科技有限公司 | Novel display mode of electrocardiographic wave pattern strip |
| CN112120689A (en) * | 2020-08-13 | 2020-12-25 | 深圳星康医疗科技有限公司 | Method and device for quickly marking long-range electrocardiogram data |
| WO2022032947A1 (en) * | 2020-08-13 | 2022-02-17 | 深圳星康医疗科技有限公司 | Method and apparatus for rapidly labeling long-term electrocardiographic data |
| CN111920399A (en) * | 2020-09-04 | 2020-11-13 | 迪姆软件(北京)有限公司 | Analysis method and device for heart rate variability |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108836317A (en) | Dynamic electrocardiogram time scatter plot conversed analysis method | |
| US5520191A (en) | Myocardial ischemia and infarction analysis and monitoring method and apparatus | |
| EP1885240B1 (en) | Analysis of electrocardiogram signals | |
| US6169919B1 (en) | System and method for quantifying alternation in an electrocardiogram signal | |
| KenKnight et al. | Regional capture of fibrillating ventricular myocardium: evidence of an excitable gap | |
| CN103399639B (en) | Brain-machine interface method and device is combined based on SSVEP and P300 | |
| CN105686826B (en) | Detecting and displaying irregular periodic waveforms | |
| EP1808125A1 (en) | Electrophysiological system for analysing an intracardiac electrocardiogram | |
| JPH01110345A (en) | Apparatus and method for detecting heart character by electric stimulation | |
| CN103829941B (en) | A kind of multidimensional electrocardiosignal imaging system and method | |
| CN102397067A (en) | Method and system for performing superposition analysis on exercise load electrocardio waveforms | |
| Ihara et al. | Atrial repolarization as observable during the PQ interval | |
| CN104586385A (en) | Electrocardiogram-based heart rate analysis method and equipment | |
| Prasad et al. | Heart rate detection using Hilbert transform | |
| Ho | Electrophysiology of arrhythmias: practical images for diagnosis and ablation | |
| US6625483B2 (en) | Method of noninvasively determining a patient's susceptibility to arrhythmia | |
| CN113384277A (en) | Electrocardiogram data classification method and classification system | |
| CN110403598A (en) | The processing method and processing device of three-dimensional electrocardio scatter plot | |
| Walcott et al. | Endocardial wave front organization during ventricular fibrillation in humans | |
| Stilli et al. | Newer data on the configuration and variability ranges of body surface maps in a sample of normal subjects | |
| CN109567794B (en) | Abnormal P wave anti-aliasing method and device based on P wave superposition | |
| CN109420254A (en) | A kind of pacemaker and cardiac pacing system for realizing CRT treatment based on single electrode | |
| Roudijk et al. | Mean temporal spatial isochrones direction as marker for activation delay in patients with arrhythmogenic cardiomyopathy | |
| Hirata et al. | Recurrent episodes of atrioventricular nodal reentrant tachycardia: Sites of ablation success, ablation endpoint, and primary culprits for recurrence | |
| US20100179396A1 (en) | Rapid method for analyzing bio-signal instantaneously by phase space complexity difference and its device |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB02 | Change of applicant information |
Address after: 510000 Second Floor, No. 604 Guangshan Second Road, Guangzhou, Guangdong Province (Room 202) Applicant after: Guangzhou Di Mao Mdt InfoTech Ltd Address before: 510000 Room 2101, Renfeng Building, 490 Tianhe Road, Tianhe District, Guangzhou City, Guangdong Province Applicant before: Guangzhou Di Mao Mdt InfoTech Ltd |
|
| CB02 | Change of applicant information | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20181120 |
|
| RJ01 | Rejection of invention patent application after publication |