CN109283121A - Pulse recognition method and apparatus, analysis instrument, storage medium - Google Patents
Pulse recognition method and apparatus, analysis instrument, storage medium Download PDFInfo
- Publication number
- CN109283121A CN109283121A CN201811183186.3A CN201811183186A CN109283121A CN 109283121 A CN109283121 A CN 109283121A CN 201811183186 A CN201811183186 A CN 201811183186A CN 109283121 A CN109283121 A CN 109283121A
- Authority
- CN
- China
- Prior art keywords
- pulse
- data
- amplitude threshold
- end point
- peak
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000012545 processing Methods 0.000 claims description 10
- 238000010586 diagram Methods 0.000 description 16
- 238000005070 sampling Methods 0.000 description 8
- 230000009977 dual effect Effects 0.000 description 6
- TVEXGJYMHHTVKP-UHFFFAOYSA-N 6-oxabicyclo[3.2.1]oct-3-en-7-one Chemical compound C1C2C(=O)OC1C=CC2 TVEXGJYMHHTVKP-UHFFFAOYSA-N 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000007792 addition Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000000977 initiatory effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000003252 repetitive effect Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/10—Investigating individual particles
- G01N15/14—Optical investigation techniques, e.g. flow cytometry
- G01N15/1456—Optical investigation techniques, e.g. flow cytometry without spatial resolution of the texture or inner structure of the particle, e.g. processing of pulse signals
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/10—Investigating individual particles
- G01N2015/1006—Investigating individual particles for cytology
Landscapes
- Chemical & Material Sciences (AREA)
- Dispersion Chemistry (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Radar Systems Or Details Thereof (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
The present invention discloses a kind of pulse recognition method and apparatus, analysis instrument, storage medium.This method comprises: pulse data to be identified is respectively compared with the first amplitude threshold and the second amplitude threshold one by one, the first amplitude threshold is greater than the second amplitude threshold；First data if it exists, the first data are more than or equal to the first amplitude threshold, then start identification first effective peak corresponding with the first amplitude threshold, and first effective peak meets the first default starting point and end-condition；Second data if it exists, second data are more than or equal to the second amplitude threshold, then start identification second effective peak corresponding with the second amplitude threshold, second effective peak meets the second default starting point and end-condition, and the peak value at second effective peak is more than or equal to default small-pulse effect height threshold less than the difference between the peak value and the second amplitude threshold at the first amplitude threshold and second effective peak.Using the technical solution in the embodiment of the present invention, the accuracy rate of pulse recognition can be improved.
Description
Technical Field
The invention relates to the technical field of medical equipment, in particular to a pulse recognition method and device, an analysis instrument and a storage medium.
Background
When medical equipment is adopted for cell classification, laser is utilized to irradiate cell samples passing through the flow chamber one by one, the optical detector converts detected characteristic optical signals into electric pulse signals and outputs the electric pulse signals to the data processing unit, the data processing unit identifies the electric pulse signals and positions each cell on a multi-dimensional coordinate system based on identified effective peaks to form a scatter diagram of a plurality of cell samples, and due to the difference of various cells on the signals, the cells of the same type are gathered together, gaps are reserved among different types, and therefore the classification purpose is achieved.
The pulse identification method in the prior art is mainly a single threshold method, namely, a threshold line is set, a pulse starting point and a pulse ending point are judged by adding a threshold to a slope, and a maximum value is obtained between the starting point and the ending point to obtain a peak value of the pulse.
However, the inventor of the present application has found that the accuracy of pulse recognition in the prior art is low, and missing recognition and erroneous recognition are likely to occur.
Disclosure of Invention
The embodiment of the invention provides a pulse identification method and device, an analysis instrument and a storage medium, which can improve the accuracy of pulse identification.
In a first aspect, an embodiment of the present invention provides a pulse identification method, where the pulse identification method includes:
comparing the pulse data to be identified with a first amplitude threshold value and a second amplitude threshold value one by one, wherein the first amplitude threshold value is larger than the second amplitude threshold value and is larger than the maximum value of a preset small pulse amplitude interval;
if the first data exist and the first data are larger than or equal to a first amplitude threshold value, a first effective peak corresponding to the first amplitude threshold value is identified, and the first effective peak meets first preset starting point and end point conditions;
and if the second data exists, the second data is greater than or equal to a second amplitude threshold value, a second effective peak corresponding to the second amplitude threshold value is identified, the second effective peak meets second preset starting point and end point conditions, the peak value of the second effective peak is smaller than the first amplitude threshold value, and the difference value between the peak value of the second effective peak and the second amplitude threshold value is greater than or equal to a preset small pulse height threshold value.
In one possible implementation manner of the first aspect, the first preset starting point and end point condition includes: a first pulse starting point and a first pulse ending point exist in the first effective peak; the first N data of the first pulse starting point are sequentially increased and are smaller than the data of the first pulse starting point, and the last M data of the first pulse starting point are sequentially increased and are larger than the data of the first pulse starting point; the first N data of the first pulse end point are sequentially reduced and are larger than the data of the first pulse end point, the last M data of the first pulse end point are sequentially reduced and are smaller than the data of the first pulse end point, and both N and M are natural numbers which are larger than or equal to 1; the first pulse start point and the first pulse end point are equal to the first amplitude threshold.
In one possible implementation of the first aspect, a peak value of the first effective peak is equal to a maximum value of a plurality of data located between a first pulse start point and a first pulse end point in the pulse data to be identified.
In one possible implementation of the first aspect, the method further comprises: and if the second data exists and the second data is greater than or equal to the second amplitude threshold, starting to identify a first effective peak corresponding to the second amplitude threshold, wherein the first effective peak meets the conditions of a second preset starting point and an end point, and the peak value is greater than the first amplitude threshold.
In one possible implementation manner of the first aspect, the second preset starting point and end point condition includes: a second pulse starting point and a second pulse ending point exist in the second effective peak; the first N data of the starting point of the second pulse are sequentially increased and are smaller than the data of the starting point of the second pulse, and the last M data of the starting point of the second pulse are sequentially increased and are larger than the data of the starting point of the second pulse; the first N data of the second pulse end point are sequentially increased and are smaller than the data of the second pulse end point, the last M data of the second pulse end point are sequentially increased and are larger than the data of the first pulse end point, and both N and M are natural numbers which are larger than or equal to 1; the second pulse start point and the second pulse end point are equal to the second amplitude threshold.
In one possible implementation of the first aspect, the peak value of the second effective peak is equal to a maximum value of a plurality of data located between the start point of the second pulse and the end point of the second pulse in the pulse data to be identified.
In a possible implementation manner of the first aspect, a width of the second effective peak is greater than or equal to a preset small pulse width threshold, and the width of the second effective peak is a distance between a second pulse starting point and a second pulse ending point.
In a possible embodiment of the first aspect, the difference between the second amplitude threshold and the baseline value of the pulse to be identified is smaller than a preset difference.
In a second aspect, an embodiment of the present invention provides a pulse recognition apparatus, including:
the data processing module is used for comparing the pulse data to be identified with a first amplitude threshold value and a second amplitude threshold value one by one, wherein the first amplitude threshold value is larger than the second amplitude threshold value, and the first amplitude threshold value is larger than the maximum value of a preset small pulse amplitude interval;
the first identification module is used for starting to identify a first effective peak corresponding to a first amplitude threshold value if first data exists and the first data is greater than or equal to the first amplitude threshold value, and the first effective peak meets first preset starting point and end point conditions;
and the second identification module is used for starting to identify a second effective peak corresponding to the second amplitude threshold if the second data exists and the second data is greater than or equal to the second amplitude threshold, wherein the second effective peak meets second preset starting and ending point conditions, the peak value of the second effective peak is smaller than the first amplitude threshold, and the difference value between the peak value of the second effective peak and the second amplitude threshold is greater than or equal to a preset small pulse height threshold.
In a possible implementation manner of the second aspect, the first identification module is further configured to: and if the second data exists and the second data is greater than or equal to the second amplitude threshold, starting to identify a first effective peak corresponding to the second amplitude threshold, wherein the first effective peak meets the conditions of a second preset starting point and an end point, and the peak value is greater than the first amplitude threshold.
In a third aspect, embodiments of the present invention provide an analysis apparatus comprising a pulse recognition apparatus according to any one of claim 8.
In a third aspect, an embodiment of the present invention provides a computer-readable storage medium, on which a program is stored, and the program, when executed by a processor, implements the pulse recognition apparatus as described above.
As described above, in one aspect, the embodiment of the present invention sets two amplitude thresholds, where the first amplitude threshold is greater than the second amplitude threshold and greater than the maximum value of the preset small pulse amplitude interval, so as to completely distinguish the large pulse from the small pulse; on the other hand, the embodiment of the invention combines the actual situation and respectively adopts different starting point and end point conditions for the large pulse and the small pulse, thereby avoiding the defect of single threshold judgment, improving the accuracy rate of pulse identification and avoiding the situations of pulse missing identification and error identification.
Drawings
The present invention will be better understood from the following description of specific embodiments thereof taken in conjunction with the accompanying drawings, in which like or similar reference characters designate like or similar features.
Fig. 1 is a schematic flow chart of a pulse identification method according to an embodiment of the present invention;
FIG. 2 is a diagram illustrating single threshold-based identification of forward-scattered fsc pulses in a cell according to an embodiment of the present invention;
FIG. 3 is a diagram of dual threshold based fsc pulse identification corresponding to FIG. 2;
FIG. 4 is a diagram illustrating single threshold based fsc pulse identification according to another embodiment of the present invention;
FIG. 5 is a diagram of dual threshold based fsc pulse identification corresponding to FIG. 4;
FIG. 6 is a diagram illustrating single threshold based fsc pulse identification according to yet another embodiment of the present invention;
FIG. 7 is a diagram of dual threshold based fsc pulse identification corresponding to FIG. 6;
fig. 8 is a schematic structural diagram of a pulse recognition device according to an embodiment of the present invention.
Detailed Description
Features and exemplary embodiments of various aspects of the present invention will be described in detail below. In the following detailed description, numerous specific details are set forth in order to provide a thorough understanding of the present invention.
The embodiment of the invention provides a pulse identification method and device, an analysis instrument and a storage medium, which are applied to medical equipment needing pulse identification.
Wherein the large pulse has a higher amplitude and the small pulse has a lower amplitude. In one example, the peak interval of the small pulse is 0-0.5, while the peak interval of the large pulse is greater than 0.5. In practical applications, the definitions of the amplitudes of the large pulse and the small pulse may vary according to the types of the objects to be detected and the types of the detection signals.
Fig. 1 is a schematic flow chart of a pulse identification method according to an embodiment of the present invention. As shown in fig. 1, the pulse recognition method includes steps 101 to 103.
In step 101, the pulse data to be identified are compared one by one with a first amplitude threshold value and a second amplitude threshold value, respectively.
The first amplitude threshold is larger than the second amplitude threshold, the first amplitude threshold is used for performing identification operation on the large pulse, the second amplitude threshold is used for performing identification operation on the small pulse, and as the amplitude definitions of the large pulse and the small pulse can change along with the difference between the type of the object to be detected and the type of the detection signal, an engineer can reasonably set the first amplitude threshold and the second amplitude threshold according to the currently input pulse waveform so as to perform distinguishing identification on the large pulse and the small pulse as much as possible.
In practical applications, a high threshold is used as a judgment value of the large pulse, and after the high threshold is set, a pulse having a peak value smaller than the high threshold is considered to be a small pulse. In one example, if the preset small pulse amplitude interval is 0-0.5, the high threshold should be greater than 0.5.
When the low threshold is used, the low threshold is required to be as close to the reference line 0 of the pulse to be recognized as possible in principle, so that the small pulse is between the high threshold and the low threshold as possible, and the recognition range of the small pulse is improved. In particular, it is possible that the difference between the low threshold and the pulse limit value is smaller than a preset difference, where the preset difference is a minimum value.
In step 102, if there is first data, where the first data is greater than or equal to a first amplitude threshold, a first effective peak corresponding to the first amplitude threshold starts to be identified, and the first effective peak satisfies a first preset start point and end point condition.
In step 103, if there is second data, where the second data is greater than or equal to a second amplitude threshold, a second effective peak corresponding to the second amplitude threshold starts to be identified, where the second effective peak meets second preset start and end point conditions, a peak value of the second effective peak is smaller than the first amplitude threshold, and a difference between the peak value of the second effective peak and the second amplitude threshold is greater than or equal to a preset small pulse height threshold.
The first preset starting point and end point condition is a pulse starting point and pulse end point judgment condition for large pulse identification, if the first pulse starting point and the first pulse end point exist in the large pulse, the large pulse is indicated to meet the first preset starting point and end point condition, the large pulse can be identified as an effective peak, and the peak value of the effective peak can be the maximum value in a plurality of data between the identified starting point and end point of the large pulse.
The second preset start point and end point condition is a pulse start point and pulse end point determination condition for identifying the small pulse, and if the small pulse has the second pulse start point and the second pulse end point, the small pulse is indicated to meet the second preset start point and end point condition, the large pulse can be identified as an effective peak, and the peak value of the effective peak can be the maximum value in the data between the start point and the end point of the identified small pulse.
In an optional embodiment, if the second data is present, the second data is greater than or equal to the second amplitude threshold, and the identification of the first valid peak corresponding to the second amplitude threshold is also started, where the first valid peak satisfies the second preset starting point and end point conditions, and the peak value is greater than the first amplitude threshold.
Because the pulse with the peak value higher than the first amplitude threshold value can simultaneously meet the first preset starting point and end point condition and the second preset starting point and end point condition, the pulse which meets the second preset starting point and end point condition and has the peak value larger than the first amplitude threshold value is identified as the first effective peak.
To understand the above steps 102 and 103, for example, suppose that two threshold lines for a large pulse and a small pulse are th _ h and th _ l respectively, data is traversed from the head or tail of the pulse data to be identified and sequentially stored in the array d, i is an integer greater than or equal to 0, d [ i ] is traversed from i ═ 0, and then d [ i ] is compared with th _ h and th _ l respectively, then:
when d [ i ] > (th _ h) is met, executing a first preset starting point and end point condition, and starting large pulse identification; when d [ i ] >, th _ l is satisfied, a second preset start point and end point condition is executed, and small pulse recognition is started.
As described above, in one aspect, the embodiment of the present invention sets two amplitude thresholds, where the first amplitude threshold is greater than the second amplitude threshold and greater than the maximum value of the preset small pulse amplitude interval, so as to completely distinguish the large pulse from the small pulse; on the other hand, the embodiment of the invention combines the actual situation and respectively adopts different starting point and end point conditions for the large pulse and the small pulse, thereby avoiding the defect of single threshold judgment, improving the accuracy rate of pulse identification and avoiding the situations of pulse missing identification and error identification.
Wherein, the conditions of the starting point of the first pulse are as follows: the first N data of the first pulse starting point are sequentially increased and are smaller than the data of the first pulse starting point, the last M data of the first pulse starting point are sequentially increased and are larger than the data of the first pulse starting point, and the first pulse starting point and the first pulse end point are equal to a first amplitude threshold value.
For example, with d [ i ] as the pulse start, d [ i ] > d [ i-1] > d [ i-2] > …. > d [ i-N ], d [ i ] < d [ i +1] > d [ i +2] > …. > d [ i + M ].
The conditions for the end of the first pulse are: the first N data of the first pulse end point are sequentially reduced and are larger than the data of the first pulse end point, and the last M data of the first pulse end point are sequentially reduced and are smaller than the data of the first pulse end point.
For example, with d [ i + x ] as the end of the pulse, d [ i ] < d [ i + x-1] > d [ i + x-2] > …. > d [ i + x-N ], [ i ] > d [ i + x +1] > d [ i + x +2] > …. > d [ i + x + M ].
The conditions for the start of the second pulse are: the first N data of the starting point of the second pulse are sequentially increased and are smaller than the data of the starting point of the second pulse, the last M data of the starting point of the second pulse are sequentially increased and are larger than the data of the starting point of the second pulse, and the starting point of the first pulse and the ending point of the first pulse are equal to the first amplitude threshold value.
For example, with d [ i ] as the pulse start, d [ i ] > d [ i-1] > d [ i-2] > …. > d [ i-N ], d [ i ] < d [ i +1] > d [ i +2] > …. > d [ i + M ].
The conditions for the end of the second pulse are: the first N data of the second pulse end point are sequentially increased and are smaller than the data of the second pulse end point, the last M data of the second pulse end point are sequentially increased and are larger than the data of the first pulse end point, and both N and M are natural numbers which are larger than or equal to 1.
For example, with d [ i + x ] as the end of the pulse, d [ i ] < d [ i + x-1] > d [ i + x-2] > …. > d [ i + x-N ], [ i ] > d [ i + x +1] > d [ i + x +2] > …. > d [ i + x + M ].
As can be seen from the above, the method for determining the first pulse start point and the first pulse end point in the first preset start point and end point condition is similar to the method for determining the second pulse start point and the second pulse end point in the first preset start point and end point condition, except that the two reference thresholds are different, the former is based on the first amplitude threshold th _ h, and the latter is based on the second amplitude threshold th _ l.
To facilitate understanding of the technical solutions in the embodiments of the present invention, the following examples are given.
FIG. 2 is a diagram illustrating single threshold-based identification of forward-scattering fsc pulses in a cell according to an embodiment of the present invention. Where the abscissa is the number of sampling points and the ordinate is the signal strength, the forward angle scattered pulses are also referred to as fsc (forwardscatter) pulses.
In the example of fig. 2, the threshold line is th2, and the fsc pulse signal is at the large peak f under the influence of the light signal and the circuit signal21May produce a small peak f22Based on a single threshold method, due to small spikes f22So that the sampling points b and c do not coincide with the large peak f21So that the large peak f is determined21An error is identified. The method specifically comprises the following steps:
b >0, the previous sampling point a of b is less than 0, but the next sampling point c of b is less than 0;
c is greater than 0, the previous sampling point b of c is greater than 0, but the next sampling point d of c is less than 0; .
FIG. 3 is a diagram illustrating dual threshold based fsc pulse identification corresponding to FIG. 2. FIG. 3 is different from FIG. 2 in that FIG. 3 employs a threshold line th-h-2 and a threshold line th-l-2, wherein the threshold value corresponding to th-h-2 is about 0.5, the threshold value corresponding to th-l-2 is smaller and close to 0, and the preset small pulse height threshold value should be smaller than 0.5.
While traversing the pulse data in FIG. 3 one by one on the basis of th-h-2, e2The data of the points is equal to th-h-2, starting a large peak f21And e is recognized, and2point sum f2All the points meet the first preset starting point and end point conditions, and the large wave crest f can be obtained21Identified as a valid peak and determined e2Point sum f2Point is large wave peak f21Starting and ending points of (e)2Point sum f2Maximum value H in a plurality of data between points21As a large peak f21Peak value of (a).
While traversing the pulse data in FIG. 3 one by one on th-l-2 basis, g2The data of the point is equal to th-l-2, so that a small spike f is started22Identification of (i) though g2Dot sum h2All points satisfy the second preset starting point and end point conditions, but are based on g2Dot sum h2The peak obtained by the point is also H21，H21Greater than a preset small pulse height threshold of 0.5, and therefore, a small spike f22Can not be identified as effective peak by mistake and can not be used for large wave peak f21The identification of (a) causes interference.
FIG. 4 is a diagram illustrating single threshold based fsc pulse identification according to another embodiment of the present invention. Wherein, the abscissa is the number of sampling points, and the ordinate is the signal intensity.
In the example of FIG. 4, where the threshold line is th4, two cells are in close proximity and passed through the flow cell, the fsc pulse signal will generate an M-wave, including the peak f41And f42，f41And f42All the peak values of (1) are more than 0.5, and all the peak values can be regarded as large peak values, because the peak values f41The value at the bottom is higher than the threshold th4, resulting in a peak f41Is unrecognized and thus is underidentified.
FIG. 5 is a diagram illustrating dual threshold based fsc pulse identification corresponding to FIG. 4. FIG. 5 is different from FIG. 4 in that FIG. 5 employs a threshold line th-h-4 and a threshold line th-l-4, wherein the threshold value corresponding to th-h-4 is about 0.5, the threshold value corresponding to th-l-4 is smaller and close to 0, and the predetermined small pulse height threshold value should be smaller than 0.5.
While traversing the pulse data in FIG. 5 one by one on the basis of th-h-4, e4The data of the point is equal to th-h-4, starting a large peak f41And e is recognized, and4point sum f4All the points meet the first preset starting point and end point conditions, and the large wave crest f can be obtained41Identified as a valid peak and determined e4Point sum f4Point is large wave peak f41Starting and ending points of (e)4Point sum f4Maximum value H in a plurality of data between points41As a large peak f41Peak value of (a).
In the same way, g4The data of the point is equal to th-h-4, starting a large peak f42And g is4Dot sum h4All the points meet the first preset starting point and end point conditions, and the large wave crest f can be obtained42Is identified asEffective peak and determine g4Dot sum h4Point is large wave peak f42Starting and ending points of (c), g4Dot sum h4Maximum value H in a plurality of data between points42As a large peak f42Thereby avoiding the peak f41Is detected.
While traversing the pulse data in FIG. 5 one by one on th-l-4 basis, p4The data of the points is equal to th-l-2, thus initiating the identification of small peaks (unknown), although p4Dot sum s4All points satisfy the second preset starting point and end point conditions, but are based on p4Dot sum s4The peak value obtained by point is H42And if the pulse height is larger than the preset small pulse height threshold value by 0.5, no small wave crest exists.
FIG. 6 is a diagram illustrating single threshold based fsc pulse identification according to a third embodiment of the present invention. Wherein, the abscissa is the number of sampling points, and the ordinate is the signal intensity.
In the example of FIG. 6, the threshold line is th6, the large peak f in the fsc pulse signal61After that, an overshoot linearity appears on the circuit processing, so that the signal value is reduced to be lower than the base line (0.0), and the wavelet peak f which is followed immediately at the moment32It is missed because the amplitude is too small to be lower than the threshold th 6.
FIG. 7 is a diagram illustrating dual threshold based fsc pulse identification corresponding to FIG. 6. FIG. 7 is different from FIG. 6 in that FIG. 7 employs a threshold line th-h-6 and a threshold line th-l-6, wherein the threshold value corresponding to th-h-6 is about 0.5, the threshold value corresponding to th-l-6 is smaller and close to 0, and the predetermined small pulse height threshold value should be smaller than 0.5.
While traversing the pulse data in FIG. 7 one by one on the basis of th-h-6, e6The data of the point is equal to th-h-6, starting a large peak f61And e is recognized, and6point sum f6All the points meet the first preset starting point and end point conditions, and the large wave crest f can be obtained61Identified as a valid peak and determined e6Point sum f6Point is large wave peak f61Starting and ending points of (e)6Point sum f6Dot of dotsMaximum value H among a plurality of data in between61As a large peak f61Peak value of (a).
While traversing the pulse data in FIG. 7 one by one on th-l-6 basis, g6The data of the point is equal to th-l-6, thus starting a small peak f62Identification of (g)6Dot sum h6All the points meet the second preset starting point and end point conditions and are based on g6Dot sum h6The peak obtained by the point is also H62，H62Less than 0.5 of the preset small pulse height threshold, therefore, the small wave crest f can be reduced62Identified as a valid peak and determined g6Dot sum h6Point is wavelet peak f62And thus will not be missed.
Furthermore, th-l-6 traverses the pulse data in FIG. 7 one by one, p6The data of the point is equal to th-l-6, so that a large peak f starts at the same time61Identification of p6Dot sum s6All the points satisfy the second preset starting point and end point conditions and are based on p6Dot sum s6Peak value H obtained by point61Is greater than th-h-6, therefore, the large peak f61Identified as a valid peak.
To avoid repetitive operations during large peak identification, p may be detected6When the data of the point is equal to th-l-6, only the identification of the small pulse is executed, and the identification of the large pulse is not executed until e6When the data of the point is equal to th-h-6, the identification of the large wave crest is switched in again.
As described above, since the shape of the large pulse at the tip is stable, the pulse is usually identified by a first large amplitude threshold and a first preset start and end point condition. In the identification process of the small pulse, in addition to the second preset starting point and end point conditions and the peak value judgment condition mentioned above, the judgment condition of the pulse width may be increased, specifically, the width of the second effective peak is greater than or equal to the preset small pulse width threshold, and the width of the second effective peak is equal to the distance between the second pulse starting point and the second pulse end point. Thus, by further limiting the pulse width and the relative height of the peak value, the interference of circuit noise and a large amount of baseline noise in the process of identifying small pulses can be solved.
Fig. 8 is a schematic structural diagram of a pulse recognition device according to an embodiment of the present invention. As shown in fig. 8, the pulse recognition apparatus includes: a data processing module 801, a first identification module 802 and a second identification module 803.
The data processing module 801 is configured to compare pulse data to be identified with a first amplitude threshold and a second amplitude threshold one by one, where the first amplitude threshold is greater than the second amplitude threshold, and the first amplitude threshold is greater than a maximum value of a preset small pulse amplitude interval.
The first identifying module 802 is configured to, if there is first data, where the first data is greater than or equal to a first amplitude threshold, start identifying a first effective peak corresponding to the first amplitude threshold, where the first effective peak meets a first preset starting point and ending point condition.
The second identifying module 803 is configured to, if there is second data, where the second data is greater than or equal to a second amplitude threshold and smaller than the first amplitude threshold, start identifying a second effective peak corresponding to the second amplitude threshold, where the second effective peak meets second preset starting and ending point conditions, a peak value of the second effective peak is smaller than the first amplitude threshold, and a difference between the peak value of the second effective peak and the second amplitude threshold is greater than or equal to a preset small pulse height threshold.
In an optional embodiment, the first identifying module 802 is further configured to, if there is second data, where the second data is greater than or equal to a second amplitude threshold, further start identifying a first valid peak corresponding to the second amplitude threshold, where the first valid peak meets a second preset start point and end point condition, and a peak value is greater than the first amplitude threshold.
Because the pulse with the peak value higher than the first amplitude threshold value can simultaneously meet the first preset starting point and end point condition and the second preset starting point and end point condition, in actual operation, the first identification module can identify the pulse with the peak value higher than the first amplitude threshold value as the first effective peak when the pulse meets the second preset starting point and end point condition.
It should be noted that, in the actual execution process, the collected pulse data may be processed in series or in parallel, and during parallel processing, parallel processing modules (modules) may be respectively set for the large pulse and the small pulse, and the collected pulse data is directly sent to the two processing modules (modules), and simultaneously, pulse recognition is performed in real time.
The embodiment of the invention also provides an analysis instrument, which comprises the pulse recognition device.
An embodiment of the present invention further provides a computer-readable storage medium, on which a program is stored, and the program, when executed by a processor, implements the pulse recognition apparatus as described above.
It should be clear that the embodiments in this specification are described in a progressive manner, and the same or similar parts in the embodiments are referred to each other, and each embodiment focuses on the differences from the other embodiments. For the device embodiments, reference may be made to the description of the method embodiments in the relevant part. Embodiments of the invention are not limited to the specific steps and structures described above and shown in the drawings. Those skilled in the art may make various changes, modifications and additions to, or change the order between the steps, after appreciating the spirit of the embodiments of the invention. Also, a detailed description of known process techniques is omitted herein for the sake of brevity.
The functional blocks shown in the above-described structural block diagrams may be implemented as hardware, software, firmware, or a combination thereof. When implemented in hardware, it may be, for example, an electronic circuit, an Application Specific Integrated Circuit (ASIC), suitable firmware, plug-in, function card, or the like. When implemented in software, the elements of an embodiment of the invention are the programs or code segments used to perform the required tasks. The program or code segments may be stored in a machine-readable medium or transmitted by a data signal carried in a carrier wave over a transmission medium or a communication link. A "machine-readable medium" may include any medium that can store or transfer information. Examples of a machine-readable medium include electronic circuits, semiconductor memory devices, ROM, flash memory, Erasable ROM (EROM), floppy disks, CD-ROMs, optical disks, hard disks, fiber optic media, Radio Frequency (RF) links, and so forth. The code segments may be downloaded via computer networks such as the internet, intranet, etc.
Embodiments of the present invention may be embodied in other specific forms without departing from its spirit or essential characteristics. For example, the algorithms described in the specific embodiments may be modified without departing from the basic spirit of the embodiments of the present invention. The present embodiments are therefore to be considered in all respects as illustrative and not restrictive, the scope of the embodiments of the invention being indicated by the appended claims rather than by the foregoing description, and all changes which come within the meaning and range of equivalency of the claims are therefore intended to be embraced therein.
Claims (12)
1. A method of pulse recognition, comprising:
comparing the pulse data to be identified with a first amplitude threshold value and a second amplitude threshold value one by one, wherein the first amplitude threshold value is larger than the second amplitude threshold value, and the first amplitude threshold value is larger than the maximum value of a preset small pulse amplitude interval;
if first data exist and the first data are larger than or equal to the first amplitude threshold value, starting to identify a first effective peak corresponding to the first amplitude threshold value, wherein the first effective peak meets first preset starting point and end point conditions;
if second data exists, the second data is larger than or equal to a second amplitude threshold value, a second effective peak corresponding to the second amplitude threshold value is identified, the second effective peak meets second preset starting point and end point conditions, the peak value of the second effective peak is smaller than the first amplitude threshold value, and the difference value between the peak value of the second effective peak and the second amplitude threshold value is larger than or equal to a preset small pulse height threshold value.
2. The method of claim 1, wherein the first predetermined start and end point conditions comprise: a first pulse start point and a first pulse end point exist in the first effective peak;
wherein,
the first N data of the first pulse starting point are sequentially increased and are smaller than the data of the first pulse starting point, and the last M data of the first pulse starting point are sequentially increased and are larger than the data of the first pulse starting point;
the first N data of the first pulse end point are sequentially reduced and are larger than the data of the first pulse end point, the last M data of the first pulse end point are sequentially reduced and are smaller than the data of the first pulse end point, and both N and M are natural numbers which are larger than or equal to 1;
the first pulse start point and the first pulse end point are equal to the first amplitude threshold.
3. The method of claim 2, wherein the peak value of the first effective peak is equal to a maximum value in the plurality of data between the first pulse start point and the first pulse end point in the data of the pulse to be identified.
4. The method of claim 1, further comprising:
and if second data exists, and the second data is greater than or equal to the second amplitude threshold, starting to identify a first effective peak corresponding to the second amplitude threshold, wherein the first effective peak meets the second preset starting point and end point conditions, and the peak value is greater than the first amplitude threshold.
5. The method according to any one of claims 1 to 4, wherein the second predetermined start and end point conditions include: a second pulse start point and a second pulse end point exist in the second effective peak;
wherein,
the first N data of the second pulse starting point are sequentially increased and are smaller than the data of the second pulse starting point, and the last M data of the second pulse starting point are sequentially increased and are larger than the data of the second pulse starting point;
the first N data of the second pulse end point are sequentially increased and are smaller than the data of the second pulse end point, the last M data of the second pulse end point are sequentially increased and are larger than the data of the first pulse end point, and both N and M are natural numbers which are larger than or equal to 1;
the second pulse start point and the second pulse end point are equal to the second amplitude threshold.
6. The method of claim 5, wherein the peak value of the second effective peak is equal to a maximum value in the plurality of data between the start point of the second pulse and the end point of the second pulse in the data of the pulse to be identified.
7. The method of claim 1, wherein a width of the second effective peak is greater than or equal to a preset small pulse width threshold, and wherein the width of the second effective peak is a distance between a start point of the second pulse and an end point of the second pulse.
8. The method according to claim 1, characterized in that the difference between the second amplitude threshold and the baseline value of the pulse to be identified is smaller than a preset difference.
9. A pulse recognition apparatus, comprising:
the data processing module is used for comparing the pulse data to be identified with a first amplitude threshold value and a second amplitude threshold value one by one, wherein the first amplitude threshold value is larger than the second amplitude threshold value, and the first amplitude threshold value is larger than the maximum value of a preset small pulse amplitude interval;
the first identification module is used for starting to identify a first effective peak corresponding to a first amplitude threshold if first data exists and the first data is greater than or equal to the first amplitude threshold, wherein the first effective peak meets first preset starting point and end point conditions;
and the second identification module is used for starting to identify a second effective peak corresponding to the second amplitude threshold if second data exists, the second effective peak meets second preset starting point and end point conditions, the peak value of the second effective peak is smaller than the first amplitude threshold, and the difference value between the peak value of the second effective peak and the second amplitude threshold is larger than or equal to a preset small pulse height threshold.
10. The apparatus of claim 9, wherein the first identification module is further configured to:
and if second data exists, and the second data is greater than or equal to the second amplitude threshold, starting to identify a first effective peak corresponding to the second amplitude threshold, wherein the first effective peak meets the second preset starting point and end point conditions, and the peak value is greater than the first amplitude threshold.
11. An analysis apparatus, comprising a pulse recognition device according to claim 9 or 10.
12. A computer-readable storage medium, on which a program is stored, which, when being executed by a processor, carries out the pulse recognition method according to any one of claims 1 to 8.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811183186.3A CN109283121B (en) | 2018-10-11 | 2018-10-11 | Pulse recognition method and device, analysis instrument and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811183186.3A CN109283121B (en) | 2018-10-11 | 2018-10-11 | Pulse recognition method and device, analysis instrument and storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109283121A true CN109283121A (en) | 2019-01-29 |
| CN109283121B CN109283121B (en) | 2021-04-09 |
Family
ID=65176779
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811183186.3A Active CN109283121B (en) | 2018-10-11 | 2018-10-11 | Pulse recognition method and device, analysis instrument and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109283121B (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111679214A (en) * | 2020-06-24 | 2020-09-18 | 东莞新能安科技有限公司 | Apparatus and method for predicting the state of charge of a battery |
| CN111860251A (en) * | 2020-07-09 | 2020-10-30 | 迈克医疗电子有限公司 | Data processing method and device |
| CN112945807A (en) * | 2021-01-29 | 2021-06-11 | 深圳市科曼医疗设备有限公司 | Automatic detection method and device based on blood cell analyzer |
| CN114674729A (en) * | 2022-03-02 | 2022-06-28 | 迈克医疗电子有限公司 | Pulse recognition method, pulse recognition device, storage medium, apparatus, and blood cell analyzer |
| WO2024187341A1 (en) * | 2023-03-13 | 2024-09-19 | 中国科学院深圳先进技术研究院 | Pulse signal detection method, apparatus and device, and storage medium |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4068169A (en) * | 1976-09-21 | 1978-01-10 | Hycel, Inc. | Method and apparatus for determining hematocrit |
| US6112160A (en) * | 1996-04-10 | 2000-08-29 | Lecroy Corporation | Optical recording measurement package |
| CN103091253A (en) * | 2011-10-31 | 2013-05-08 | 深圳迈瑞生物医疗电子股份有限公司 | Method for eliminating interference side lobe signal, system thereof and blood cell analyzer |
| CN104678423A (en) * | 2015-03-10 | 2015-06-03 | 四川中测辐射科技有限公司 | Double-channel counting system and measurement method of dose equivalent in high dose condition |
| CN104736049A (en) * | 2012-09-11 | 2015-06-24 | 德尔格医疗系统有限公司 | A system and method for detecting a characteristic in an ECG waveform |
| US20150301086A1 (en) * | 2012-11-12 | 2015-10-22 | Rohde & Schwarz Gmbh & Co. Kg | A method and a device for determining a trigger condition for a rare signal event |
| CN106574891A (en) * | 2014-08-06 | 2017-04-19 | 贝克曼考尔特公司 | Evaluation of multi-peak events using a flow cytometer |
-
2018
- 2018-10-11 CN CN201811183186.3A patent/CN109283121B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4068169A (en) * | 1976-09-21 | 1978-01-10 | Hycel, Inc. | Method and apparatus for determining hematocrit |
| US6112160A (en) * | 1996-04-10 | 2000-08-29 | Lecroy Corporation | Optical recording measurement package |
| CN103091253A (en) * | 2011-10-31 | 2013-05-08 | 深圳迈瑞生物医疗电子股份有限公司 | Method for eliminating interference side lobe signal, system thereof and blood cell analyzer |
| CN104736049A (en) * | 2012-09-11 | 2015-06-24 | 德尔格医疗系统有限公司 | A system and method for detecting a characteristic in an ECG waveform |
| US20150301086A1 (en) * | 2012-11-12 | 2015-10-22 | Rohde & Schwarz Gmbh & Co. Kg | A method and a device for determining a trigger condition for a rare signal event |
| CN106574891A (en) * | 2014-08-06 | 2017-04-19 | 贝克曼考尔特公司 | Evaluation of multi-peak events using a flow cytometer |
| CN104678423A (en) * | 2015-03-10 | 2015-06-03 | 四川中测辐射科技有限公司 | Double-channel counting system and measurement method of dose equivalent in high dose condition |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111679214A (en) * | 2020-06-24 | 2020-09-18 | 东莞新能安科技有限公司 | Apparatus and method for predicting the state of charge of a battery |
| CN111860251A (en) * | 2020-07-09 | 2020-10-30 | 迈克医疗电子有限公司 | Data processing method and device |
| CN111860251B (en) * | 2020-07-09 | 2023-09-15 | 迈克医疗电子有限公司 | Data processing method and device |
| CN112945807A (en) * | 2021-01-29 | 2021-06-11 | 深圳市科曼医疗设备有限公司 | Automatic detection method and device based on blood cell analyzer |
| CN114674729A (en) * | 2022-03-02 | 2022-06-28 | 迈克医疗电子有限公司 | Pulse recognition method, pulse recognition device, storage medium, apparatus, and blood cell analyzer |
| CN114674729B (en) * | 2022-03-02 | 2023-11-21 | 迈克医疗电子有限公司 | Pulse identification method, pulse identification device, pulse identification storage medium, pulse identification equipment and blood cell analyzer |
| WO2024187341A1 (en) * | 2023-03-13 | 2024-09-19 | 中国科学院深圳先进技术研究院 | Pulse signal detection method, apparatus and device, and storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109283121B (en) | 2021-04-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109283121B (en) | Pulse recognition method and device, analysis instrument and storage medium | |
| US7860688B2 (en) | Signal baseline processing device and processing method thereof | |
| CN106353548B (en) | Method and apparatus for classifying signal waveforms to be measured | |
| CN107101984A (en) | Signal waveform characteristic detection method, device, storage medium and computer equipment | |
| CN1968228A (en) | Waveform recognition method in signal processing | |
| CN106908655A (en) | Photosignal peak-value detection method and system | |
| CN114185030B (en) | A radar signal intra-pulse modulation identification method and system | |
| CN1988383A (en) | Method for identifying wave form in signal processing | |
| CN104270331A (en) | BPSK signal blind identification result effectiveness evaluation method based on HOUGH transformation | |
| CN115078616A (en) | Multi-window spectral peak identification method, device, medium and product based on signal-to-noise ratio | |
| CN113848391A (en) | Pulse signal detection and extraction method | |
| KR101300650B1 (en) | A method of recognizing pri modulation type based on support vector machine, and an apparatus of recognizing pri modulation type | |
| CN118759491A (en) | A smoke and dust interference-resistant detection method based on narrow pulse width laser | |
| CN118707442A (en) | Signal arrival time estimation method, medium and device based on envelope amplitude judgment | |
| CN117491947A (en) | Object trajectory tracking method, device, equipment and storage medium | |
| CN110632563B (en) | A method for measuring parameters of intrapulse frequency encoded signal based on short-time Fourier transform | |
| CN102004177B (en) | Oscilloscope and method for identifying time sequence of universal serial bus signals by using same | |
| CN107179303A (en) | Droplet fluorescence detection method, device, system, storage medium and computer equipment | |
| KR102733677B1 (en) | A method for determining a single sampling frequency of the classification target signal in order to predict a modulation scheme of the classification target signal, and a method and apparatus for predicting the modulation scheme using the classification target signal sampled with the single sampling frequency | |
| CN109412901B (en) | Method and system for detecting continuity of acquired data based on time domain processing | |
| CN114675273A (en) | Method, device, electronic device and vehicle for suppressing interference signal | |
| JPH0843539A (en) | Received signal processing method used for buried object detection | |
| US7006925B2 (en) | Method of determining a base sequence for nucleic acid | |
| US12222415B2 (en) | Mobile device for detecting object and operation method therefor | |
| CN106963386B (en) | eye blinking recognition method and device |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |