CN111449649B - Electrocardiogram data processing method and device and computer readable storage medium - Google Patents
Electrocardiogram data processing method and device and computer readable storage medium Download PDFInfo
- Publication number
- CN111449649B CN111449649B CN201910048588.0A CN201910048588A CN111449649B CN 111449649 B CN111449649 B CN 111449649B CN 201910048588 A CN201910048588 A CN 201910048588A CN 111449649 B CN111449649 B CN 111449649B
- Authority
- CN
- China
- Prior art keywords
- segment
- atrial fibrillation
- heart rate
- preset
- rate value
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
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/361—Detecting fibrillation
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 invention discloses a method for processing electrocardiographic data, which comprises the following steps: generating an interval trend graph and an electrocardiographic waterfall graph according to the electrocardiographic data when the electrocardiographic data are acquired; detecting a to-be-selected atrial fibrillation fragment in the interval trend chart by taking a preset detection parameter as a detection condition; and analyzing the atrial fibrillation fragments to be selected according to the electrocardio waterfall graph to obtain atrial fibrillation fragments, and marking the atrial fibrillation fragments. The invention also discloses a processing device of the electrocardio data and a computer readable storage medium. The invention improves the accuracy of electrocardiograph data detection.
Description
Technical Field
The present invention relates to the field of medical signal processing, and in particular, to a method for processing electrocardiographic data, an apparatus for processing electrocardiographic data, and a computer readable storage medium.
Background
Among all current electrocardiographic examinations, dynamic electrocardiographic examination is an important examination method, and atrial fibrillation is an important arrhythmia screening in dynamic electrocardiographic examination. Atrial fibrillation is divided into whole-course atrial fibrillation, long Cheng Fang atrial fibrillation and paroxysmal atrial fibrillation, wherein the detection of paroxysmal atrial fibrillation is very difficult because the duration of paroxysmal atrial fibrillation is generally only a few seconds or a few minutes, and the f wave, which is an important sign, has the characteristics of small amplitude and high frequency of fibrillation. In the prior art, doctors generally perform overview analysis on long-range dynamic data, human eyes are used for identifying short-range atrial fibrillation, and atrial fibrillation events are manually defined, so that the workload of the doctors is increased, errors or missed detection are likely to occur due to inaccuracy of human factors, and clinical risks are easily caused.
Disclosure of Invention
The invention mainly aims to provide a processing method of electrocardiographic data, a processing device of electrocardiographic data and a computer readable storage medium, and aims to improve the accuracy of electrocardiographic data detection.
In order to achieve the above object, the present invention provides a method for processing electrocardiographic data, the method for processing electrocardiographic data comprising the steps of:
generating an interval trend graph and an electrocardiographic waterfall graph according to the electrocardiographic data when the electrocardiographic data are acquired;
detecting a to-be-selected atrial fibrillation fragment in the interval trend chart by taking a preset detection parameter as a detection condition;
and analyzing the atrial fibrillation fragments to be selected according to the electrocardio waterfall graph to obtain atrial fibrillation fragments, and marking the atrial fibrillation fragments.
Optionally, the step of determining the segments of atrial fibrillation to be selected in the interval trend graph includes:
dividing the interval trend graph into a plurality of segments with a preset number of heartbeats;
and sequentially determining whether each segment is a segment of atrial fibrillation to be selected according to detection parameters, wherein the detection parameters comprise at least one of dispersion, duration and heart rate value.
Optionally, the step of sequentially determining whether each segment is a segment of atrial fibrillation to be selected according to the detected parameters includes:
acquiring the dispersion degree of each segment in the interval trend graph, and judging whether the dispersion degree of each segment in the interval trend graph is larger than a preset dispersion degree or not;
and when the dispersion is larger than the preset dispersion, taking the segment as the atrial fibrillation segment to be selected.
Optionally, after the step of obtaining the dispersion of each segment in the interval trend graph and determining whether the dispersion of each segment in the interval trend graph is greater than a preset dispersion, the method for processing electrocardiographic data further includes:
when the dispersion is larger than the preset dispersion, judging whether the duration of the segment is larger than the preset duration;
and when the duration time is longer than the preset duration time, taking the segment as the atrial fibrillation segment to be selected.
Optionally, after the step of determining whether the duration of the segment is greater than the preset duration when the dispersion is greater than the preset dispersion, the method for processing electrocardiograph data further includes:
judging whether the heart rate value in the segment is larger than a preset heart rate value or not when the duration time is longer than the preset duration time;
and when the heart rate value is larger than the preset heart rate value, taking the segment as the atrial fibrillation segment to be selected.
Optionally, after the step of determining whether the heart rate value in the segment is greater than the preset heart rate value when the duration is greater than the preset duration, the method for processing electrocardiograph data further includes:
when the heart rate value is smaller than or equal to the preset heart rate value, dividing the segment into a plurality of sub-segments by taking a heart beat corresponding to the heart rate value as a demarcation point;
and when the duration time of the sub-segment is longer than the preset duration time, taking the sub-segment as the atrial fibrillation segment to be selected.
Optionally, the step of determining the atrial fibrillation segment in the electrocardiographic waterfall map according to the atrial fibrillation segment to be selected includes:
acquiring a position range of occurrence of P waves in the electrocardio-fall map, and acquiring an electrocardio-fall picture segment corresponding to the atrial fibrillation segment to be selected in the electrocardio-fall map;
and determining an electrocardio-fall picture segment in which the P wave does not appear according to the position range and the electrocardio-fall picture segment, and taking the electrocardio-fall picture segment in which the P wave does not appear as the atrial fibrillation segment.
Optionally, after the step of analyzing the atrial fibrillation segment according to the electrocardiograph waterfall graph to obtain an atrial fibrillation segment and marking the atrial fibrillation segment, the method for processing electrocardiograph data further includes:
acquiring parameter information corresponding to preset parameters of the atrial fibrillation fragments from the electrocardiographic waterfall map, wherein the parameter information comprises at least one of starting time, duration time, average dispersion and average heart rate value of the atrial fibrillation fragments;
and outputting the parameter information.
In order to achieve the above object, the present invention further provides an apparatus for processing electrocardiographic data, the apparatus for processing electrocardiographic data comprising:
the processing method comprises the steps of a memory, a processor and a processing program of electrocardiographic data, wherein the processing program of electrocardiographic data is stored in the memory and can run on the processor, and the processing method of electrocardiographic data is realized when the processing program of electrocardiographic data is executed by the processor.
In order to achieve the above object, the present invention further provides a computer readable storage medium, on which a processing program of electrocardiographic data is stored, the processing program of electrocardiographic data implementing the steps of the above-described processing method of electrocardiographic data when executed by a processor.
According to the electrocardio data processing method, the electrocardio data processing device and the computer readable storage medium, when the electrocardio data is acquired, an interval trend chart and an electrocardio waterfall chart are generated according to the electrocardio data, a to-be-selected atrial fibrillation fragment is detected in the interval trend chart by taking preset detection parameters as detection conditions, the atrial fibrillation fragment is obtained by analyzing the to-be-selected atrial fibrillation fragment according to the electrocardio waterfall chart, and the atrial fibrillation fragment is marked. The invention improves the accuracy of electrocardiograph data detection.
Drawings
FIG. 1 is a schematic diagram of a hardware operating environment of a terminal according to an embodiment of the present invention;
FIG. 2 is a flowchart of a first embodiment of a method for processing electrocardiographic data according to the present invention;
FIG. 3 is a flowchart of a second embodiment of a method for processing electrocardiographic data according to the present invention;
fig. 4 is a schematic diagram of an intermittent trend chart and an electrocardiographic waterfall chart according to the present invention.
The achievement of the objects, functional features and advantages of the present invention will be further described with reference to the accompanying drawings, in conjunction with the embodiments.
Detailed Description
It should be understood that the specific embodiments described herein are for purposes of illustration only and are not intended to limit the scope of the invention.
Atrial fibrillation is divided into whole-course atrial fibrillation, long Cheng Fang atrial fibrillation and paroxysmal atrial fibrillation, wherein the detection of paroxysmal atrial fibrillation is very difficult because the duration of paroxysmal atrial fibrillation is generally only a few seconds or a few minutes, and the f wave, which is an important sign, has the characteristics of small amplitude and high frequency of fibrillation. In the prior art, doctors generally perform overview analysis on long-range dynamic data, human eyes are used for identifying short-range atrial fibrillation, and atrial fibrillation events are manually defined, so that the workload of the doctors is increased, errors or missed detection are likely to occur due to inaccuracy of human factors, and clinical risks are easily caused. The invention provides a processing method of electrocardio data, which comprises the steps of generating an interval trend graph and an electrocardio waterfall graph according to electrocardio data when the electrocardio data are acquired, detecting a to-be-selected atrial fibrillation segment in the interval trend graph by taking a preset detection parameter as a detection condition, analyzing the to-be-selected atrial fibrillation segment according to the electrocardio waterfall graph to obtain the atrial fibrillation segment, and marking the atrial fibrillation segment, so that the accuracy of electrocardio data detection is improved.
As shown in fig. 1, fig. 1 is a schematic diagram of a hardware operating environment of a terminal according to an embodiment of the present invention.
The terminal of the embodiment of the invention comprises a processing device of electrocardiographic data.
As shown in fig. 1, fig. 1 is a schematic diagram of a hardware operating environment of a terminal according to an embodiment of the present invention.
Referring to fig. 1, the terminal may include: a processor 1001, such as a CPU, memory 1002, a communications bus 1003. The communication bus 1003 is used to implement connection communication between the constituent elements in the terminal. The memory 1002 may be a high-speed RAM memory or a stable memory (non-volatile memory), such as a disk memory. The memory 1002 may alternatively be a storage device separate from the processor 1001 described above.
It will be appreciated by those skilled in the art that the structure of the terminal shown in fig. 1 is not limiting of the terminal and may include more or fewer components than shown, or may combine certain components, or a different arrangement of components.
As shown in fig. 1, a processing program of electrocardiographic data may be included in a memory 1002 as a computer storage medium.
In the server shown in fig. 1, the processor 1001 may be configured to call a processing program of electrocardiographic data stored in the memory 1002, and perform the following operations:
generating an interval trend graph and an electrocardiographic waterfall graph according to the electrocardiographic data when the electrocardiographic data are acquired;
detecting a to-be-selected atrial fibrillation fragment in the interval trend chart by taking a preset detection parameter as a detection condition;
and analyzing the atrial fibrillation fragments to be selected according to the electrocardio waterfall graph to obtain atrial fibrillation fragments, and marking the atrial fibrillation fragments.
Further, the processor 1001 may be configured to call a processing program of the electrocardiographic data stored in the memory 1002, and perform the following operations:
dividing the interval trend graph into a plurality of segments with a preset number of heartbeats;
and sequentially determining whether each segment is a segment of atrial fibrillation to be selected according to detection parameters, wherein the detection parameters comprise at least one of dispersion, duration and heart rate value.
Further, the processor 1001 may be configured to call a processing program of the electrocardiographic data stored in the memory 1002, and perform the following operations:
acquiring the dispersion degree of each segment in the interval trend graph, and judging whether the dispersion degree of each segment in the interval trend graph is larger than a preset dispersion degree or not;
and when the dispersion is larger than the preset dispersion, taking the segment as the atrial fibrillation segment to be selected.
Further, the processor 1001 may be configured to call a processing program of the electrocardiographic data stored in the memory 1002, and perform the following operations:
when the dispersion is larger than the preset dispersion, judging whether the duration of the segment is larger than the preset duration;
and when the duration time is longer than the preset duration time, taking the segment as the atrial fibrillation segment to be selected.
Further, the processor 1001 may be configured to call a processing program of the electrocardiographic data stored in the memory 1002, and perform the following operations:
judging whether the heart rate value in the segment is larger than a preset heart rate value or not when the duration time is longer than the preset duration time;
and when the heart rate value is larger than the preset heart rate value, taking the segment as the atrial fibrillation segment to be selected.
Further, the processor 1001 may be configured to call a processing program of the electrocardiographic data stored in the memory 1002, and perform the following operations:
when the heart rate value is smaller than or equal to the preset heart rate value, dividing the segment into a plurality of sub-segments by taking a heart beat corresponding to the heart rate value as a demarcation point;
and when the duration time of the sub-segment is longer than the preset duration time, taking the sub-segment as the atrial fibrillation segment to be selected.
Further, the processor 1001 may be configured to call a processing program of the electrocardiographic data stored in the memory 1002, and perform the following operations:
acquiring a position range of occurrence of P waves in the electrocardio-fall map, and acquiring an electrocardio-fall picture segment corresponding to the atrial fibrillation segment to be selected in the electrocardio-fall map;
and determining an electrocardio-fall picture segment in which the P wave does not appear according to the position range and the electrocardio-fall picture segment, and taking the electrocardio-fall picture segment in which the P wave does not appear as the atrial fibrillation segment.
Further, the processor 1001 may be configured to call a processing program of the electrocardiographic data stored in the memory 1002, and perform the following operations:
acquiring parameter information corresponding to preset parameters of the atrial fibrillation fragments from the electrocardiographic waterfall map, wherein the parameter information comprises at least one of starting time, duration time, average dispersion and average heart rate value of the atrial fibrillation fragments;
and outputting the parameter information.
Referring to fig. 2, in an embodiment, the method for processing electrocardiographic data includes:
step S10, generating an interval trend graph and an electrocardiographic waterfall graph according to electrocardiographic data when the electrocardiographic data are acquired;
step S20, detecting a to-be-selected atrial fibrillation fragment in the interval trend chart by taking a preset detection parameter as a detection condition;
and step S30, analyzing the atrial fibrillation fragments to be selected according to the electrocardio waterfall diagram to obtain atrial fibrillation fragments, and marking the atrial fibrillation fragments.
In this embodiment, the execution body is an electrocardiographic data processing device, and the electrocardiographic data processing device may be an independent terminal or a terminal including an electrocardiographic data processing program.
After the processing program of the electrocardio data is started, initializing operation is carried out, when a card reading instruction is received, a path corresponding to the card reading instruction is obtained, and electrocardio data in the dynamic electrocardio recording box is read according to the path. And carrying out whole-course atrial fibrillation analysis on the electrocardiograph data to obtain an overall analysis result of the whole-course atrial fibrillation, wherein the overall analysis result at least comprises an interval trend graph and an electrocardiograph waterfall graph.
As shown in fig. 4, the upper graph in fig. 4 is an interval trend graph. The interval trend graph refers to an RR interval trend graph, the horizontal axis of the interval trend graph is time, the vertical axis of the interval trend graph is RR interval, the interval trend graph reflects the variation trend of the RR interval of whole-course heart beat, and the dispersion of the RR interval is intuitively displayed. RR interval is the time between R waves in two QRS waves, which is used to calculate the ventricular rate. The QRS complex reflects the changes in left and right ventricular depolarization potentials and time, with the first downward wave being the Q wave, then the upward wave being the R wave, and then the downward wave being the S wave. The P wave is an atrial depolarization wave, represents the excitation of the left atrium and the right atrium, and has important significance for diagnosing and distinguishing diagnosis of arrhythmia by analyzing the P wave.
As shown in fig. 4, the lower graph in fig. 4 is an electrocardiographic waterfall graph. The electrocardio waterfall diagram is a graphical presentation form of waveform characteristics, all cardiac waves are longitudinally overlapped in overlook, and different waveforms are represented by different colors and different brightnesses, so that whether a certain segment has a P wave or not can be judged through the electrocardio waterfall diagram, and whether the segment is an atrial fibrillation segment (one of the characteristics of the atrial fibrillation segment is that the P wave does not exist) is determined.
In this embodiment, the atrial fibrillation segment to be selected is initially determined according to the interval trend chart, and the atrial fibrillation segment is finally determined from the atrial fibrillation segments to be selected through the electrocardiographic waterfall chart. The step of determining the selected atrial fibrillation fragments according to the interval trend graph may be: and (3) scanning and analyzing the interval trend graph, and screening out the atrial fibrillation fragments to be selected by taking preset detection parameters as detection conditions. Specifically, the interval trend graph is divided into a plurality of segments by a preset number of heartbeats, and whether each segment is an atrial fibrillation segment to be selected or not is sequentially determined according to detection parameters, wherein the detection parameters comprise at least one of dispersion, duration and heart rate value. Optionally, the preset number is three.
The step of obtaining the dispersion of the segments comprises: acquiring the time difference between adjacent intervals in the segment; calculating the dispersion of the fragments according to the time difference. Specifically, taking 3 heart beats as an example, dividing the interval trend graph into a plurality of segments, wherein R peaks of 3 heart beats in time sequence are respectively R1, R2 and R3, then two RR intervals are respectively R1R2 and R2R3, and then the RR interval dispersion of the preset segment is s1= (r2r3—r1r2)/R1R 2, wherein R1R2 is a time period between R1 and R2, and R2R3 is a time period between R2 and R3.
Optionally, acquiring the dispersion of each segment in the interval trend graph, judging whether the dispersion of each segment in the interval trend graph is larger than a preset dispersion, and if so, taking the segment as the segment to be selected. Taking 3 heart beats as an example, the interval trend graph is divided into a plurality of segments. The R peaks of 3 heartbeats in time sequence are R1, R2 and R3 respectively, then the two RR intervals are R1R2 and R2R3 respectively, then the RR interval dispersion of the fragments is s1=
(R2R 3-R1R 2)/R1R 2, wherein R1R2 is the time period between R1 and R2R3 is the time period between R2 and R3. And when S1 is less than or equal to S, the preset fragment is not considered.
Optionally, if the duration of the segment is longer than the preset duration, the segment is used as the segment of the atrial fibrillation to be selected. Taking 3 heart beats as an example, dividing an interval trend graph into a plurality of segments, judging the size relation between the duration time T1 of the segments and the preset time T, taking the segments as the segments to be selected for atrial fibrillation when T1 is more than T, and taking the preset segments not into consideration when T1 is less than or equal to T.
Optionally, if the heart rate value in the segment is greater than the preset heart rate value, the segment is used as the atrial fibrillation segment to be selected. Taking 3 heart beats as an example, dividing the interval trend graph into a plurality of segments, judging the magnitude relation between the heart rate value HRx of the segments and the preset heart rate value HR, and taking the segments as the segments to be selected for atrial fibrillation when the heart rate values HRx of the segments are all larger than HR. And if the heart rate value in the segment is smaller than or equal to the preset heart rate value, taking a heart beat corresponding to the heart rate value as a demarcation point, dividing the segment into a plurality of sub-segments, and taking the sub-segments as the atrial fibrillation segments to be selected when the duration of the sub-segments is longer than the preset duration. Taking 3 beats as an example, dividing the interval trend graph into a plurality of segments, when the heart rate value HRx of the segments is smaller than or equal to HR, dividing the segments into an A1 segment and an A2 segment by taking the heart rate value as a dividing point (a plurality of sub-segments can be divided by taking the heart rate value as the dividing point, HRx of each sub-segment is larger than HR), judging whether the duration T1 of the A1 segment is larger than the preset duration T again, if yes, the A1 segment is a meeting condition segment, judging whether the duration T1 of the A2 segment is larger than the preset duration T, if yes, the A2 segment is also a meeting condition segment, and then both the A1 segment and the A2 segment are taken as the atrial fibrillation segments to be selected.
When at least two detection parameters are used, the detection order of the detection parameters is not limited. For example, if the dispersion and the duration are used for detection, it may be firstly determined whether the dispersion of each segment in the interval trend graph is greater than a preset dispersion, then determine whether the duration of the segment with the dispersion greater than the preset dispersion is greater than a preset duration, or firstly determine whether the duration of each segment in the interval trend graph is greater than a preset duration, and then determine whether the dispersion of the segment with the duration greater than the preset duration is greater than the preset dispersion, that is, the detection order may be sequentially: dispersion-duration, or duration-dispersion. Alternatively, taking 3 heart beats as an example, the interval trend graph is divided into a plurality of segments. The R peaks of 3 heartbeats in time sequence are R1, R2 and R3, respectively, then the two RR intervals are R1R2 and R2R3, respectively, and then the RR interval dispersion of the fragment is s1= (r2r3—r1r2)/R1R 2, where R1R2 is the period between R1 and R2, and R2R3 is the period between R2 and R3. And when S1 is more than S, judging the size relation between the duration time T1 of the segment and the preset time T, when T1 is more than T, taking the segment as the atrial fibrillation segment to be selected, and when T1 is less than or equal to T, taking the preset segment into consideration.
For example, if the dispersion and the heart rate value are used for detection, whether the dispersion of each segment in the interval trend graph is larger than the preset dispersion or not can be judged, and then whether the heart rate value of the segment with the dispersion larger than the preset dispersion is larger than the preset heart rate value or not can be judged, or whether the heart rate value of each segment in the interval trend graph is larger than the preset heart rate value or not can be judged, and then whether the dispersion of the segment with the heart rate value larger than the preset heart rate value is larger than the preset dispersion or not can be judged.
For example, the duration and the heart rate value are used for detection, whether the duration of each segment in the interval trend graph is longer than the preset duration can be judged first, whether the heart rate value of the segment with the duration longer than the preset duration is longer than the preset heart rate value can be judged, whether the heart rate value of each segment in the interval trend graph is longer than the preset heart rate value can be judged first, and whether the duration of the segment with the heart rate value longer than the preset heart rate value is longer than the preset duration can be judged.
For example, using the dispersion, duration and heart rate values, the detection sequence may be: dispersion-duration-heart rate value, dispersion-heart rate value-duration, duration-dispersion-heart rate value, duration-heart rate value-dispersion, heart rate value-dispersion-duration, heart rate value-duration-dispersion.
Optionally, if the dispersion of the segment is greater than a preset dispersion, the duration of the segment is longer than a preset duration, and the heart rate value in the segment is greater than the preset heart rate value, the segment is used as the segment to be selected for atrial fibrillation. Taking 3 heart beats as an example, the interval trend graph is divided into a plurality of segments. The R peaks of 3 heartbeats in time sequence are R1, R2 and R3, respectively, then the two RR intervals are R1R2 and R2R3, respectively, and then the RR interval dispersion of the fragment is s1= (r2r3—r1r2)/R1R 2, where R1R2 is the period between R1 and R2, and R2R3 is the period between R2 and R3. And when S1 is more than S, judging the size relation between the duration time T1 of the segment and the preset time T, when T1 is more than T, judging the size relation between the heart rate value HRx of the segment and the preset heart rate value HR, and when the heart rate values HRx of the segment are all more than HR, taking the segment as the atrial fibrillation segment to be selected.
If the dispersion of the segment is larger than the preset dispersion, the duration of the segment is longer than the preset duration, when the heart rate value in the segment is smaller than or equal to the preset heart rate value, the segment is divided into a plurality of sub-segments by taking the heart beat corresponding to the heart rate value as a demarcation point, and when the duration of the sub-segment is longer than the preset duration, the sub-segment is taken as the atrial fibrillation segment to be selected. Taking 3 heart beats as an example, the interval trend graph is divided into a plurality of segments. The R peaks of 3 heartbeats in time sequence are R1, R2 and R3, respectively, then the two RR intervals are R1R2 and R2R3, respectively, and then the RR interval dispersion of the fragment is s1= (r2r3—r1r2)/R1R 2, where R1R2 is the period between R1 and R2, and R2R3 is the period between R2 and R3. When S1 > S, judging the size relation between the duration T1 of the segment and the preset time T, when T1 > T, judging the size relation between the heart rate value HRx of the segment and the preset heart rate value HR, when the heart rate value HRx of the segment is smaller than or equal to HR, taking the heart rate value as a demarcation point, splitting the segment into an A1 segment and an A2 segment (a plurality of sub-segments can be split by taking the heart rate value as the demarcation point, and HRx of each sub-segment is larger than HR), judging whether the duration T1 of the A1 segment is larger than the preset duration T again, if yes, the A1 segment is a meeting condition segment, judging whether the duration T1 of the A2 segment is larger than the preset duration T, if yes, the A2 segment is also a meeting condition segment, and the A1 segment and the A2 segment are taken as the atrial fibrillation segments to be selected.
In this embodiment, the step of determining the atrial fibrillation segment in the electrocardiographic waterfall map according to the atrial fibrillation segment to be selected includes: acquiring a position range of occurrence of P waves in the electrocardio-fall map, and acquiring an electrocardio-fall picture segment corresponding to the atrial fibrillation segment to be selected in the electrocardio-fall map; and determining an electrocardio-fall picture segment in which the P wave does not appear according to the position range and the electrocardio-fall picture segment, and taking the electrocardio-fall picture segment in which the P wave does not appear as the atrial fibrillation segment. Or analyzing the electrocardio-fall picture segments one by one, and taking the electrocardio-fall picture segments as atrial fibrillation segments if the electrocardio-fall picture segments do not have the graphic characteristic values of the P waves. The graphics characteristic values include, but are not limited to, pixel point distribution, color values, brightness, saturation and the like of the P-wave position area.
It should be noted that, the atrial fibrillation segment to be selected determined according to the interval trend chart may also be directly used as an atrial fibrillation segment, without performing secondary screening of the electrocardiographic waterfall chart.
It should be noted that, the method described in this embodiment is suitable for detecting short Cheng Fang fibrillation, where the atrial fibrillation segment has the characteristics of large RR interval dispersion and large heart rate value, and the short atrial fibrillation has the characteristic of short duration, so that the dispersion, heart rate value and duration are used to screen the atrial fibrillation segment to be selected. The duration of the segment is set to be longer than the preset time in order to avoid that the segment with too short duration causes errors to the detection result.
In the technical scheme disclosed in the embodiment, when the electrocardiographic data is acquired, an interval trend chart and an electrocardiographic waterfall chart are generated according to the electrocardiographic data, a to-be-selected atrial fibrillation fragment is detected in the interval trend chart by taking a preset detection parameter as a detection condition, the atrial fibrillation fragment is analyzed according to the electrocardiographic waterfall chart to obtain the atrial fibrillation fragment, and the atrial fibrillation fragment is marked. Thus, the accuracy of short-range atrial fibrillation detection is improved.
In an embodiment, on the basis of the embodiment shown in fig. 2, the electrocardiographic data processing method further includes:
generating an interval trend graph and an electrocardiographic waterfall graph according to the electrocardiographic data when the electrocardiographic data are acquired;
detecting a to-be-selected atrial fibrillation fragment in the interval trend chart by taking a preset detection parameter as a detection condition;
and analyzing the atrial fibrillation fragments to be selected according to the electrocardio waterfall graph to obtain atrial fibrillation fragments, and marking the atrial fibrillation fragments.
The method described in this embodiment is applicable to the detection of long Cheng Fang fibrillation. And the atrial fibrillation fragments to be selected are preliminarily determined according to the interval trend graph, and the atrial fibrillation fragments are finally determined from the atrial fibrillation fragments to be selected through an electrocardio waterfall graph. The step of determining the selected atrial fibrillation fragments according to the interval trend graph may be: and (3) scanning and analyzing the interval trend graph, and screening out the atrial fibrillation fragments to be selected by taking preset detection parameters as detection conditions. Specifically, the interval trend graph is divided into a plurality of segments by a preset number of heartbeats, whether each segment is a candidate atrial fibrillation segment is sequentially determined according to detection parameters, and the detection parameters comprise at least one of dispersion and heart rate values. Optionally, the preset number is three.
Optionally, acquiring the dispersion of each segment in the interval trend graph, judging whether the dispersion of each segment in the interval trend graph is larger than a preset dispersion, and if so, taking the segment as the segment to be selected. Taking 3 heart beats as an example, the interval trend graph is divided into a plurality of segments. The R peaks of 3 heartbeats in time sequence are R1, R2 and R3, respectively, then the two RR intervals are R1R2 and R2R3, respectively, and then the RR interval dispersion of the fragment is s1= (r2r3—r1r2)/R1R 2, where R1R2 is the period between R1 and R2, and R2R3 is the period between R2 and R3. And when S1 is less than or equal to S, the preset fragment is not considered.
Optionally, if the heart rate value in the segment is greater than the preset heart rate value, the segment is used as the atrial fibrillation segment to be selected. Taking 3 heart beats as an example, dividing the interval trend graph into a plurality of segments, judging the magnitude relation between the heart rate value HRx of the segments and the preset heart rate value HR, and taking the segments as the segments to be selected for atrial fibrillation when the heart rate values HRx of the segments are all larger than HR. And if the heart rate value in the segment is smaller than or equal to the preset heart rate value, taking a heart beat corresponding to the heart rate value as a demarcation point, dividing the segment into a plurality of sub-segments, and taking the sub-segments as the atrial fibrillation segments to be selected when the duration of the sub-segments is longer than the preset duration. Taking 3 beats as an example, dividing the interval trend graph into a plurality of segments, when the heart rate value HRx of the segments is smaller than or equal to HR, dividing the segments into an A1 segment and an A2 segment by taking the heart rate value as a dividing point (a plurality of sub-segments can be divided by taking the heart rate value as the dividing point, HRx of each sub-segment is larger than HR), judging whether the duration T1 of the A1 segment is larger than the preset duration T again, if yes, the A1 segment is a meeting condition segment, judging whether the duration T1 of the A2 segment is larger than the preset duration T, if yes, the A2 segment is also a meeting condition segment, and then both the A1 segment and the A2 segment are taken as the atrial fibrillation segments to be selected.
When at least two detection parameters are used, the detection order of the detection parameters is not limited. For example, if the dispersion and the heart rate value are used for detection, it may be firstly determined whether the dispersion of each segment in the interval trend graph is greater than a preset dispersion, then it may be determined whether the heart rate value of the segment with the dispersion greater than the preset dispersion is greater than a preset heart rate value, or it may be firstly determined whether the heart rate value of each segment in the interval trend graph is greater than the preset heart rate value, and then it may be determined whether the dispersion of the segment with the heart rate value greater than the preset heart rate value is greater than the preset dispersion, that is, the detection order may be sequentially: dispersion-heart rate value, also can be heart rate value-dispersion.
Alternatively, taking 3 heart beats as an example, the interval trend graph is divided into a plurality of segments. The R peaks of 3 heartbeats in time sequence are R1, R2 and R3, respectively, then the two RR intervals are R1R2 and R2R3, respectively, and then the RR interval dispersion of the fragment is s1= (r2r3—r1r2)/R1R 2, where R1R2 is the period between R1 and R2, and R2R3 is the period between R2 and R3. And when S1 is more than S, judging the magnitude relation between the heart rate value HRx of the segment and the preset heart rate value HR, and when the heart rate values HRx of the segments are all larger than HR, taking the segment as the segment of the atrial fibrillation to be selected.
If the dispersion of the segments is larger than the preset dispersion, dividing the segments into a plurality of sub-segments by taking heart beats corresponding to the heart rate values as demarcation points when the heart rate values in the segments are smaller than or equal to the preset heart rate values, and taking the sub-segments as the segments to be selected for atrial fibrillation when the duration of the sub-segments is longer than the preset duration. Taking 3 heart beats as an example, the interval trend graph is divided into a plurality of segments. The R peaks of 3 heartbeats in time sequence are R1, R2 and R3, respectively, then the two RR intervals are R1R2 and R2R3, respectively, and then the RR interval dispersion of the fragment is s1= (r2r3—r1r2)/R1R 2, where R1R2 is the period between R1 and R2, and R2R3 is the period between R2 and R3. When S1 > S, judging the magnitude relation between the heart rate value HRx of the segment and the preset heart rate value HR, when the heart rate value HRx of the segment is smaller than or equal to HR, taking the heart rate value as a demarcation point, splitting the segment into an A1 segment and an A2 segment (a plurality of sub-segments can be split by taking the heart rate value as the demarcation point, and HRx of each sub-segment is larger than HR), judging whether the duration T1 of the A1 segment is larger than the preset duration T, if yes, the A1 segment is a coincidence segment, and if not, the duration T1 of the A2 segment is larger than the preset duration T, if yes, the A2 segment is also a coincidence segment, and then the A1 segment and the A2 segment are taken as the atrial fibrillation segments to be selected.
In this embodiment, the step of determining the atrial fibrillation segment in the electrocardiographic waterfall map according to the atrial fibrillation segment to be selected includes: acquiring a position range of occurrence of P waves in the electrocardio-fall map, and acquiring an electrocardio-fall picture segment corresponding to the atrial fibrillation segment to be selected in the electrocardio-fall map; and determining an electrocardio-fall picture segment in which the P wave does not appear according to the position range and the electrocardio-fall picture segment, and taking the electrocardio-fall picture segment in which the P wave does not appear as the atrial fibrillation segment. Or analyzing the electrocardio-fall picture segments one by one, and taking the electrocardio-fall picture segments as atrial fibrillation segments if the electrocardio-fall picture segments do not have the graphic characteristic values of the P waves. The graphics characteristic values include, but are not limited to, pixel point distribution, color values, brightness, saturation and the like of the P-wave position area.
It should be noted that, the atrial fibrillation segment to be selected determined according to the interval trend chart may also be directly used as an atrial fibrillation segment, without performing secondary screening of the electrocardiographic waterfall chart.
In the technical scheme disclosed in the embodiment, when the electrocardiographic data is acquired, an interval trend chart and an electrocardiographic waterfall chart are generated according to the electrocardiographic data, a to-be-selected atrial fibrillation fragment is detected in the interval trend chart by taking a preset detection parameter as a detection condition, the atrial fibrillation fragment is analyzed according to the electrocardiographic waterfall chart to obtain the atrial fibrillation fragment, and the atrial fibrillation fragment is marked. In this way, the accuracy of the long Cheng Fang fibrillation detection is improved.
In an embodiment, as shown in fig. 3, based on the embodiment shown in fig. 2, the method for processing electrocardiographic data further includes, after the steps of analyzing the atrial fibrillation segment according to the electrocardiographic waterfall diagram to obtain an atrial fibrillation segment and marking the atrial fibrillation segment:
step S40, acquiring parameter information corresponding to preset parameters of the atrial fibrillation fragments in the electrocardiographic waterfall, wherein the parameter information comprises at least one of starting time, duration time, average dispersion and average heart rate value of the atrial fibrillation fragments;
and S50, outputting the parameter information.
In this embodiment, the atrial fibrillation fragments are marked, saved, and recorded, and at the same time, parameter information corresponding to preset parameters of the atrial fibrillation fragments is presented on an interface in a list form, where the parameter information includes at least one of a start time, a duration time, an average dispersion, and an average heart rate value of the atrial fibrillation fragments. The system generates and outputs a scatter plot, a histogram, a trend plot, etc., and further generates and outputs an analysis report.
In the technical scheme disclosed in the embodiment, after the atrial fibrillation segment is acquired, the parameter information of the atrial fibrillation segment is output, so that the analysis by a doctor is facilitated.
The invention also provides an electrocardio data processing device, which comprises an electrocardio data processing program, wherein the electrocardio data processing program is configured to realize the steps of the electrocardio data processing method under the main execution body.
The embodiment of the invention also provides a computer readable storage medium, wherein the computer readable storage medium stores a processing program of the electrocardio data, and the processing program of the electrocardio data is executed by a processor to realize the steps of the method for processing the electrocardio data under the main execution body by the air conditioner or the server.
The foregoing embodiment numbers of the present invention are merely for the purpose of description, and do not represent the advantages or disadvantages of the embodiments.
From the above description of the embodiments, it will be clear to those skilled in the art that the above-described embodiment method may be implemented by means of software plus a necessary general hardware platform, but of course may also be implemented by means of hardware, but in many cases the former is a preferred embodiment. Based on such understanding, the technical solution of the present invention may be embodied essentially or in a part contributing to the prior art in the form of a software product stored in a storage medium (e.g. ROM/RAM, magnetic disk, optical disk) as described above, comprising instructions for causing a terminal device (which may be a television, a mobile phone, a computer, a server, an air conditioner, or a network device, etc.) to perform the method according to the embodiments of the present invention.
The foregoing description is only of the preferred embodiments of the present invention, and is not intended to limit the scope of the invention, but rather is intended to cover any equivalents of the structures or equivalent processes disclosed herein or in the alternative, which may be employed directly or indirectly in other related arts.
Claims (6)
1. An electrocardiographic data processing device, comprising a memory, a processor and an electrocardiographic data processing program stored in the memory and capable of running on the processor, wherein the electrocardiographic data processing program, when executed by the processor, realizes an electrocardiographic data processing method, and the electrocardiographic data processing method comprises the following steps:
generating an interval trend graph and an electrocardiographic waterfall graph according to the electrocardiographic data when the electrocardiographic data are acquired;
detecting a to-be-selected atrial fibrillation fragment in the interval trend chart by taking a preset detection parameter as a detection condition;
analyzing the atrial fibrillation fragments to be selected according to the electrocardio waterfall graph to obtain atrial fibrillation fragments, and marking the atrial fibrillation fragments;
wherein the detection parameters comprise duration and heart rate values, and the step of detecting the segments of atrial fibrillation to be selected in the interval trend graph comprises:
dividing the interval trend graph into a plurality of segments with a preset number of heartbeats;
judging whether the heart rate value in the segment is larger than a preset heart rate value or not when the duration time of the segment is longer than a preset time;
when the heart rate value of the segment is larger than the preset heart rate value, the segment is used as the atrial fibrillation segment to be selected;
when the heart rate value of the segment is smaller than or equal to a preset heart rate value, dividing the segment into a plurality of sub-segments by taking a heart beat corresponding to the heart rate value as a demarcation point;
and when the duration time of the sub-segment is longer than the preset duration time, taking the sub-segment as the atrial fibrillation segment to be selected.
2. The apparatus for processing electrocardiographic data according to claim 1, wherein the step of analyzing the segments of atrial fibrillation to be selected based on the electrocardiographic waterfall map to obtain segments of atrial fibrillation comprises:
acquiring a position range of occurrence of P waves in the electrocardio-fall map, and acquiring an electrocardio-fall picture segment corresponding to the atrial fibrillation segment to be selected in the electrocardio-fall map;
and determining an electrocardio-fall picture segment in which the P wave does not appear according to the position range and the electrocardio-fall picture segment, and taking the electrocardio-fall picture segment in which the P wave does not appear as the atrial fibrillation segment.
3. The apparatus for processing electrocardiographic data according to any one of claims 1-2, wherein after the step of analyzing the selected atrial fibrillation fragments according to the electrocardiographic waterfall map to obtain atrial fibrillation fragments and marking the atrial fibrillation fragments, the method for processing electrocardiographic data further comprises:
acquiring parameter information corresponding to preset parameters of the atrial fibrillation fragments from the electrocardiographic waterfall map, wherein the parameter information comprises at least one of starting time, duration time, average dispersion and average heart rate value of the atrial fibrillation fragments;
and outputting the parameter information.
4. A computer-readable storage medium, wherein a processing program of electrocardiographic data is stored on the computer-readable storage medium, and the processing program of electrocardiographic data implements an electrocardiographic data processing method when executed by a processor, the method comprising the steps of:
generating an interval trend graph and an electrocardiographic waterfall graph according to the electrocardiographic data when the electrocardiographic data are acquired;
detecting a to-be-selected atrial fibrillation fragment in the interval trend chart by taking a preset detection parameter as a detection condition;
analyzing the atrial fibrillation fragments to be selected according to the electrocardio waterfall graph to obtain atrial fibrillation fragments, and marking the atrial fibrillation fragments;
wherein the detection parameters comprise duration and heart rate values, and the step of detecting the segments of atrial fibrillation to be selected in the interval trend graph comprises:
dividing the interval trend graph into a plurality of segments with a preset number of heartbeats;
judging whether the heart rate value in the segment is larger than a preset heart rate value or not when the duration time of the segment is longer than a preset time;
when the heart rate value of the segment is larger than the preset heart rate value, the segment is used as the atrial fibrillation segment to be selected;
when the heart rate value of the segment is smaller than or equal to a preset heart rate value, dividing the segment into a plurality of sub-segments by taking a heart beat corresponding to the heart rate value as a demarcation point;
and when the duration time of the sub-segment is longer than the preset duration time, taking the sub-segment as the atrial fibrillation segment to be selected.
5. The storage medium of claim 4, wherein the step of analyzing the selected atrial fibrillation fragments based on the electrocardiographic waterfall map to obtain atrial fibrillation fragments comprises:
acquiring a position range of occurrence of P waves in the electrocardio-fall map, and acquiring an electrocardio-fall picture segment corresponding to the atrial fibrillation segment to be selected in the electrocardio-fall map;
and determining an electrocardio-fall picture segment in which the P wave does not appear according to the position range and the electrocardio-fall picture segment, and taking the electrocardio-fall picture segment in which the P wave does not appear as the atrial fibrillation segment.
6. The storage medium of any one of claims 4-5, wherein after the step of analyzing the candidate atrial fibrillation fragments according to the electrocardiographic waterfall map to obtain atrial fibrillation fragments and marking the atrial fibrillation fragments, the method for processing electrocardiographic data further comprises:
acquiring parameter information corresponding to preset parameters of the atrial fibrillation fragments from the electrocardiographic waterfall map, wherein the parameter information comprises at least one of starting time, duration time, average dispersion and average heart rate value of the atrial fibrillation fragments;
and outputting the parameter information.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910048588.0A CN111449649B (en) | 2019-01-18 | 2019-01-18 | Electrocardiogram data processing method and device and computer readable storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910048588.0A CN111449649B (en) | 2019-01-18 | 2019-01-18 | Electrocardiogram data processing method and device and computer readable storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111449649A CN111449649A (en) | 2020-07-28 |
| CN111449649B true CN111449649B (en) | 2023-07-11 |
Family
ID=71671505
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910048588.0A Active CN111449649B (en) | 2019-01-18 | 2019-01-18 | Electrocardiogram data processing method and device and computer readable storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111449649B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116685268A (en) * | 2020-11-30 | 2023-09-01 | 深圳迈瑞生物医疗电子股份有限公司 | Monitoring system, method for comprehensive management of atrial fibrillation and method for displaying monitored data |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2008137539A1 (en) * | 2007-05-01 | 2008-11-13 | Medtronic, Inc. | Method and apparatus for detecting arrhythmias in a medical device |
| CN104367318A (en) * | 2014-11-07 | 2015-02-25 | 北京海思敏医疗技术有限公司 | Detection method and device for atrial fibrillation |
| CN107296604A (en) * | 2017-08-29 | 2017-10-27 | 心云（北京）医疗器械有限公司 | A kind of atrial fibrillation determination methods |
| CN108577823A (en) * | 2018-04-27 | 2018-09-28 | 京东方科技集团股份有限公司 | A kind of arrhythmia detection device and arrhythmia detection system |
| CN109044338A (en) * | 2018-08-08 | 2018-12-21 | 广州视源电子科技股份有限公司 | Atrial fibrillation detection apparatus and storage medium |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8326407B2 (en) * | 2004-10-19 | 2012-12-04 | University Of Washington | Long-term monitoring for discrimination of different heart rhythms |
| US8233973B2 (en) * | 2008-05-02 | 2012-07-31 | Spacelabs Healthcare, Llc | Methods for detection of cardiac arrhythmias |
| US9901276B2 (en) * | 2015-02-18 | 2018-02-27 | Medtronic, Inc. | Method and apparatus for identifying sick sinus syndrome in an implantable cardiac monitoring device |
| US10631744B2 (en) * | 2016-04-13 | 2020-04-28 | Cardiac Pacemakers, Inc. | AF monitor and offline processing |
| CN106264518A (en) * | 2016-08-04 | 2017-01-04 | 东莞市第三人民医院 | A kind of based on finger tip Electrocardiographic atrial fibrillation detection method and device |
| US10744334B2 (en) * | 2016-10-18 | 2020-08-18 | Cardiac Pacemakers, Inc. | Systems and methods for arrhythmia detection |
-
2019
- 2019-01-18 CN CN201910048588.0A patent/CN111449649B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2008137539A1 (en) * | 2007-05-01 | 2008-11-13 | Medtronic, Inc. | Method and apparatus for detecting arrhythmias in a medical device |
| CN104367318A (en) * | 2014-11-07 | 2015-02-25 | 北京海思敏医疗技术有限公司 | Detection method and device for atrial fibrillation |
| CN107296604A (en) * | 2017-08-29 | 2017-10-27 | 心云（北京）医疗器械有限公司 | A kind of atrial fibrillation determination methods |
| CN108577823A (en) * | 2018-04-27 | 2018-09-28 | 京东方科技集团股份有限公司 | A kind of arrhythmia detection device and arrhythmia detection system |
| CN109044338A (en) * | 2018-08-08 | 2018-12-21 | 广州视源电子科技股份有限公司 | Atrial fibrillation detection apparatus and storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111449649A (en) | 2020-07-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108186011B (en) | Atrial fibrillation detection method, atrial fibrillation detection device and readable storage medium | |
| AU2009243424B2 (en) | ECG signal analysis tool | |
| US10682118B2 (en) | Ultrasound system and method for analyzing cardiac periodicity | |
| CN111433860B (en) | User interface for analyzing an electrocardiogram | |
| CN109310354B (en) | Method and device for presenting ST events in electrocardiogram | |
| CN110464333B (en) | Method and device for storing electrocardiogram data | |
| CN110393520B (en) | Electrocardiogram data analysis method and device, electrocardiogram analyzer and storage medium | |
| US20090275850A1 (en) | Electrocardiographic (ECG) Data Analysis Systems and Methods | |
| CN110495872A (en) | Ecg analysis method, apparatus, equipment and medium based on picture and heart-beat messages | |
| US20160183832A1 (en) | Heart Electrophysiological Signal Analysis System | |
| CN110448295B (en) | Electrocardiographic waveform editing and displaying method and device and computer equipment | |
| CN111449649B (en) | Electrocardiogram data processing method and device and computer readable storage medium | |
| CN111150386A (en) | Method, device and storage medium for electrocardio analysis | |
| CN118021316B (en) | Dynamic electrocardiographic data processing method and device, electronic equipment and storage medium | |
| US9386935B2 (en) | Electrocardiogram analyzer | |
| US11452476B2 (en) | Detection report data generation method | |
| CN110545713B (en) | Physiological information waveform processing method, program, computer-readable storage medium, and physiological information waveform processing device | |
| CN113647960B (en) | Arrhythmia rapid identification method, terminal equipment and storage medium | |
| CN109431488B (en) | Method for identifying typical data in electrocardiogram dynamic real-time analysis data | |
| CN110916649B (en) | Processing device, processing method and detection device for long-range electrocardiogram scatter diagram | |
| CN112120688B (en) | Electrocardiosignal processing method, equipment and computer readable storage medium | |
| CN116807490A (en) | Electrocardiogram processing method, electrocardiogram processing device, medical equipment and storage medium | |
| CN111134656B (en) | Electrocardiogram data detection method and device and computer readable storage medium | |
| CN111493856B (en) | Correction method and device for heart beat classification, electrocardiogram checking system and storage medium | |
| CN112353400A (en) | Heart beat classification method and device, computer equipment and storage medium |
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 |