CN110477889B - Method and device for identifying reflection point of radial artery pressure waveform - Google Patents

Method and device for identifying reflection point of radial artery pressure waveform Download PDF

Info

Publication number
CN110477889B
CN110477889B CN201811168960.3A CN201811168960A CN110477889B CN 110477889 B CN110477889 B CN 110477889B CN 201811168960 A CN201811168960 A CN 201811168960A CN 110477889 B CN110477889 B CN 110477889B
Authority
CN
China
Prior art keywords
waveform
point
differential signal
starting point
value
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201811168960.3A
Other languages
Chinese (zh)
Other versions
CN110477889A (en
Inventor
张永亮
叶骏
李道清
张启莲
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
BEIJING DONGLIANG HEALTH TECHNOLOGY Co.,Ltd.
Hefei Yiyang Health Technology Co.,Ltd.
Original Assignee
Beijing Dongliang Health Technology Co ltd
Hefei Yiyang Health Technology Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing Dongliang Health Technology Co ltd, Hefei Yiyang Health Technology Co ltd filed Critical Beijing Dongliang Health Technology Co ltd
Priority to CN201811168960.3A priority Critical patent/CN110477889B/en
Publication of CN110477889A publication Critical patent/CN110477889A/en
Application granted granted Critical
Publication of CN110477889B publication Critical patent/CN110477889B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, 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/021Measuring pressure in heart or blood vessels
    • A61B5/02108Measuring pressure in heart or blood vessels from analysis of pulse wave characteristics
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7225Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/725Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Veterinary Medicine (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • Public Health (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Physiology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Psychiatry (AREA)
  • Cardiology (AREA)
  • Power Engineering (AREA)
  • Vascular Medicine (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Abstract

The application discloses a method and a device for identifying reflection points of a radial artery pressure waveform. The method comprises the following steps: filtering the radial pulse waveform; generating a differential signal of the radial pulse waveform, and determining a maximum value of the differential signal; determining a starting point of the radial pulse waveform based on a maximum of the differential signal; determining a period of the radial pulse waveform based on the starting point; performing baseline wandering processing on the radial pulse waveform based on the starting point and calibrating the waveform; and identifying reflection points of the calibrated waveform. By adopting the method, analysis can be completely carried out based on signal waveforms without adopting additional hardware, and compared with the complicated and fussy analysis process in the prior art, the method reduces the complexity and difficulty of analysis and improves the accuracy of analysis results.

Description

Method and device for identifying reflection point of radial artery pressure waveform
Technical Field
The present application relates to the field of waveform detection, and in particular, to a method and apparatus for identifying reflection points of a radial artery pressure waveform.
Background
The pulse is an important biodynamic signal for representing the heart function and the blood circulation system of the human body, and integrates the pumping activity of the heart and various physiological and pathological information carried in the transmission process of pressure waves along the whole body arterial tree. The existing clinical experiment results show that the waveform characteristics of the pulse wave are closely related to the cardiovascular diseases. The pulse wave signals are detected by the pulse wave sensor, and various cardiovascular function parameters are analyzed by combining the waveform characteristics of the pulse waves and the blood pressure information. The method has the characteristics of simple and convenient operation, no wound, comprehensive detection parameters and the like, and is widely used in practical application.
Currently, the pulse waveform characteristics of major superficial arteries of a human body, such as the vagus artery, the radial artery, the carotid artery, etc., are commonly used in clinic, wherein the radial artery measurement of the wrist is the most common one. The method is simple to apply in clinical environment, low in measurement cost, free of sequelae, capable of continuously detecting for a long time and accurate in obtained measurement result. Corresponding characteristic parameters can be extracted from the measured pulse wave waveform through a series of signal processing methods, so that relevant index parameters for evaluating the cardiovascular function are obtained. Therefore, a practical method for accurately extracting waveform feature parameters from a measured pulse waveform is needed.
The method comprises the steps of collecting human body radial artery pulse wave signals through a photoelectric sensor, carrying out filtering amplification and AD conversion on the signals by using a microprocessor to reduce noise, then carrying out frequency domain transformation, wavelet processing and second-order difference to determine characteristic points of human body radial artery pulse wave waveforms, calculating accurate numerical values of human body radial artery pulse wave conduction time, establishing a regression equation for measuring blood pressure, and finally obtaining blood pressure values. The processing method adopted when the waveform characteristics are obtained by the method is complex and tedious, and more hardware resources are consumed, so that a method which can be used for extracting the characteristic point parameters more easily is needed, and the subsequent analysis and processing process is facilitated.
Disclosure of Invention
It is an object of the present application to overcome the above problems or to at least partially solve or mitigate the above problems.
According to one aspect of the present application, there is provided a method of identifying a reflection point of a radial artery pressure waveform, comprising:
a filtering step: filtering the radial pulse waveform;
a differential signal generation step: generating a differential signal of the radial pulse waveform, and determining a maximum value of the differential signal;
a starting point searching step: determining a starting point of the radial pulse waveform using the maxima and the differential signal;
a period determining step: determining a period of the radial pulse waveform based on the starting point;
waveform calibration: performing baseline wandering processing on the radial pulse waveform based on the starting point and calibrating the waveform; and
a characteristic point identification step: reflection points of the calibrated waveform are identified.
By adopting the method, analysis can be completely carried out based on signal waveforms without adopting additional hardware, and compared with the complicated and fussy analysis process in the prior art, the method reduces the complexity and difficulty of analysis and improves the accuracy of analysis results. The method utilizes algorithm operations such as convolution filtering, first-order difference, convolution and the like, reduces the complex process of data processing to a certain extent, reduces the hardware computing resources required to be occupied, enables the feature points of the radial pulse waveform to be easier to extract, and is beneficial to the subsequent analysis processing process.
Optionally, the smoothing filtering step includes: and performing convolution filtering processing on the radial pulse waveform and a preset two-dimensional array.
Optionally, in the differential signal generating step, the calculation formula for generating the differential signal is:
vtCpDifSig[i]=vtSmoothedData[i+1]-vtSmoothedData[i-1]+2*(vtSmoothedData[i+2]-vtSmoothedData[i-2])
wherein vtCpDifSig represents a differential signal; vtsmootheneddata represents the filtered radial pulse waveform; i is an outer loop variable representing the ith bit signal value.
The method adopts a mathematical calculation mode to calculate the differential signal, can realize operation without adopting hardware of a differential circuit, and is convenient, quick, simple and realizable.
Optionally, in the differential signal generating step, the method for determining the waveform maximum value includes:
and (3) a step of segmented searching: dividing the differential signal waveform into eight sections, and respectively searching a maximum value point of each section of waveform;
and an average value calculation step: calculating an average value AvgSegmentMax of maximum values of the eight sections of waveforms; and
a maximum value determining step: and searching a maximum value point in the whole differential signal waveform, and determining the value of the maximum value point as the maximum value point of the differential signal waveform under the condition that the amplitude of the maximum value point is greater than 0.6 AvgSegmentMax.
By adopting the method, the maximum value is searched in a segmented mode, the operation speed can be improved in a parallel processing mode, the obtained average value of the maximum value can be used as a threshold value to judge the maximum value of the whole waveform, and the judgment accuracy is improved.
Optionally, the starting point searching step includes: determining a search area based on the maximum value point of the differential signal, searching a forward zero crossing point in the search area, and taking the zero crossing point as a waveform starting point under the condition that the forward zero crossing point exists; otherwise, the minimum value point is used as the starting point of the waveform.
Optionally, in the period determining step, the waveform period is determined by calculating a distance difference between two starting points, and the waveform period is calculated by:
m_vtPeriod[i]=m_vtFootPos[i]-m_vtFootPos[i-1]
wherein m _ vtPeriod is a waveform period; m _ vtFootPos is used as a starting point; i is an outer loop variable representing the ith bit signal value.
Due to the adoption of the technical scheme, compared with the prior art, the method and the device reduce the processing process of pulse waveform data to a certain extent, reduce the consumption of computing resources, and are easier to extract waveform characteristic points, so that the subsequent analysis processing process of utilizing pulse waveforms is more convenient.
Optionally, in the waveform calibration step, a baseline wander process is performed according to a slope of the waveform start point value and a waveform offset, where the slope is calculated by:
Figure BDA0001821919510000031
the waveform offset calculation formula is as follows:
ΔvtSlope=vtSlope[i]*(j-m_vtFootPos[i])
wherein vtSlope represents the slope; Δ vtSlope represents the waveform offset; vtsmootheneddata represents the filtered radial pulse waveform; m _ vtFootPos is used as a starting point; i is an outer loop variable and represents the ith signal value; j is an inner loop variable.
Optionally, the feature point identifying step includes:
and (3) convolution step: acquiring the fourth derivative of the calibrated waveform through convolution based on a set scale;
a reflection point determining step: searching the first reverse zero crossing point of the fourth derivative in the region from the starting point to 0.67 waveform periods, wherein the zero crossing point is a reflection point.
According to another aspect of the present application, there is also provided an apparatus for identifying a reflection point of a radial artery pressure waveform, including:
a filtering module configured to perform filtering processing on a radial artery pulse waveform;
a differential signal generation module configured to generate a differential signal of the radial pulse waveform, determine a maximum of the differential signal;
a starting point search module configured to determine a starting point of the radial pulse waveform using the maximum value and the differential signal;
a period determination module configured to determine a period of the radial pulse waveform based on the starting point;
a waveform calibration module configured to perform de-baseline wander processing on the radial pulse waveform and calibrate a waveform based on the starting point;
a feature point identification module configured to identify reflection points of the calibrated waveform.
By adopting the device, analysis can be carried out completely based on signal waveforms, and can be realized without adopting extra hardware, compared with the complicated and fussy analysis process in the prior art, the complexity and difficulty of analysis are reduced, and the accuracy of an analysis result is also improved. The method utilizes algorithm operations such as convolution filtering, first-order difference, convolution and the like, reduces the complex process of data processing to a certain extent, reduces the hardware computing resources required to be occupied, enables the feature points of the radial pulse waveform to be easier to extract, and is beneficial to the subsequent analysis processing process.
Optionally, the differential signal generating module includes:
a segment searching module configured to divide the radial artery pulse waveform into eight segments, and search a maximum value point of each segment of the waveform respectively;
an average value calculating module configured to calculate an average value AvgSegmentMax of maxima of the eight-segment waveform; and
a maximum determination module configured to search for a maximum point within an entire radial pulse waveform, and in the event that the magnitude of the maximum point is greater than 0.6 AvgSegmentMax, determine the value of the maximum point as the maximum point of the radial pulse waveform.
The above and other objects, advantages and features of the present application will become more apparent to those skilled in the art from the following detailed description of specific embodiments thereof, taken in conjunction with the accompanying drawings.
Drawings
The accompanying drawings, which are included to provide a further understanding of the application and are incorporated in and constitute a part of this application, illustrate embodiment(s) of the application and together with the description serve to explain the application and not to limit the application. Some specific embodiments of the present application will be described in detail hereinafter by way of illustration and not limitation with reference to the accompanying drawings. The same reference numbers in the drawings identify the same or similar elements or components. Those skilled in the art will appreciate that the drawings are not necessarily drawn to scale. In the drawings:
FIG. 1 is a diagram illustrating a hardware configuration of a computer device for performing a method for identifying reflection points of a radial artery pressure waveform according to an embodiment of the present application;
FIG. 2 is a schematic flow chart diagram of a method of identifying reflection points of a radial artery pressure waveform according to one embodiment of the present application;
FIG. 3 is a schematic overview flow chart of a method of identifying reflection points of a radial artery pressure waveform in accordance with an embodiment of the present application;
FIG. 4 is a schematic flow chart diagram of a method of searching for a reflection point according to an embodiment of the present application;
FIG. 5 is a schematic illustration of a location of a reflection point on a radial pulse waveform according to an embodiment of the present application;
FIG. 6 is a schematic block diagram of a device for identifying reflection points of a radial artery pressure waveform according to one embodiment of the present application;
FIG. 7 is a block diagram of one embodiment of a computing device of the present application;
FIG. 8 is a block diagram of one embodiment of a computer-readable storage medium of the present application.
Detailed Description
In order to make the technical solutions better understood by those skilled in the art, the technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application, and it is obvious that the described embodiments are only partial embodiments of the present application, but 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 application.
It should be noted that the terms "first," "second," and the like in the description and claims of this application and in the drawings described above are used for distinguishing between similar elements and not necessarily for describing a particular sequential or chronological order. It is to be understood that the data so used is interchangeable under appropriate circumstances such that the embodiments of the application described herein are capable of operation in sequences other than those illustrated or described herein. Furthermore, the terms "comprises," "comprising," and "having," and any variations thereof, are intended to cover a non-exclusive inclusion, such that a process, method, system, article, or apparatus that comprises a list of steps or elements is not necessarily limited to those steps or elements expressly listed, but may include other steps or elements not expressly listed or inherent to such process, method, article, or apparatus.
Example 1
There is also provided, in accordance with an embodiment of the present application, an embodiment of a method for identifying characteristic points of a radial artery pressure waveform, where it is noted that the steps illustrated in the flowchart of the drawings may be performed in a computer system, such as a set of computer-executable instructions, and that while a logical order is illustrated in the flowchart, in some cases the steps illustrated or described may be performed in an order different than here.
The method provided by the first embodiment of the present application may be executed in a mobile terminal, a computer terminal, or a similar computing device. Fig. 1 shows a block diagram of a hardware structure of a computer apparatus (or mobile device) for the method of the present application. As shown in fig. 1, computer apparatus 10 (or mobile device 10) may include one or more processors (shown as 102a, 102b, … …, 102n, which may include, but are not limited to, a processing device such as a microprocessor MCU or a programmable logic device FPGA), a memory 104 for storing data, and a transmission module for communication functions. Besides, the method can also comprise the following steps: a display, an input/output interface (I/O interface), a Universal Serial Bus (USB) port (which may be included as one of the ports of the I/O interface), a network interface, a power source, and/or a camera. It will be understood by those skilled in the art that the structure shown in fig. 1 is only an illustration and is not intended to limit the structure of the electronic device. For example, computer device 10 may also include more or fewer components than shown in FIG. 1, or have a different configuration than shown in FIG. 1.
It should be noted that the one or more processors and/or other data processing circuitry described above may be referred to generally herein as "data processing circuitry". The data processing circuitry may be embodied in whole or in part in software, hardware, firmware, or any combination thereof. Further, the data processing circuit may be a single stand-alone processing module, or incorporated in whole or in part into any of the other elements in the computer apparatus 10 (or mobile device). As referred to in the embodiments of the application, the data processing circuit acts as a processor control (e.g. selection of a variable resistance termination path connected to the interface).
The memory 104 can be used for storing software programs and modules of application software, such as program instructions/data storage devices corresponding to the methods in the embodiments of the present application, and the processor executes various functional applications and data processing by executing the software programs and modules stored in the memory 104, that is, implementing the methods of the application programs described above. The memory 104 may include high speed random access memory, and may also include non-volatile memory, such as one or more magnetic storage devices, flash memory, or other non-volatile solid-state memory. In some examples, memory 104 may further include memory located remotely from the processor, which may be connected to computer device 10 via a network. Examples of such networks include, but are not limited to, the internet, intranets, local area networks, mobile communication networks, and combinations thereof.
The transmission device is used for receiving or transmitting data via a network. Specific examples of the network described above may include a wireless network provided by a communication provider of computer device 10. In one example, the transmission device includes a Network adapter (NIC) that can be connected to other Network devices through a base station to communicate with the internet. In one example, the transmission device may be a Radio Frequency (RF) module, which is used for communicating with the internet in a wireless manner.
The display may be, for example, a touch screen type Liquid Crystal Display (LCD) that may enable a user to interact with a user interface of the computer device 10 (or mobile device).
Under the operating environment, the application provides a method for identifying characteristic points of a radial artery pressure waveform, particularly reflection points (E points). Fig. 2 is a schematic flow diagram of a method of identifying reflection points of a radial artery pressure waveform according to one embodiment of the present application. The method may comprise the steps of:
s100, a filtering step: filtering the radial pulse waveform;
s200 differential signal generation step: determining a maximum of the radial pulse waveform and generating a differential signal of the radial pulse waveform;
s300, a starting point searching step: determining a starting point of the radial pulse waveform based on a maximum of the differential signal;
s400, period determination: determining a period of the radial pulse waveform based on the starting point;
s500, waveform calibration: performing baseline wandering processing on the radial pulse waveform based on the starting point and calibrating the waveform; and
s600, feature point identification: reflection points of the calibrated waveform are identified.
By adopting the method, analysis can be completely carried out based on signal waveforms without adopting additional hardware, and compared with the complicated and fussy analysis process in the prior art, the method reduces the complexity and difficulty of analysis and improves the accuracy of analysis results. The method utilizes algorithm operations such as convolution filtering, first-order difference, convolution and the like, reduces the complex process of data processing to a certain extent, reduces the hardware computing resources required to be occupied, enables the feature points of the radial pulse waveform to be easier to extract, and is beneficial to the subsequent analysis processing process.
Optionally, the smoothing filtering step includes: and performing convolution filtering processing on the radial pulse waveform and a preset two-dimensional array. Fig. 3 is a general flow diagram of a method of identifying reflection points of a radial artery pressure waveform according to one embodiment of the present application. For example, the array may be an 11 x 14 two dimensional array: QG [11] [14 ]. And performing convolution filtering processing on the pulse waveform and the filter coefficient to obtain a filtered pulse waveform.
Wherein, each element in the specific two-dimensional array QG [11] [14] can be a specific numerical value in the following table one:
watch 1
Figure BDA0001821919510000071
Figure BDA0001821919510000081
Optionally, in the differential signal generating step, the calculation formula for generating the differential signal is:
vtCpDifSig[i]=vtSmoothedData[i+1]-vtSmoothedData[i-1]+2*(vtSmoothedData[i+2]-vtSmoothedData[i-2])
wherein vtCpDifSig represents a differential signal; vtsmootheneddata represents the filtered radial pulse waveform; i is an outer loop variable representing the ith bit signal value.
The method adopts a mathematical calculation mode to calculate the differential signal, can realize operation without adopting hardware of a differential circuit, and is convenient, quick, simple and realizable.
Optionally, in the differential signal generating step, the method for determining the waveform maximum value includes:
and (3) a step of segmented searching: dividing the differential signal waveform into eight sections, and respectively searching a maximum value point of each section of waveform;
and an average value calculation step: calculating an average value AvgSegmentMax of maximum values of the eight sections of waveforms;
a maximum value determining step: and searching a maximum value point in the whole differential signal waveform, and determining the value of the maximum value point as the maximum value point of the differential signal waveform under the condition that the amplitude of the maximum value point is greater than 0.6 AvgSegmentMax.
By adopting the method, the maximum value is searched in a segmented mode, the operation speed can be improved in a parallel processing mode, the obtained average value of the maximum value can be used as a threshold value to judge the maximum value of the whole waveform, and the judgment accuracy is improved.
Optionally, the starting point searching step includes: determining a search area based on the maximum value point of the differential signal, searching a forward zero crossing point in the search area, and taking the zero crossing point as a waveform starting point under the condition that the forward zero crossing point exists; otherwise, the minimum value point is used as the starting point of the waveform.
The process of determining the starting point by the difference signal maximum value point is as follows: searching forward for 200ms through a maximum value point (generally occurring after the starting point) of the differential signal to determine a search area, wherein the search area is used as a waveform starting point when a forward zero crossing point of the differential signal occurs; when there is no forward zero crossing point, the minimum point of the search area is taken as a starting point.
Optionally, in the period determining step, the waveform period is determined by calculating a distance difference between two starting points, and the waveform period is calculated by:
m_vtPeriod[i]=m_vtFootPos[i]-m_vtFootPos[i-1]
wherein m _ vtPeriod is a waveform period; m _ vtFootPos is used as a starting point; i is an outer loop variable representing the ith bit signal value.
Due to the adoption of the technical scheme, compared with the prior art, the method and the device reduce the processing process of pulse waveform data to a certain extent, reduce the consumption of computing resources, and are easier to extract waveform characteristic points, so that the subsequent analysis processing process of utilizing pulse waveforms is more convenient.
Referring to fig. 2 and 3, optionally, after the period determining step, the method further includes the peak point determining step, which includes searching for a waveform maximum point in an interval from the starting point to 2/3 of the period, and using the waveform maximum point as the peak point.
The method searches the maximum point of the waveform in the interval from the starting point to 2/3 of the period, so that the number of data processing can be reduced, and the calculation speed can be improved.
Optionally, in the waveform calibration step, a baseline wander process is performed according to a slope of the waveform start point value and a waveform offset, where the slope is calculated by:
Figure BDA0001821919510000091
the waveform offset calculation formula is as follows:
ΔvtSlope=vtSlope[i]*(j-m_vtFootPos[i])
wherein vtSlope represents the slope; Δ vtSlope represents the waveform offset; vtsmootheneddata represents the filtered radial pulse waveform; m _ vtFootPos is used as a starting point; i is an outer loop variable and represents the ith signal value; j is an inner loop variable.
Optionally, in one embodiment, the waveform scaling is scaling the amplitude of the radial artery waveform to have blood pressure units based on the systolic and diastolic pressures of the brachial artery.
Referring to fig. 2 and 3, optionally, the feature point identifying step includes:
and (3) convolution step: based on a set scale, obtaining the calibrated waveform fourth-order derivative vtScale2ConvSigNtoP through convolution; for example, a scale1 of 50, a scale2 of 100;
a reflection point determining step: searching the first reverse zero crossing point of the fourth derivative in the region from the starting point to 0.67 waveform periods, wherein the zero crossing point is a reflection point.
Fig. 4 is a schematic flow chart diagram of a method of searching for a reflection point according to an embodiment of the present application. By setting the scale, the fourth derivative of the waveform is obtained by convolution operation, so that the obtained waveform has less burrs and does not need subsequent filtering treatment. The reverse zero crossing of the fourth derivative is then searched for in the first 2/3 region of the entire waveform period to obtain the reflection point E of the waveform.
Wherein, the flow of obtaining fourth derivative by the double-scale convolution comprises:
(1) extracting waveform data after the first convolution; (2) and extracting waveform data after secondary convolution.
Extracting a convolution, and setting a double loop, wherein an outer loop variable is i, an inner loop variable is j, and when the absolute value of i-j is less than 0.5 Scale1, the value of the convolution is (j + initial point) amplitude value ((-8.0/3.0 (i-j) absolute value/Scale 1) + (5.0/6.0)); when the absolute value of Scale1> (i-j) is >0.5 Scale1, the first convolution value is (j + starting point) × (absolute value of 7.0/3.0 (i-j)/Scale 1) + (13.0/12.0)); when Scale1< ═ absolute value of (i-j), the first convolution value is (j + starting point) amplitude value ((-absolute value of 1.0/6.0 (i-j)/Scale 1) + (1.0/4.0)).
Extracting secondary convolution, wherein by setting a double loop, an outer loop variable is i, an inner loop variable is j, and when the absolute value of i-j is less than 0.5 Scale2, the value of the secondary convolution is the amplitude value of the primary convolution (j + starting point) (-8.0/3.0 (i-j) absolute value/Scale 2) + (5.0/6.0)); when the absolute value of Scale2> (i-j) is > Scale1/2, the value of the second convolution is the amplitude value of the first convolution (j + starting point) ((7.0/3.0: (i-j) absolute value/Scale 2) + (13.0/12.0)); when Scale2 ═ absolute value (i-j), quadratic convolution value ═ amplitude value of one convolution (j + starting point) (-absolute value of 1.0/6.0 ═ i-j)/Scale 2) + (1.0/4.0)); the second convolution value is the fourth derivative vtScale2convSigNtoP of the waveform.
FIG. 5 is a schematic illustration of a location of a reflection point on a radial pulse waveform according to an embodiment of the present application. The figure contains a complete single cycle radial pulse waveform consisting primarily of ascending and descending branches. Ascending branches and descending branches form a main wave 1, an incision on the descending branches is called as the descending isthmus 3, and a severe pre-pulsation wave 2, also called as an tidal wave, often appears between the main wave 1 and the descending isthmus 3. The severe beat wave 4, also called the descending central wave, appears next to the descending central isthmus 3. Both the above-mentioned waves and isthmus are the main components constituting the pulse wave.
Specifically, the point B of the dominant wave 1 trough is an aortic opening point, which represents a pressure turning point at which the aortic valve begins to open and marks the beginning of the rapid ejection period of the heart; point C is the peak of the main wave 1, which is the highest point of the aortic pressure and represents the highest point of the systolic arterial pressure; point E is the peak of the dicrotic wavefront 2, the left ventricular ejection stop, and the arterial pressure rise due to blood backflow. When the smooth muscle of the arterial wall of the human body is less and the elastic fiber is more, the main wave C forming the pulse wave is high and sharp, and the tidal wave E is not obvious. On the contrary, the smooth muscle of the artery wall is increased, the elastic fiber of the artery wall is less, the speed of the reflected wave is increased, and the tide wave E becomes obvious gradually. Therefore, the fluctuation change of the pulse wave characteristic points can directly reflect the change of the blood vessel resistance and the blood vessel wall elasticity of the human body, and is an objective important mark for measuring different physiological and pathological changes of the human body.
Example 2
There is also provided, in accordance with an embodiment of the present application, a device for identifying reflection points of a radial artery pressure waveform, the device corresponding to the method described in embodiment 1. Fig. 6 is a schematic block diagram of a device for identifying reflection points of a radial artery pressure waveform according to one embodiment of the present application. The apparatus may include:
a filtering module 100 configured to perform filtering processing on a radial artery pulse waveform;
a differential signal generation module 200 configured to generate a differential signal of the radial pulse waveform, determine a maximum of the differential signal;
a starting point search module 300 configured to determine a starting point of the radial pulse waveform based on a maximum value of the differential signal;
a period determination module 400 configured to determine a period of the radial pulse waveform based on the starting point;
a waveform calibration module 500 configured for de-baseline wander processing and waveform calibration of the radial pulse waveform based on the starting point; and
a feature point identification module 600 configured to identify reflection points of the calibrated waveform.
By adopting the device, analysis can be carried out completely based on signal waveforms, and can be realized without adopting extra hardware, compared with the complicated and fussy analysis process in the prior art, the complexity and difficulty of analysis are reduced, and the accuracy of an analysis result is also improved. The device utilizes algorithm operations such as convolution filtering, first-order difference, convolution and the like, reduces the complex process of data processing to a certain extent, reduces hardware computing resources required to be occupied, enables the feature points of the radial pulse waveform to be easier to extract, and is beneficial to the subsequent analysis and processing process.
Optionally, in the differential signal generating module, the calculation formula for generating the differential signal is:
vtCpDifSig[i]=vtSmoothedData[i+1]-vtSmoothedData[i-1]+2*(vtSmoothedData[i+2]-vtSmoothedData[i-2])
wherein vtCpDifSig represents a differential signal; vtsmootheneddata represents the filtered radial pulse waveform; i is an outer loop variable representing the ith bit signal value.
The device adopts the mode of mathematical computation to calculate the differential signal, and can realize the operation without adopting the hardware of a differential circuit, thereby being convenient, rapid and simple to realize.
Optionally, the differential signal generating module includes:
a segment searching module configured to divide the differential signal waveform into eight segments, and search for a maximum value point of each segment waveform respectively;
an average value calculating module configured to calculate an average value AvgSegmentMax of maxima of the eight-segment waveform;
a maximum determination module configured to search for a maximum point within the entire differential signal waveform, determine a value of the maximum point as the maximum point of the radial pulse waveform if the magnitude of the maximum point is greater than 0.6 at AvgSegmentMax.
By adopting the device, the maximum value is searched in a segmented mode, the operation speed can be improved in a parallel processing mode, the obtained average value of the maximum value can be used as a threshold value to judge the maximum value of the whole waveform, and the judgment accuracy is improved.
Optionally, the starting point searching module is further configured to: determining a search area based on the maximum value point of the differential signal, searching a forward zero crossing point in the search area, and taking the zero crossing point as a waveform starting point under the condition that the forward zero crossing point exists; otherwise, the minimum value point is used as the starting point of the waveform.
Optionally, in the period determining module, the waveform period is determined by calculating a distance difference between two starting points, and the waveform period is calculated by:
m_vtPeriod[i]=m_vtFootPos[i]-m_vtFootPos[i-1]
wherein m _ vtPeriod is a waveform period; m _ vtFootPos is used as a starting point; i is an outer loop variable representing the ith bit signal value.
Due to the adoption of the technical scheme, compared with the prior art, the method and the device reduce the processing process of pulse waveform data to a certain extent, reduce the consumption of computing resources, and are easier to extract waveform characteristic points, so that the subsequent analysis processing process of utilizing pulse waveforms is more convenient.
Optionally, after the period determining module, the apparatus may further include the peak point determining module further configured to search for a waveform maximum point in an interval from the starting point to 2/3 of the period, and use the waveform maximum point as the peak point.
The device can reduce the number of data processing and improve the calculation speed by searching the maximum value point of the waveform in the interval from the starting point to 2/3 of the period.
Optionally, in the waveform calibration module, the baseline wander is removed according to a slope of the waveform start point value and a waveform offset, where the slope is calculated by:
Figure BDA0001821919510000121
the waveform offset calculation formula is as follows:
ΔvtSlope=vtSlope[i]*(j-m_vtFootPos[i])
wherein vtSlope represents the slope; Δ vtSlope represents the waveform offset; vtsmootheneddata represents the filtered radial pulse waveform; m _ vtFootPos is used as a starting point; i is an outer loop variable and represents the ith signal value; j is an inner loop variable.
Optionally, in one embodiment, the waveform scaling is scaling the amplitude of the radial artery waveform to have blood pressure units based on the systolic and diastolic pressures of the brachial artery.
Optionally, the feature point identification module includes:
a convolution module configured to obtain the calibrated waveform fourth-order derivative vtScale2convSigNtoP through convolution based on a set scale; for example, a scale1 of 50, a scale2 of 100;
a reflection point determination module configured to search for a first inverted zero crossing of the fourth derivative from the starting point to 0.67 waveform period, the zero crossing being a reflection point.
The above-mentioned serial numbers of the embodiments of the present application are merely for description and do not represent the merits of the embodiments.
In the above embodiments of the present application, the descriptions of the respective embodiments have respective emphasis, and for parts that are not described in detail in a certain embodiment, reference may be made to related descriptions of other embodiments.
In the embodiments provided in the present application, it should be understood that the disclosed technology can be implemented in other ways. The above-described embodiments of the apparatus are merely illustrative, and for example, the division of the units is only one type of division of logical functions, and there may be other divisions when actually implemented, for example, a plurality of units or components may be combined or may be integrated into another system, or some features may be omitted, or not executed. In addition, the shown or discussed mutual coupling or direct coupling or communication connection may be an indirect coupling or communication connection through some interfaces, units or modules, and may be in an electrical or other form.
The units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the units can be selected according to actual needs to achieve the purpose of the solution of the embodiment.
In addition, functional units in the embodiments of the present application may be integrated into one processing unit, or each unit may exist alone physically, or two or more units are integrated into one unit. The integrated unit can be realized in a form of hardware, and can also be realized in a form of a software functional unit.
The integrated unit, if implemented in the form of a software functional unit and sold or used as a stand-alone product, may be stored in a computer readable storage medium. Based on such understanding, the technical solution of the present application may be substantially implemented or contributed to by the prior art, or all or part of the technical solution may be embodied in a software product, which is stored in a storage medium and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present application. And the aforementioned storage medium includes: a U-disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a removable hard disk, a magnetic or optical disk, and other various media capable of storing program codes.
Example 3
An aspect of embodiments of the present application provides a computing device, referring to fig. 7, comprising a memory 1120, a processor 1110 and a computer program stored in said memory 1120 and executable by said processor 1110, the computer program being stored in a space 1130 for program code in the memory 1120, the computer program, when executed by the processor 1110, implementing a method for performing any one of the above-mentioned method steps 1131 according to the present application.
An aspect of embodiments of the present application also provides a computer-readable storage medium. Referring to fig. 8, the computer readable storage medium comprises a storage unit for program code provided with a program 1131' for performing one of the above-mentioned method steps according to the present application, which program is executed by a processor.
An aspect of an embodiment of the present application also provides a computer program product containing instructions, including computer readable code, which when executed by a computing device, causes the computing device to perform the method as described above.
In the above embodiments, the implementation may be wholly or partially realized by software, hardware, firmware, or any combination thereof. When implemented in software, may be implemented in whole or in part in the form of a computer program product. The computer program product includes one or more computer instructions. When loaded and executed by a computer, cause the computer to perform, in whole or in part, the procedures or functions described in accordance with the embodiments of the application. The computer may be a general purpose computer, a special purpose computer, a network of computers, or other programmable device. The computer instructions may be stored in a computer readable storage medium or transmitted from one computer readable storage medium to another, for example, from one website site, computer, server, or data center to another website site, computer, server, or data center via wired (e.g., coaxial cable, fiber optic, Digital Subscriber Line (DSL)) or wireless (e.g., infrared, wireless, microwave, etc.). The computer-readable storage medium can be any available medium that can be accessed by a computer or a data storage device, such as a server, a data center, etc., that incorporates one or more of the available media. The usable medium may be a magnetic medium (e.g., floppy Disk, hard Disk, magnetic tape), an optical medium (e.g., DVD), or a semiconductor medium (e.g., Solid State Disk (SSD)), among others.
Those of skill would further appreciate that the various illustrative components and algorithm steps described in connection with the embodiments disclosed herein may be implemented as electronic hardware, computer software, or combinations of both, and that the various illustrative components and steps have been described above generally in terms of their functionality in order to clearly illustrate this interchangeability of hardware and software. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the implementation. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present application.
It will be understood by those skilled in the art that all or part of the steps in the method for implementing the above embodiments may be implemented by a program, and the program may be stored in a computer-readable storage medium, where the storage medium is a non-transitory medium, such as a random access memory, a read only memory, a flash memory, a hard disk, a solid state disk, a magnetic tape (magnetic tape), a floppy disk (floppy disk), an optical disk (optical disk), and any combination thereof.
The above description is only for the preferred embodiment of the present application, but the scope of the present application is not limited thereto, and any changes or substitutions that can be easily conceived by those skilled in the art within the technical scope of the present application should be covered within the scope of the present application. Therefore, the protection scope of the present application shall be subject to the protection scope of the claims.

