Baseline Processing Method, Identification Method and Processing Device of Pulse Signal

Abstract

The present application discloses a baseline processing method, an identification method, and a processing device for pulse signals. The baseline processing method includes: sequentially obtaining data of each sampling point in the pulse signal to be processed; obtaining the average value of the data of the current sampling point and multiple sampling points within its predetermined range; determining whether the data of the current sampling point is greater than the corresponding average value; if so, using the baseline value at the position of the previous sampling point as the baseline value at the position of the current sampling point; otherwise, obtaining the baseline value at the position of the current sampling point by using the baseline value of the previous sampling point, the data of the previous sampling point, and the hysteresis; wherein, the hysteresis is from the same source as the pulse signal to be processed and is related to the pulse signal after the baseline has been removed; subtracting the corresponding baseline value from the data of the current sampling point. By the above method, the baseline value obtained at the current sampling point can better follow the baseline, and at the same time, the computational amount in the process of processing the baseline of the pulse signal can be effectively reduced, releasing more FPGA computing space.

Description

Technical Field

[0001] This application relates to the technical field of signal processing, and particularly to a baseline processing method, an identification method, and a processing device for pulse signals. Background Art

[0002] During the pulse identification process, since the original signal is easily affected by hardware background elevation, circuit noise, optical background, optical noise, etc. and there are many interferences, it is necessary to first calculate the baseline of the original data before pulse identification, and this baseline will calibrate the background value of the original data and filter out noise.

[0003] Currently, common methods for calculating the baseline value include mean filtering algorithm, median algorithm, Gaussian filtering algorithm, etc. However, among them, the mean filtering algorithm is easily affected by the particle concentration, resulting in excessive excision of valid pulse signals and filtering out smaller signal values, leading to generally low peak values, and further resulting in inaccurate pulse identification process. In addition, since the multi-channel pulse signals are detected by multiple sensors from different angles, using the existing technology to identify each single channel separately has a large operation dimension and occupies too much operation space of the FPGA chip. Summary of the Invention

[0004] This application provides a baseline processing method, an identification method, and a processing device for pulse signals, which can make the baseline value at the current sampling point more follow the baseline, effectively improve the accuracy of the pulse signal identification process, and at the same time can effectively reduce the calculation amount in the pulse signal baseline processing process and release more FPGA calculation space.

[0005] To solve the above technical problems, a technical solution adopted in this application is: providing a baseline processing method for pulse signals, including: sequentially obtaining the data of each sampling point in the pulse signal to be processed; obtaining the average value of the data of the current sampling point and multiple sampling points within a predetermined range around the current sampling point; determining whether the data of the current sampling point is greater than the corresponding average value; if so, using the baseline value at the previous sampling point position as the baseline value at the current sampling point position; otherwise, obtaining the baseline value at the current sampling point position by using the baseline value at the previous sampling point position, the data of the previous sampling point, and the hysteresis; wherein, the hysteresis is of the same source as the pulse signal to be processed and is related to the pulse signal after the baseline has been removed; subtracting the corresponding baseline value from the data of the current sampling point.

[0006] Among them, the step of obtaining the baseline value at the current sampling point position by using the baseline value at the previous sampling point position, the data of the previous sampling point, and the hysteresis includes: obtaining the following degree of the current sampling point by using the hysteresis and the baseline value of the previous sampling point, and referring to the data of the previous sampling point to obtain an average value; using the average value as the baseline value at the current sampling point position.

[0007] Among them, before the step of using the average value as the baseline value at the current sampling point position, it further includes: performing a rounding operation on the average value.

[0008] Among them, the predetermined range is related to the smoothness, and the smoothness is related to the pulse signal that has the same source as the pulse signal to be processed and has had the baseline removed.

[0009] Among them, the step of obtaining the average value of the data of the current sampling point and multiple sampling points within a predetermined range around the current sampling point includes: obtaining the average value of the data of multiple sampling points within one smoothness range before and after the current sampling point.

[0010] Among them, before the step of sequentially obtaining the data of each sampling point in the pulse signal to be processed, it includes: obtaining the median value of all peak widths in at least one pulse signal that has the same source as the pulse signal to be processed and has had the baseline removed; setting the hysteresis and the smoothness according to the median value.

[0011] Among them, the hysteresis is between half of the median value and the median value.

[0012] Among them, the smoothness is between twice the median value and eight times the median value.

[0013] To solve the above technical problems, another technical solution adopted by this application is: to provide a method for identifying multi-channel pulse signals, including: performing the baseline processing method mentioned in any of the above embodiments on the pulse signals of each channel; determining the peak information of the main pulse signal; using the peak information of the main pulse signal to determine the peak information of other sub-pulse signals with the same sampling time.

[0014] To solve the above technical problems, another technical solution adopted by this application is: to provide a pulse signal processing device, including a processor and a memory that are coupled to each other, and the processor and the memory cooperate with each other to implement the baseline processing method mentioned in any of the above embodiments or the method for identifying multi-channel pulse signals mentioned in any of the above embodiments.

[0015] The beneficial effects of the present application are as follows: The present application provides a method for processing the baseline of a pulse signal, which includes successively obtaining the data of each sampling point in the pulse signal to be processed; obtaining the average value of the data of the current sampling point and a plurality of sampling points within a predetermined range around the current sampling point; determining whether the data of the current sampling point is greater than the corresponding average value; if so, taking the baseline value at the position of the previous sampling point as the baseline value at the position of the current sampling point; otherwise, obtaining the baseline value at the position of the current sampling point by using the baseline value at the position of the previous sampling point, the data of the previous sampling point, and the hysteresis; wherein, the hysteresis has the same source as the pulse signal to be processed and is related to the pulse signal after the baseline has been removed; subtracting the corresponding baseline value from the data of the current sampling point. Through the above method, during the transmission of the pulse signal, based on the data of consecutive sampling points, the hysteresis, and the baseline value at the position of the previous sampling point, the baseline value at the position of the current sampling point can be obtained. This baseline value more closely follows the true baseline, can weaken the height of the original pulse signal, and can accommodate more small pulse signals, thereby further improving the accuracy of the pulse signal recognition process; compared with the mean filtering algorithm for the baseline value in the prior art, this method for processing the baseline of the pulse signal can effectively reduce the amount of calculation and release more FPGA computing space. Description of the Drawings

[0016] In order to more clearly illustrate the technical solutions in the embodiments of the present application, the following will briefly introduce the drawings required for the description of the embodiments. Obviously, the following drawings are only some embodiments of the present application. For those of ordinary skill in the art, without creative efforts, other drawings can also be obtained based on these drawings, where:

[0017] Figure 1 is a schematic flowchart of an embodiment of the method for processing the baseline of the pulse signal of the present application;

[0018] Figure 2 is Figure 1 a schematic flowchart of an embodiment before step S101 in

[0019] Figure 3 is Figure 1 a schematic flowchart of an embodiment of step S105;

[0020] Figure 4 is a schematic flowchart of an embodiment of the method for identifying multi-channel pulse signals in the present application;

[0021] Figure 5 is a schematic framework diagram of an embodiment of the device for processing the baseline of the pulse signal of the present application;

[0022] Figure 6It is a schematic structural diagram of an embodiment of the baseline processing method for the pulse signal of the present application;

[0023] Figure 7 It is a schematic framework diagram of an embodiment of the device with the function of baseline processing or recognition for the pulse signal of the present application. Specific embodiments

[0024] Next, the technical solutions in the embodiments of the present application will be clearly and completely described in conjunction with the accompanying drawings in the embodiments of the present application. Obviously, the described embodiments are only a part of the embodiments of the present application, rather than all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those of ordinary skill in the art without creative efforts shall fall within the protection scope of the present application.

[0025] Please refer to Figure 1 , Figure 1 It is a schematic flowchart of an embodiment of the baseline processing method for the pulse signal of the present application. This baseline processing method includes:

[0026] S101: Sequentially obtain the data of each sampling point in the pulse signal to be processed.

[0027] Specifically, the pulse signal to be processed can be formed by multiple sampling points. The abscissa of the sampling point is the sampling time, and the unit can be seconds, etc. The ordinate of the sampling point is the sampling data, which can be the light intensity value, etc. The data of each sampling point in the above step S101 refers to the ordinate data of each sampling point; and in the above step S101, the sampling points can be sequentially obtained in the order of decreasing or increasing sampling time.

