CN110123297B - Heart rate measuring method and device, computer equipment and storage medium - Google Patents
Heart rate measuring method and device, computer equipment and storage medium Download PDFInfo
- Publication number
- CN110123297B CN110123297B CN201910268887.5A CN201910268887A CN110123297B CN 110123297 B CN110123297 B CN 110123297B CN 201910268887 A CN201910268887 A CN 201910268887A CN 110123297 B CN110123297 B CN 110123297B
- Authority
- CN
- China
- Prior art keywords
- signal point
- value
- maximum
- heart rate
- observed
- 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/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/024—Detecting, measuring or recording pulse rate or heart rate
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Cardiology (AREA)
- Biomedical Technology (AREA)
- Medical Informatics (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Physiology (AREA)
- Heart & Thoracic Surgery (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
The invention discloses a heart rate measuring method, a heart rate measuring device, computer equipment and a storage medium, and relates to the field of artificial intelligence. The heart rate measuring method comprises the following steps: acquiring a target amplitude spectrum according to the heart rate measurement signal, wherein the target amplitude spectrum comprises signal points, and coordinate values of the signal points comprise amplitude values and frequency values; acquiring a maximum signal point in the signal points according to the amplitude value and the frequency value of the signal points; acquiring a preset number of maximum signal points to be observed from the maximum signal points, wherein the amplitude value of the maximum signal points to be observed is larger than the amplitude values of other residual maximum signal points; searching a target maximum signal point for dominating heart rate measurement in the maximum signal point to be observed according to the amplitude value of the maximum signal point to be observed; and obtaining a target heart rate value according to the target maximum signal point. By adopting the heart rate measuring method, the heart rate can be effectively measured in actual processing, and the accuracy of heart rate measurement is ensured.
Description
[ technical field ] A method for producing a semiconductor device
The invention relates to the field of artificial intelligence, in particular to a heart rate measuring method and device, computer equipment and a storage medium.
[ background of the invention ]
The heart rate measurement can be realized by filtering an original signal by a filter, performing Fourier transform on the filtered signal to obtain an amplitude spectrum, and then using the frequency corresponding to the amplitude peak value of the amplitude spectrum as the heart rate. However, in actual processing, there are still some abnormal peak points in the amplitude spectrum due to insufficient quality of the acquired signal and excessive noise, so that the measurer cannot accurately measure the heart rate, and a heart rate measurement error may occur.
[ summary of the invention ]
In view of this, embodiments of the present invention provide a method and an apparatus for heart rate measurement, a computer device, and a storage medium, so as to solve the problem that heart rate measurement cannot be accurately performed in actual processing at present.
In a first aspect, an embodiment of the present invention provides a heart rate measurement method, including:
acquiring a target amplitude spectrum according to the heart rate measurement signal, wherein the target amplitude spectrum comprises signal points, and the coordinate values of the signal points comprise an amplitude value and a frequency value;
acquiring a maximum signal point in the signal points according to the amplitude value and the frequency value of the signal points;
acquiring a preset number of maximum signal points to be observed from the maximum signal points, wherein the amplitude value of the maximum signal points to be observed is larger than the amplitude values of other residual maximum signal points;
searching a target maximum signal point for dominating heart rate measurement in the maximum signal point to be observed according to the amplitude value of the maximum signal point to be observed;
and obtaining a target heart rate value according to the target maximum signal point.
The above aspect and any possible implementation manner further provide an implementation manner, where acquiring a maximum signal point from the signal points according to the amplitude value and the frequency value of the signal point includes:
searching a target signal point in a preset frequency value interval, wherein the amplitude value of the target signal point is larger than the amplitude values of other residual signal points;
if the frequency value is found, judging whether a first signal point and a second signal point exist in the preset frequency value interval, wherein the frequency value of the first signal point is larger than the frequency value of the target signal point, and the frequency value of the second signal point is smaller than the frequency value of the target signal point;
and if the first signal point and the second signal point exist in the preset frequency value interval, determining the target signal point as a maximum signal point and acquiring the maximum signal point.
The above aspect and any possible implementation further provide an implementation in which, according to the amplitude value of the maximum signal point to be observed, finding a target maximum signal point that dominates a heart rate measurement among the maximum signal points to be observed includes:
acquiring a maximum value signal point to be observed with the largest amplitude value as a candidate maximum value signal point according to the amplitude value of the maximum value signal point to be observed;
respectively calculating the difference value between the candidate maximum signal point and the amplitude value of each of the other remaining maximum signal points to be observed, and calculating the absolute value of each difference value;
and if the absolute value of each difference value is larger than a first preset value, determining the candidate maximum signal point as a target maximum signal point for leading heart rate measurement.
The above aspect and any possible implementation further provide an implementation in which the finding, from the amplitude values of the maximum signal points to be observed, a target maximum signal point that dominates a heart rate measurement among the maximum signal points to be observed includes:
screening a third signal point set from the maximum signal points to be observed according to the amplitude values of the maximum signal points to be observed, wherein the third signal point set comprises a plurality of third signal points, and the absolute values of the difference values of the third signal points and the amplitude values of other remaining maximum signal points to be observed are all larger than a second preset value;
determining a fourth signal point set in the third signal point set, wherein the fourth signal point set comprises a plurality of fourth signal points, and the frequency values of the fourth signal points are in a multiple relation, and the fourth signal point with the smallest frequency value is taken as a reference in the multiple relation;
and if the amplitude value of the fourth signal point is reduced along with the increase of the frequency value or increased along with the increase of the frequency value, determining the fourth signal point with the maximum amplitude value as a target maximum value signal point for leading the heart rate measurement.
The above aspect and any possible implementation further provide an implementation in which the finding, from the amplitude values of the maximum signal points to be observed, a target maximum signal point that dominates a heart rate measurement among the maximum signal points to be observed includes:
screening a fifth signal point set from the maximum signal points to be observed according to the amplitude values of the maximum signal points to be observed, wherein the fifth signal point set comprises a plurality of fifth signal points, and the absolute values of the difference values of the fifth signal points and the amplitude values of other remaining maximum signal points to be observed are all larger than a third preset value;
determining a sixth signal point set in the fifth signal points, wherein the sixth signal point set comprises a plurality of sixth signal points, and the absolute value of the difference between the frequency values of any two adjacent sixth signal points is smaller than a fourth preset value;
and if the number of the sixth signal points is within a preset numerical value interval, determining the sixth signal points as target maximum signal points.
The above aspects and any possible implementations further provide an implementation in which the deriving a target heart rate value from the target maximum signal points includes:
and calculating the mean value of the amplitude values of the sixth signal points in the sixth signal point set, and determining the mean value of the amplitude values as the target heart rate value.
In a second aspect, an embodiment of the present invention provides a heart rate measuring apparatus, including:
the target amplitude spectrum acquisition module is used for acquiring a target amplitude spectrum according to the heart rate measurement signal, the target amplitude spectrum comprises signal points, and the coordinate values of the signal points comprise amplitude values and frequency values;
the maximum signal point acquisition module is used for acquiring a maximum signal point in the signal points according to the amplitude value and the frequency value of the signal points;
the device comprises a to-be-observed maximum signal point acquisition module, a to-be-observed maximum signal point acquisition module and a to-be-observed maximum signal point acquisition module, wherein the to-be-observed maximum signal points are acquired from a preset number of the maximum signal points, and the amplitude value of each to-be-observed maximum signal point is larger than the amplitude values of other residual maximum signal points;
the searching module is used for searching a target maximum signal point which is used for leading heart rate measurement in the maximum signal point to be observed according to the amplitude value of the maximum signal point to be observed;
and the target heart rate value acquisition module is used for acquiring a target heart rate value according to the target maximum signal point.
In a third aspect, a computer device comprises a memory, a processor and a computer program stored in the memory and executable on the processor, the processor implementing the steps of the heart rate measurement method described above when executing the computer program.
In a fourth aspect, the present invention provides a computer-readable storage medium, which includes a computer program, and when the computer program is executed by a processor, the computer program implements the steps of the heart rate measurement method described above.
In the embodiment of the invention, firstly, a target amplitude spectrum is obtained according to a heart rate measurement signal, the target amplitude spectrum comprises signal points, coordinate values of the signal points comprise an amplitude value and a frequency value, then, a maximum signal point in the signal points is obtained according to the amplitude value and the frequency value of the signal points, and the maximum signal point on the target amplitude spectrum is analyzed to provide a technical premise for realizing the subsequent obtaining of a target heart rate value; then, a preset number of maximum signal points to be observed are obtained from the maximum signal points, wherein the amplitude value of the maximum signal point to be observed is larger than the amplitude values of other remaining maximum signal points, so that the analysis of the target amplitude spectrum can be simplified to be the analysis of the maximum preset number of maximum signal points, and the efficiency and the feasibility of implementation of heart rate measurement can be improved; and then searching a target maximum signal point for leading heart rate measurement in the maximum signal points to be observed according to the amplitude value of the maximum signal point to be observed, analyzing from the potential relation between the maximum signal points to be observed based on the amplitude value, searching the target maximum signal point for leading heart rate measurement, and obtaining a target heart rate value according to the target maximum signal point, thereby realizing accurate heart rate measurement in actual processing.
[ description of the drawings ]
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings needed to be used in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings based on these drawings without inventive labor.
FIG. 1 is a flow chart of a method of measuring heart rate according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of searching for a target maximum signal point according to an embodiment of the present invention;
FIG. 3 is another schematic diagram of finding a target maximum signal point in accordance with an embodiment of the present invention;
FIG. 4 is a schematic diagram of finding a target maximum signal point according to an embodiment of the present invention;
FIG. 5 is a schematic view of a heart rate measuring device in accordance with an embodiment of the present invention;
FIG. 6 is a schematic diagram of a computing device according to an embodiment of the invention.
[ detailed description ] A
In order to better understand the technical scheme of the invention, the following detailed description of the embodiments of the invention is made with reference to the accompanying drawings.
It should be understood that the described embodiments are only some embodiments of the invention, and not all embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The terminology used in the embodiments of the invention is for the purpose of describing particular embodiments only and is not intended to be limiting of the invention. As used in the examples of the present invention and the appended claims, the singular forms "a", "an", and "the" are intended to include the plural forms as well, unless the context clearly indicates otherwise.
It should be understood that the term "and/or" as used herein is merely a field that describes the same of an associated object, meaning that three relationships may exist, e.g., A and/or B, may represent: a exists alone, A and B exist simultaneously, and B exists alone. In addition, the character "/" herein generally indicates that the former and latter related objects are in an "or" relationship.
It should be understood that although the terms first, second, third, etc. may be used to describe preset ranges, etc. in embodiments of the present invention, these preset ranges should not be limited to these terms. These terms are only used to distinguish preset ranges from each other. For example, the first preset range may also be referred to as a second preset range, and similarly, the second preset range may also be referred to as the first preset range, without departing from the scope of the embodiments of the present invention.
The word "if," as used herein, may be interpreted as "at \8230; \8230when" or "when 8230; \823030when" or "in response to a determination" or "in response to a detection", depending on the context. Similarly, the phrases "if determined" or "if detected (a stated condition or event)" may be interpreted as "when determined" or "in response to a determination" or "when detected (a stated condition or event)" or "in response to a detection (a stated condition or event)", depending on the context.
Fig. 1 shows a flowchart of the heart rate measuring method of the present embodiment. The heart rate measuring method can be applied to a heart rate measuring system, and can be realized by adopting the heart rate measuring system when heart rate measurement is carried out. The heart rate measurement system can be applied to computer equipment, wherein the computer equipment can perform human-computer interaction with a user, and the computer equipment comprises, but is not limited to, computers, smart phones, tablets and other equipment. As shown in fig. 1, the heart rate measuring method includes the steps of:
s10: and acquiring a target amplitude spectrum according to the heart rate measurement signal, wherein the target amplitude spectrum comprises signal points, and the coordinate values of the signal points comprise amplitude values and frequency values.
In the frequency domain description of the signal, frequency values are used as independent variables, and amplitude values of frequency value components constituting the signal are used as dependent variables, and such a frequency function is called an amplitude spectrum (also called an amplitude spectrum) and characterizes the distribution of the amplitude values of the signal along with the frequency. The target amplitude spectrum is the amplitude spectrum to be analyzed, which is related to the heart rate measurement. It is understood that the target amplitude spectrum may be obtained by acquiring an original signal obtained by measuring a heart rate, filtering the original signal by using a filter, and performing fourier transform on the filtered signal. The target amplitude spectrum comprises signal points on a frequency function, each signal point is provided with a corresponding frequency value and an amplitude value, and the distribution condition of the amplitude value along with the frequency value at a specific moment is represented. In this embodiment, the target amplitude spectrum is obtained according to the heart rate measurement signal, so that the target amplitude spectrum is analyzed based on the distribution of the amplitude value along with the frequency value in the target amplitude spectrum, and the accuracy of heart rate measurement is improved.
S20: and acquiring a maximum signal point in the signal points according to the amplitude value and the frequency value of the signal points.
The maximum value is that in a certain minimum interval, an independent variable value x exists in the function, and independent variables larger than the independent variable value x and smaller than the independent variable value x exist in the function, and function values corresponding to the independent variables are all smaller than the function values corresponding to the independent variables x, so that the function values are called as the maximum value. "maximum" is a concept of locality. In the present embodiment, the maximum value signal point indicates a signal point having a maximum value property.
It can be understood that, when measuring the heart rate, the measuring method is to use the frequency value corresponding to the amplitude peak of the amplitude spectrum (the maximum signal point is the signal point with the amplitude peak) as the heart rate, which is not problematic in that only one amplitude peak in the amplitude spectrum is measured in an ideal state, whereas in an actual situation, a plurality of amplitude peaks may occur in the amplitude spectrum at the same time. In this embodiment, the plurality of amplitude peaks may be represented by acquiring a maximum signal point on the target amplitude spectrum according to the amplitude value and the frequency value of the signal point, and further analyzing the target amplitude spectrum based on the maximum signal point.
Further, in step S20, acquiring a maximum signal point in the signal points according to the amplitude value and the frequency value of the signal point, specifically including:
s21: and searching a target signal point in a preset frequency value interval, wherein the amplitude value of the target signal point is larger than the amplitude values of other residual signal points.
In one embodiment, the target signal point can be found by the magnitude of the amplitude value of the signal point within the preset frequency value interval, and when the amplitude value of one signal point is larger than the amplitude values of other residual signal points, the signal point is taken as the target signal point.
S22: if the frequency value is found, whether a first signal point and a second signal point exist in the preset frequency value interval is judged, the frequency value of the first signal point is larger than that of the target signal point, and the frequency value of the second signal point is smaller than that of the target signal point.
In an embodiment, after the target signal point is found, the frequency function is removed from the monotonic function, so as to determine whether the target signal point is a maximum signal point.
S23: and if the first signal point and the second signal point exist in the preset frequency value interval, determining the target signal point as a maximum signal point and acquiring the maximum signal point.
In steps S21-S23, an embodiment of obtaining a maximum signal point is provided, which can quickly obtain a maximum signal point in the target amplitude spectrum. In addition, the maximum signal point can be obtained by performing primary derivation, secondary derivation and the like on the target amplitude spectrum.
S30: and acquiring a preset number of maximum signal points to be observed from the maximum signal points, wherein the amplitude value of the maximum signal points to be observed is larger than the amplitude values of other residual maximum signal points.
The maximum signal points to be observed are maximum signal points for observation, and the number of the maximum signal points to be observed can be preset and determined. Particularly, the maximum signal point to be observed has a limiting condition, that is, the amplitude value of the maximum signal point to be observed is larger than the amplitude values of the other remaining maximum signal points, that is, the obtained maximum signal points to be observed are all maximum values of the first several (preset number) rows of amplitude values, and if the preset number is 5, the obtained amplitude value of the maximum signal point to be observed is the maximum signal point of the first 5 rows of amplitude values in the target amplitude spectrum.
In an embodiment, the preset number may be specifically 5, the preset number of maximum signal points to be observed is obtained from the maximum signal points, which indicates that no matter how many amplitude peak values exist in the target amplitude spectrum, only the maximum signal points of the first 5 amplitude peak values with the highest amplitude values are taken for analysis, and the rest removal is not considered, so that the analysis of the target amplitude spectrum can be simplified to analyze the maximum signal points of the first 5 amplitude peak values with the highest amplitude values, and the efficiency of heart rate measurement and the feasibility of implementation are effectively improved.
S40: and searching a target maximum value signal point for leading heart rate measurement in the maximum value signal points to be observed according to the amplitude value of the maximum value signal points to be observed.
It can be understood that although there are a plurality of maximum signal points to be observed on the target amplitude spectrum, the heart rate value cannot be accurately detected from the ideal heart rate detection perspective, the heart rate detection can be implemented from the perspective of "whether there is a target maximum signal point that dominates the heart rate measurement in the maximum signal points to be observed", and the target maximum signal point that dominates the heart rate measurement. The angle of the target maximum signal point of the dominant heart rate measurement is analyzed from the potential relation between the maximum signal points to be observed based on the amplitude value, the target maximum signal point of the dominant heart rate measurement is searched in the maximum signal points to be observed, and the heart rate measurement can be realized on the premise that the accuracy is guaranteed.
Further, before step S40, that is, before searching for a target maximum signal point of the dominant heart rate measurement in the maximum signal points to be observed according to the amplitude values of the maximum signal points to be observed, the method further includes:
s411: and acquiring the maximum amplitude value and the minimum amplitude value in the maximum signal points to be observed.
S412: and judging whether the absolute value of the difference value between the maximum amplitude value and the minimum amplitude value in the maximum signal point to be observed is smaller than a preset reference value or not.
S413: and if the maximum value signal points are smaller than the target maximum value signal points, determining that the target maximum value signal points which are used for leading the heart rate measurement do not exist in the maximum value signal points to be observed.
It can be understood that, if the difference between the maximum amplitude value and the minimum amplitude value in the maximum signal points to be observed is smaller than the preset reference value, a plurality of maximum signal points to be observed with similar amplitude values are displayed on the target amplitude spectrum, the amplitude peak values corresponding to the maximum signal points to be observed on the target amplitude spectrum are relatively close, and for this situation, it is considered that none of the maximum signal points to be observed is in the position of the dominant heart rate measurement, and the user is directly reminded to measure the heart rate again for an invalid measurement.
In steps S411 to S413, a determination is performed in advance before searching for the target maximum signal point for the dominant heart rate measurement, and whether the measurement is invalid is determined by comparing whether the absolute value of the difference between the maximum amplitude value and the minimum amplitude value in the maximum signal point to be observed is smaller than a preset reference value. The embodiment provides a specific implementation method for determining that no target maximum signal point for dominant heart rate measurement exists in maximum signal points to be observed, and the efficiency of heart rate detection can be improved.
Further, in step S40, according to the amplitude value of the maximum signal point to be observed, a target maximum signal point that dominates the heart rate measurement is searched for in the maximum signal point to be observed, which specifically includes:
s421: and acquiring the maximum value signal point to be observed with the largest amplitude value as a candidate maximum value signal point according to the amplitude value of the maximum value signal point to be observed.
S422: and respectively calculating the difference value between the candidate maximum signal point and the amplitude value of each of the other residual maximum signal points to be observed, and calculating the absolute value of each difference value.
S423: and if the absolute value of each difference value is larger than the first preset value, determining the candidate maximum signal point as a target maximum signal point for leading the heart rate measurement.
Specifically, as shown in fig. 2, it is assumed that there are 5 maximum signal points to be observed (maximum points where black dots are located), and when absolute values of differences between the maximum amplitude value (i.e., the amplitude value of the candidate maximum signal point) in the maximum signal points to be observed and the amplitude values of the other remaining maximum signal points to be observed are all greater than a first preset value, it indicates that there is a maximum signal point in the maximum signal points to be observed whose amplitude value is much higher than the amplitude values of the other remaining maximum signal points to be observed, and at this time, it is considered that there is a target maximum signal point which is dominant heart rate measurement in the target amplitude spectrum, and the target maximum signal point is the maximum signal point in the maximum signal points to be observed. It is understood that when one amplitude peak exists in the target amplitude spectrum much higher than other amplitude peaks, the other amplitude peaks can be considered as the amplitude peaks where the error occurs approximately, and the target amplitude spectrum can be considered as having only one amplitude peak. The judgment method can accurately measure the heart rate.
In steps S421 to S423, a specific implementation manner of finding a target maximum signal point of the maximum signal points to be observed, which is dominant in heart rate measurement, is provided, and by calculating absolute values of differences between the candidate maximum signal points and amplitude values of other remaining maximum signal points to be observed, and determining the candidate maximum signal point as the target maximum signal point of the dominant heart rate measurement when the absolute values of the differences are greater than a first preset value, accurate measurement of the heart rate can be achieved.
Further, in step S40, according to the amplitude value of the maximum signal point to be observed, a target maximum signal point of the dominant heart rate measurement is searched for in the maximum signal point to be observed, which specifically includes:
s431: and screening a third signal point set from the maximum signal points to be observed according to the amplitude values of the maximum signal points to be observed, wherein the third signal point set comprises a plurality of third signal points, and the absolute values of the difference values of the third signal points and the amplitude values of other residual maximum signal points to be observed are all larger than a second preset value.
S432: and in the third signal point set, determining a fourth signal point set, wherein the fourth signal point set comprises a plurality of fourth signal points, and the frequency values of the fourth signal points are in a multiple relation, and the fourth signal point with the smallest frequency value is taken as a reference in the multiple relation of the frequency values.
Wherein, when the frequency values are in a multiple relation, the third signal point set can be determined as a fourth signal point set, i.e. the fourth signal point set is more than the third signal point set by a defined condition.
S433: if the amplitude value of the fourth signal point decreases with increasing frequency value or increases with increasing frequency value, the fourth signal point with the largest amplitude value is determined as the target maximum signal point that dominates the heart rate measurement.
It is to be understood that the relationship that the fourth signal point belongs to the frequency doubling peak can be determined when the amplitude value of the fourth signal point decreases with increasing frequency value (signal is weaker and weaker) or increases with increasing frequency value (signal is stronger and stronger).
The second predetermined value may be equal to the first predetermined value.
It can be understood that if the absolute values of the differences between the amplitude values of the plurality of maximum signal points to be observed and the amplitude values of the other remaining maximum signal points to be observed are greater than the second preset value, it is determined that a plurality of similar amplitude values exist in the target amplitude spectrum, and this condition cannot obtain an accurate heart rate measurement result from the perspective of heart rate measurement, but the analysis is performed from the perspective of frequency values, and when the frequency values of the plurality of maximum signal points to be observed are in a multiple relationship, a relationship of frequency doubling peaks exists between the plurality of similar amplitude peak values. The frequency doubling peak refers to a peak appearing at a position multiplied by a number on a frequency value axis, such as one peak appearing at 1.1Hz, one peak appearing at 2.2Hz, one peak appearing at 3.3Hz, \8230, and the frequency value is in a multiple relation based on a signal point with the minimum frequency value, such as two times and three times of 1.1Hz for 2.2Hz and 3.3Hz, respectively. In this embodiment, the frequency doubling peak represents that the signal satisfies a regular signal attenuation or enhancement of the frequency value, and the maximum value signal point to be observed with the largest amplitude value among the maximum value signal points to be observed can be used as the target maximum value signal point for dominating the heart rate measurement.
As shown in fig. 3, absolute values of differences between the amplitude values of 3 maximum signal points to be observed and the amplitude values of other remaining maximum signal points to be observed are all greater than a second preset value, which are corresponding maximum signal points to be observed of 1Hz, 2Hz, and 3Hz, respectively, frequency values of the 3 maximum signal points to be observed are in a multiple relationship and belong to a frequency doubling peak, and for this case, the maximum signal point to be observed with the largest amplitude value among the maximum signal points to be observed can be used as a target maximum signal point for dominant heart rate measurement.
In steps S431 to S433, a method for searching for a target maximum signal point for dominant heart rate measurement is provided, in which whether absolute values of differences between amplitude values of a plurality of maximum signal points to be observed and amplitude values of other remaining maximum signal points to be observed in the amplitude values of the maximum signal points to be observed are greater than a second preset value is determined, and when frequency values are in a multiple relation, whether a frequency doubling peak exists in a target amplitude spectrum is determined, and if a frequency doubling peak exists, a maximum signal point to be observed (i.e., a fourth signal point with a maximum amplitude value) of the maximum signal points to be observed can be used as a target maximum signal point for dominant heart rate measurement.
Further, in step S40, according to the amplitude value of the maximum signal point to be observed, a target maximum signal point that dominates the heart rate measurement is searched for in the maximum signal point to be observed, which specifically includes:
s441: and screening a fifth signal point set from the maximum signal points to be observed according to the amplitude values of the maximum signal points to be observed, wherein the fifth signal point set comprises a plurality of fifth signal points, and the absolute values of the difference values of the fifth signal points and the amplitude values of other residual maximum signal points to be observed are all larger than a third preset value.
S442: and determining a sixth signal point set in the fifth signal points, wherein the sixth signal point set comprises a plurality of sixth signal points, and the absolute value of the difference value of the frequency values of any two adjacent sixth signal points is smaller than a fourth preset value.
And when the absolute value of the difference value of the frequency values of any two adjacent signal points is smaller than a fourth preset value, determining that the fifth signal point set is a sixth signal point set, namely the sixth signal point set is larger than the fifth signal point set by a limiting condition.
S443: and if the number of the sixth signal points is within the preset numerical value interval, determining the sixth signal points as the target maximum signal points.
The third predetermined value may be equal to the first predetermined value.
It can be understood that if the absolute values of the differences between the amplitude values of the plurality of maximum signal points to be observed and the amplitude values of the other remaining maximum signal points to be observed are all greater than the third preset value, this situation may not result in an accurate heart rate measurement from the viewpoint of heart rate measurement, but if the absolute values of the differences between any two adjacent frequency values of the plurality of maximum signal points to be observed are all less than the third preset value, that is, the frequency values are very close to each other, then the plurality of maximum signal points to be observed may be approximately regarded as one maximum signal point, and accordingly, the plurality of amplitude peaks may be approximately regarded as one peak. (different from the frequency doubling peak, the frequency values of the maximum signal points to be observed of the frequency doubling peak generally have larger difference (larger than a fourth preset value), when a plurality of amplitude peak values are approximated to one peak value, because the frequency values of the corresponding maximum signal points to be observed are very close, the existing error can be ignored at the moment, the effective heart rate measurement is realized, and the repeated heart rate detection is reduced.
As shown in fig. 4, absolute values of differences between the amplitude values of the 3 maximum signal points to be observed and the amplitude values of the other remaining maximum signal points to be observed are all greater than a third preset value, and the frequency values of the 3 maximum signal points to be observed have no multiple relationship, but the absolute value of the difference between any two adjacent frequency values of the 3 maximum signal points to be observed is less than the third preset value, and at this time, the 3 maximum signal points to be observed can all be used as the target maximum signal point for dominant heart rate measurement.
Preferably, the sixth signal points included in the sixth signal point set may be two or three, and the obtained heart rate measurement result is more accurate.
In steps S441-S443, a method for searching for a target maximum signal point for dominant heart rate measurement is provided, in which by determining whether absolute values of differences between amplitude values of a plurality of maximum signal points to be observed and amplitude values of other remaining maximum signal points to be observed in a target amplitude spectrum are all greater than a third preset value, and absolute values of differences between frequency values of the plurality of maximum signal points to be observed are all less than a fourth preset value, a target maximum signal point for dominant heart rate measurement is found, thereby enabling effective heart rate measurement and reducing repeated heart rate detection.
S50: and obtaining a target heart rate value according to the target maximum value signal point.
Further, in step S50, obtaining a target heart rate value according to the target maximum value signal point includes: and calculating the mean value of the amplitude values of the sixth signal points in the sixth signal point set, and determining the mean value of the amplitude values as the target heart rate value.
It can be understood that, for the case that a plurality of sixth signal points in the sixth signal point set all exist as target maximum signal points for dominating heart rate measurement, an averaging processing mode may be adopted, and since frequency values of the sixth signal points are very close, an existing error may be ignored at this time, and accurate heart rate measurement is achieved.
It can be understood that, for the case that only one target maximum signal point dominates the target maximum signal point of the heart rate measurement, the method of obtaining the target heart rate value is to directly use the frequency value corresponding to the target maximum signal point as the target heart rate value.
Further, steps S421-S423, S431-S433 and S441-S443 are in parallel.
Further, if the target amplitude spectrum does not belong to the above steps S421-S423, S431-S433, and S441-S443, it can be considered that there is no maximum signal point to be observed on the target amplitude spectrum that dominates the heart rate measurement.
Further, if the maximum signal point to be observed, which is not dominant in heart rate measurement, does not exist on the target amplitude spectrum, reminding information is sent to remind the user to measure the heart rate again.
It can be understood that if there is no maximum signal point to be observed on the target amplitude spectrum, which is dominant in heart rate measurement, the obtained target amplitude spectrum is considered to be chaotic, and the heart rate measurement needs to be performed again for the situation, and the user is reminded by sending reminding information, so that the user performs the heart rate measurement after eliminating possible interference factors, and the accurate measurement of the heart rate is ensured.
In the embodiment of the invention, firstly, a target amplitude spectrum is obtained according to a heart rate measurement signal, the target amplitude spectrum comprises signal points, and coordinate values of the signal points comprise amplitude values and frequency values; then, a preset number of maximum signal points to be observed are obtained from the maximum signal points, wherein the amplitude value of the maximum signal point to be observed is larger than the amplitude values of other residual maximum signal points, so that the analysis of the target amplitude spectrum can be simplified to be carried out on the maximum preset number of maximum signal points, and the efficiency and implementation feasibility of heart rate measurement can be improved; and then searching a target maximum signal point for leading heart rate measurement in the maximum signal points to be observed according to the amplitude value of the maximum signal point to be observed, analyzing from the potential relation between the maximum signal points to be observed based on the amplitude value, searching the target maximum signal point for leading heart rate measurement, and obtaining a target heart rate value according to the target maximum signal point, thereby realizing accurate heart rate measurement in actual processing.
It should be understood that, the sequence numbers of the steps in the foregoing embodiments do not imply an execution sequence, and the execution sequence of each process should be determined by its function and inherent logic, and should not constitute any limitation to the implementation process of the embodiments of the present invention.
Based on the heart rate measurement method provided in the embodiment, the embodiment of the invention further provides an embodiment of a device for realizing the steps and the method in the embodiment of the method.
Fig. 5 shows a schematic block diagram of a heart rate measuring device in one-to-one correspondence with the heart rate measuring method in the embodiment. As shown in fig. 5, the heart rate measuring apparatus includes a target amplitude spectrum acquiring module 10, a maximum signal point acquiring module 20, a maximum signal point acquiring module 30 to be observed, a searching module 40, and a target heart rate value acquiring module 50. The implementation functions of the target amplitude spectrum acquisition module 10, the maximum signal point acquisition module 20, the maximum signal point to be observed acquisition module 30, the search module 40, and the target heart rate value acquisition module 50 correspond to the steps corresponding to the heart rate measurement method in the embodiment one to one, and in order to avoid redundancy, this embodiment is not described in detail.
And the target amplitude spectrum acquisition module 10 is configured to acquire a target amplitude spectrum according to the heart rate measurement signal, where the target amplitude spectrum includes signal points, and the coordinate values of the signal points include amplitude values and frequency values.
And a maximum signal point obtaining module 20, configured to obtain a maximum signal point in the signal points according to the amplitude value and the frequency value of the signal point.
The maximum signal point to be observed acquiring module 30 is configured to acquire a preset number of maximum signal points to be observed from the maximum signal points, where the amplitude value of the maximum signal point to be observed is greater than the amplitude values of the other remaining maximum signal points.
And the searching module 40 is configured to search, according to the amplitude value of the maximum signal point to be observed, a target maximum signal point which dominates the heart rate measurement from the maximum signal point to be observed.
And a target heart rate value obtaining module 50, configured to obtain a target heart rate value according to the target maximum signal point.
Alternatively, the maximum signal point acquisition module 20 includes a target signal point search unit, a signal point judgment unit, and a maximum signal point acquisition unit.
And the target signal point searching unit is used for searching a target signal point in a preset frequency value interval, and the amplitude value of the target signal point is larger than the amplitude values of other residual signal points.
And the signal point judging unit is used for judging whether a first signal point and a second signal point exist in the preset frequency value interval if the first signal point is found, wherein the frequency value of the first signal point is larger than the frequency value of the target signal point, and the frequency value of the second signal point is smaller than the frequency value of the target signal point.
And the maximum signal point acquisition unit is used for determining the target signal point as a maximum signal point and acquiring the maximum signal point if the first signal point and the second signal point exist in the preset frequency value interval.
Optionally, the finding module 40 includes a candidate maximum signal point acquisition unit, a difference absolute value calculation unit, and a first determination unit.
And the candidate maximum signal point acquisition unit is used for acquiring the maximum signal point to be observed with the largest amplitude value as the candidate maximum signal point according to the amplitude value of the maximum signal point to be observed.
And the difference absolute value calculating unit is used for calculating the difference between the candidate maximum signal point and the amplitude value of each of the other residual maximum signal points to be observed respectively and calculating the absolute value of each difference.
And the first determining unit is used for determining the candidate maximum signal point as a target maximum signal point for leading the heart rate measurement if the absolute value of each difference value is greater than a first preset value.
Optionally, the finding module 40 further includes a third signal point set obtaining unit, a fourth signal point set obtaining unit and a second determining unit.
And the third signal point set acquisition unit is used for screening out a third signal point set from the maximum signal points to be observed according to the amplitude values of the maximum signal points to be observed, the third signal point set comprises a plurality of third signal points, and the absolute values of the difference values of the third signal points and the amplitude values of other remaining maximum signal points to be observed are all larger than a second preset value.
And the fourth signal point set acquisition unit is used for determining a fourth signal point set in the third signal point set, wherein the fourth signal point set comprises a plurality of fourth signal points, and the frequency values of the fourth signal points are in a multiple relation, and the fourth signal point with the smallest frequency value is taken as a reference in the multiple relation.
And the second determining unit is used for determining the fourth signal point with the maximum amplitude value as the target maximum value signal point for dominating the heart rate measurement if the amplitude value of the fourth signal point is reduced along with the increase of the frequency value.
Optionally, the finding module 40 further includes a fifth signal point set obtaining unit, a sixth signal point set obtaining unit and a third determining unit.
And the fifth signal point set acquisition unit is used for screening a fifth signal point set from the maximum signal points to be observed according to the amplitude values of the maximum signal points to be observed, the fifth signal point set comprises a plurality of fifth signal points, and the absolute values of the difference values of the fifth signal points and the amplitude values of other remaining maximum signal points to be observed are all larger than a third preset value.
And the sixth signal point set acquisition unit is used for determining a sixth signal point set in the fifth signal points, wherein the sixth signal point set comprises a plurality of sixth signal points, and the absolute value of the frequency value difference value of any two adjacent sixth signal points is smaller than a fourth preset value.
And the third determining unit is used for determining the sixth signal point as the target maximum signal point if the number of the sixth signal points is within the preset numerical value interval.
Optionally, the target heart rate value obtaining module 50 is specifically configured to calculate an amplitude value mean value of the sixth signal point in the sixth signal point set, and determine the amplitude value mean value as the target heart rate value.
In the embodiment of the invention, firstly, a target amplitude spectrum is obtained according to a heart rate measurement signal, the target amplitude spectrum comprises signal points, coordinate values of the signal points comprise an amplitude value and a frequency value, then, a maximum signal point in the signal points is obtained according to the amplitude value and the frequency value of the signal points, and the maximum signal point on the target amplitude spectrum is analyzed to provide a technical premise for realizing the subsequent obtaining of a target heart rate value; then, a preset number of maximum signal points to be observed are obtained from the maximum signal points, wherein the amplitude value of the maximum signal point to be observed is larger than the amplitude values of other remaining maximum signal points, so that the analysis of the target amplitude spectrum can be simplified to be the analysis of the maximum preset number of maximum signal points, and the efficiency and the feasibility of implementation of heart rate measurement can be improved; and then searching a target maximum signal point for leading heart rate measurement in the maximum signal points to be observed according to the amplitude value of the maximum signal point to be observed, analyzing from the potential relation between the maximum signal points to be observed based on the amplitude value, searching the target maximum signal point for leading heart rate measurement, and obtaining a target heart rate value according to the target maximum signal point, thereby realizing accurate heart rate measurement in actual processing.
The present embodiment provides a computer-readable storage medium, where a computer program is stored, and when the computer program is executed by a processor, the heart rate measurement method in the embodiments is implemented, which is not described herein again to avoid repetition. Alternatively, the computer program is executed by the processor to implement the functions of the modules/units in the heart rate measuring device in the embodiments, which are not repeated herein to avoid repetition.
Fig. 6 is a schematic diagram of a computer device provided by an embodiment of the present invention. As shown in fig. 6, the computer device 60 of this embodiment includes: a processor 61, a memory 62, and a computer program 63 stored in the memory 62 and capable of running on the processor 61, wherein the computer program 63 when executed by the processor 61 implements the heart rate measurement method in the embodiment, which is not repeated herein. Alternatively, the computer program 63, when executed by the processor 61, implements the functions of each model/unit in the heart rate measuring device in the embodiments, which are not described herein again to avoid repetition.
The computing device 60 may be a desktop computer, a notebook, a palm top computer, a cloud server, or other computing device. The computer device 60 may include, but is not limited to, a processor 61, a memory 62. Those skilled in the art will appreciate that fig. 6 is merely an example of a computer device 60 and is not intended to limit the computer device 60 and that it may include more or less components than those shown, or some of the components may be combined, or different components, e.g., the computer device may also include input output devices, network access devices, buses, etc.
The Processor 61 may be a Central Processing Unit (CPU), other general purpose Processor, a Digital Signal Processor (DSP), an Application Specific Integrated Circuit (ASIC), a Field Programmable Gate Array (FPGA) or other Programmable logic device, discrete Gate or transistor logic, discrete hardware components, etc. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like.
The storage 62 may be an internal storage unit of the computer device 60, such as a hard disk or a memory of the computer device 60. The memory 62 may also be an external storage device of the computer device 60, such as a plug-in hard disk provided on the computer device 60, a Smart Media Card (SMC), a Secure Digital (SD) Card, a Flash memory Card (Flash Card), and the like. Further, memory 62 may also include both internal and external storage devices for computer device 60. The memory 62 is used to store computer programs and other programs and data required by the computer device. The memory 62 may also be used to temporarily store data that has been output or is to be output.
It will be apparent to those skilled in the art that, for convenience and brevity of description, only the above-mentioned division of the functional units and modules is illustrated, and in practical applications, the above-mentioned function distribution may be performed by different functional units and modules according to needs, that is, the internal structure of the apparatus is divided into different functional units or modules, so as to perform all or part of the functions described above.
The above examples are only intended to illustrate the technical solution of the present invention, but not to limit it; although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; such modifications and substitutions do not substantially depart from the spirit and scope of the embodiments of the present invention, and are intended to be included within the scope of the present invention.
Claims (4)
1. A method of heart rate measurement, the method comprising:
acquiring a target amplitude spectrum according to the heart rate measurement signal, wherein the target amplitude spectrum comprises signal points, and coordinate values of the signal points comprise amplitude values and frequency values;
acquiring a maximum signal point in the signal points according to the amplitude value and the frequency value of the signal points;
acquiring a preset number of maximum signal points to be observed from the maximum signal points, wherein the amplitude value of the maximum signal points to be observed is larger than the amplitude values of other residual maximum signal points;
searching a target maximum signal point for dominating heart rate measurement in the maximum signal point to be observed according to the amplitude value of the maximum signal point to be observed;
obtaining a target heart rate value according to the target maximum signal point;
wherein determining the target maximum signal point specifically comprises:
screening a third signal point set from the maximum signal points to be observed according to the amplitude values of the maximum signal points to be observed, wherein the third signal point set comprises a plurality of third signal points, and the absolute values of the difference values of the third signal points and the amplitude values of other remaining maximum signal points to be observed are all larger than a second preset value;
determining a fourth signal point set in the third signal point set, wherein the fourth signal point set comprises a plurality of fourth signal points, and the frequency values of the fourth signal points are in a multiple relation, and the fourth signal point with the smallest frequency value is taken as a reference in the multiple relation of the frequency values;
and if the amplitude value of the fourth signal point is increased along with the increase of the frequency value, determining the fourth signal point with the maximum amplitude value as a target maximum signal point for dominating the heart rate measurement.
2. A heart rate measuring device, the device comprising:
the target amplitude spectrum acquisition module is used for acquiring a target amplitude spectrum according to the heart rate measurement signal, wherein the target amplitude spectrum comprises signal points, and coordinate values of the signal points comprise amplitude values and frequency values;
the maximum signal point acquisition module is used for acquiring a maximum signal point in the signal points according to the amplitude value and the frequency value of the signal points;
the device comprises a to-be-observed maximum signal point acquisition module, a to-be-observed maximum signal point acquisition module and a to-be-observed maximum signal point acquisition module, wherein the to-be-observed maximum signal points are acquired from a preset number of the maximum signal points, and the amplitude value of each to-be-observed maximum signal point is larger than the amplitude values of other residual maximum signal points;
the searching module is used for searching a target maximum value signal point for leading heart rate measurement in the maximum value signal point to be observed according to the amplitude value of the maximum value signal point to be observed;
the target heart rate value acquisition module is used for acquiring a target heart rate value according to the target maximum signal point;
wherein the search module is specifically configured to:
screening a third signal point set from the maximum signal points to be observed according to the amplitude values of the maximum signal points to be observed, wherein the third signal point set comprises a plurality of third signal points, and the absolute values of the difference values of the third signal points and the amplitude values of other remaining maximum signal points to be observed are all larger than a second preset value;
determining a fourth signal point set in the third signal point set, wherein the fourth signal point set comprises a plurality of fourth signal points, and the frequency values of the fourth signal points are in a multiple relation, and the fourth signal point with the smallest frequency value is taken as a reference in the multiple relation;
and if the amplitude value of the fourth signal point is increased along with the increase of the frequency value, determining the fourth signal point with the maximum amplitude value as a target maximum signal point for dominating the heart rate measurement.
3. A computer device comprising a memory, a processor and a computer program stored in the memory and executable on the processor, characterized in that the processor, when executing the computer program, carries out the steps of the heart rate measurement method according to claim 1.
4. A computer-readable storage medium, in which a computer program is stored which, when being executed by a processor, carries out the steps of the heart rate measurement method according to claim 1.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910268887.5A CN110123297B (en) | 2019-04-04 | 2019-04-04 | Heart rate measuring method and device, computer equipment and storage medium |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910268887.5A CN110123297B (en) | 2019-04-04 | 2019-04-04 | Heart rate measuring method and device, computer equipment and storage medium |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110123297A CN110123297A (en) | 2019-08-16 |
CN110123297B true CN110123297B (en) | 2022-11-25 |
Family
ID=67569345
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910268887.5A Active CN110123297B (en) | 2019-04-04 | 2019-04-04 | Heart rate measuring method and device, computer equipment and storage medium |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110123297B (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110755055A (en) * | 2019-10-18 | 2020-02-07 | 上海掌门科技有限公司 | Method and equipment for determining waveform evaluation information of pulse waveform |
CN110935159B (en) * | 2019-12-25 | 2021-08-03 | 青岛英派斯健康科技股份有限公司 | Fitness equipment heart rate display method, device, equipment and storage medium |
CN113251945B (en) * | 2021-05-17 | 2022-07-12 | 东北大学秦皇岛分校 | Demodulation method of line profile imaging device and imaging device |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103230267B (en) * | 2013-05-14 | 2015-06-03 | 北京理工大学 | Anti-movement-interference extraction method for pulse rates |
CN104873184B (en) * | 2015-05-08 | 2017-09-29 | 上海与德通讯技术有限公司 | A kind of method and system for measuring heart rate |
US10123746B2 (en) * | 2015-05-08 | 2018-11-13 | Texas Instruments Incorporated | Accuracy of heart rate estimation from photoplethysmographic (PPG) signals |
CN106264505A (en) * | 2016-07-21 | 2017-01-04 | 浙江师范大学 | A kind of heart rate spectral peak system of selection based on support vector machine |
EP3305180A1 (en) * | 2016-10-05 | 2018-04-11 | Murata Manufacturing Co., Ltd. | Method and apparatus for monitoring heartbeats |
CN109222948B (en) * | 2018-09-17 | 2021-07-13 | 歌尔科技有限公司 | Method, device, electronic equipment and storage medium for eliminating motion interference noise |
-
2019
- 2019-04-04 CN CN201910268887.5A patent/CN110123297B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN110123297A (en) | 2019-08-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110123297B (en) | Heart rate measuring method and device, computer equipment and storage medium | |
CN110363076B (en) | Personnel information association method and device and terminal equipment | |
US10832032B2 (en) | Facial recognition method, facial recognition system, and non-transitory recording medium | |
CN108389195B (en) | Image detection method and device | |
CN105354721B (en) | Method and device for identifying machine operation behavior | |
CN109165657A (en) | A kind of image feature detection method and device based on improvement SIFT | |
WO2022100028A1 (en) | Interface traffic anomaly detection method and apparatus, terminal device, and storage medium | |
CN111626360A (en) | Method, device, equipment and storage medium for detecting fault type of boiler | |
CN110348717B (en) | Base station value scoring method and device based on grid granularity | |
US20180288121A1 (en) | Streaming media play mode determination method and apparatus | |
US20150293598A1 (en) | Method for processing information and electronic device | |
CN114297665A (en) | Intelligent contract vulnerability detection method and device based on deep learning | |
CN112985867A (en) | Steering engine test method, device, equipment and storage medium | |
CN112463564A (en) | Method and device for determining correlation index influencing host state | |
CN115756173A (en) | Eye tracking method, system, storage medium and computing equipment | |
CN108629219B (en) | Method and device for identifying one-dimensional code | |
US20170185831A1 (en) | Method and device for distinguishing finger and wrist | |
CN114202804A (en) | Behavior action recognition method and device, processing equipment and storage medium | |
CN107229621B (en) | Method and device for cleaning difference data | |
CN109286798B (en) | Frame position identification method and system of projection screen and terminal equipment | |
EP2998760A1 (en) | Method, user equipment, system and computer readable medium for localizing an user equipment | |
CN105786614A (en) | Loading calculation method and loading calculation system for processor in electronic device | |
CN112826460B (en) | Physiological signal frequency extraction method, device, physiological signal acquisition equipment and medium | |
KR20200117582A (en) | Method and apparatus for mental workload assessment of multitasking user using physiological data and simulated data | |
EP4239542A1 (en) | Assistance system, assistance method, 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 |