Claims (8)

1. A method of identifying a reflection point of a radial artery pressure waveform, comprising:
a filtering step: filtering the radial pulse waveform;
a differential signal generation step: generating a differential signal of the radial pulse waveform, and determining a maximum value of the differential signal;
a starting point searching step: determining a starting point of the radial pulse waveform based on a maximum of the differential signal;
a period determining step: determining a period of the radial pulse waveform based on the starting point;
waveform calibration: performing baseline-removing drift processing on the radial pulse waveform based on the starting point and calibrating the waveform; and
a characteristic point identification step: identifying reflection points of the calibrated waveform;
in the waveform calibration step, baseline drift removal processing is performed according to the slope of the starting point value and the waveform offset, and the calculation formula of the slope is as follows:
Figure FDA0003257419510000011
the waveform offset calculation formula is as follows:
ΔvtSlope=vtSlope[i]*(j-m_vtFootPos[i])
wherein vtSlope represents the slope; Δ vtSlope represents the waveform offset; vtsmootheneddata represents the filtered radial pulse waveform; m _ vtFootPos is used as a starting point; m _ vtPeriod is a waveform period; i is an outer loop variable and represents the ith signal value; j is an inner loop variable;
the characteristic point identification step includes:
and (3) convolution step: acquiring the fourth derivative of the calibrated waveform through convolution based on a set scale; and
a reflection point determining step: searching the first reverse zero crossing point of the fourth derivative in the region from the starting point to 0.67 waveform periods, wherein the zero crossing point is a reflection point.
2. The method of claim 1, wherein the filtering step is a smoothing filtering step, the smoothing filtering step comprising: and performing convolution filtering processing on the radial pulse waveform and a preset two-dimensional array.
3. The method according to claim 1, wherein in the differential signal generating step, the calculation formula for generating the differential signal is:
vtCpDifSig[i]=vtSmoothedData[i+1]-vtSmoothedData[i-1]+2*(vtSmoothedData[i+2]-vtSmoothedData[i-2])
wherein vtCpDifSig represents a differential signal; vtsmootheneddata represents the filtered radial pulse waveform; i is an outer loop variable representing the ith bit signal value.
4. The method according to claim 1, wherein in the differential signal generation step, the method of determining the maximum value of the differential signal includes:
and (3) a step of segmented searching: dividing the differential signal waveform into eight sections, and respectively searching a maximum value point of each section of waveform;
and an average value calculation step: calculating an average value AvgSegmentMax of maximum values of the eight sections of waveforms; and
a maximum value determining step: and searching a maximum value point in the whole differential signal waveform, and determining the value of the maximum value point as the maximum value point of the differential signal under the condition that the amplitude of the maximum value point is greater than 0.6 AvgSegmentMax.
5. The method of claim 1, wherein the starting point searching step comprises: determining a search area based on the maximum value point of the differential signal, searching a forward zero crossing point in the search area, and taking the zero crossing point as a waveform starting point under the condition that the forward zero crossing point exists; otherwise, the minimum value point is used as the starting point of the waveform.
6. The method according to claim 1, wherein in the period determining step, the waveform period is determined by calculating a distance difference between two starting points, and the waveform period is calculated by:
m_vtPeriod[i]=m_vtFootPos[i]-m_vtFootPos[i-1]
wherein m _ vtPeriod is a waveform period; m _ vtFootPos is used as a starting point; i is an outer loop variable representing the ith bit signal value.
7. An apparatus for identifying reflection points of a radial artery pressure waveform, comprising:
a filtering module configured to perform filtering processing on a radial artery pulse waveform;
a differential signal generation module configured to generate a differential signal of the radial pulse waveform, determine a maximum of the differential signal;
a starting point search module configured to determine a starting point of the radial pulse waveform based on a maximum value of the differential signal;
a period determination module configured to determine a period of the radial pulse waveform based on the starting point;
a waveform calibration module configured to perform de-baseline wander processing on the radial pulse waveform based on the starting point and calibrate the waveform; and
a feature point identification module configured to identify reflection points of the calibrated waveform;
and performing baseline wandering removal processing according to the slope of the starting point value and the waveform offset, wherein the calculation formula of the slope is as follows:
Figure FDA0003257419510000031
the waveform offset calculation formula is as follows:
ΔvtSlope=vtSlope[i]*(j-m_vtFootPos[i])
wherein vtSlope represents the slope; Δ vtSlope represents the waveform offset; vtsmootheneddata represents the filtered radial pulse waveform; m _ vtFootPos is used as a starting point; m _ vtPeriod is a waveform period; i is an outer loop variable and represents the ith signal value; j is an inner loop variable;
the step of feature point identification comprises:
and (3) convolution step: acquiring the fourth derivative of the calibrated waveform through convolution based on a set scale; and
a reflection point determining step: searching the first reverse zero crossing point of the fourth derivative in the region from the starting point to 0.67 waveform periods, wherein the zero crossing point is a reflection point.
8. The apparatus of claim 7, wherein the differential signal generating module comprises:
a segment searching module configured to divide the differential signal waveform into eight segments, and search for a maximum value point of each segment waveform respectively;
an average value calculating module configured to calculate an average value AvgSegmentMax of maxima of the eight-segment waveform; and
a maximum determination module configured to search for a maximum point within the entire differential signal waveform, determine a value of the maximum point as the maximum point of the differential signal waveform if the magnitude of the maximum point is greater than 0.6 at AvgSegmentMax.
CN201811168960.3A 2018-10-08 2018-10-08 Method and device for identifying reflection point of radial artery pressure waveform Active CN110477889B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811168960.3A CN110477889B (en) 2018-10-08 2018-10-08 Method and device for identifying reflection point of radial artery pressure waveform

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811168960.3A CN110477889B (en) 2018-10-08 2018-10-08 Method and device for identifying reflection point of radial artery pressure waveform