[0028] S102: Obtain the average value of the data of the current sampling point and multiple sampling points within a predetermined range around the current sampling point.

[0029] Specifically, in the calculation process of the FPGA, since the sampling points are processed one by one, it is necessary to cache the data of the sampling point and the data of multiple sampling points within a predetermined range around it, so that the FPGA can only perform subsequent operations on this sampling point during the calculation.

[0030] In this embodiment, the above-mentioned predetermined range is related to the smoothness, and the smoothness is related to the pulse signal with the same source as the pulse signal to be processed and the baseline removed. Among them, the same source means that it is collected by the same sensor and the corresponding collection conditions (such as sheath flow velocity, etc.) are roughly the same. The smoothness here refers to the flatness and smoothness of the baseline. In order to ensure that the calculated baseline value follows the previous baseline and finally forms a complete and smooth baseline.

[0031] In this embodiment, the above step S102 further includes: obtaining the average value of the data of a plurality of sampling points within a smoothness range before and after the current sampling point. The smoothness range can be set according to needs. However, if the smoothness range is too large, it will cause a large amount of calculation; if the smoothness range is too small, the calculation result may not be accurate enough. The specific setting process of the smoothness will be introduced in detail in the subsequent embodiments and will not be elaborated here. Through the above embodiment, the average value of the sampling points within a suitable range can be obtained.

[0032] Of course, in other embodiments, it is also possible to obtain the average value of the data of a plurality of sampling points within two or three smoothness ranges before and after the sampling point. There is no specific limitation here, as long as the average value of the data of the sampling points within at least one smoothness range can be obtained.

[0033] S103: Determine whether the data of the current sampling point is greater than the corresponding average value.

[0034] S104: If so, use the baseline value at the position of the previous sampling point as the baseline value at the position of the current sampling point.

[0035] Specifically, when the data of the sampling point is greater than the corresponding average value, it indicates that the current sampling point is at the pulse position. Therefore, the baseline value of the non-pulse position before the current position can be directly used as the baseline value of the current position.

[0036] S105: If not, obtain the baseline value at the position of the current sampling point by using the baseline value at the position of the previous sampling point, the data of the previous sampling point, and the hysteresis; wherein, the hysteresis is from the same source as the pulse signal to be processed and is related to the pulse signal after the baseline has been removed.

[0037] Specifically, the hysteresis refers to the maximum deviation between the output values of the positive and negative strokes corresponding to the same input quantity during the same calibration when performing a full measurement range calibration under the same working conditions. Its value is expressed as the percentage of the maximum deviation or half of the maximum deviation to the full-scale output value. The numerical value of the hysteresis can be set according to the peak value and peak width of the pulse signal. The specific principle of calculating the baseline value based on the baseline value at the position of the previous sampling point, the data of the previous sampling point, and the hysteresis will be introduced in detail in the subsequent embodiments and will not be elaborated here.

[0038] S106: Subtract the corresponding baseline value from the data of the current sampling point.

[0039] Through the above-described embodiments, during the transmission of the pulse signal, based on the data of consecutive sampling points, the hysteresis, and the baseline value at the position of the previous sampling point, the baseline value at the current sampling point can be obtained. This baseline value better follows the true baseline, can weaken the height of the original pulse signal, and can accommodate more small pulse signals, thereby further improving the accuracy of the pulse signal recognition process. Compared with the mean filtering algorithm for the baseline value in the prior art, this baseline processing method for the pulse signal can effectively reduce the amount of calculation and release more FPGA computing space.

[0040] In this embodiment, please refer to Figure 2 , Figure 2 is Figure 1 a schematic flowchart of the previous embodiment before step S101. Before the above step S101, it further includes:

[0041] S201: Obtain the median value of all peak widths in at least one pulse signal that has the same source as the pulse signal to be processed and has had its baseline removed.

[0042] Specifically, for the signals obtained from different channels, their hysteresis and smoothness should be set correspondingly. The process of obtaining the median value can be as follows: obtain all the peak widths in at least one pulse signal that has the same source as the pulse signal to be processed and has had its baseline removed; plot a peak width distribution histogram; and obtain the median value according to this peak width distribution histogram.

[0043] S202: Set the hysteresis and smoothness according to the median value.

[0044] In this embodiment, the hysteresis is between half of the median value and the median value. The size of the hysteresis is specifically adjusted according to the peak value and peak width of the pulse signal. Generally speaking, when the median value is 21, the above hysteresis should be greater than or equal to 4. Preferably, the lower the hysteresis, the stronger the followability, and the more the peak value will be reduced; the higher the hysteresis, the worse the followability, and the more the peak value will be retained, but it will follow better during the blank signal and the noise suppression performance is worse. Considering the above factors comprehensively, the above hysteresis can be selected as 32. Through the above embodiments, the appropriate hysteresis can be adjusted as needed.

[0045] In this embodiment, the smoothness is between twice the median value and eight times the median value. Through the above embodiments, the appropriate smoothness value can be adjusted as needed. And when there are more pulses in the pulse signal to be processed, the smoothness can be set smaller.

[0046] In this embodiment, please refer to Figure 3 , Figure 3 is Figure 1 a schematic flowchart of one embodiment of step S105. The above step S105 includes:

[0047] S301: Obtain the following degree of the current sampling point using the hysteresis and the baseline value of the previous sampling point, and refer to the data of the previous sampling point to obtain the mean value.

[0048] Specifically, the above mean value is calculated using the following degree of the current sampling point and the data of the previous sampling point. Calculating the mean value based on the hysteresis and the baseline value of the previous sampling point can ensure that the current sampling point follows the previous baseline and there will be no numerical mutation; referring to the background value of the previous sampling point can ensure that the obtained mean value does not perturb the previous baseline.

[0049] S302: Take the mean value as the baseline value at the position of the current sampling point.

[0050] Through the above implementation manner, a filtering process with hysteresis smoothing can be achieved. The baseline value at the current pulse signal is calculated based on the data, baseline value, and hysteresis of the previous sampling point, without deviating from the original trajectory, without perturbing and mutating the original baseline, making the newly obtained baseline closer to the true baseline.

[0051] In another implementation manner, before the above step S302, it further includes: performing a rounding-up operation on the mean value. Specifically, the rounding-up operation here refers to rounding up. As long as there is a decimal, the integer value in front is incremented by one. This operation step is to be compatible with the calculation method of the FPGA. The way of rounding up the mean value is not specifically limited here, as long as the effect of rounding up can be achieved. For example, the hysteresis can be used to perform a rounding-up operation on the data of the current sampling point. Through the above implementation manner, the calculated mean value can achieve the effect of rounding up by the FPGA. In this implementation manner, please refer to Figure 4 , Figure 4 is a schematic flowchart of an implementation manner of the multi-channel pulse signal recognition method in this application. The steps of the above multi-channel pulse signal recognition method include:

[0052] S401: Perform the above baseline processing method on the pulse signals of each channel.

[0053] S402: Determine the peak information of the main pulse signal.

[0054] Specifically, first determine the starting position of the main pulse signal. For example, a threshold can be set, and the amplitude of each sampling point of the main pulse signal is obtained in sequence. When the amplitude of the main pulse signal is greater than the set threshold, this point is taken as the starting position of the main pulse signal. Of course, in other embodiments, the amplitudes of multiple sampling points can also be continuously obtained. When the amplitudes of the continuous multiple sampling points all meet the above requirements, one of the sampling points can be taken as the starting position. For example, when the amplitudes of 3 consecutive sampling points are greater than the set threshold, any one of the three sampling points can be taken as the starting position.

[0055] Determine the main pulse peak point of the main pulse signal based on the starting position. For example, after obtaining the starting position of the main pulse signal, take three consecutive sampling points. If the value of the middle sampling point is greater than the values of the sampling points on both sides, then use it as the main pulse peak point.

[0056] Determine the main pulse valley point of the main pulse signal based on the starting position. For example, after obtaining the starting position of the main pulse signal, take three consecutive sampling points. If the value of the middle sampling point is less than the values of the sampling points on both sides, then use it as the main pulse valley point. Among them, if the valley value is less than 90% of the peak value or there is no valley point in this main pulse signal, directly enter the step of determining the end point.

[0057] Find the highest peak point based on multiple peak points in the pulse signal. For example, take the maximum value among all peak points as the highest peak point corresponding to this pulse signal.