Publications (2)

Publication Number Publication Date
CN110477889A CN110477889A (en) 2019-11-22
CN110477889B true CN110477889B (en) 2021-11-02

Family

ID=68545749

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811168960.3A Active CN110477889B (en) 2018-10-08 2018-10-08 Method and device for identifying reflection point of radial artery pressure waveform

Country Status (1)

Country Link
CN (1) CN110477889B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112932475B (en) * 2021-02-01 2023-02-21 武汉泰利美信医疗科技有限公司 Method and device for calculating blood oxygen saturation, electronic equipment and storage medium

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5027641A (en) * 1989-02-23 1991-07-02 Costello Jr Leo F Oscillometric non-invasive blood pressure simulator
CN102307520A (en) * 2009-02-03 2012-01-04 欧姆龙健康医疗事业株式会社 Pulse wave analyzing device and pulse wave analyzing method
CN103169456A (en) * 2013-03-29 2013-06-26 深圳职业技术学院 Processing method and processing system for pulse wave signals
CN106691406A (en) * 2017-01-05 2017-05-24 大连理工大学 Detection method of vascular elasticity and blood pressure based on single probe photoplethysmography pulse wave
CN107898443A (en) * 2017-11-21 2018-04-13 深圳先进技术研究院 Dicrotic pulse wave detecting method, device and computer-readable storage medium

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5027641A (en) * 1989-02-23 1991-07-02 Costello Jr Leo F Oscillometric non-invasive blood pressure simulator
CN102307520A (en) * 2009-02-03 2012-01-04 欧姆龙健康医疗事业株式会社 Pulse wave analyzing device and pulse wave analyzing method
CN103169456A (en) * 2013-03-29 2013-06-26 深圳职业技术学院 Processing method and processing system for pulse wave signals
CN106691406A (en) * 2017-01-05 2017-05-24 大连理工大学 Detection method of vascular elasticity and blood pressure based on single probe photoplethysmography pulse wave
CN107898443A (en) * 2017-11-21 2018-04-13 深圳先进技术研究院 Dicrotic pulse wave detecting method, device and computer-readable storage medium