[0058] After determining the data of the highest peak point of the main pulse, the end point position can also be determined. For example, a threshold can be set. When the data of the current sampling point is less than the set threshold, determine the current sampling point as the end position.

[0059] Determine the pulse width data of the main pulse signal based on the starting position and the end position of the main pulse signal. Among them, determine the front peak width data of the main pulse signal based on the number of sampling points between the starting position and the sampling point corresponding to the main pulse peak data; determine the rear peak width data of the main pulse signal based on the number of sampling points between the sampling point corresponding to the main pulse peak data and the end position; determine the full pulse width data of the main pulse signal based on the number of sampling points obtained between the starting position and the end position.

[0060] S403: Determine the peak information of other sub-pulse signals with the same sampling time by using the peak information of the main pulse signal.

[0061] Specifically, use the peak time and peak width of the pre-scattered signal to find the peak information of other sub-pulse signals.

[0062] In another embodiment, after the above step S403, it includes performing normalization calculation on all obtained peaks, which can achieve a better display effect. It can also include judging whether the pulse peak data and pulse width data are valid. For example, it can be judged by whether the pulse peak of the target pulse signal is greater than the set peak threshold and whether the pulse width of the target pulse signal is greater than the set pulse width threshold. The above method can judge whether the current pulse is a valid pulse under the conditions of effective pulse amplitude, effective front peak width, effective rear peak width, and effective full peak width.

[0063] Please refer to Figure 5 , Figure 5It is a schematic framework diagram of an embodiment of the baseline processing device for pulse signals in the present application. The device 100 includes an acquisition module 10, a judgment module 12, and a processing module 14. Among them, the acquisition module 10 is used to sequentially acquire the data of each sampling point in the pulse signal to be processed and the average value of the data of the current sampling point and a plurality of sampling points within a predetermined range around the current sampling point. The judgment module 12 is used to judge whether the data of the current sampling point is greater than the corresponding average value. The processing module 14 is used to, in response to the data of the current sampling point being greater than the corresponding average value, take the baseline value at the position of the previous sampling point as the baseline value at the position of the current sampling point. The processing module 24 is further used to, in response to the data of the current sampling point being less than or equal to the corresponding average value, obtain the baseline value at the position of the current sampling point by using the baseline value at the position of the previous sampling point, the data of the previous sampling point, and the hysteresis, wherein the hysteresis has the same source as the pulse signal to be processed and is related to the pulse signal after the baseline has been removed. The processing module 24 is further used to subtract the corresponding baseline value from the data of the current sampling point. Through the above embodiment, during the transmission of the pulse signal, based on the data of continuous sampling points, the hysteresis, and the baseline value at the position of the previous sampling point, the baseline value at the current sampling point can be obtained. This baseline value follows the true baseline more closely, can weaken the height of the original pulse signal, and can accommodate more small pulse signals, thereby further improving the accuracy of the pulse signal recognition process. Compared with the mean filtering algorithm for the baseline value in the prior art, this baseline processing method for pulse signals can effectively reduce the amount of calculation and release more FPGA computing space.

[0064] Please refer to Figure 6 , Figure 6 is a schematic structural diagram of an embodiment of the baseline processing method for pulse signals in the present application. The device 200 includes a memory 20 and a processor 22 that are coupled to each other. The memory 20 stores program instructions, and the processor 22 is used to execute the program instructions to implement the baseline processing and recognition methods for pulse signals mentioned in any of the above embodiments.

[0065] Specifically, the processor 22 can also be referred to as a CPU (Central Processing Unit). The processor 22 may be an integrated circuit chip with signal processing capabilities. The processor 22 can also be a general-purpose processor, a digital signal processor (DSP), an application-specific integrated circuit (ASIC), a field-programmable gate array (FPGA), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components. The general-purpose processor can be a microprocessor or the processor can also be any conventional processor, etc. Additionally, the processor 22 can be implemented jointly by multiple integrated circuit chips.

[0066] Please refer to Figure 7 , Figure 7 which is a schematic framework diagram of an embodiment of the device with the function of baseline processing or recognition of pulse signals in the present application. The device 300 has program data 30 that can be read by a computer, and the program data 30 can also be executed by the processor. The program data 30 is used to implement the negotiation method mentioned in any of the above embodiments. Among them, the program data 30 can be stored in the above-mentioned device 300 with storage function in the form of a software product, including several instructions to enable a computer device (which can be a personal computer, a server, or a network device, etc.) or a processor to execute all or part of the steps of the methods described in various embodiments of the present application. The aforementioned storage device includes: various media that can store program codes, such as USB flash drives, mobile hard disks, read-only memories (ROM, Read-Only Memory), random access memories (RAM, Random Access Memory), magnetic disks, or optical discs, or terminal devices such as computers, servers, mobile phones, and tablets.

[0067] In summary, different from the prior art, the present application can obtain the baseline value at the current sampling point based on the data of continuous sampling points, the hysteresis, and the baseline value at the position of the previous sampling point during the transmission of pulse signals. This baseline value follows the true baseline more closely, can weaken the height of the original pulse signal, and can accommodate more small pulse signals, thereby further improving the accuracy of the pulse signal recognition process; compared with the mean filtering algorithm of the baseline value in the prior art, this baseline processing method for pulse signals can effectively reduce the amount of calculation and release more FPGA computing space.

[0068] The above are only embodiments of the present application, and do not thereby limit the patent scope of the present application. Any equivalent structure or equivalent process transformation made by using the content of the specification and drawings of the present application, or directly or indirectly applied in other related technical fields, shall similarly be included within the patent protection scope of the present application.

Claims

1. A method for processing the baseline of a pulse signal, characterized in that, Comprising: Successively obtaining data of each sampling point in a pulse signal to be processed; Obtaining an average value of data of the current sampling point and a plurality of sampling points within a predetermined range around the current sampling point; Determining whether the data of the current sampling point is greater than the corresponding average value; If so, using the baseline value at the position of the previous sampling point as the baseline value at the position of the current sampling point; otherwise, obtaining the baseline value at the position of the current sampling point by using the baseline value at the position of the previous sampling point, the data of the previous sampling point, and a hysteresis; wherein, the hysteresis is from the same source as the pulse signal to be processed and is related to the pulse signal after baseline removal; Subtracting the corresponding baseline value from the data of the current sampling point; The step of obtaining the baseline value at the position of the current sampling point by using the baseline value at the position of the previous sampling point, the data of the previous sampling point, and a hysteresis includes: Obtaining a following degree of the current sampling point by using the hysteresis and the baseline value of the previous sampling point, and referring to the data of the previous sampling point to obtain an average value; Using the average value as the baseline value at the position of the current sampling point.

2. The baseline processing method according to claim 1, wherein Before the step of using the average value as the baseline value at the position of the current sampling point, further including: Performing a rounding operation on the average value.

3. The baseline processing method according to claim 1, wherein The predetermined range is related to smoothness, and the smoothness is from the same source as the pulse signal to be processed and is related to the pulse signal after baseline removal.

4. The baseline processing method according to claim 3, characterized in that, The step of obtaining an average value of data of the current sampling point and a plurality of sampling points within a predetermined range around the current sampling point includes: Obtaining an average value of data of a plurality of sampling points within one smoothness range before and after the current sampling point.

5. The baseline processing method according to claim 3, wherein Before the step of successively obtaining data of each sampling point in the pulse signal to be processed, including: Obtaining a median value of all peak widths in at least one pulse signal that is from the same source as the pulse signal to be processed and after baseline removal; Setting the hysteresis and the smoothness according to the median value.

6. The baseline processing method according to claim 5, wherein The hysteresis is between 1 / 2 of the median value and the median value.

7. The baseline processing method according to claim 5, wherein The smoothness is between 2 times of the median value and 8 times of the median value.

8. A method for identifying multi-channel pulse signals, characterized in that, Comprising: Performing the baseline processing method according to any one of claims 1-7 on the pulse signal of each channel; Determining peak information of the main pulse signal; Using the peak information of the main pulse signal to determine peak information of other sub-pulse signals with the same sampling time.

9. A processing device for pulse signals, characterized in that, Comprising a processor and a memory that are mutually coupled, and the processor and the memory cooperate with each other to implement the baseline processing method according to any one of claims 1-7 or the identification method of multi-channel pulse signals according to any one of claims 8.

Images