Also Published As

Publication number Publication date
CN110477889A (en) 2019-11-22

Similar Documents

Publication Publication Date Title
EP2992820B1 (en) Measuring blood pressure
CN107928643B (en) Apparatus and method for estimating biometric information
CN107530005B (en) Method and device for deriving mean arterial pressure of a subject
US10285651B2 (en) On-demand heart rate estimation based on optical measurements
CN109872820B (en) Method, device, equipment and storage medium for measuring blood pressure without cuff
CN106659404B (en) Continuous blood pressure measuring method, device and equipment
WO2018072212A1 (en) Calibration method for blood pressure measuring device, and blood pressure measuring device
CN108992054B (en) Pulse signal peak point detection method and device
EP2825977A2 (en) System and method for analysis and reconstruction of variable pulse-width signals having low sampling rates
CN110755059A (en) Blood pressure waveform monitoring method and device
KR20200021208A (en) Apparatus and method for estimating blood pressure
CN104068841B (en) A kind of measuring method and device measuring Indices of Systolic Time parameter
CN110477889B (en) Method and device for identifying reflection point of radial artery pressure waveform
CN103549951A (en) P wave information measuring method based on electrocardiogram 12 lead correlation computing
CN110495863B (en) Method and device for identifying characteristic points of radial artery pressure waveform central isthmus
CN110840427A (en) Continuous blood pressure measuring method, device and equipment based on volume pulse wave signals
CN110200642A (en) A kind of measurement method and terminal of cognitive load and psychological pressure
CN110693474B (en) Method and device for identifying characteristic points of radial artery pressure waveform dicrotic wave
CN110840426A (en) Device for estimating blood pressure
WO2024041149A1 (en) Method and device for calibrating blood pressure measurement device, electronic device, and storage medium
CN108652604A (en) It is a kind of based on electrocardiosignal without air bag blood pressure detecting method and system
CN111000543A (en) Device for estimating blood pressure
EP3435853B1 (en) A method of deriving systolic blood pressure and/or diastolic blood pressure of a subject
CN103720462A (en) Pulse wave signal analyzing method and device
CN114469018B (en) Pulse diagnosis fingering recommendation method and device, electronic equipment and storage medium

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
CB03 Change of inventor or designer information

Inventor after: Zhang Yongliang

Inventor after: Ye Jun

Inventor after: Li Daoqing

Inventor after: Zhang Qilian

Inventor before: Li Daoqing

Inventor before: Zhang Qilian

CB03 Change of inventor or designer information
TA01 Transfer of patent application right

Effective date of registration: 20210414

Address after: 230000 1707, building 1, science and Technology Industrial Park (Zhongwang sub Park), 198 Mingzhu Avenue, high tech Zone, Hefei City, Anhui Province

Applicant after: Hefei Yiyang Health Technology Co.,Ltd.

Applicant after: BEIJING DONGLIANG HEALTH TECHNOLOGY Co.,Ltd.

Address before: 100038 room 2033, 2 / F, 21 YANGFANGDIAN Road, Haidian District, Beijing

Applicant before: Heyunli Health Technology (Beijing) Co.,Ltd.

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant